1515// /
1616// / \author Adrian Fereydon Nassirpour <adrian.fereydon.nassirpour@cern.ch>
1717
18- #include < CCDB/BasicCCDBManager.h>
19-
20- #include < iostream>
21- #include < vector>
22- #include < string>
23- #include < optional>
24-
25- #include < TLorentzVector.h>
26- #include < TVector2.h>
27-
28- #include " Framework/ASoA.h"
29- #include " Framework/AnalysisDataModel.h"
30- #include " Framework/AnalysisTask.h"
31- #include " Framework/HistogramRegistry.h"
32- #include " Framework/runDataProcessing.h"
33- #include " ReconstructionDataFormats/Track.h"
18+ #include " PWGJE/Core/FastJetUtilities.h"
19+ #include " PWGJE/Core/JetDerivedDataUtilities.h"
20+ #include " PWGJE/DataModel/EMCALClusters.h"
21+ #include " PWGJE/DataModel/Jet.h"
3422
3523#include " Common/Core/RecoDecay.h"
3624#include " Common/Core/TrackSelection.h"
3725#include " Common/Core/TrackSelectionDefaults.h"
3826#include " Common/Core/trackUtilities.h"
3927#include " Common/DataModel/EventSelection.h"
40- #include " Common/DataModel/TrackSelectionTables.h"
4128#include " Common/DataModel/Multiplicity.h"
4229#include " Common/DataModel/PIDResponse.h"
43- #include " CommonConstants/PhysicsConstants.h"
44-
45- #include " PWGJE/Core/FastJetUtilities.h"
46- #include " PWGJE/Core/JetDerivedDataUtilities.h"
47- #include " PWGJE/DataModel/Jet.h"
48- #include " PWGJE/DataModel/EMCALClusters.h"
49- #include " EMCALBase/Geometry.h"
50- #include " EMCALCalib/BadChannelMap.h"
30+ #include " Common/DataModel/TrackSelectionTables.h"
5131
32+ #include " CommonConstants/PhysicsConstants.h"
33+ #include " CommonDataFormat/InteractionRecord.h"
34+ #include " DataFormatsEMCAL/AnalysisCluster.h"
5235#include " DataFormatsEMCAL/Cell.h"
5336#include " DataFormatsEMCAL/Constants.h"
54- #include " DataFormatsEMCAL/AnalysisCluster.h"
55- #include " DataFormatsParameters/GRPObject.h"
5637#include " DataFormatsParameters/GRPMagField.h"
57-
38+ # include " DataFormatsParameters/GRPObject.h "
5839#include " DetectorsBase/Propagator.h"
40+ #include " EMCALBase/Geometry.h"
41+ #include " EMCALCalib/BadChannelMap.h"
42+ #include " Framework/ASoA.h"
43+ #include " Framework/AnalysisDataModel.h"
44+ #include " Framework/AnalysisTask.h"
45+ #include " Framework/HistogramRegistry.h"
46+ #include " Framework/runDataProcessing.h"
47+ #include " ReconstructionDataFormats/Track.h"
48+ #include < CCDB/BasicCCDBManager.h>
5949
60- #include " CommonDataFormat/InteractionRecord.h"
50+ #include < TLorentzVector.h>
51+ #include < TVector2.h>
52+
53+ #include < iostream>
54+ #include < optional>
55+ #include < string>
56+ #include < vector>
6157
6258using namespace o2 ;
6359using namespace o2 ::framework;
@@ -106,7 +102,7 @@ struct statPromptPhoton {
106102 Configurable<bool > cfgSkimmedTrigger{" cfgSkimmedTrigger" , false , " Enables trigger for skimmied datasets (2023 onwards)" };
107103 Configurable<std::string> cfgTriggerMasks{" cfgTriggerMasks" , " " , " possible JE Trigger masks: fJetChLowPt,fJetChHighPt,fTrackLowPt,fTrackHighPt,fJetD0ChLowPt,fJetD0ChHighPt,fJetLcChLowPt,fJetLcChHighPt,fEMCALReadout,fJetFullHighPt,fJetFullLowPt,fJetNeutralHighPt,fJetNeutralLowPt,fGammaVeryHighPtEMCAL,fGammaVeryHighPtDCAL,fGammaHighPtEMCAL,fGammaHighPtDCAL,fGammaLowPtEMCAL,fGammaLowPtDCAL,fGammaVeryLowPtEMCAL,fGammaVeryLowPtDCAL" };
108104 Configurable<bool > cfgDebug{" cfgDebug" , false , " Enables debug information for local running" };
109-
105+
110106 int trackFilter = -1 ;
111107 std::vector<int > triggerMaskBits;
112108
@@ -259,10 +255,10 @@ struct statPromptPhoton {
259255 using jMCClusters = o2::soa::Join<o2::aod::JMcClusterLbs, o2::aod::JClusters, o2::aod::JClusterTracks>;
260256 using jClusters = o2::soa::Join<o2::aod::JClusters, o2::aod::JClusterTracks>;
261257 using jselectedCollisions = soa::Join<aod::JCollisions, aod::JCollisionBCs, aod::JCollisionPIs, aod::EvSels, aod::JEMCCollisionLbs, aod::JMcCollisionLbs>;
262- using jselectedDataCollisions = soa::Join<aod::JCollisions, aod::JCollisionBCs, aod::JCollisionPIs, aod::EvSels, aod::JEMCCollisionLbs>;
263- // using jselectedDataCollisions = soa::Join<aod::JCollisions, aod::JCollisionBCs, aod::JCollisionPIs, aod::JCollisionMcInfos, aod::EvSels, aod::JEMCCollisionLbs>;
258+ using jselectedDataCollisions = soa::Join<aod::JCollisions, aod::JCollisionBCs, aod::JCollisionPIs, aod::EvSels, aod::JEMCCollisionLbs>;
259+ // using jselectedDataCollisions = soa::Join<aod::JCollisions, aod::JCollisionBCs, aod::JCollisionPIs, aod::JCollisionMcInfos, aod::EvSels, aod::JEMCCollisionLbs>;
264260 using jfilteredCollisions = soa::Filtered<jselectedCollisions>;
265- using jfilteredDataCollisions = soa::Filtered<jselectedDataCollisions>;
261+ using jfilteredDataCollisions = soa::Filtered<jselectedDataCollisions>;
266262 using jfilteredMCClusters = soa::Filtered<jMCClusters>;
267263 using jfilteredClusters = soa::Filtered<jClusters>;
268264
@@ -639,20 +635,19 @@ struct statPromptPhoton {
639635 }
640636 histos.fill (HIST (" REC_nEvents" ), 2.5 );
641637
642-
643- if (cfgSkimmedTrigger){
638+ if (cfgSkimmedTrigger) {
644639 if (!jetderiveddatautilities::selectTrigger (collision, triggerMaskBits)) {
645- return ;
640+ return ;
646641 }
647- }// JE Software Triggers
648-
642+ } // JE Software Triggers
643+
649644 histos.fill (HIST (" REC_nEvents" ), 3.5 );
650-
645+
651646 double weight = 1 ;
652647 if (collision.has_mcCollision ()) {
653648 weight = collision.mcCollision ().weight ();
654649 }
655-
650+
656651 bool noTrk = true ;
657652 for (auto & track : tracks) {
658653 if (cfgJETracks) {
@@ -1071,12 +1066,12 @@ struct statPromptPhoton {
10711066 }
10721067 }
10731068 double pthadsum = GetPtHadSum (tracks, track, cfgMinR, cfgMaxR, true , false , true );
1074- histos.fill (HIST (" REC_Trigger_V_PtHadSum_Nch" ), sternPt, pthadsum,weight);
1069+ histos.fill (HIST (" REC_Trigger_V_PtHadSum_Nch" ), sternPt, pthadsum, weight);
10751070 if (sterntrigger) {
10761071 bool doStern = true ;
10771072 double sterncount = 1.0 ;
10781073 while (doStern) {
1079- histos.fill (HIST (" REC_Trigger_V_PtHadSum_Stern" ), sterncount, pthadsum, (2.0 / sternPt)* weight);
1074+ histos.fill (HIST (" REC_Trigger_V_PtHadSum_Stern" ), sterncount, pthadsum, (2.0 / sternPt) * weight);
10801075 if (sterncount < sternPt) {
10811076 sterncount++;
10821077 } else {
@@ -1100,9 +1095,9 @@ struct statPromptPhoton {
11001095 }
11011096 if ((nEventsData + 1 ) % 10000 == 0 ) {
11021097 std::cout << " Processed Data Events: " << nEventsData << std::endl;
1103- std::cout << " Events Trigger Bit: " << collision.triggerSel () << std::endl;
1104- std::cout << " Trigger Mask Bit: " << triggerMaskBits[0 ] << std::endl;
1105- std::cout << " Trigger Mask Cfg Line: " << cfgTriggerMasks << std::endl;
1098+ std::cout << " Events Trigger Bit: " << collision.triggerSel () << std::endl;
1099+ std::cout << " Trigger Mask Bit: " << triggerMaskBits[0 ] << std::endl;
1100+ std::cout << " Trigger Mask Cfg Line: " << cfgTriggerMasks << std::endl;
11061101 }
11071102 }
11081103
@@ -1122,14 +1117,14 @@ struct statPromptPhoton {
11221117
11231118 histos.fill (HIST (" DATA_nEvents" ), 2.5 );
11241119
1125- if (cfgSkimmedTrigger){
1120+ if (cfgSkimmedTrigger) {
11261121 if (!jetderiveddatautilities::selectTrigger (collision, triggerMaskBits)) {
1127- return ;
1122+ return ;
11281123 }
1129- }// JE Software Triggers
1130-
1124+ } // JE Software Triggers
1125+
11311126 histos.fill (HIST (" DATA_nEvents" ), 3.5 );
1132-
1127+
11331128 bool noTrk = true ;
11341129 for (auto & track : tracks) {
11351130
0 commit comments