2424// Framework and other headers after
2525#include " PWGJE/Core/FastJetUtilities.h"
2626#include " PWGJE/Core/JetDerivedDataUtilities.h"
27+ #include " PWGJE/Core/JetFinder.h"
28+ #include " PWGJE/Core/JetFindingUtilities.h"
29+ #include " PWGJE/Core/JetSubstructureUtilities.h"
2730#include " PWGJE/Core/JetUtilities.h"
2831#include " PWGJE/DataModel/EMCALClusters.h"
2932#include " PWGJE/DataModel/GammaJetAnalysisTree.h"
3033#include " PWGJE/DataModel/Jet.h"
31- #include " PWGJE/Core/JetFinder.h"
32- #include " PWGJE/Core/JetFindingUtilities.h"
33- #include " PWGJE/Core/JetSubstructureUtilities.h"
3434
3535#include " Common/Core/RecoDecay.h"
3636#include " Common/Core/TrackSelection.h"
@@ -71,14 +71,14 @@ struct GammaJetTreeProducer {
7171 // analysis tree
7272 // charged jets
7373 // photon candidates
74- Produces<aod::GjChargedJets> chargedJetsTable; // detector level jets
75- Produces<aod::GjEvents> eventsTable; // rec events
76- Produces<aod::GjGammas> gammasTable; // detector level clusters
77- Produces<aod::GjMCEvents> mcEventsTable; // mc collisions information
78- Produces<aod::GjMCParticles> mcParticlesTable; // gen level particles (photons and pi0)
79- Produces<aod::GjGammaMCInfos> gammaMCInfosTable; // detector level clusters MC information
80- Produces<aod::GjChJetMCInfos> chJetMCInfosTable; // detector level charged jets MC information
81- Produces<aod::GjMCJets> mcJetsTable; // gen level jets
74+ Produces<aod::GjChargedJets> chargedJetsTable; // detector level jets
75+ Produces<aod::GjEvents> eventsTable; // rec events
76+ Produces<aod::GjGammas> gammasTable; // detector level clusters
77+ Produces<aod::GjMCEvents> mcEventsTable; // mc collisions information
78+ Produces<aod::GjMCParticles> mcParticlesTable; // gen level particles (photons and pi0)
79+ Produces<aod::GjGammaMCInfos> gammaMCInfosTable; // detector level clusters MC information
80+ Produces<aod::GjChJetMCInfos> chJetMCInfosTable; // detector level charged jets MC information
81+ Produces<aod::GjMCJets> mcJetsTable; // gen level jets
8282 Produces<aod::GjJetSubstructures> jetSubstructuresTable; // jet substructure observables
8383
8484 HistogramRegistry mHistograms {" GammaJetTreeProducerHisto" };
@@ -932,7 +932,7 @@ struct GammaJetTreeProducer {
932932 if (collision.posZ () > mVertexCut ) {
933933 continue ;
934934 }
935- if (!jetderiveddatautilities::selectCollision (collision, eventSelectionBits,true ,true ,rctLabel)) {
935+ if (!jetderiveddatautilities::selectCollision (collision, eventSelectionBits, true , true , rctLabel)) {
936936 continue ;
937937 }
938938 if (!jetderiveddatautilities::selectTrigger (collision, triggerMaskBits)) {
@@ -1208,11 +1208,11 @@ struct GammaJetTreeProducer {
12081208
12091209 fastjet::ClusterSequenceArea clusterSeq (jetReclusterer.findJets (jetConstituents, jetReclustered));
12101210 jetReclustered = sorted_by_pt (jetReclustered);
1211-
1211+
12121212 if (jetReclustered.size () == 0 ) {
12131213 return ;
12141214 }
1215-
1215+
12161216 fastjet::PseudoJet daughterSubJet = jetReclustered[0 ];
12171217 fastjet::PseudoJet parentSubJet1;
12181218 fastjet::PseudoJet parentSubJet2;
@@ -1231,7 +1231,7 @@ struct GammaJetTreeProducer {
12311231 // Continue with harder parent
12321232 daughterSubJet = parentSubJet1;
12331233 }
1234-
1234+
12351235 // Fill one row per jet with all splittings stored as vectors
12361236 // Pass the jet's global index to associate this substructure entry with the jet
12371237 jetSubstructuresTable (jetGlobalIndex, energyMotherVec, ptLeadingVec, ptSubLeadingVec, thetaVec);
@@ -1248,15 +1248,15 @@ struct GammaJetTreeProducer {
12481248 int32_t storedColIndex = getStoredColIndex (collision);
12491249 if (storedColIndex == -1 )
12501250 return ;
1251-
1251+
12521252 // build kd tree for tracks (needed for perpendicular cone rho calculation)
12531253 buildKdTree (collision, tracks);
1254-
1254+
12551255 // loop over charged jets
12561256 for (const auto & jet : chargedJets) {
12571257 // Fill jet table and get the stored jet's global index
12581258 int64_t jetGlobalIndex = fillChargedJetTable (storedColIndex, jet, tracks);
1259-
1259+
12601260 // Fill substructure table if enabled and jet was stored
12611261 if (calculateJetSubstructure && jetGlobalIndex >= 0 ) {
12621262 fillSubstructureTable (jetGlobalIndex, jet, tracks);
@@ -1265,9 +1265,6 @@ struct GammaJetTreeProducer {
12651265 }
12661266 PROCESS_SWITCH (GammaJetTreeProducer, processChargedJetsData, " Process charged jets" , true );
12671267
1268-
1269-
1270-
12711268 Preslice<aod::JetParticles> ParticlesPerMCCollisions = aod::jmcparticle::mcCollisionId;
12721269 // / \brief Processes MC particles and fills MC particle table
12731270 // / \param collision The collision to process
@@ -1362,11 +1359,11 @@ struct GammaJetTreeProducer {
13621359 if (!collision.has_mcCollision ()) {
13631360 return ;
13641361 }
1365-
1362+
13661363 // build kd tree for mc particles (needed for perpendicular cone rho calculation)
13671364 auto particlesPerMcCollision = mcgenparticles.sliceBy (ParticlesPerMCCollisions, collision.mcCollisionId ());
13681365 buildKdTree (collision, particlesPerMcCollision);
1369-
1366+
13701367 int localIndex = 0 ;
13711368 auto pjetsPerMcCollision = chargedJets.sliceBy (PJetsPerMCCollisions, collision.mcCollisionId ());
13721369 for (const auto & pjet : pjetsPerMcCollision) {
@@ -1394,10 +1391,10 @@ struct GammaJetTreeProducer {
13941391 int32_t storedColIndex = getStoredColIndex (collision);
13951392 if (storedColIndex == -1 )
13961393 return ;
1397-
1394+
13981395 // build kd tree for tracks (needed for perpendicular cone rho calculation)
13991396 buildKdTree (collision, tracks);
1400-
1397+
14011398 // loop over charged jets
14021399 for (const auto & jet : chargedJets) {
14031400 // Fill jet table and get the stored jet's global index
0 commit comments