@@ -99,15 +99,15 @@ struct HfTaskOmegac0ToOmegapi {
9999 Preslice<Omegac0CandsMlKF> candOmegacKFMlPerCollision = aod::hf_cand_xic0_omegac0::collisionId;
100100 PresliceUnsorted<CollisionsWithMcLabels> colPerMcCollision = aod::mccollisionlabel::mcCollisionId;
101101
102- ConfigurableAxis thnConfigAxisPromptScore{" thnConfigAxisPromptScore" , {100 , 0 , 1 }, " Prompt score" };
103102 ConfigurableAxis thnConfigAxisMass{" thnConfigAxisMass" , {700 , 2.4 , 3.1 }, " Cand. inv. mass" };
104103 ConfigurableAxis thnConfigAxisPt{" thnConfigAxisPt" , {500 , 0 , 50 }, " Cand. pT" };
105104 ConfigurableAxis thnConfigAxisPtB{" thnConfigAxisPtB" , {500 , 0 , 50 }, " Cand. beauty mother pT" };
106105 ConfigurableAxis thnConfigAxisY{" thnConfigAxisY" , {20 , -1 , 1 }, " Cand. rapidity" };
107- ConfigurableAxis thnConfigAxisCent{" thnConfigAxisCent" , {100 , 0 , 100 }, " Centrality" };
108106 ConfigurableAxis thnConfigAxisOrigin{" thnConfigAxisOrigin" , {3 , -0.5 , 2.5 }, " Cand. origin" };
109107 ConfigurableAxis thnConfigAxisMatchFlag{" thnConfigAxisMatchFlag" , {15 , -7.5 , 7.5 }, " Cand. MC match flag" };
110- ConfigurableAxis thnConfigAxisNumPvContr{" thnConfigAxisNumPvContr" , {200 , -0.5 , 199.5 }, " PV contributors" };
108+ ConfigurableAxis thnConfigAxisNumPvContr{" thnConfigAxisNumPvContr" , {200 , -0.5 , 199.5 }, " Coll. num. PV contributors" };
109+ ConfigurableAxis thnConfigAxisCent{" thnConfigAxisCent" , {100 , 0 , 100 }, " Coll. centrality precentile" };
110+ ConfigurableAxis thnConfigAxisPromptScore{" thnConfigAxisPromptScore" , {100 , 0 , 1 }, " Prompt score" };
111111 HistogramRegistry registry{" registry" , {}};
112112
113113 void init (InitContext&)
@@ -119,23 +119,29 @@ struct HfTaskOmegac0ToOmegapi {
119119 LOGP (fatal, " One and only one process function should be enabled at a time." );
120120 }
121121
122- const AxisSpec thnAxisMass{thnConfigAxisMass, " inv . mass (#Omega#pi) (GeV/#it{c}^{2})" };
122+ const AxisSpec thnAxisMass{thnConfigAxisMass, " Inv . mass (#Omega#pi) (GeV/#it{c}^{2})" };
123123 const AxisSpec thnAxisPt{thnConfigAxisPt, " #it{p}_{T} (GeV/#it{c})" };
124124 const AxisSpec thnAxisPtB{thnConfigAxisPtB, " #it{p}_{T}^{B} (GeV/#it{c})" };
125125 const AxisSpec thnAxisY{thnConfigAxisY, " y" };
126126 const AxisSpec thnAxisOrigin{thnConfigAxisOrigin, " Origin" };
127127 const AxisSpec thnAxisMatchFlag{thnConfigAxisMatchFlag, " MC match flag" };
128- const AxisSpec thnAxisNumPvContr{thnConfigAxisNumPvContr, " PV contributors" };
128+ const AxisSpec thnAxisNumPvContr{thnConfigAxisNumPvContr, " Number of primary vtx. contributors" };
129+ const AxisSpec thnAxisCent{thnConfigAxisCent, " Centrality percentile" };
130+ const AxisSpec thnAxisCentMc{thnConfigAxisCent, " Centrality percentile (from gen. MC info)" };
129131 const AxisSpec thnAxisPromptScore{thnConfigAxisPromptScore, " BDT score prompt" };
130- const AxisSpec thnAxisCent{thnConfigAxisCent, " Centrality" };
131132
132133 std::vector<AxisSpec> axes = {thnAxisMass, thnAxisPt, thnAxisY};
133134 std::vector<AxisSpec> axesMcGen = {thnAxisPt, thnAxisPtB, thnAxisY, thnAxisOrigin};
134135
135- if (doprocessDataKFParticleFT0C || doprocessDataKFParticleMlFT0C || doprocessDataKFParticleFT0M || doprocessDataKFParticleMlFT0M || doprocessMcKFParticleFT0M || doprocessMcKFParticleMlFT0M ) {
136+ if (doprocessDataKFParticleFT0C || doprocessDataKFParticleMlFT0C || doprocessDataKFParticleFT0M || doprocessDataKFParticleMlFT0M) {
136137 axes.push_back (thnAxisCent);
137138 axes.push_back (thnConfigAxisNumPvContr);
138- axesMcGen.push_back (thnAxisCent);
139+ }
140+
141+ if (doprocessMcKFParticleFT0M || doprocessMcKFParticleMlFT0M) {
142+ axes.push_back (thnAxisCentMc);
143+ axes.push_back (thnConfigAxisNumPvContr);
144+ axesMcGen.push_back (thnAxisCentMc);
139145 axesMcGen.push_back (thnConfigAxisNumPvContr);
140146 }
141147
@@ -145,7 +151,7 @@ struct HfTaskOmegac0ToOmegapi {
145151
146152 if (doprocessMcKFParticleFT0M || doprocessMcKFParticleMlFT0M) {
147153 registry.add (" hMcGenWithRecoColl" , " Gen. #Omega_{c}^{0} from charm and beauty (associated to a reco collision)" , HistType::kTHnSparseD , axesMcGen);
148- registry.add (" hNumRecoCollPerMcColl" , " Number of reco collisions associated to a mc collision;Num. reco. coll. per Mc coll.;" , {HistType::kTH1D , {{10 , -1 .5 , 8 .5 }}});
154+ registry.add (" hNumRecoCollPerMcColl" , " Number of reco collisions associated to a mc collision;Num. reco. coll. per Mc coll.;" , {HistType::kTH1D , {{10 , -0 .5 , 9 .5 }}});
149155 registry.get <THnSparse>(HIST (" hMcGenWithRecoColl" ))->Sumw2 ();
150156 }
151157
0 commit comments