@@ -326,7 +326,8 @@ struct sigma0builder {
326326 doprocessMonteCarlo +
327327 doprocessMonteCarloWithTOF +
328328 doprocessV0QA +
329- doprocessV0MCQA > 1 ) {
329+ doprocessV0MCQA >
330+ 1 ) {
330331 LOGF (fatal, " You have enabled more than one process function. Please check your configuration! Aborting now." );
331332 }
332333
@@ -485,7 +486,6 @@ struct sigma0builder {
485486 histos.get <TH1>(HIST (" KStarSel/hSelectionStatistics" ))->GetXaxis ()->SetBinLabel (3 , " KStar Y Window" );
486487
487488 histos.add (" KStarSel/hKStarMassSelected" , " hKStarMassSelected" , kTH1F , {axisKStarMass});
488-
489489 }
490490
491491 if (doAssocStudy && (doprocessMonteCarlo || doprocessMonteCarloWithTOF)) {
@@ -507,7 +507,7 @@ struct sigma0builder {
507507 histos.add (" MCQA/hPhotonMotherSize" , " hPhotonMotherSize" , kTH1D , {{10 , -0 .5f , +9 .5f }});
508508 histos.add (" MCQA/hPhotonMCProcess" , " hPhotonMCProcess" , kTH1D , {{50 , -0 .5f , 49 .5f }});
509509 histos.add (" MCQA/hPhotonMotherMCProcess" , " hPhotonMotherMCProcess" , kTH1D , {{50 , -0 .5f , 49 .5f }});
510-
510+
511511 // photons for kstar analysis
512512 histos.add (" MCQA/h2dPhoton_KShortNMothersVsPDG" , " h2dPhoton_KShortNMothersVsPDG" , kTHnSparseD , {{10 , -0 .5f , +9 .5f }, {10001 , -5000 .5f , +5000 .5f }});
513513 histos.add (" MCQA/h2dPhoton_KShortNMothersVsMCProcess" , " h2dPhoton_KShortNMothersVsMCProcess" , kTH2D , {{10 , -0 .5f , +9 .5f }, {50 , -0 .5f , 49 .5f }});
@@ -648,7 +648,6 @@ struct sigma0builder {
648648 histos.add (" V0QA/h3dYKShortMass" , " h3dYKShortMass" , kTH3D , {axisRapidity, axisPt, axisK0SMass});
649649 histos.add (" V0QA/h3dYRKShortMass" , " h3dYRKShortMass" , kTH3D , {axisRapidity, axisRadius, axisK0SMass});
650650 histos.add (" V0QA/h2dTrueKStarKShortMass" , " h2dTrueKStarKShortMass" , kTH2D , {axisPt, axisK0SMass});
651-
652651 }
653652
654653 if (doprocessV0Generated) {
@@ -678,7 +677,7 @@ struct sigma0builder {
678677 histos.add (" V0QA/h2dGenAntiLambdaVsMultMC" , " h2dGenAntiLambdaVsMultMC" , kTH2D , {axisNch, axisPt});
679678
680679 histos.add (" V0QA/h2dGenKShort" , " h2dGenKShort" , kTH2D , {axisCentrality, axisPt});
681-
680+
682681 histos.add (" V0QA/h2dGenKShortVsMultMC_RecoedEvt" , " h2dGenKShortVsMultMC_RecoedEvt" , kTH2D , {axisNch, axisPt});
683682
684683 histos.add (" V0QA/h2dGenKShortVsMultMC" , " h2dGenKShortVsMultMC" , kTH2D , {axisNch, axisPt});
@@ -810,7 +809,7 @@ struct sigma0builder {
810809 if ((v01MC.pdgCode () == PDG_t::kGamma ) && (v01MC.pdgCodeMother () == PDG_t::kSigma0 ) && (v02MC.pdgCode () == PDG_t::kLambda0 ) && (v02MC.pdgCodeMother () == PDG_t::kSigma0 ) && (v01.motherMCPartId () == v02.motherMCPartId ()))
811810 fIsSigma0 = true ;
812811
813- bool fIsKStar = false ;
812+ bool fIsKStar = false ;
814813 if ((v01MC.pdgCode () == PDG_t::kGamma ) && (v01MC.pdgCodeMother () == o2::constants::physics::Pdg::kK0Star892 ) && (v02MC.pdgCode () == PDG_t::kK0Short ) && (v02MC.pdgCodeMother () == o2::constants::physics::Pdg::kK0Star892 ) && (v01.motherMCPartId () == v02.motherMCPartId ()))
815814 fIsKStar = true ;
816815
@@ -846,8 +845,8 @@ struct sigma0builder {
846845 auto const & MCMothersList_v02 = MCParticle_v02.template mothers_as <aod::McParticles>();
847846
848847 if (!MCMothersList_v01.empty () && !MCMothersList_v02.empty ()) { // Are there mothers?
849- auto const & MCMother_v01 = MCMothersList_v01.front (); // First mother
850- auto const & MCMother_v02 = MCMothersList_v02.front (); // First mother
848+ auto const & MCMother_v01 = MCMothersList_v01.front (); // First mother
849+ auto const & MCMother_v02 = MCMothersList_v02.front (); // First mother
851850
852851 if (MCMother_v01.globalIndex () == MCMother_v02.globalIndex ()) { // Is it the same mother?
853852
@@ -922,11 +921,11 @@ struct sigma0builder {
922921 histos.fill (HIST (" MCQA/hSigma0MCCheck" ), 3 ); // mismatch
923922
924923 // Check association correctness
925- if (fIsKStar && (MCinfo.V0PairPDGCode == o2::constants::physics::Pdg::kK0Star892 ))
924+ if (fIsKStar && (MCinfo.V0PairPDGCode == o2::constants::physics::Pdg::kK0Star892 ))
926925 histos.fill (HIST (" MCQA/hKStarMCCheck" ), 1 ); // match
927926 if (fIsKStar && !(MCinfo.V0PairPDGCode == o2::constants::physics::Pdg::kK0Star892 ))
928927 histos.fill (HIST (" MCQA/hKStarMCCheck" ), 2 ); // mismatch
929- if (!fIsKStar && (MCinfo.V0PairPDGCode == o2::constants::physics::Pdg::kK0Star892 ))
928+ if (!fIsKStar && (MCinfo.V0PairPDGCode == o2::constants::physics::Pdg::kK0Star892 ))
930929 histos.fill (HIST (" MCQA/hKStarMCCheck" ), 3 ); // mismatch
931930 }
932931 }
@@ -1478,7 +1477,7 @@ struct sigma0builder {
14781477 // Pi0
14791478 if (fillPi0Tables && MCGenInfo.IsPi0 ) {
14801479 pi0Gens (MCGenInfo.IsProducedByGenerator , MCGenInfo.MCPt , mcParticle.y ()); // optional table to store generated pi0 candidates. Be careful, this is a large table!
1481- pi0GenCollRefs (MCGenInfo.MCCollId ); // link to stramccollision table
1480+ pi0GenCollRefs (MCGenInfo.MCCollId ); // link to stramccollision table
14821481 }
14831482
14841483 // Sigma0/ASigma0
@@ -1546,15 +1545,15 @@ struct sigma0builder {
15461545
15471546 // Lambda/ALambda-specific
15481547 histos.fill (HIST (MainDir[mode]) + HIST (" /hLambdaLifeTime" ), fLambdaLifeTime );
1549- histos.fill (HIST (MainDir[mode]) + HIST (" /hLambdaY" ), v0.yLambda ());
1548+ histos.fill (HIST (MainDir[mode]) + HIST (" /hLambdaY" ), v0.yLambda ());
15501549 histos.fill (HIST (MainDir[mode]) + HIST (" /hLambdaMass" ), v0.mLambda ());
15511550 histos.fill (HIST (MainDir[mode]) + HIST (" /hALambdaMass" ), v0.mAntiLambda ());
15521551 histos.fill (HIST (MainDir[mode]) + HIST (" /h2dMassLambdaVsK0S" ), v0.mLambda (), v0.mK0Short ());
15531552 histos.fill (HIST (MainDir[mode]) + HIST (" /h2dMassLambdaVsGamma" ), v0.mLambda (), v0.mGamma ());
15541553
15551554 // KShort-specific
15561555 histos.fill (HIST (MainDir[mode]) + HIST (" /hKShortLifeTime" ), fKShortLifeTime );
1557- histos.fill (HIST (MainDir[mode]) + HIST (" /hKShortY" ), v0.yK0Short ());
1556+ histos.fill (HIST (MainDir[mode]) + HIST (" /hKShortY" ), v0.yK0Short ());
15581557 histos.fill (HIST (MainDir[mode]) + HIST (" /hKShortMass" ), v0.mK0Short ());
15591558 histos.fill (HIST (MainDir[mode]) + HIST (" /h2dMassK0SvsLambda" ), v0.mK0Short (), v0.mK0Short ());
15601559 histos.fill (HIST (MainDir[mode]) + HIST (" /h2dMassK0SVsGamma" ), v0.mK0Short (), v0.mGamma ());
@@ -1568,12 +1567,11 @@ struct sigma0builder {
15681567 // Optional MC selection
15691568 if (photonSelections.doMCAssociation ) {
15701569 if constexpr (requires { gamma.motherMCPartId (); }) {
1571- if (!gamma.has_v0MCCore ())
1570+ if (!gamma.has_v0MCCore ())
15721571 return false ;
15731572 auto gammaMC = gamma.template v0MCCore_as <soa::Join<aod::V0MCCores, aod::V0MCCollRefs>>();
15741573 if (gammaMC.pdgCode () != PDG_t::kGamma )
15751574 return false ;
1576-
15771575 }
15781576 }
15791577
@@ -1666,12 +1664,12 @@ struct sigma0builder {
16661664 // Optional MC selection
16671665 if (lambdaSelections.doMCAssociation ) {
16681666 if constexpr (requires { lambda.motherMCPartId (); }) {
1669- if (!lambda.has_v0MCCore ())
1667+ if (!lambda.has_v0MCCore ())
16701668 return false ;
16711669
16721670 auto lambdaMC = lambda.template v0MCCore_as <soa::Join<aod::V0MCCores, aod::V0MCCollRefs>>();
1673- if (TMath::Abs (lambdaMC.pdgCode ()) != PDG_t::kLambda0 )
1674- return false ;
1671+ if (TMath::Abs (lambdaMC.pdgCode ()) != PDG_t::kLambda0 )
1672+ return false ;
16751673 }
16761674 }
16771675
@@ -1761,22 +1759,22 @@ struct sigma0builder {
17611759 // Optional MC selection
17621760 if (kshortSelections.doMCAssociation ) {
17631761 if constexpr (requires { kshort.motherMCPartId (); }) {
1764- if (!kshort.has_v0MCCore ())
1762+ if (!kshort.has_v0MCCore ())
17651763 return false ;
17661764
17671765 auto kshortMC = kshort.template v0MCCore_as <soa::Join<aod::V0MCCores, aod::V0MCCollRefs>>();
17681766 if (TMath::Abs (kshortMC.pdgCode ()) != PDG_t::kK0Short )
17691767 return false ;
1770- }
17711768 }
1772-
1769+ }
1770+
17731771 // V0 type selection
17741772 if (kshort.v0Type () != 1 )
17751773 return false ;
17761774
17771775 if (useMLScores) {
17781776 // if (kshort.k0ShortBDTScore() <= kshortSelections.KShort_MLThreshold)
1779- return false ;
1777+ return false ;
17801778
17811779 } else {
17821780 // KShort basic selection criteria:
@@ -1814,7 +1812,7 @@ struct sigma0builder {
18141812
18151813 histos.fill (HIST (" KShortSel/hSelectionStatistics" ), 9 .);
18161814
1817- if (kshort.qtarm () < kshortSelections.KShortArmenterosCoefficient * TMath::Abs (kshort.alpha ()))
1815+ if (kshort.qtarm () < kshortSelections.KShortArmenterosCoefficient * TMath::Abs (kshort.alpha ()))
18181816 return false ;
18191817
18201818 histos.fill (HIST (" KShortSel/hSelectionStatistics" ), 10 .);
@@ -2353,10 +2351,8 @@ struct sigma0builder {
23532351 histos.fill (HIST (" V0QA/h3dKShortMass" ), centrality, v0.pt (), v0.mK0Short ());
23542352 histos.fill (HIST (" V0QA/h3dYKShortMass" ), v0.yK0Short (), v0.pt (), v0.mK0Short ());
23552353 histos.fill (HIST (" V0QA/h3dYRKShortMass" ), v0.yK0Short (), v0.v0radius (), v0.mK0Short ());
2356-
23572354 }
23582355
2359-
23602356 // MC part:
23612357 if constexpr (requires { v0.motherMCPartId (); }) {
23622358 if (v0.has_v0MCCore ()) {
@@ -2394,7 +2390,7 @@ struct sigma0builder {
23942390 if (v0MC.pdgCodeMother () == o2::constants::physics::Pdg::kK0Star892 ) { // If from kstar
23952391 histos.fill (HIST (" V0QA/h2dTrueKStarKShortMass" ), v0MC.ptMC (), v0.mK0Short ());
23962392 }
2397- }
2393+ }
23982394 }
23992395 }
24002396 }
@@ -2455,4 +2451,5 @@ struct sigma0builder {
24552451WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
24562452{
24572453 return WorkflowSpec{adaptAnalysisTask<sigma0builder>(cfgc)};
2458- }
2454+ }
2455+
0 commit comments