1616#include " PWGLF/DataModel/LFKinkDecayTables.h"
1717#include " PWGLF/Utils/svPoolCreator.h"
1818
19- #include " Common/DataModel/Centrality.h"
2019#include " Common/Core/RecoDecay.h"
2120#include " Common/Core/trackUtilities.h"
21+ #include " Common/DataModel/Centrality.h"
2222
2323#include < CCDB/BasicCCDBManager.h>
2424#include < CommonConstants/MathConstants.h>
@@ -322,7 +322,7 @@ struct kinkBuilder {
322322 hSelKinkedTrackQA->GetXaxis ()->SetBinLabel (11 , " MothLastLayerRadiusCheck" );
323323 hMothDaughSignsInit = qaRegistry.add <TH2>(" hMothDaughSignsInit" , " hMothDaughSignsInit; Sign Mother; Sign Daughter" , {HistType::kTH2F , {{3 , -1.5 , 1.5 }, {3 , -1.5 , 1.5 }}});
324324 hMothDaughSignsFinal = qaRegistry.add <TH2>(" hMothDaughSignsFinal" , " hMothDaughSignsFinal; Sign Mother; Sign Daughter" , {HistType::kTH2F , {{3 , -1.5 , 1.5 }, {3 , -1.5 , 1.5 }}});
325- hZDiff = qaRegistry.add <TH2>(" hZDiff" , " hZDiff; #Delta z (#mu m);(Q_{Mother}, Q_{Daughter})" , HistType::kTH2F , {zDiffBins, {4 , -0.5 , 3.5 }});
325+ hZDiff = qaRegistry.add <TH2>(" hZDiff" , " hZDiff; #Delta z (#mu m);(Q_{Mother}, Q_{Daughter})" , HistType::kTH2F , {zDiffBins, {4 , -0.5 , 3.5 }});
326326 hZDiff->GetYaxis ()->SetBinLabel (1 , " (+,+)" );
327327 hZDiff->GetYaxis ()->SetBinLabel (2 , " (-,-)" );
328328 hZDiff->GetYaxis ()->SetBinLabel (3 , " (+,-)" );
@@ -355,12 +355,12 @@ struct kinkBuilder {
355355
356356 if (doprocessMc || doprocessMcWCent) {
357357 if (skipBkgCands) {
358- hRecCandidates = qaRegistry.add <TH2>(" hRecCandidates" , " hRecCandidates;Counts;" , {HistType::kTH2F , {{kNMatchedDecays , -0.5 , static_cast <float >(kNMatchedDecays )- 0.5 }, absPtAxis}});
358+ hRecCandidates = qaRegistry.add <TH2>(" hRecCandidates" , " hRecCandidates;Counts;" , {HistType::kTH2F , {{kNMatchedDecays , -0.5 , static_cast <float >(kNMatchedDecays ) - 0.5 }, absPtAxis}});
359359 hRecCandidates->GetXaxis ()->SetBinLabel (1 , " #Sigma^{-} #rightarrow n#pi^{-}" );
360360 hRecCandidates->GetXaxis ()->SetBinLabel (2 , " #Sigma^{+} #rightarrow n#pi^{+}" );
361361 hRecCandidates->GetXaxis ()->SetBinLabel (3 , " #Sigma^{+} #rightarrow p#pi^{0}" );
362362 }
363- hGenCandidates = qaRegistry.add <TH2>(" hGenCandidates" , " hGenCandidates;Counts;" , {HistType::kTH2F , {{kNMatchedDecays , -0.5 , static_cast <float >(kNMatchedDecays )- 0.5 }, absPtAxis}});
363+ hGenCandidates = qaRegistry.add <TH2>(" hGenCandidates" , " hGenCandidates;Counts;" , {HistType::kTH2F , {{kNMatchedDecays , -0.5 , static_cast <float >(kNMatchedDecays ) - 0.5 }, absPtAxis}});
364364 hGenCandidates->GetXaxis ()->SetBinLabel (1 , " #Sigma^{-} #rightarrow n#pi^{-}" );
365365 hGenCandidates->GetXaxis ()->SetBinLabel (2 , " #Sigma^{+} #rightarrow n#pi^{+}" );
366366 hGenCandidates->GetXaxis ()->SetBinLabel (3 , " #Sigma^{+} #rightarrow p#pi^{0}" );
@@ -437,15 +437,15 @@ struct kinkBuilder {
437437 if (candidate.itsNClsInnerBarrel () != 0 )
438438 return false ;
439439 hSelDaugQA->Fill (4 .f , isPositive);
440-
440+
441441 if (candidate.itsNCls () >= 4 )
442442 return false ;
443443 hSelDaugQA->Fill (5 .f , isPositive);
444-
444+
445445 if (candidate.tpcNClsCrossedRows () <= 0.8 * candidate.tpcNClsFindable ())
446446 return false ;
447447 hSelDaugQA->Fill (6 .f , isPositive);
448-
448+
449449 if (candidate.tpcNClsFound () <= nTPCClusMinDaug)
450450 return false ;
451451 hSelDaugQA->Fill (7 .f , isPositive);
@@ -463,7 +463,7 @@ struct kinkBuilder {
463463 auto trackDaug = tracks.rawIteratorAt (svCand.tr1Idx );
464464
465465 // Fill Selections QA histo
466- int chargeCombSvCand = 2 * unlikeSignBkg + (trackMoth.sign () == -1 ? 1 : 0 );
466+ int chargeCombSvCand = 2 * unlikeSignBkg + (trackMoth.sign () == -1 ? 1 : 0 );
467467 hSelKinkedTrackQA->Fill (0 .f , chargeCombSvCand); // all candidates bin
468468 hMothDaughSignsInit->Fill (trackMoth.sign (), trackDaug.sign ());
469469
@@ -641,8 +641,9 @@ struct kinkBuilder {
641641 LOG (info) << " Task initialized for run " << mRunNumber << " with magnetic field " << mBz << " kZG" ;
642642 }
643643
644- template <typename TColls, typename TTracks, typename TAmbiTracks>
645- void buildSvPool (const TColls& collisions, const TTracks& tracks, const TAmbiTracks& ambiguousTracks, const aod::BCs& bcs) {
644+ template <typename TColls, typename TTracks, typename TAmbiTracks>
645+ void buildSvPool (const TColls& collisions, const TTracks& tracks, const TAmbiTracks& ambiguousTracks, const aod::BCs& bcs)
646+ {
646647 svCreator.clearPools ();
647648 svCreator.fillBC2Coll (collisions, bcs);
648649
@@ -661,8 +662,9 @@ struct kinkBuilder {
661662 }
662663 }
663664
664- template <typename TColls, typename TTracks>
665- void fillOutputsData (const TColls& collisions, const TTracks& tracks, const aod::AmbiguousTracks& ambiTracks, const aod::BCs& bcs) {
665+ template <typename TColls, typename TTracks>
666+ void fillOutputsData (const TColls& collisions, const TTracks& tracks, const aod::AmbiguousTracks& ambiTracks, const aod::BCs& bcs)
667+ {
666668 kinkCandidates.clear ();
667669
668670 buildSvPool (collisions, tracks, ambiTracks, bcs);
@@ -703,15 +705,16 @@ struct kinkBuilder {
703705 }
704706 PROCESS_SWITCH (kinkBuilder, processDataWCentSel, " Data processing with centrality selection" , false );
705707
706- template <bool checkKinkDaugPdg, typename TMother>
707- int matchKinkDecay (const TMother& motherPart, const aod::McParticles& mcParticles) {
708+ template <bool checkKinkDaugPdg, typename TMother>
709+ int matchKinkDecay (const TMother& motherPart, const aod::McParticles& mcParticles)
710+ {
708711 int pdgMother = motherPart.pdgCode ();
709712 int8_t sign = 0 ;
710713 int pdgCodeNeutralDaug{-1 }, pdgCodeChargedDaug{-1 };
711714 std::array<int , 2 > finState = {-1 , -1 };
712715 switch (std::abs (pdgMother)) {
713716 case PDG_t::kSigmaMinus : {
714- // Swap the sign of the neutral decay products in case of anti-particles
717+ // Swap the sign of the neutral decay products in case of anti-particles
715718 pdgCodeNeutralDaug = (pdgMother > 0 ) ? +PDG_t::kNeutron : -PDG_t::kNeutron ;
716719 pdgCodeChargedDaug = (pdgMother > 0 ) ? +PDG_t::kPiMinus : +PDG_t::kPiPlus ;
717720 finState = {pdgCodeChargedDaug, pdgCodeNeutralDaug}; // Both decay channels have the same neutral daughter
@@ -723,7 +726,7 @@ struct kinkBuilder {
723726 case PDG_t::kSigmaPlus : {
724727 // Swap the sign of the neutral decay products in case of anti-particles
725728 pdgCodeNeutralDaug = (pdgMother > 0 ) ? +PDG_t::kNeutron : -PDG_t::kNeutron ;
726- pdgCodeChargedDaug = (pdgMother > 0 ) ? +PDG_t::kPiPlus : +PDG_t::kPiMinus ;
729+ pdgCodeChargedDaug = (pdgMother > 0 ) ? +PDG_t::kPiPlus : +PDG_t::kPiMinus ;
727730 finState = {pdgCodeChargedDaug, pdgCodeNeutralDaug};
728731 if (RecoDecay::isMatchedMCGen<false , true >(mcParticles, motherPart, pdgMother, finState, true , &sign, DepthMcMatchMax)) {
729732 return kSigmaPlusToPiPlusNeutron ;
@@ -746,8 +749,9 @@ struct kinkBuilder {
746749 return -1 ;
747750 }
748751
749- template <typename TColls, typename TTracks>
750- void fillOutputsMc (const TColls& mcRecoCollisions, const TTracks& tracksMc, const aod::AmbiguousTracks& ambiTracksMc, const aod::BCs& bcs, const aod::McParticles& mcParticles) {
752+ template <typename TColls, typename TTracks>
753+ void fillOutputsMc (const TColls& mcRecoCollisions, const TTracks& tracksMc, const aod::AmbiguousTracks& ambiTracksMc, const aod::BCs& bcs, const aod::McParticles& mcParticles)
754+ {
751755 kinkCandidates.clear ();
752756
753757 buildSvPool (mcRecoCollisions, tracksMc, ambiTracksMc, bcs);
@@ -828,7 +832,7 @@ struct kinkBuilder {
828832 McRecoCollisions const & mcRecoCollisions,
829833 aod::McParticles const & mcParticles,
830834 TracksFullMc const & tracksMc,
831- aod::AmbiguousTracks const & ambiTracksMc,
835+ aod::AmbiguousTracks const & ambiTracksMc,
832836 aod::BCs const & bcs)
833837 {
834838 fillOutputsMc (mcRecoCollisions, tracksMc, ambiTracksMc, bcs, mcParticles);
@@ -839,7 +843,7 @@ struct kinkBuilder {
839843 McRecoCollisionsCentSel const & mcRecoCollisions,
840844 aod::McParticles const & mcParticles,
841845 TracksFullMc const & tracksMc,
842- aod::AmbiguousTracks const & ambiTracksMc,
846+ aod::AmbiguousTracks const & ambiTracksMc,
843847 aod::BCs const & bcs)
844848 {
845849 fillOutputsMc (mcRecoCollisions, tracksMc, ambiTracksMc, bcs, mcParticles);
0 commit comments