@@ -225,7 +225,6 @@ struct NonPromptCascadeTask {
225225 AxisSpec multAxisZoom = {7000 , 3000 , 10000 , " Multiplicity FT0M" };
226226 AxisSpec nTracksAxis = {100 , 0 ., 100 ., " NTracksGlobal" };
227227
228-
229228 void initCCDB (aod::BCsWithTimestamps::iterator const & bc)
230229 {
231230 if (mRunNumber == bc.runNumber ()) {
@@ -262,7 +261,7 @@ struct NonPromptCascadeTask {
262261
263262 std::vector<double > ptBinning = {0.4 , 0.8 , 1.2 , 1.6 , 2.0 , 2.4 , 2.8 , 3.2 , 3.6 , 4.0 , 4.4 , 4.8 , 5.2 , 5.6 , 6.0 };
264263 // AxisSpec ptAxis = {ptBinning, "#it{p}_{T} (GeV/#it{c})"};
265-
264+
266265 std::array<std::string, 7 > cutsNames{" # candidates" , " hasTOF" , " nClusTPC" , " nSigmaTPCbach" , " nSigmaTPCprotontrack" , " nSigmaTPCpiontrack" , " cosPA" };
267266 auto cutsOmega{std::get<std::shared_ptr<TH2>>(mRegistry .add (" h_PIDcutsOmega" , " ;;Invariant mass (GeV/#it{c}^{2})" , HistType::kTH2D , {{cutsNames.size (), -0.5 , -0.5 + cutsNames.size ()}, {125 , 1.650 , 1.700 }}))};
268267 auto cutsXi{std::get<std::shared_ptr<TH2>>(mRegistry .add (" h_PIDcutsXi" , " ;;Invariant mass (GeV/#it{c}^{2})" , HistType::kTH2D , {{6 , -0.5 , 5.5 }, {125 , 1.296 , 1.346 }}))};
@@ -275,7 +274,6 @@ struct NonPromptCascadeTask {
275274 cutsOmega->GetYaxis ()->SetBinLabel (iBin + 1 , cutsNames[iBin].c_str ());
276275 cutsXi->GetYaxis ()->SetBinLabel (iBin + 1 , cutsNames[iBin].c_str ());
277276 }
278-
279277 }
280278
281279 template <typename CollisionType, typename TrackType>
@@ -341,24 +339,23 @@ struct NonPromptCascadeTask {
341339 {
342340 // std::cout << "Filling mult histos" << std::endl;
343341 for (const auto & coll : collisions) {
344- if (!mRuns .count (mRunNumber )) {
342+ if (!mRuns .count (mRunNumber )) {
345343 std::string histName = " mult/hMultVsCent_run" + std::to_string (mRunNumber );
346- mHistPointer1 = std::get<std::shared_ptr<TH2>>(mRegistry .add (histName.c_str (), histName.c_str (), HistType::kTH2F , {centAxis, multAxis}));
344+ mHistPointer1 = std::get<std::shared_ptr<TH2>>(mRegistry .add (histName.c_str (), histName.c_str (), HistType::kTH2F , {centAxis, multAxis}));
347345 histName = " mult/hMultVsCentZoom_run" + std::to_string (mRunNumber );
348- mHistPointer2 = std::get<std::shared_ptr<TH2>>(mRegistry .add (histName.c_str (), histName.c_str (), HistType::kTH2F , {centAxisZoom, multAxisZoom}));
346+ mHistPointer2 = std::get<std::shared_ptr<TH2>>(mRegistry .add (histName.c_str (), histName.c_str (), HistType::kTH2F , {centAxisZoom, multAxisZoom}));
349347 histName = " mult/hNTracksVsCent_run" + std::to_string (mRunNumber );
350- mHistPointer3 = std::get<std::shared_ptr<TH2>>(mRegistry .add (histName.c_str (), histName.c_str (), HistType::kTH2F , {centAxis, nTracksAxis}));
348+ mHistPointer3 = std::get<std::shared_ptr<TH2>>(mRegistry .add (histName.c_str (), histName.c_str (), HistType::kTH2F , {centAxis, nTracksAxis}));
351349 histName = " mult/hNTracksVsCentZoom_run" + std::to_string (mRunNumber );
352- mHistPointer4 = std::get<std::shared_ptr<TH2>>(mRegistry .add (histName.c_str (), histName.c_str (), HistType::kTH2F , {centAxisZoom, nTracksAxis}));
350+ mHistPointer4 = std::get<std::shared_ptr<TH2>>(mRegistry .add (histName.c_str (), histName.c_str (), HistType::kTH2F , {centAxisZoom, nTracksAxis}));
353351 mRuns .insert (mRunNumber );
354352 }
355- // std::cout << coll.selection_bit(aod::evsel::kNoSameBunchPileup) << std::endl;
353+ // std::cout << coll.selection_bit(aod::evsel::kNoSameBunchPileup) << std::endl;
356354 mHistPointer1 ->Fill (coll.centFT0M (), coll.multFT0M ());
357355 mHistPointer2 ->Fill (coll.centFT0M (), coll.multFT0M ());
358356 mHistPointer3 ->Fill (coll.centFT0M (), coll.multNTracksGlobal ());
359357 mHistPointer4 ->Fill (coll.centFT0M (), coll.multNTracksGlobal ());
360358
361-
362359 mRegistry .fill (HIST (" hMultVsCent" ), coll.centFT0M (), coll.multFT0M ());
363360 mRegistry .fill (HIST (" hMultVsCentZoom" ), coll.centFT0M (), coll.multFT0M ());
364361 mRegistry .fill (HIST (" hNTracksVsCent" ), coll.centFT0M (), (float )coll.multNTracksGlobal ());
@@ -575,7 +572,7 @@ struct NonPromptCascadeTask {
575572 if (mToiMap .count (bc.globalBC ())) {
576573 toiMask = mToiMap [bc.globalBC ()];
577574 }
578- candidates.emplace_back (NPCascCandidate{mRunNumber , mcParticleID, trackedCascGlobalIndex, itsTrackGlobalIndex, candidate.collisionId (), matchingChi2, deltaPtITSCascade, deltaPtCascade, cascITSclsSize, hasReassociatedClusters, hasFakeReassociation, isGoodMatch, isGoodCascade, pdgCodeMom, itsTrackPDG, fromHF[0 ], fromHF[1 ],
575+ candidates.emplace_back (NPCascCandidate{mRunNumber , mcParticleID, trackedCascGlobalIndex, itsTrackGlobalIndex, candidate.collisionId (), matchingChi2, deltaPtITSCascade, deltaPtCascade, cascITSclsSize, hasReassociatedClusters, hasFakeReassociation, isGoodMatch, isGoodCascade, pdgCodeMom, itsTrackPDG, fromHF[0 ], fromHF[1 ],
579576 collision.numContrib (), cascPVContribs, collision.collisionTimeRes (), primaryVertex.getX (), primaryVertex.getY (), primaryVertex.getZ (),
580577 cascadeLvector.pt (), cascadeLvector.eta (), cascadeLvector.phi (),
581578 protonTrack.pt (), protonTrack.eta (), pionTrack.pt (), pionTrack.eta (), bachelor.pt (), bachelor.eta (),
0 commit comments