@@ -147,6 +147,7 @@ struct HfDataCreatorHiddenCharmReduced {
147147
148148 if (config.fillHistograms ) {
149149 const AxisSpec axisPt{360 , 0 ., 36 ., " #it{p}_{T}^{proton} (GeV/#it{c})" };
150+ const AxisSpec axisPtEtaC{100 , 0 ., 10 ., " #it{p}_{T}^{#eta_{c}} (GeV/#it{c})" };
150151 const AxisSpec axisEta{100 , -1 ., 1 ., " #eta" };
151152 const AxisSpec axisDca{400 , -2 ., 2 ., " DCA_{xy} to primary vertex (cm)" };
152153 const AxisSpec axisNSigmaTPC{100 , -5 ., 5 ., " n#sigma_{TPC}" };
@@ -169,10 +170,9 @@ struct HfDataCreatorHiddenCharmReduced {
169170 registry.add (" hTpcChi2NCl" , " Selected proton tracks;#it{p}_{T}^{track} (GeV/#it{c}); TpcChi2NCl" , {HistType::kTH2D , {axisPt, {10 , 0 ., 10 , " TpcChi2NCl" }}});
170171 registry.add (" hItsNCls" , " Selected proton tracks;#it{p}_{T}^{track} (GeV/#it{c});ItsNCls" , {HistType::kTH2D , {axisPt, {7 , 0 , 7 , " ItsNCls" }}});
171172 registry.add (" hItsChi2NCl" , " Selected proton tracks;#it{p}_{T}^{track} (GeV/#it{c}); ItsChi2NCl" , {HistType::kTH2D , {axisPt, {10 , 0 ., 10 , " ItsChi2NCl" }}});
172- registry.add (" hInvMass" , " Invariant mass of selected proton with all other tracks in the event;#it{p}_{T}^{proton} (GeV/#it{c});" , {HistType::kTH2D , {axisPt, axisMass}});
173- registry.add (" hInvMassSparse" , " Invariant mass of selected proton with all other tracks in the event;#it{p}_{T}^{proton} (GeV/#it{c});" , {HistType::kTH2D , {axisPt, axisMass}});
173+ registry.add (" hInvMass" , " Invariant mass of selected protons pairs;#it{p}_{T}^{#eta_{c}} (GeV/#it{c});" , {HistType::kTH2D , {axisPtEtaC, axisMass}});
174174 registry.add (" hDeDxTPCProton" , " Selected proton tracks;#it{p}_{T}^{track} (GeV/#it{c});TPC dE/dx (a.u.)" , {HistType::kTH2D , {axisPt, AxisSpec{100 , 0 ., 200 ., " TPC dE/dx (a.u.)" }}});
175- registry.add (" hMassPtCutVars" , " Charmonia candidates;#it{M} (p#overline{p}) (GeV/#it{c}^{2});#it{p}_{T}(cc) (GeV/#it{c}); n#sigma_{TPC}^{p1 }; n#sigma_{TOF}^{p1 }; dca_{xy}^{p1 }; n#sigma_{TPC}^{p2} ; n#sigma_{TOF}^{p2} ; dca_{xy}^{p2} " , {HistType::kTHnSparseF , {axisMass, axisPt, axisNSigmaTPC, axisNSigmaTOF, axisDca, axisNSigmaTPC, axisNSigmaTOF, axisDca}});
175+ registry.add (" hMassPtCutVars" , " Charmonia candidates;#it{M} (p#overline{p}) (GeV/#it{c}^{2});#it{p}_{T}(cc) (GeV/#it{c}); p_{T}^{p}; n#sigma_{TPC}^{p }; n#sigma_{TOF}^{p }; dca_{xy}^{p }; p_{T}^{#overline{p}}; n#sigma_{TPC}^{#overline{p}} ; n#sigma_{TOF}^{#overline{p}} ; dca_{xy}^{#overline{p}} " , {HistType::kTHnSparseF , {axisMass, axisPtEtaC, axisPt, axisNSigmaTPC, axisNSigmaTOF, axisDca, axisPt , axisNSigmaTPC, axisNSigmaTOF, axisDca}});
176176 }
177177
178178 // init HF event selection helper
@@ -340,6 +340,7 @@ struct HfDataCreatorHiddenCharmReduced {
340340 registry.fill (HIST (" hTpcChi2NCl" ), trk.pt (), trk.tpcChi2NCl ());
341341 registry.fill (HIST (" hItsNCls" ), trk.pt (), trk.itsNCls ());
342342 registry.fill (HIST (" hItsChi2NCl" ), trk.pt (), trk.itsChi2NCl ());
343+ registry.fill (HIST (" hDeDxTPCProton" ), trk.pt (), trk.tpcSignal ());
343344 }
344345 }
345346 // skip event if not at least a pair of tracks skimming selection
@@ -368,7 +369,11 @@ struct HfDataCreatorHiddenCharmReduced {
368369 float ptEtac = RecoDecay::pt (RecoDecay::sumOfVec (pVec1, pVec2));
369370 registry.fill (HIST (" hInvMass" ), ptEtac, invMass);
370371 if (config.fillSparses .value ) {
371- registry.fill (HIST (" hMassPtCutVars" ), invMass, ptEtac, t1.tpcNSigmaPr (), t1.tofNSigmaPr (), t1.dcaXY (), t2.tpcNSigmaPr (), t2.tofNSigmaPr (), t2.dcaXY ());
372+ if (t1.sign () > 0 ) {
373+ registry.fill (HIST (" hMassPtCutVars" ), invMass, ptEtac, t1.pt (), t1.tpcNSigmaPr (), t1.tofNSigmaPr (), t1.dcaXY (), t2.pt (), t2.tpcNSigmaPr (), t2.tofNSigmaPr (), t2.dcaXY ());
374+ } else {
375+ registry.fill (HIST (" hMassPtCutVars" ), invMass, ptEtac, t2.pt (), t2.tpcNSigmaPr (), t2.tofNSigmaPr (), t2.dcaXY (), t1.pt (), t1.tpcNSigmaPr (), t1.tofNSigmaPr (), t1.dcaXY ());
376+ }
372377 }
373378 }
374379 }
0 commit comments