1818#include " PWGUD/DataModel/UDTables.h"
1919
2020#include " Common/Core/RecoDecay.h"
21+
2122#include " CCDB/BasicCCDBManager.h"
2223#include " CommonConstants/PhysicsConstants.h"
2324#include " DataFormatsParameters/GRPECSObject.h"
@@ -40,43 +41,42 @@ using namespace ROOT::Math;
4041namespace dimu
4142{
4243// dimuon
43- DECLARE_SOA_COLUMN (RunNumber, runNumber, int );
44- DECLARE_SOA_COLUMN (M, m, float );
45- DECLARE_SOA_COLUMN (Energy, energy, float );
46- DECLARE_SOA_COLUMN (Px, px, float );
47- DECLARE_SOA_COLUMN (Py, py, float );
48- DECLARE_SOA_COLUMN (Pz, pz, float );
49- DECLARE_SOA_COLUMN (Pt, pt, float );
50- DECLARE_SOA_COLUMN (Rap, rap, float );
51- DECLARE_SOA_COLUMN (Phi, phi, float );
52-
53- // tracks positive (p) and negative (n) stored
54- DECLARE_SOA_COLUMN (EnergyP, energyP, float );
55- DECLARE_SOA_COLUMN (Pxp, pxp, float );
56- DECLARE_SOA_COLUMN (Pyp, pyp, float );
57- DECLARE_SOA_COLUMN (Pzp, pzp, float );
58- DECLARE_SOA_COLUMN (Ptp, ptp, float );
59- DECLARE_SOA_COLUMN (Etap, etap, float );
60- DECLARE_SOA_COLUMN (Phip, phip, float );
61-
62- DECLARE_SOA_COLUMN (EnergyN, energyN, float );
63- DECLARE_SOA_COLUMN (Pxn, pxn, float );
64- DECLARE_SOA_COLUMN (Pyn, pyn, float );
65- DECLARE_SOA_COLUMN (Pzn, pzn, float );
66- DECLARE_SOA_COLUMN (Ptn, ptn, float );
67- DECLARE_SOA_COLUMN (Etan, etan, float );
68- DECLARE_SOA_COLUMN (Phin, phin, float );
69-
44+ DECLARE_SOA_COLUMN (RunNumber, runNumber, int );
45+ DECLARE_SOA_COLUMN (M, m, float );
46+ DECLARE_SOA_COLUMN (Energy, energy, float );
47+ DECLARE_SOA_COLUMN (Px, px, float );
48+ DECLARE_SOA_COLUMN (Py, py, float );
49+ DECLARE_SOA_COLUMN (Pz, pz, float );
50+ DECLARE_SOA_COLUMN (Pt, pt, float );
51+ DECLARE_SOA_COLUMN (Rap, rap, float );
52+ DECLARE_SOA_COLUMN (Phi, phi, float );
53+
54+ // tracks positive (p) and negative (n) stored
55+ DECLARE_SOA_COLUMN (EnergyP, energyP, float );
56+ DECLARE_SOA_COLUMN (Pxp, pxp, float );
57+ DECLARE_SOA_COLUMN (Pyp, pyp, float );
58+ DECLARE_SOA_COLUMN (Pzp, pzp, float );
59+ DECLARE_SOA_COLUMN (Ptp, ptp, float );
60+ DECLARE_SOA_COLUMN (Etap, etap, float );
61+ DECLARE_SOA_COLUMN (Phip, phip, float );
62+
63+ DECLARE_SOA_COLUMN (EnergyN, energyN, float );
64+ DECLARE_SOA_COLUMN (Pxn, pxn, float );
65+ DECLARE_SOA_COLUMN (Pyn, pyn, float );
66+ DECLARE_SOA_COLUMN (Pzn, pzn, float );
67+ DECLARE_SOA_COLUMN (Ptn, ptn, float );
68+ DECLARE_SOA_COLUMN (Etan, etan, float );
69+ DECLARE_SOA_COLUMN (Phin, phin, float );
70+
7071} // namespace dimu
7172
7273namespace o2 ::aod
7374{
74- DECLARE_SOA_TABLE (DiMu, " AOD" , " DIMU" ,
75+ DECLARE_SOA_TABLE (DiMu, " AOD" , " DIMU" ,
7576 dimu::RunNumber,
7677 dimu::M, dimu::Energy, dimu::Px, dimu::Py, dimu::Pz, dimu::Pt, dimu::Rap, dimu::Phi,
77- dimu::EnergyP, dimu::Pxp, dimu::Pyp, dimu::Pzp, dimu::Ptp, dimu::Etap, dimu::Phip,
78- dimu::EnergyN, dimu::Pxn, dimu::Pyn, dimu::Pzn, dimu::Ptn, dimu::Etan, dimu::Phin);
79-
78+ dimu::EnergyP, dimu::Pxp, dimu::Pyp, dimu::Pzp, dimu::Ptp, dimu::Etap, dimu::Phip,
79+ dimu::EnergyN, dimu::Pxn, dimu::Pyn, dimu::Pzn, dimu::Ptn, dimu::Etan, dimu::Phin);
8080
8181} // namespace o2::aod
8282using namespace o2 ;
@@ -97,8 +97,6 @@ const int kMaxChi2MFTMatch = 30;
9797const float kMaxZDCTime = 2 .;
9898const float kMaxZDCTimeHisto = 10 .;
9999
100-
101-
102100struct UpcPolarisationJpsiIncoh {
103101
104102 using CandidatesFwd = soa::Join<o2::aod::UDCollisions, o2::aod::UDCollisionsSelsFwd>;
@@ -130,19 +128,19 @@ struct UpcPolarisationJpsiIncoh {
130128 Configurable<int > nBinsPhi{" nBinsPhi" , 600 , " N bins in phi histo" };
131129 Configurable<float > lowPhi{" lowPhi" , -Pi, " lower limit in phi histo" };
132130 Configurable<float > highPhi{" highPhi" , Pi, " upper limit in phi histo" };
133- // pT of single muons
131+ // pT of single muons
134132 Configurable<int > nBinsPtSingle{" nBinsPtSingle" , 500 , " N bins in pT histo single muon" };
135133 Configurable<float > lowPtSingle{" lowPtSingle" , 0 ., " lower limit in pT histo single muon" };
136134 Configurable<float > highPtSingle{" highPtSingle" , 2 ., " upper limit in pT histo single muon" };
137- // eta of single muons
135+ // eta of single muons
138136 Configurable<int > nBinsEtaSingle{" nBinsEtaSingle" , 250 , " N bins in eta histo single muon" };
139137 Configurable<float > lowEtaSingle{" lowEtaSingle" , -4.5 , " lower limit in eta histo single muon" };
140138 Configurable<float > highEtaSingle{" highEtaSingle" , -2 ., " upper limit in eta histo single muon" };
141- // phi of single muons
139+ // phi of single muons
142140 Configurable<int > nBinsPhiSingle{" nBinsPhiSingle" , 600 , " N bins in phi histo single muon" };
143141 Configurable<float > lowPhiSingle{" lowPhiSingle" , -Pi, " lower limit in phi histo single muon" };
144142 Configurable<float > highPhiSingle{" highPhiSingle" , Pi, " upper limit in phi histo single muon" };
145- // ZDC
143+ // ZDC
146144 Configurable<int > nBinsZDCen{" nBinsZDCen" , 200 , " N bins in ZN energy" };
147145 Configurable<float > lowEnZN{" lowEnZN" , -50 ., " lower limit in ZN energy histo" };
148146 Configurable<float > highEnZN{" highEnZN" , 250 ., " upper limit in ZN energy histo" };
@@ -180,8 +178,6 @@ struct UpcPolarisationJpsiIncoh {
180178 registry.add (" hEtaTrkNeg" , " #eta of negative muons;;#counts" , kTH1D , {axisEtaSingle});
181179 registry.add (" hPhiTrkPos" , " #varphi of positive muons;;#counts" , kTH1D , {axisPhiSingle});
182180 registry.add (" hPhiTrkNeg" , " #varphi of negative muons;;#counts" , kTH1D , {axisPhiSingle});
183-
184-
185181 }
186182
187183 // template function that fills a map with the collision id of each udcollision as key
@@ -287,11 +283,11 @@ struct UpcPolarisationJpsiIncoh {
287283 }
288284 }
289285
290- // select opposite charge events only
286+ // select opposite charge events only
291287 if (cand.netCharge () != 0 ) {
292288 return ;
293289 }
294-
290+
295291 // MCH-MID match selection
296292 int nMIDs = 0 ;
297293 if (tr1.chi2MatchMCHMID () > 0 )
@@ -355,20 +351,18 @@ struct UpcPolarisationJpsiIncoh {
355351 registry.fill (HIST (" hPhiTrkPos" ), p1.Phi ());
356352 registry.fill (HIST (" hPhiTrkNeg" ), p2.Phi ());
357353
358-
359- // store the event to save it into a tree
354+ // store the event to save it into a tree
360355 if (tr1.sign () > 0 ) {
361356 dimuSel (cand.runNumber (),
362357 p.M (), p.E (), p.Px (), p.Py (), p.Pz (), p.Pt (), p.Rapidity (), p.Phi (),
363- p1.E (), p1.Px (), p1.Py (), p1.Pz (), p1.Pt (), p1.Eta (), p1.Phi (),
358+ p1.E (), p1.Px (), p1.Py (), p1.Pz (), p1.Pt (), p1.Eta (), p1.Phi (),
364359 p2.E (), p2.Px (), p2.Py (), p2.Pz (), p2.Pt (), p2.Eta (), p2.Phi ());
365360 } else {
366361 dimuSel (cand.runNumber (),
367362 p.M (), p.E (), p.Px (), p.Py (), p.Pz (), p.Pt (), p.Rapidity (), p.Phi (),
368363 p2.E (), p2.Px (), p2.Py (), p2.Pz (), p2.Pt (), p2.Eta (), p2.Phi (),
369364 p1.E (), p1.Px (), p1.Py (), p1.Pz (), p1.Pt (), p1.Eta (), p1.Phi ());
370365 }
371-
372366 }
373367 // PROCESS FUNCTION
374368 void processData (CandidatesFwd const & eventCandidates,
0 commit comments