|
11 | 11 |
|
12 | 12 | // author: Arvind Khuntia (arvind.khuntia@cern.ch) INFN Bologna, Italy |
13 | 13 |
|
| 14 | +#include <string> |
| 15 | +#include <vector> |
14 | 16 | #include <TLorentzVector.h> |
15 | 17 | #include <TVector2.h> |
16 | 18 |
|
|
31 | 33 | #include "CommonConstants/PhysicsConstants.h" |
32 | 34 | #include "ReconstructionDataFormats/Track.h" |
33 | 35 |
|
| 36 | +#include "PWGLF/DataModel/LFParticleIdentification.h" |
| 37 | + |
34 | 38 | #include "PWGJE/Core/FastJetUtilities.h" |
35 | 39 | #include "PWGJE/Core/JetDerivedDataUtilities.h" |
36 | 40 | #include "PWGJE/DataModel/Jet.h" |
@@ -92,6 +96,7 @@ struct nucleiInJets { |
92 | 96 | Configurable<bool> isMC{"isMC", false, "flag for the MC"}; |
93 | 97 | Configurable<bool> isWithJetEvents{"isWithJetEvents", true, "Events with at least one jet"}; |
94 | 98 | Configurable<bool> isWithLeadingJet{"isWithLeadingJet", true, "Events with leading jet"}; |
| 99 | + Configurable<bool> useLfTpcPid{"useLfTpcPid", true, "Events with custom TPC parameters"}; |
95 | 100 |
|
96 | 101 | Configurable<double> cfgtrkMinPt{"cfgtrkMinPt", 0.15, "set track min pT"}; |
97 | 102 | Configurable<double> cfgtrkMaxEta{"cfgtrkMaxEta", 0.8, "set track max Eta"}; |
@@ -149,7 +154,9 @@ struct nucleiInJets { |
149 | 154 | using TrackCandidates = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTPCFullKa, |
150 | 155 | aod::pidTPCFullPr, aod::pidTPCFullDe, aod::pidTPCFullHe, aod::pidTPCFullTr, aod::pidTOFFullPi, aod::pidTOFFullKa, |
151 | 156 | aod::pidTOFFullPr, aod::pidTOFFullDe, aod::pidTOFFullTr, aod::pidTOFFullHe, aod::pidTOFbeta, aod::TOFSignal>; |
152 | | - |
| 157 | + using TrackCandidatesLfPid = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTPCFullKa, |
| 158 | + aod::pidTPCLfFullPr, aod::pidTPCLfFullDe, aod::pidTPCLfFullHe, aod::pidTPCLfFullTr, aod::pidTOFFullPi, aod::pidTOFFullKa, |
| 159 | + aod::pidTOFFullPr, aod::pidTOFFullDe, aod::pidTOFFullTr, aod::pidTOFFullHe, aod::pidTOFbeta, aod::TOFSignal>; |
153 | 160 | using TrackCandidatesMC = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTPCFullKa, |
154 | 161 | aod::pidTPCFullPr, aod::pidTPCFullDe, aod::pidTPCFullHe, aod::pidTPCFullTr, aod::pidTOFFullPi, |
155 | 162 | aod::pidTOFFullKa, aod::pidTOFFullPr, aod::pidTOFFullDe, aod::pidTOFFullTr, aod::pidTOFFullHe, |
@@ -963,8 +970,13 @@ struct nucleiInJets { |
963 | 970 | jetHist.fill(HIST("jet/h1JetEvents"), 0.5); |
964 | 971 |
|
965 | 972 | for (auto& track : tracks) { |
966 | | - auto trk = track.track_as<TrackCandidates>(); |
967 | | - fillTrackInfo<false>(trk, chargedjets, leadingJetWithPtEtaPhi); |
| 973 | + if (useLfTpcPid) { |
| 974 | + auto trk = track.track_as<TrackCandidatesLfPid>(); |
| 975 | + fillTrackInfo<false>(trk, chargedjets, leadingJetWithPtEtaPhi); |
| 976 | + } else { |
| 977 | + auto trk = track.track_as<TrackCandidates>(); |
| 978 | + fillTrackInfo<false>(trk, chargedjets, leadingJetWithPtEtaPhi); |
| 979 | + } |
968 | 980 | } |
969 | 981 | } |
970 | 982 |
|
|
0 commit comments