@@ -147,7 +147,7 @@ struct GammaJetTreeProducer {
147147 return true ;
148148 }
149149
150- double ch_iso_in_cone (const auto & cluster, aod:: JetTracks const & tracks, float radius = 0.4 )
150+ double ch_iso_in_cone (const auto & cluster, JetTracks const & tracks, float radius = 0.4 )
151151 {
152152 double iso = 0 ;
153153 for (auto track : tracks) {
@@ -162,7 +162,20 @@ struct GammaJetTreeProducer {
162162 }
163163 return iso;
164164 }
165- double ch_perp_cone_rho (const auto & cluster, aod::JetTracks const & tracks, float radius = 0.4 )
165+
166+ void runTrackQA (const auto & collision, JetTracks const & tracks)
167+ {
168+ for (auto track : tracks) {
169+ if (!isTrackSelected (track)) {
170+ continue ;
171+ }
172+ mHistograms .fill (HIST (" trackPt" ), track.pt ());
173+ mHistograms .fill (HIST (" trackPtEtaPhi" ), track.pt (), track.eta (), track.phi ());
174+ mHistograms .fill (HIST (" trackPtEtaOccupancy" ), track.pt (), track.eta (), collision.trackOccupancyInTimeRange ());
175+ }
176+ }
177+
178+ double ch_perp_cone_rho (const auto & object, JetTracks const & tracks, float radius = 0.4 )
166179 {
167180 double ptSumLeft = 0 ;
168181 double ptSumRight = 0 ;
@@ -197,7 +210,7 @@ struct GammaJetTreeProducer {
197210 // ---------------------
198211 // Processing functions
199212 // ---------------------
200- void processClearMaps (aod:: JetCollisions const &)
213+ void processClearMaps (JetCollisions const &)
201214 {
202215 collisionMapping.clear ();
203216 }
@@ -208,7 +221,7 @@ struct GammaJetTreeProducer {
208221 // an integer instead
209222 Filter clusterDefinitionSelection = (o2::aod::jcluster::definition == mClusterDefinition );
210223 // Process clusters
211- void processClusters (soa::Join<aod:: JetCollisions, aod::BkgChargedRhos, aod::JCollisionBCs>::iterator const & collision, selectedClusters const & clusters, aod:: JetTracks const & tracks)
224+ void processClusters (soa::Join<JetCollisions, aod::BkgChargedRhos, aod::JCollisionBCs>::iterator const & collision, selectedClusters const & clusters, JetTracks const & tracks)
212225 {
213226 if (!isEventAccepted (collision)) {
214227 return ;
@@ -266,7 +279,7 @@ struct GammaJetTreeProducer {
266279
267280 Filter jetCuts = aod::jet::pt > jetPtMin;
268281 // Process charged jets
269- void processChargedJets (soa::Join<aod:: JetCollisions, aod::BkgChargedRhos, aod::JCollisionBCs>::iterator const & collision, soa::Filtered<soa::Join<aod::ChargedJets, aod::ChargedJetConstituents>> const & chargedJets, aod:: JetTracks const &)
282+ void processChargedJets (soa::Join<JetCollisions, aod::BkgChargedRhos, aod::JCollisionBCs>::iterator const & collision, soa::Filtered<soa::Join<aod::ChargedJets, aod::ChargedJetConstituents>> const & chargedJets, JetTracks const & tracks )
270283 {
271284 // event selection
272285 if (!isEventAccepted (collision)) {
@@ -281,7 +294,7 @@ struct GammaJetTreeProducer {
281294 nconst = 0 ;
282295 leadingTrackPt = 0 ;
283296 // loop over constituents
284- for (auto & constituent : jet.template tracks_as <aod:: JetTracks>()) {
297+ for (auto & constituent : jet.template tracks_as <JetTracks>()) {
285298 mHistograms .fill (HIST (" chjetpt_vs_constpt" ), jet.pt (), constituent.pt ());
286299 nconst++;
287300 if (constituent.pt () > leadingTrackPt) {
0 commit comments