@@ -288,19 +288,7 @@ struct HfTagSelTracks {
288288 HistogramRegistry registry{
289289 " registry" ,
290290 {{" hRejTracks" , " Tracks;;entries" , {HistType::kTH1F , {{15 , 0.5 , 15.5 }}}},
291- {" hPtNoCuts" , " all tracks;#it{p}_{T}^{track} (GeV/#it{c});entries" , {HistType::kTH1F , {{100 , 0 ., 10 .}}}},
292- // 2-prong histograms
293- {" hPtCuts2Prong" , " tracks selected for 2-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});entries" , {HistType::kTH1F , {{100 , 0 ., 10 .}}}},
294- {" hDCAToPrimXYVsPtCuts2Prong" , " tracks selected for 2-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});DCAxy to prim. vtx. (cm);entries" , {HistType::kTH2F , {{100 , 0 ., 10 .}, {400 , -2 ., 2 .}}}},
295- {" hEtaCuts2Prong" , " tracks selected for 2-prong vertexing;#it{#eta};entries" , {HistType::kTH1F , {{static_cast <int >(1.2 * etaMax2Prong * 100 ), -1.2 * etaMax2Prong, 1.2 * etaMax2Prong}}}},
296- // 3-prong histograms
297- {" hPtCuts3Prong" , " tracks selected for 3-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});entries" , {HistType::kTH1F , {{100 , 0 ., 10 .}}}},
298- {" hDCAToPrimXYVsPtCuts3Prong" , " tracks selected for 3-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});DCAxy to prim. vtx. (cm);entries" , {HistType::kTH2F , {{100 , 0 ., 10 .}, {400 , -2 ., 2 .}}}},
299- {" hEtaCuts3Prong" , " tracks selected for 3-prong vertexing;#it{#eta};entries" , {HistType::kTH1F , {{static_cast <int >(1.2 * etaMax3Prong * 100 ), -1.2 * etaMax3Prong, 1.2 * etaMax3Prong}}}},
300- // bachelor (for cascades) histograms
301- {" hPtCutsV0bachelor" , " tracks selected for 3-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});entries" , {HistType::kTH1F , {{100 , 0 ., 10 .}}}},
302- {" hDCAToPrimXYVsPtCutsV0bachelor" , " tracks selected for V0-bachelor vertexing;#it{p}_{T}^{track} (GeV/#it{c});DCAxy to prim. vtx. (cm);entries" , {HistType::kTH2F , {{100 , 0 ., 10 .}, {400 , -2 ., 2 .}}}},
303- {" hEtaCutsV0bachelor" , " tracks selected for 3-prong vertexing;#it{#eta};entries" , {HistType::kTH1F , {{static_cast <int >(1.2 * etaMaxBach * 100 ), -1.2 * etaMaxBach, 1.2 * etaMaxBach}}}}}};
291+ {" hPtNoCuts" , " all tracks;#it{p}_{T}^{track} (GeV/#it{c});entries" , {HistType::kTH1F , {{360 , 0 ., 36 .}}}}}};
304292
305293 static const int nCuts = 4 ;
306294
@@ -309,6 +297,19 @@ struct HfTagSelTracks {
309297
310298 void init (InitContext const &)
311299 {
300+ // 2-prong histograms
301+ registry.add (" hPtCuts2Prong" , " tracks selected for 2-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});entries" , {HistType::kTH1F , {{360 , 0 ., 36 .}}});
302+ registry.add (" hDCAToPrimXYVsPtCuts2Prong" , " tracks selected for 2-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});DCAxy to prim. vtx. (cm);entries" , {HistType::kTH2F , {{360 , 0 ., 36 .}, {400 , -2 ., 2 .}}});
303+ registry.add (" hEtaCuts2Prong" , " tracks selected for 2-prong vertexing;#it{#eta};entries" , {HistType::kTH1F , {{static_cast <int >(1.2 * etaMax2Prong * 100 ), -1.2 * etaMax2Prong, 1.2 * etaMax2Prong}}});
304+ // 3-prong histograms
305+ registry.add (" hPtCuts3Prong" , " tracks selected for 3-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});entries" , {HistType::kTH1F , {{360 , 0 ., 36 .}}});
306+ registry.add (" hDCAToPrimXYVsPtCuts3Prong" , " tracks selected for 3-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});DCAxy to prim. vtx. (cm);entries" , {HistType::kTH2F , {{360 , 0 ., 36 .}, {400 , -2 ., 2 .}}});
307+ registry.add (" hEtaCuts3Prong" , " tracks selected for 3-prong vertexing;#it{#eta};entries" , {HistType::kTH1F , {{static_cast <int >(1.2 * etaMax3Prong * 100 ), -1.2 * etaMax3Prong, 1.2 * etaMax3Prong}}});
308+ // bachelor (for cascades) histograms
309+ registry.add (" hPtCutsV0bachelor" , " tracks selected for 3-prong vertexing;#it{p}_{T}^{track} (GeV/#it{c});entries" , {HistType::kTH1F , {{360 , 0 ., 36 .}}});
310+ registry.add (" hDCAToPrimXYVsPtCutsV0bachelor" , " tracks selected for V0-bachelor vertexing;#it{p}_{T}^{track} (GeV/#it{c});DCAxy to prim. vtx. (cm);entries" , {HistType::kTH2F , {{360 , 0 ., 36 .}, {400 , -2 ., 2 .}}});
311+ registry.add (" hEtaCutsV0bachelor" , " tracks selected for 3-prong vertexing;#it{#eta};entries" , {HistType::kTH1F , {{static_cast <int >(1.2 * etaMaxBach * 100 ), -1.2 * etaMaxBach, 1.2 * etaMaxBach}}});
312+
312313 cutsSingleTrack = {cutsTrack2Prong, cutsTrack3Prong, cutsTrackBach};
313314 std::string cutNames[nCuts + 1 ] = {" selected" , " rej pT" , " rej eta" , " rej track quality" , " rej dca" };
314315 std::string candNames[CandidateType::NCandidateTypes] = {" 2-prong" , " 3-prong" , " bachelor" };
@@ -864,19 +865,19 @@ struct HfTrackIndexSkimsCreator {
864865 {" hVtx2ProngZ" , " 2-prong candidates;#it{z}_{sec. vtx.} (cm);entries" , {HistType::kTH1F , {{1000 , -20 ., 20 .}}}},
865866 {" hNCand2Prong" , " 2-prong candidates preselected;# of candidates;entries" , {HistType::kTH1F , {{2000 , 0 ., 200000 .}}}},
866867 {" hNCand2ProngVsNTracks" , " 2-prong candidates preselected;# of selected tracks;# of candidates;entries" , {HistType::kTH2F , {{2500 , 0 ., 25000 .}, {2000 , 0 ., 200000 .}}}},
867- {" hmassD0ToPiK " , " D^{0} candidates;inv. mass (#pi K) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
868- {" hmassJpsiToEE " , " J/#psi candidates;inv. mass (e^{#plus} e^{#minus}) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
869- {" hmassJpsiToMuMu " , " J/#psi candidates;inv. mass (#mu^{#plus} #mu^{#minus}) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
868+ {" hMassD0ToPiK " , " D^{0} candidates;inv. mass (#pi K) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
869+ {" hMassJpsiToEE " , " J/#psi candidates;inv. mass (e^{#plus} e^{#minus}) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
870+ {" hMassJpsiToMuMu " , " J/#psi candidates;inv. mass (#mu^{#plus} #mu^{#minus}) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
870871 // 3-prong histograms
871872 {" hVtx3ProngX" , " 3-prong candidates;#it{x}_{sec. vtx.} (cm);entries" , {HistType::kTH1F , {{1000 , -2 ., 2 .}}}},
872873 {" hVtx3ProngY" , " 3-prong candidates;#it{y}_{sec. vtx.} (cm);entries" , {HistType::kTH1F , {{1000 , -2 ., 2 .}}}},
873874 {" hVtx3ProngZ" , " 3-prong candidates;#it{z}_{sec. vtx.} (cm);entries" , {HistType::kTH1F , {{1000 , -20 ., 20 .}}}},
874875 {" hNCand3Prong" , " 3-prong candidates preselected;# of candidates;entries" , {HistType::kTH1F , {{5000 , 0 ., 500000 .}}}},
875876 {" hNCand3ProngVsNTracks" , " 3-prong candidates preselected;# of selected tracks;# of candidates;entries" , {HistType::kTH2F , {{2500 , 0 ., 25000 .}, {5000 , 0 ., 500000 .}}}},
876- {" hmassDPlusToPiKPi " , " D^{#plus} candidates;inv. mass (#pi K #pi) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
877- {" hmassLcToPKPi " , " #Lambda_{c} candidates;inv. mass (p K #pi) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
878- {" hmassDsToPiKK " , " D_{s} candidates;inv. mass (K K #pi) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
879- {" hmassXicToPKPi " , " #Xi_{c} candidates;inv. mass (p K #pi) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}}}};
877+ {" hMassDPlusToPiKPi " , " D^{#plus} candidates;inv. mass (#pi K #pi) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
878+ {" hMassLcToPKPi " , " #Lambda_{c} candidates;inv. mass (p K #pi) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
879+ {" hMassDsToPiKK " , " D_{s} candidates;inv. mass (K K #pi) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}},
880+ {" hMassXicToPKPi " , " #Xi_{c} candidates;inv. mass (p K #pi) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}}}};
880881
881882 static const int n2ProngDecays = hf_cand_prong2::DecayType::N2ProngDecays; // number of 2-prong hadron types
882883 static const int n3ProngDecays = hf_cand_prong3::DecayType::N3ProngDecays; // number of 3-prong hadron types
@@ -1626,20 +1627,20 @@ struct HfTrackIndexSkimsCreator {
16261627 auto mass2Prong = RecoDecay::m (arrMom, arrMass2Prong[iDecay2P][0 ]);
16271628 switch (iDecay2P) {
16281629 case hf_cand_prong2::DecayType::D0ToPiK:
1629- registry.fill (HIST (" hmassD0ToPiK " ), mass2Prong);
1630+ registry.fill (HIST (" hMassD0ToPiK " ), mass2Prong);
16301631 break ;
16311632 case hf_cand_prong2::DecayType::JpsiToEE:
1632- registry.fill (HIST (" hmassJpsiToEE " ), mass2Prong);
1633+ registry.fill (HIST (" hMassJpsiToEE " ), mass2Prong);
16331634 break ;
16341635 case hf_cand_prong2::DecayType::JpsiToMuMu:
1635- registry.fill (HIST (" hmassJpsiToMuMu " ), mass2Prong);
1636+ registry.fill (HIST (" hMassJpsiToMuMu " ), mass2Prong);
16361637 break ;
16371638 }
16381639 }
16391640 if (whichHypo2Prong[iDecay2P] >= 2 ) {
16401641 auto mass2Prong = RecoDecay::m (arrMom, arrMass2Prong[iDecay2P][1 ]);
16411642 if (iDecay2P == hf_cand_prong2::DecayType::D0ToPiK) {
1642- registry.fill (HIST (" hmassD0ToPiK " ), mass2Prong);
1643+ registry.fill (HIST (" hMassD0ToPiK " ), mass2Prong);
16431644 }
16441645 }
16451646 }
@@ -1828,30 +1829,30 @@ struct HfTrackIndexSkimsCreator {
18281829 auto mass3Prong = RecoDecay::m (arr3Mom, arrMass3Prong[iDecay3P][0 ]);
18291830 switch (iDecay3P) {
18301831 case hf_cand_prong3::DecayType::DPlusToPiKPi:
1831- registry.fill (HIST (" hmassDPlusToPiKPi " ), mass3Prong);
1832+ registry.fill (HIST (" hMassDPlusToPiKPi " ), mass3Prong);
18321833 break ;
18331834 case hf_cand_prong3::DecayType::DsToPiKK:
1834- registry.fill (HIST (" hmassDsToPiKK " ), mass3Prong);
1835+ registry.fill (HIST (" hMassDsToPiKK " ), mass3Prong);
18351836 break ;
18361837 case hf_cand_prong3::DecayType::LcToPKPi:
1837- registry.fill (HIST (" hmassLcToPKPi " ), mass3Prong);
1838+ registry.fill (HIST (" hMassLcToPKPi " ), mass3Prong);
18381839 break ;
18391840 case hf_cand_prong3::DecayType::XicToPKPi:
1840- registry.fill (HIST (" hmassXicToPKPi " ), mass3Prong);
1841+ registry.fill (HIST (" hMassXicToPKPi " ), mass3Prong);
18411842 break ;
18421843 }
18431844 }
18441845 if (whichHypo3Prong[iDecay3P] >= 2 ) {
18451846 auto mass3Prong = RecoDecay::m (arr3Mom, arrMass3Prong[iDecay3P][1 ]);
18461847 switch (iDecay3P) {
18471848 case hf_cand_prong3::DecayType::DsToPiKK:
1848- registry.fill (HIST (" hmassDsToPiKK " ), mass3Prong);
1849+ registry.fill (HIST (" hMassDsToPiKK " ), mass3Prong);
18491850 break ;
18501851 case hf_cand_prong3::DecayType::LcToPKPi:
1851- registry.fill (HIST (" hmassLcToPKPi " ), mass3Prong);
1852+ registry.fill (HIST (" hMassLcToPKPi " ), mass3Prong);
18521853 break ;
18531854 case hf_cand_prong3::DecayType::XicToPKPi:
1854- registry.fill (HIST (" hmassXicToPKPi " ), mass3Prong);
1855+ registry.fill (HIST (" hMassXicToPKPi " ), mass3Prong);
18551856 break ;
18561857 }
18571858 }
@@ -2031,30 +2032,30 @@ struct HfTrackIndexSkimsCreator {
20312032 auto mass3Prong = RecoDecay::m (arr3Mom, arrMass3Prong[iDecay3P][0 ]);
20322033 switch (iDecay3P) {
20332034 case hf_cand_prong3::DecayType::DPlusToPiKPi:
2034- registry.fill (HIST (" hmassDPlusToPiKPi " ), mass3Prong);
2035+ registry.fill (HIST (" hMassDPlusToPiKPi " ), mass3Prong);
20352036 break ;
20362037 case hf_cand_prong3::DecayType::DsToPiKK:
2037- registry.fill (HIST (" hmassDsToPiKK " ), mass3Prong);
2038+ registry.fill (HIST (" hMassDsToPiKK " ), mass3Prong);
20382039 break ;
20392040 case hf_cand_prong3::DecayType::LcToPKPi:
2040- registry.fill (HIST (" hmassLcToPKPi " ), mass3Prong);
2041+ registry.fill (HIST (" hMassLcToPKPi " ), mass3Prong);
20412042 break ;
20422043 case hf_cand_prong3::DecayType::XicToPKPi:
2043- registry.fill (HIST (" hmassXicToPKPi " ), mass3Prong);
2044+ registry.fill (HIST (" hMassXicToPKPi " ), mass3Prong);
20442045 break ;
20452046 }
20462047 }
20472048 if (whichHypo3Prong[iDecay3P] >= 2 ) {
20482049 auto mass3Prong = RecoDecay::m (arr3Mom, arrMass3Prong[iDecay3P][1 ]);
20492050 switch (iDecay3P) {
20502051 case hf_cand_prong3::DecayType::DsToPiKK:
2051- registry.fill (HIST (" hmassDsToPiKK " ), mass3Prong);
2052+ registry.fill (HIST (" hMassDsToPiKK " ), mass3Prong);
20522053 break ;
20532054 case hf_cand_prong3::DecayType::LcToPKPi:
2054- registry.fill (HIST (" hmassLcToPKPi " ), mass3Prong);
2055+ registry.fill (HIST (" hMassLcToPKPi " ), mass3Prong);
20552056 break ;
20562057 case hf_cand_prong3::DecayType::XicToPKPi:
2057- registry.fill (HIST (" hmassXicToPKPi " ), mass3Prong);
2058+ registry.fill (HIST (" hMassXicToPKPi " ), mass3Prong);
20582059 break ;
20592060 }
20602061 }
@@ -2148,7 +2149,7 @@ struct HfTrackIndexSkimsCreatorCascades {
21482149 {{" hVtx2ProngX" , " 2-prong candidates;#it{x}_{sec. vtx.} (cm);entries" , {HistType::kTH1F , {{1000 , -2 ., 2 .}}}},
21492150 {" hVtx2ProngY" , " 2-prong candidates;#it{y}_{sec. vtx.} (cm);entries" , {HistType::kTH1F , {{1000 , -2 ., 2 .}}}},
21502151 {" hVtx2ProngZ" , " 2-prong candidates;#it{z}_{sec. vtx.} (cm);entries" , {HistType::kTH1F , {{1000 , -20 ., 20 .}}}},
2151- {" hmass2 " , " 2-prong candidates;inv. mass (K0s p) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}}}};
2152+ {" hMass2 " , " 2-prong candidates;inv. mass (K0s p) (GeV/#it{c}^{2});entries" , {HistType::kTH1F , {{500 , 0 ., 5 .}}}}}};
21522153
21532154 // using MyTracks = aod::BigTracksMC;
21542155 // Partition<MyTracks> selectedTracks = aod::hf_seltrack::isSelProng >= 4;
@@ -2341,7 +2342,7 @@ struct HfTrackIndexSkimsCreatorCascades {
23412342 registry.fill (HIST (" hVtx2ProngX" ), posCasc[0 ]);
23422343 registry.fill (HIST (" hVtx2ProngY" ), posCasc[1 ]);
23432344 registry.fill (HIST (" hVtx2ProngZ" ), posCasc[2 ]);
2344- registry.fill (HIST (" hmass2 " ), mass2K0sP);
2345+ registry.fill (HIST (" hMass2 " ), mass2K0sP);
23452346 }
23462347
23472348 } // loop over V0s
0 commit comments