@@ -154,6 +154,7 @@ struct nuclei_in_jets {
154154 registryQC.add (" sumPtUE" , " sumPtUE" , HistType::kTH1F , {{500 , 0 , 50 , " #it{p}_{T} (GeV/#it{c})" }});
155155 registryQC.add (" nJets_found" , " nJets_found" , HistType::kTH1F , {{10 , 0 , 10 , " #it{n}_{Jet}" }});
156156 registryQC.add (" nJets_selected" , " nJets_selected" , HistType::kTH1F , {{10 , 0 , 10 , " #it{n}_{Jet}" }});
157+ registryQC.add (" event_selection_jets" , " event_selection_jets" , HistType::kTH1F , {{10 , 0 , 10 , " counter" }});
157158 registryQC.add (" dcaxy_vs_pt" , " dcaxy_vs_pt" , HistType::kTH2F , {{100 , 0.0 , 5.0 , " #it{p}_{T} (GeV/#it{c})" }, {2000 , -0.05 , 0.05 , " DCA_{xy} (cm)" }});
158159 registryQC.add (" dcaz_vs_pt" , " dcaz_vs_pt" , HistType::kTH2F , {{100 , 0.0 , 5.0 , " #it{p}_{T} (GeV/#it{c})" }, {2000 , -0.05 , 0.05 , " DCA_{z} (cm)" }});
159160 registryQC.add (" jet_ue_overlaps" , " jet_ue_overlaps" , HistType::kTH2F , {{20 , 0.0 , 20.0 , " #it{n}_{jet}" }, {200 , 0.0 , 200.0 , " #it{n}_{overlaps}" }});
@@ -440,6 +441,7 @@ struct nuclei_in_jets {
440441 {
441442 // Event Counter: before event selection
442443 registryData.fill (HIST (" number_of_events_data" ), 0.5 );
444+ registryQC.fill (HIST (" event_selection_jets" ), 0.5 ); // all events before jet selection
443445
444446 // Event Selection
445447 if (!collision.sel8 ())
@@ -591,8 +593,27 @@ struct nuclei_in_jets {
591593 if (n_jets_selected == 0 )
592594 return ;
593595 registryData.fill (HIST (" number_of_events_data" ), 3.5 );
596+ registryQC.fill (HIST (" event_selection_jets" ), 1.5 ); // events with pTjet>10 GeV/c selected
594597 // ************************************************************************************************************************************
595598
599+ // Leading Track
600+ double pt_max (0 );
601+
602+ // Loop over Reconstructed Tracks
603+ for (auto const & track : tracks) {
604+
605+ if (!passedTrackSelectionForJetReconstruction (track))
606+ continue ;
607+
608+ if (track.pt () > pt_max) {
609+ pt_max = track.pt ();
610+ }
611+ }
612+ // Event Counter: Skip Events with pt<pt_leading_min
613+ if (pt_max < 5.0 )
614+ return ;
615+ registryQC.fill (HIST (" event_selection_jets" ), 2.5 ); // events with pTleading > 5 GeV/c selected
616+
596617 // Overlaps
597618 int nOverlaps (0 );
598619 for (int i = 0 ; i < static_cast <int >(jet.size ()); i++) {
0 commit comments