@@ -1155,56 +1155,56 @@ struct AntinucleiInJets {
11551155 if (requireIsVertexTOFmatched && !collision.selection_bit (o2::aod::evsel::kIsVertexTOFmatched ))
11561156 return ;
11571157 registryData.fill (HIST (" number_of_events_data" ), 8.5 );
1158-
1158+
11591159 // Loop over all tracks for Full Event histograms
11601160 for (auto const & track : tracks) {
1161- // Apply standard track selection
1162- if (!passedTrackSelection (track))
1163- continue ;
1164-
1165- // Define variables
1166- double nsigmaTPCPr = track.tpcNSigmaPr ();
1167- double nsigmaTOFPr = track.tofNSigmaPr ();
1168- double nsigmaTPCDe = track.tpcNSigmaDe ();
1169- double nsigmaTOFDe = track.tofNSigmaDe ();
1170- double pt = track.pt ();
1161+ // Apply standard track selection
1162+ if (!passedTrackSelection (track))
1163+ continue ;
11711164
1172- // ITS PID logic
1173- bool passedItsPidProt (true ), passedItsPidDeut (true );
1174- double nSigmaITSprot = static_cast <double >(itsResponse.nSigmaITS <o2::track::PID::Proton>(track));
1175- double nSigmaITSdeut = static_cast <double >(itsResponse.nSigmaITS <o2::track::PID::Deuteron>(track));
1165+ // Define variables
1166+ double nsigmaTPCPr = track.tpcNSigmaPr ();
1167+ double nsigmaTOFPr = track.tofNSigmaPr ();
1168+ double nsigmaTPCDe = track.tpcNSigmaDe ();
1169+ double nsigmaTOFDe = track.tofNSigmaDe ();
1170+ double pt = track.pt ();
11761171
1177- if (applyItsPid && pt < ptMaxItsPidProt && (nSigmaITSprot < nSigmaItsMin || nSigmaITSprot > nSigmaItsMax)) {
1178- passedItsPidProt = false ;
1179- }
1180- if (applyItsPid && pt < ptMaxItsPidDeut && (nSigmaITSdeut < nSigmaItsMin || nSigmaITSdeut > nSigmaItsMax)) {
1181- passedItsPidDeut = false ;
1172+ // ITS PID logic
1173+ bool passedItsPidProt (true ), passedItsPidDeut (true );
1174+ double nSigmaITSprot = static_cast <double >(itsResponse.nSigmaITS <o2::track::PID::Proton>(track));
1175+ double nSigmaITSdeut = static_cast <double >(itsResponse.nSigmaITS <o2::track::PID::Deuteron>(track));
1176+
1177+ if (applyItsPid && pt < ptMaxItsPidProt && (nSigmaITSprot < nSigmaItsMin || nSigmaITSprot > nSigmaItsMax)) {
1178+ passedItsPidProt = false ;
1179+ }
1180+ if (applyItsPid && pt < ptMaxItsPidDeut && (nSigmaITSdeut < nSigmaItsMin || nSigmaITSdeut > nSigmaItsMax)) {
1181+ passedItsPidDeut = false ;
1182+ }
1183+
1184+ // Fill histograms for antimatter
1185+ if (track.sign () < 0 ) {
1186+
1187+ // Antiprotons Full Event
1188+ if (passedItsPidProt) {
1189+ registryData.fill (HIST (" antiproton_fullEvent_tpc" ), pt, nsigmaTPCPr);
1190+
1191+ // Require TOF matching and preliminary TPC cut
1192+ if (nsigmaTPCPr > minNsigmaTpc && nsigmaTPCPr < maxNsigmaTpc && track.hasTOF ()) {
1193+ registryData.fill (HIST (" antiproton_fullEvent_tof" ), pt, nsigmaTOFPr);
11821194 }
1195+ }
11831196
1184- // Fill histograms for antimatter
1185- if (track.sign () < 0 ) {
1186-
1187- // Antiprotons Full Event
1188- if (passedItsPidProt) {
1189- registryData.fill (HIST (" antiproton_fullEvent_tpc" ), pt, nsigmaTPCPr);
1190-
1191- // Require TOF matching and preliminary TPC cut
1192- if (nsigmaTPCPr > minNsigmaTpc && nsigmaTPCPr < maxNsigmaTpc && track.hasTOF ()) {
1193- registryData.fill (HIST (" antiproton_fullEvent_tof" ), pt, nsigmaTOFPr);
1194- }
1195- }
1196-
1197- // Antideuterons Full Event
1198- if (passedItsPidDeut) {
1199- registryData.fill (HIST (" antideuteron_fullEvent_tpc" ), pt, nsigmaTPCDe);
1200-
1201- // Require TOF matching and preliminary TPC cut
1202- if (nsigmaTPCDe > minNsigmaTpc && nsigmaTPCDe < maxNsigmaTpc && track.hasTOF ()) {
1203- registryData.fill (HIST (" antideuteron_fullEvent_tof" ), pt, nsigmaTOFDe);
1204- }
1205- }
1206- }
1207- }
1197+ // Antideuterons Full Event
1198+ if (passedItsPidDeut) {
1199+ registryData.fill (HIST (" antideuteron_fullEvent_tpc" ), pt, nsigmaTPCDe);
1200+
1201+ // Require TOF matching and preliminary TPC cut
1202+ if (nsigmaTPCDe > minNsigmaTpc && nsigmaTPCDe < maxNsigmaTpc && track.hasTOF ()) {
1203+ registryData.fill (HIST (" antideuteron_fullEvent_tof" ), pt, nsigmaTOFDe);
1204+ }
1205+ }
1206+ }
1207+ }
12081208
12091209 // Loop over reconstructed tracks
12101210 int id (-1 );
@@ -1236,7 +1236,7 @@ struct AntinucleiInJets {
12361236 fastjet::ClusterSequenceArea cs (fjParticles, jetDef, areaDef);
12371237 std::vector<fastjet::PseudoJet> jets = fastjet::sorted_by_pt (cs.inclusive_jets ());
12381238 auto [rhoPerp, rhoMPerp] = jetutilities::estimateRhoPerpCone (fjParticles, jets[0 ], rJet);
1239-
1239+
12401240 // Loop over reconstructed jets
12411241 bool isAtLeastOneJetSelected = false ;
12421242 for (const auto & jet : jets) {
0 commit comments