1515
1616// / \author Niveditha Ram , IP2I
1717
18- #include < vector>
19- #include < unordered_map>
18+ #include " PWGUD/DataModel/UDTables.h"
2019
21- #include " Framework/runDataProcessing.h"
22- #include " Framework/O2DatabasePDGPlugin.h"
23- #include " Framework/AnalysisTask.h"
24- #include " Framework/AnalysisDataModel.h"
2520#include " CCDB/BasicCCDBManager.h"
26- #include " DataFormatsParameters/GRPLHCIFData.h"
2721#include " DataFormatsParameters/GRPECSObject.h"
28- #include " PWGUD/DataModel/UDTables.h"
22+ #include " DataFormatsParameters/GRPLHCIFData.h"
23+ #include " Framework/AnalysisDataModel.h"
24+ #include " Framework/AnalysisTask.h"
25+ #include " Framework/O2DatabasePDGPlugin.h"
26+ #include " Framework/runDataProcessing.h"
2927
3028#include " TLorentzVector.h"
31- #include " TSystem.h"
3229#include " TMath.h"
3330#include " TRandom3.h"
31+ #include " TSystem.h"
32+
33+ #include < unordered_map>
34+ #include < vector>
3435
35- // table for saving tree with info on data
36+ // table for saving tree with info on data
3637namespace dimu
3738{
38- // dimuon
39+ // dimuon
3940DECLARE_SOA_COLUMN (RunNumber, runNumber, int );
4041DECLARE_SOA_COLUMN (M, m, float );
4142DECLARE_SOA_COLUMN (Pt, pt, float );
4243DECLARE_SOA_COLUMN (Rap, rap, float );
4344DECLARE_SOA_COLUMN (Phi, phi, float );
44- }
45+ } // namespace dimu
4546
4647namespace o2 ::aod
4748{
4849DECLARE_SOA_TABLE (DiMu, " AOD" , " DIMU" ,
4950 dimu::RunNumber,
5051 dimu::M, dimu::Pt, dimu::Rap, dimu::Phi);
51- } // namespace o2::aod
52-
53-
54-
52+ } // namespace o2::aod
5553
5654using namespace o2 ;
5755using namespace o2 ::framework;
@@ -87,7 +85,7 @@ struct upcPolarisationJPsiIncorr {
8785 Produces<o2::aod::DiMu> dimuSel;
8886 // defining histograms using histogram registry: different histos for the different process functions
8987 HistogramRegistry registry{" registry" , {}, OutputObjHandlingPolicy::AnalysisObject, true , true };
90-
88+
9189 // CONFIGURABLES
9290 static constexpr double Pi = o2::constants::math::PI;
9391 // pT of muon pairs
@@ -111,7 +109,7 @@ struct upcPolarisationJPsiIncorr {
111109 Configurable<float > lowPhi{" lowPhi" , -Pi, " lower limit in phi histo" };
112110 Configurable<float > highPhi{" highPhi" , Pi, " upper limit in phi histo" };
113111
114- // Analysis cuts
112+ // Analysis cuts
115113 Configurable<float > maxJpsiMass{" maxJpsiMass" , 3.18 , " Maximum of the jpsi peak for peak cut" };
116114 Configurable<float > minJpsiMass{" minJpsiMass" , 3.0 , " Minimum of the jpsi peak for peak cut" };
117115
@@ -128,7 +126,7 @@ struct upcPolarisationJPsiIncorr {
128126 const AxisSpec axisEta{nBinsEta, lowEta, highEta, " #eta" };
129127 const AxisSpec axisRapidity{nBinsRapidity, lowRapidity, highRapidity, " Rapidity" };
130128 const AxisSpec axisPhi{nBinsPhi, lowPhi, highPhi, " #varphi" };
131-
129+
132130 // histos
133131 // data and reco MC
134132 registry.add (" hMass" , " Invariant mass of muon pairs;;#counts" , kTH1D , {axisMass});
@@ -138,7 +136,6 @@ struct upcPolarisationJPsiIncorr {
138136 registry.add (" hPhi" , " #varphi of muon pairs;;#counts" , kTH1D , {axisPhi});
139137 }
140138
141-
142139 // retrieve particle mass (GeV/c^2) from TDatabasePDG
143140 float particleMass (int pid)
144141 {
@@ -251,7 +248,7 @@ struct upcPolarisationJPsiIncorr {
251248 for (unsigned int i = 0 ; i < ampsV0A.size (); ++i) {
252249 if (std::abs (ampsRelBCsV0A[i]) <= 1 ) {
253250 if (ampsV0A[i] > kMaxAmpV0A )
254- return ;
251+ return ;
255252 }
256253 }
257254
@@ -308,20 +305,16 @@ struct upcPolarisationJPsiIncorr {
308305 if (p.Rapidity () > highRapidity)
309306 return ;
310307
311-
312308 // fill the histos without looking at neutron emission
313309 registry.fill (HIST (" hMass" ), p.M ());
314310 registry.fill (HIST (" hPt" ), p.Pt ());
315311 registry.fill (HIST (" hEta" ), p.Eta ());
316312 registry.fill (HIST (" hRapidity" ), p.Rapidity ());
317313 registry.fill (HIST (" hPhi" ), p.Phi ());
318314
319- dimuSel (cand.runNumber (),p.M (), p.Pt (), p.Rapidity (), p.Phi ());
320-
321-
315+ dimuSel (cand.runNumber (), p.M (), p.Pt (), p.Rapidity (), p.Phi ());
322316 }
323317
324-
325318 // PROCESS FUNCTION
326319 void processData (CandidatesFwd const & eventCandidates,
327320 o2::aod::UDZdcsReduced const & ZDCs,
@@ -330,7 +323,7 @@ struct upcPolarisationJPsiIncorr {
330323
331324 // map with the tracks
332325 std::unordered_map<int32_t , std::vector<int32_t >> tracksPerCand;
333- // takes a tracks table with a coloumn of collision ID and makes it into a map of collision ID to each track.
326+ // takes a tracks table with a coloumn of collision ID and makes it into a map of collision ID to each track.
334327 collectCandIDs (tracksPerCand, fwdTracks);
335328
336329 // map with the ZDC info
@@ -361,7 +354,6 @@ struct upcPolarisationJPsiIncorr {
361354 }
362355
363356 PROCESS_SWITCH (upcPolarisationJPsiIncorr, processData, " " , true );
364-
365357};
366358
367359WorkflowSpec defineDataProcessing (ConfigContext const & cfgc)
0 commit comments