@@ -208,6 +208,8 @@ struct TrHeAnalysis {
208208 Configurable<float > cfgCutMinItsClusterSizeHe{" cfgCutMinItsClusterSizeHe" , 1 .f , " Minimum ITS Cluster Size for He" };
209209 Configurable<float > cfgCutMaxItsClusterSizeH3{" cfgCutMaxItsClusterSizeH3" , 4 .f , " Maximum ITS Cluster Size for Tr" };
210210 Configurable<float > cfgCutMinItsClusterSizeH3{" cfgCutMinItsClusterSizeH3" , 1 .f , " Minimum ITS Cluster Size for Tr" };
211+ Configurable<float > cfgCutMinTofMassH3{" cfgCutMinTofMassH3" , 2 .24f , " Minimum Tof mass H3" };
212+ Configurable<float > cfgCutMaxTofMassH3{" cfgCutMaxTofMassH3" , 3 .32f , " Maximum TOF mass H3" };
211213 // Set the kinematic and PID cuts for tracks
212214 struct : ConfigurableGroup {
213215 Configurable<float > pCut{" pCut" , 0 .3f , " Value of the p selection for spectra (default 0.3)" };
@@ -273,7 +275,7 @@ struct TrHeAnalysis {
273275 h->GetXaxis ()->SetBinLabel (6 , " Sel8 cut" );
274276 h->GetXaxis ()->SetBinLabel (7 , " Z-vert Cut" );
275277 histos.add <TH1>(" histogram/cuts" , " cuts" , HistType::kTH1D ,
276- {{12 , -0.5 , 11 .5 }});
278+ {{13 , -0.5 , 12 .5 }});
277279 auto hCuts = histos.get <TH1>(HIST (" histogram/cuts" ));
278280 hCuts->GetXaxis ()->SetBinLabel (1 , " total" );
279281 hCuts->GetXaxis ()->SetBinLabel (2 , " p cut" );
@@ -287,6 +289,7 @@ struct TrHeAnalysis {
287289 hCuts->GetXaxis ()->SetBinLabel (10 , " TPC refit" );
288290 hCuts->GetXaxis ()->SetBinLabel (11 , " ITS refit" );
289291 hCuts->GetXaxis ()->SetBinLabel (12 , " ITS cluster size" );
292+ hCuts->GetXaxis ()->SetBinLabel (13 , " TOF mass cut" );
290293 for (int i = 0 ; i < nParticles; i++) {
291294 particles.push_back (Particle (particleNames.at (i), particlePdgCodes.at (i),
292295 particleMasses.at (i), particleCharge.at (i),
@@ -382,6 +385,10 @@ struct TrHeAnalysis {
382385 histos.fill (HIST (" histogram/cuts" ), 12 );
383386 continue ;
384387 }
388+ if (track.mass () < cfgCutMinTofMassH3 || track.mass () > cfgCutMaxTofMassH3) {
389+ histos.fill (HIST (" histogram/cuts" ), 13 );
390+ continue ;
391+ }
385392 histos.fill (HIST (" histogram/H3/H3-TPCsignVsTPCmomentum" ),
386393 track.tpcInnerParam () / (1 .f * track.sign ()),
387394 track.tpcSignal ());
@@ -531,6 +538,10 @@ struct TrHeAnalysis {
531538 histos.fill (HIST (" histogram/cuts" ), 12 );
532539 continue ;
533540 }
541+ if (track.mass () < cfgCutMinTofMassH3 || track.mass () > cfgCutMaxTofMassH3) {
542+ histos.fill (HIST (" histogram/cuts" ), 13 );
543+ continue ;
544+ }
534545 histos.fill (HIST (" histogram/H3/H3-TPCsignVsTPCmomentum" ),
535546 track.tpcInnerParam () / (1 .f * track.sign ()),
536547 track.tpcSignal ());
0 commit comments