@@ -63,17 +63,17 @@ struct LfTreeCreatorNuclei {
6363
6464 hEvents.add (" eventSelection" , " eventSelection" , kTH1D , {{6 , -0.5 , 5.5 }});
6565 auto h = hEvents.get <TH1>(HIST (" eventSelection" ));
66- h->GetXaxis ()->SetBinLabel (1 , " total " );
66+ h->GetXaxis ()->SetBinLabel (1 , " z-vertex cut " );
6767 h->GetXaxis ()->SetBinLabel (2 , " sel8" );
6868 h->GetXaxis ()->SetBinLabel (3 , " TFborder" );
69- h->GetXaxis ()->SetBinLabel (4 , " z-vertex " );
70- h->GetXaxis ()->SetBinLabel (5 , " not empty " );
69+ h->GetXaxis ()->SetBinLabel (4 , " not empty " );
70+ h->GetXaxis ()->SetBinLabel (5 , " |z|<10 normalization " );
7171 h->GetXaxis ()->SetBinLabel (6 , " With a good track" );
7272 customTrackSelection = myTrackSelection ();
7373 }
7474
7575 // track
76- Configurable<float > yCut{" yCut" , 1 .f , " Rapidity cut" };
76+ // Configurable<float> yCut{"yCut", 1.f, "Rapidity cut"};
7777 Configurable<float > cfgCutDCAxy{" cfgCutDCAxy" , 2 .0f , " DCAxy range for tracks" };
7878 Configurable<float > cfgCutDCAz{" cfgCutDCAz" , 2 .0f , " DCAz range for tracks" };
7979 Configurable<float > cfgCutEta{" cfgCutEta" , 0 .8f , " Eta range for tracks" };
@@ -106,7 +106,7 @@ struct LfTreeCreatorNuclei {
106106 Configurable<float > minNCrossedRowsOverFindableClustersTPC{" minNCrossedRowsOverFindableClustersTPC" , 0 .8f , " Additional cut on the minimum value of the ratio between crossed rows and findable clusters in the TPC" };
107107
108108 Filter collisionFilter = (aod::collision::posZ < cfgHighCutVertex && aod::collision::posZ > cfgLowCutVertex);
109- // Filter trackFilter = (nabs(aod::track::eta) < cfgCutEta) && (requireGlobalTrackInFilter());
109+
110110 Filter etaFilter = (nabs(aod::track::eta) < cfgCutEta);
111111 Filter trackFilter = (trackSelType.value == 0 && requireGlobalTrackWoDCAInFilter()) ||
112112 (trackSelType.value == 1 && requireGlobalTrackInFilter()) ||
@@ -301,20 +301,26 @@ struct LfTreeCreatorNuclei {
301301 {
302302 for (const auto & collision : collisions) {
303303 hEvents.fill (HIST (" eventSelection" ), 0 );
304+
304305 if (useEvsel && !collision.sel8 ()) {
305306 continue ;
306307 }
307308 hEvents.fill (HIST (" eventSelection" ), 1 );
309+
308310 if (removeTFBorder && !collision.selection_bit (aod::evsel::kNoTimeFrameBorder ))
309311 continue ;
310312 hEvents.fill (HIST (" eventSelection" ), 2 );
311- if (collision.posZ () >= cfgHighCutVertex && collision.posZ () <= cfgLowCutVertex)
312- continue ;
313- hEvents.fill (HIST (" eventSelection" ), 3 );
313+
314314 const auto & tracksInCollision = tracks.sliceBy (perCollision, collision.globalIndex ());
315315 if (doSkim && tracksInCollision.size () == 0 )
316316 continue ;
317- hEvents.fill (HIST (" eventSelection" ), 4 );
317+ hEvents.fill (HIST (" eventSelection" ), 3 );
318+
319+ // Fill the norm. column with good events with |z| < 10 cm before skimming
320+ if (collision.posZ () < 10 && collision.posZ () > -10 ) {
321+ hEvents.fill (HIST (" eventSelection" ), 4 );
322+ }
323+
318324 if (doSkim && (trackSelType.value == 3 ) && !checkQuality<false >(collision, tracksInCollision))
319325 continue ;
320326 fillForOneEvent<false >(collision, tracksInCollision);
@@ -329,18 +335,25 @@ struct LfTreeCreatorNuclei {
329335 aod::BCs const &, aod::McCollisions const & mcCollisions, aod::McParticles const & mcParticles)
330336 {
331337 for (const auto & collision : collisions) {
338+ hEvents.fill (HIST (" eventSelection" ), 0 );
339+
332340 if (useEvsel && !collision.sel8 ()) {
333341 continue ;
334342 }
335343 hEvents.fill (HIST (" eventSelection" ), 1 );
344+
336345 if (removeTFBorder && !collision.selection_bit (aod::evsel::kNoTimeFrameBorder ))
337346 continue ;
338347 hEvents.fill (HIST (" eventSelection" ), 2 );
339- if (collision.posZ () >= cfgHighCutVertex && collision.posZ () <= cfgLowCutVertex)
340- continue ;
341- hEvents.fill (HIST (" eventSelection" ), 3 );
348+
349+ // Fill the norm. column with good events with |z| < 10 cm before skimming
350+ if (collision.posZ () < 10 && collision.posZ () > -10 ) {
351+ hEvents.fill (HIST (" eventSelection" ), 4 );
352+ }
353+
342354 const auto & tracksInCollision = tracks.sliceBy (perCollision, collision.globalIndex ());
343355 fillForOneEvent<true >(collision, tracksInCollision);
356+ hEvents.fill (HIST (" eventSelection" ), 5 );
344357 }
345358 }
346359
0 commit comments