@@ -640,8 +640,8 @@ struct HfTaskCharmPolarisation {
640640 std::vector<AxisSpec> hgenPromptAxes = {thnAxisPt, thnAxisNumPvContributors, thnAxisY, thnAxisCosThetaStarRandom, thnAxisDausAcc, thnAxisResoChannelLc, thnAxisCharge};
641641 std::vector<AxisSpec> hgenNonPromptAxes = {thnAxisPt, thnAxisNumPvContributors, thnAxisY, thnAxisCosThetaStarRandom, thnAxisPtB, thnAxisDausAcc, thnAxisResoChannelLc, thnAxisCharge};
642642 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
643- hgenPromptAxes.push_back (thnAxisCentrality);
644- hgenNonPromptAxes.push_back (thnAxisCentrality);
643+ hgenPromptAxes.push_back (thnAxisCentrality);
644+ hgenNonPromptAxes.push_back (thnAxisCentrality);
645645 }
646646 registry.add (" hGenPromptRandom" , " THn for polarisation studies with cosThStar w.r.t. random axis and BDT scores for generated prompt D*+ candidates" , HistType::kTHnSparseF , hgenPromptAxes);
647647 registry.add (" hGenNonPromptRandom" , " THn for polarisation studies with cosThStar w.r.t. random axis and BDT scores for generated non-prompt D*+ candidates" , HistType::kTHnSparseF , hgenNonPromptAxes);
@@ -1207,7 +1207,7 @@ struct HfTaskCharmPolarisation {
12071207 } else {
12081208 registry.fill (HIST (" hRecoPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, nMuons);
12091209 }
1210- } else {
1210+ } else {
12111211 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
12121212 registry.fill (HIST (" hPartRecoPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, nMuons, centrality);
12131213 } else {
@@ -1219,13 +1219,13 @@ struct HfTaskCharmPolarisation {
12191219 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
12201220 registry.fill (HIST (" hRecoPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons, centrality);
12211221 } else {
1222- registry.fill (HIST (" hRecoPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons);
1222+ registry.fill (HIST (" hRecoPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons);
12231223 }
12241224 } else {
12251225 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
12261226 registry.fill (HIST (" hPartRecoPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons, centrality);
12271227 } else {
1228- registry.fill (HIST (" hPartRecoPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons);
1228+ registry.fill (HIST (" hPartRecoPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons);
12291229 }
12301230 }
12311231 }
@@ -1237,27 +1237,27 @@ struct HfTaskCharmPolarisation {
12371237 if (activateTrackingSys) {
12381238 if (!isPartRecoDstar) {
12391239 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
1240- registry.fill (HIST (" hRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, nMuons, ptBhadMother, centrality);
1240+ registry.fill (HIST (" hRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, nMuons, ptBhadMother, centrality);
12411241 } else {
12421242 registry.fill (HIST (" hRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, nMuons, ptBhadMother);
12431243 }
12441244 } else {
12451245 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
1246- registry.fill (HIST (" hPartRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, nMuons, ptBhadMother, centrality);
1246+ registry.fill (HIST (" hPartRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, nMuons, ptBhadMother, centrality);
12471247 } else {
12481248 registry.fill (HIST (" hPartRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], absEtaMin, numItsClsMin, numTpcClsMin, nMuons, ptBhadMother);
12491249 }
12501250 }
12511251 } else {
12521252 if (!isPartRecoDstar) {
12531253 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
1254- registry.fill (HIST (" hRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons, ptBhadMother, centrality);
1254+ registry.fill (HIST (" hRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons, ptBhadMother, centrality);
12551255 } else {
12561256 registry.fill (HIST (" hRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons, ptBhadMother);
12571257 }
12581258 } else {
12591259 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
1260- registry.fill (HIST (" hPartRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons, ptBhadMother, centrality);
1260+ registry.fill (HIST (" hPartRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons, ptBhadMother, centrality);
12611261 } else {
12621262 registry.fill (HIST (" hPartRecoNonPromptRandom" ), invMassCharmHad, ptCharmHad, numPvContributors, std::abs (rapCharmHad), invMassD0, cosThetaStar, outputMl[0 ], /* outputMl[1],*/ outputMl[2 ], nMuons, ptBhadMother);
12631263 }
@@ -1442,13 +1442,13 @@ struct HfTaskCharmPolarisation {
14421442 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
14431443 registry.fill (HIST (" hGenPromptRandom" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, areDausInAcc, resoChannelLc, charge, centrality);
14441444 } else {
1445- registry.fill (HIST (" hGenPromptRandom" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, areDausInAcc, resoChannelLc, charge);
1445+ registry.fill (HIST (" hGenPromptRandom" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, areDausInAcc, resoChannelLc, charge);
14461446 }
14471447 } else {
14481448 if (doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
14491449 registry.fill (HIST (" hGenPartRecoPromptRandom" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, areDausInAcc, resoChannelLc, charge, centrality);
14501450 } else {
1451- registry.fill (HIST (" hGenPartRecoPromptRandom" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, areDausInAcc, resoChannelLc, charge);
1451+ registry.fill (HIST (" hGenPartRecoPromptRandom" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, areDausInAcc, resoChannelLc, charge);
14521452 }
14531453 }
14541454 } else { // non-prompt
@@ -1463,7 +1463,7 @@ struct HfTaskCharmPolarisation {
14631463 registry.fill (HIST (" hGenPartRecoNonPromptRandom" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, ptBhadMother, areDausInAcc, resoChannelLc, charge, centrality);
14641464 } else {
14651465 registry.fill (HIST (" hGenPartRecoNonPromptRandom" ), ptCharmHad, numPvContributors, std::abs (rapCharmHad), cosThetaStar, ptBhadMother, areDausInAcc, resoChannelLc, charge);
1466- }
1466+ }
14671467 }
14681468 }
14691469 }
0 commit comments