@@ -218,6 +218,9 @@ struct correlateStrangeness {
218218 if (std::abs (track.dcaXY ()) > systCuts.dcaXYconstant + systCuts.dcaXYpTdep * std::abs (track.signed1Pt ())) {
219219 return false ;
220220 }
221+ if (track.pt () > axisRanges[3 ][1 ] || track.pt () < axisRanges[3 ][0 ]) {
222+ return false ;
223+ }
221224 return true ;
222225 }
223226 void fillCorrelationsV0 (aod::TriggerTracks const & triggers, aod::AssocV0s const & assocs, bool mixing, float pvz, float mult)
@@ -266,8 +269,7 @@ struct correlateStrangeness {
266269 continue ;
267270 if (ptassoc < axisRanges[2 ][0 ] || ptassoc > axisRanges[2 ][1 ])
268271 continue ;
269- if (pttrigger < axisRanges[3 ][0 ] || pttrigger > axisRanges[3 ][1 ])
270- continue ;
272+
271273 TH2F * hEfficiencyV0[3 ];
272274 hEfficiencyV0[0 ] = hEfficiencyK0Short;
273275 hEfficiencyV0[1 ] = hEfficiencyLambda;
@@ -359,8 +361,7 @@ struct correlateStrangeness {
359361 continue ;
360362 if (ptassoc < axisRanges[2 ][0 ] || ptassoc > axisRanges[2 ][1 ])
361363 continue ;
362- if (pttrigger < axisRanges[3 ][0 ] || pttrigger > axisRanges[3 ][1 ])
363- continue ;
364+
364365 TH2F * hEfficiencyCascade[4 ];
365366 hEfficiencyCascade[0 ] = hEfficiencyXiMinus;
366367 hEfficiencyCascade[1 ] = hEfficiencyXiPlus;
@@ -431,8 +432,6 @@ struct correlateStrangeness {
431432 continue ;
432433 if (ptassoc < axisRanges[2 ][0 ] || ptassoc > axisRanges[2 ][1 ])
433434 continue ;
434- if (pttrigger < axisRanges[3 ][0 ] || pttrigger > axisRanges[3 ][1 ])
435- continue ;
436435
437436 if (!mixing)
438437 histos.fill (HIST (" sameEvent/Pion" ), deltaphi, deltaeta, ptassoc, pttrigger, pvz, mult);
@@ -1274,6 +1273,9 @@ struct correlateStrangeness {
12741273 return ;
12751274 if (!bestCollisionINELgtZERO)
12761275 return ;
1276+ if (bestCollisionFT0Mpercentile > axisRanges[5 ][1 ] || bestCollisionFT0Mpercentile < axisRanges[5 ][0 ]) {
1277+ return ;
1278+ }
12771279 }
12781280
12791281 int iteratorNum = -1 ;
@@ -1337,9 +1339,10 @@ struct correlateStrangeness {
13371339
13381340 for (Int_t iTrigger = 0 ; iTrigger < triggerIndices.size (); iTrigger++) {
13391341 auto triggerParticle = mcParticles.iteratorAt (triggerIndices[iTrigger]);
1340- // if (!mcParticle) {
1341- // continue;
1342- // }
1342+ // check range of trigger particle
1343+ if (triggerParticle.pt () > axisRanges[3 ][1 ] || triggerParticle.pt () < axisRanges[3 ][0 ]) {
1344+ continue ;
1345+ }
13431346 Double_t getatrigger = triggerParticle.eta ();
13441347 Double_t gphitrigger = triggerParticle.phi ();
13451348 Double_t pttrigger = triggerParticle.pt ();
@@ -1353,6 +1356,14 @@ struct correlateStrangeness {
13531356 Double_t getaassoc = assocParticle.eta ();
13541357 Double_t gphiassoc = assocParticle.phi ();
13551358 Double_t ptassoc = assocParticle.pt ();
1359+
1360+ // skip if basic ranges not met
1361+ if (gphiassoc < axisRanges[0 ][0 ] || gphiassoc > axisRanges[0 ][1 ])
1362+ continue ;
1363+ if (getaassoc < axisRanges[1 ][0 ] || getaassoc > axisRanges[1 ][1 ])
1364+ continue ;
1365+ if (ptassoc < axisRanges[2 ][0 ] || ptassoc > axisRanges[2 ][1 ])
1366+ continue ;
13561367 histos.fill (HIST (" ClosureTest/sameEvent/" ) + HIST (particlenames[index]), ComputeDeltaPhi (gphitrigger, gphiassoc), getatrigger - getaassoc, ptassoc, pttrigger, collision.posZ (), bestCollisionFT0Mpercentile);
13571368 }
13581369 }
0 commit comments