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 // LOG(info) << "[kink] ";
668670 // LOG(info) << "[kink] *****************************************";
@@ -707,15 +709,16 @@ struct kinkBuilder {
707709 }
708710 PROCESS_SWITCH (kinkBuilder, processDataWCentSel, " Data processing with centrality selection" , false );
709711
710- template <bool checkKinkDaugPdg, typename TMother>
711- int matchKinkDecay (const TMother& motherPart, const aod::McParticles& mcParticles) {
712+ template <bool checkKinkDaugPdg, typename TMother>
713+ int matchKinkDecay (const TMother& motherPart, const aod::McParticles& mcParticles)
714+ {
712715 int pdgMother = motherPart.pdgCode ();
713716 int8_t sign = 0 ;
714717 int pdgCodeNeutralDaug{-1 }, pdgCodeChargedDaug{-1 };
715718 std::array<int , 2 > finState = {-1 , -1 };
716719 switch (std::abs (pdgMother)) {
717720 case PDG_t::kSigmaMinus : {
718- // Swap the sign of the neutral decay products in case of anti-particles
721+ // Swap the sign of the neutral decay products in case of anti-particles
719722 pdgCodeNeutralDaug = (pdgMother > 0 ) ? +PDG_t::kNeutron : -PDG_t::kNeutron ;
720723 pdgCodeChargedDaug = (pdgMother > 0 ) ? +PDG_t::kPiMinus : +PDG_t::kPiPlus ;
721724 finState = {pdgCodeChargedDaug, pdgCodeNeutralDaug}; // Both decay channels have the same neutral daughter
@@ -727,7 +730,7 @@ struct kinkBuilder {
727730 case PDG_t::kSigmaPlus : {
728731 // Swap the sign of the neutral decay products in case of anti-particles
729732 pdgCodeNeutralDaug = (pdgMother > 0 ) ? +PDG_t::kNeutron : -PDG_t::kNeutron ;
730- pdgCodeChargedDaug = (pdgMother > 0 ) ? +PDG_t::kPiPlus : +PDG_t::kPiMinus ;
733+ pdgCodeChargedDaug = (pdgMother > 0 ) ? +PDG_t::kPiPlus : +PDG_t::kPiMinus ;
731734 finState = {pdgCodeChargedDaug, pdgCodeNeutralDaug};
732735 if (RecoDecay::isMatchedMCGen<false , true >(mcParticles, motherPart, pdgMother, finState, true , &sign, DepthMcMatchMax)) {
733736 return kSigmaPlusToPiPlusNeutron ;
@@ -750,8 +753,9 @@ struct kinkBuilder {
750753 return -1 ;
751754 }
752755
753- template <typename TColls, typename TTracks>
754- void fillOutputsMc (const TColls& mcRecoCollisions, const TTracks& tracksMc, const aod::AmbiguousTracks& ambiTracksMc, const aod::BCs& bcs, const aod::McParticles& mcParticles) {
756+ template <typename TColls, typename TTracks>
757+ void fillOutputsMc (const TColls& mcRecoCollisions, const TTracks& tracksMc, const aod::AmbiguousTracks& ambiTracksMc, const aod::BCs& bcs, const aod::McParticles& mcParticles)
758+ {
755759 kinkCandidates.clear ();
756760
757761 buildSvPool (mcRecoCollisions, tracksMc, ambiTracksMc, bcs);
@@ -832,7 +836,7 @@ struct kinkBuilder {
832836 McRecoCollisions const & mcRecoCollisions,
833837 aod::McParticles const & mcParticles,
834838 TracksFullMc const & tracksMc,
835- aod::AmbiguousTracks const & ambiTracksMc,
839+ aod::AmbiguousTracks const & ambiTracksMc,
836840 aod::BCs const & bcs)
837841 {
838842 fillOutputsMc (mcRecoCollisions, tracksMc, ambiTracksMc, bcs, mcParticles);
@@ -843,7 +847,7 @@ struct kinkBuilder {
843847 McRecoCollisionsCentSel const & mcRecoCollisions,
844848 aod::McParticles const & mcParticles,
845849 TracksFullMc const & tracksMc,
846- aod::AmbiguousTracks const & ambiTracksMc,
850+ aod::AmbiguousTracks const & ambiTracksMc,
847851 aod::BCs const & bcs)
848852 {
849853 fillOutputsMc (mcRecoCollisions, tracksMc, ambiTracksMc, bcs, mcParticles);
0 commit comments