@@ -170,6 +170,7 @@ struct PidDiHadron {
170170 ConfigurableAxis axisNsigmaTPC{" axisNsigmaTPC" , {80 , -5 , 5 }, " nsigmaTPC axis" };
171171 ConfigurableAxis axisNsigmaTOF{" axisNsigmaTOF" , {80 , -5 , 5 }, " nsigmaTOF axis" };
172172 ConfigurableAxis axisNsigmaITS{" axisNsigmaITS" , {80 , -5 , 5 }, " nsigmaITS axis" };
173+ ConfigurableAxis axisTpcSignal{" axisTpcSignal" , {250 , 0 , 250 }, " dEdx axis for TPC" };
173174
174175 Configurable<LabeledArray<int >> cfgUseEventCuts{" cfgUseEventCuts" , {LongArrayInt[0 ], 15 , 1 , {" Filtered Events" , " Sel8" , " kNoTimeFrameBorder" , " kNoITSROFrameBorder" , " kNoSameBunchPileup" , " kIsGoodZvtxFT0vsPV" , " kNoCollInTimeRangeStandard" , " kIsGoodITSLayersAll" , " kNoCollInRofStandard" , " kNoHighMultCollInPrevRof" , " Occupancy" , " Multcorrelation" , " T0AV0ACut" , " kIsVertexITSTPC" , " kTVXinTRD" }, {" EvCuts" }}, " Labeled array (int) for various cuts on resonances" };
175176 Configurable<LabeledArray<float >> nSigmas{" nSigmas" , {LongArrayFloat[0 ], 6 , 3 , {" UpCut_pi" , " UpCut_ka" , " UpCut_pr" , " LowCut_pi" , " LowCut_ka" , " LowCut_pr" }, {" TPC" , " TOF" , " ITS" }}, " Labeled array for n-sigma values for TPC, TOF, ITS for pions, kaons, protons (positive and negative)" };
@@ -498,6 +499,9 @@ struct PidDiHadron {
498499 histos.add (" TofItsNsigma_before" , " " , {HistType::kTHnSparseD , {{axisNsigmaITS, axisNsigmaTOF, axisPt}}});
499500 histos.add (" TofItsNsigma_after" , " " , {HistType::kTHnSparseD , {{axisNsigmaITS, axisNsigmaTOF, axisPt}}});
500501 }
502+
503+ histos.add (" TpcdEdx_ptwise" , " " , {HistType::kTH2D , {{axisTpcSignal, axisPt}}});
504+ histos.add (" TpcdEdx_ptwise_afterCut" , " " , {HistType::kTH2D , {{axisTpcSignal, axisPt}}});
501505 }
502506 }
503507
@@ -835,12 +839,18 @@ struct PidDiHadron {
835839 histos.fill (HIST (" TofItsNsigma_before" ), itsResponse.nSigmaITS <o2::track::PID::Proton>(track1), track1.tofNSigmaPr (), track1.pt ());
836840 }
837841
842+ if (cfgGetNsigmaQA)
843+ histos.fill (HIST (" TpcdEdx_ptwise" ), track1.tpcSignal (), track1.pt ());
844+
838845 if (cfgPIDParticle && getNsigmaPID (track1) != cfgPIDParticle)
839846 continue ; // if PID is selected, check if the track has the right PID
840847
841848 if (cfgGetNsigmaQA)
842849 fillNsigmaAfterCut (track1, cfgPIDParticle);
843850
851+ if (cfgGetNsigmaQA)
852+ histos.fill (HIST (" TpcdEdx_ptwise_afterCut" ), track1.tpcSignal (), track1.pt ());
853+
844854 if (!getEfficiencyCorrection (triggerWeight, track1.eta (), track1.pt (), posZ))
845855 continue ;
846856 if (system == SameEvent) {
0 commit comments