3232#include " TMath.h"
3333#include " TRandom3.h"
3434
35- // table for saving tree with info on data
35+ // table for saving tree with info on data
3636namespace dimu
3737{
38- // dimuon
38+ // dimuon
3939DECLARE_SOA_COLUMN (RunNumber, runNumber, int );
4040DECLARE_SOA_COLUMN (M, m, float );
4141DECLARE_SOA_COLUMN (Pt, pt, float );
@@ -48,11 +48,7 @@ namespace o2::aod
4848DECLARE_SOA_TABLE (DiMu, " AOD" , " DIMU" ,
4949 dimu::RunNumber,
5050 dimu::M, dimu::Pt, dimu::Rap, dimu::Phi);
51- } // namespace o2::aod
52-
53-
54-
55-
51+ } // namespace o2::aod
5652using namespace o2 ;
5753using namespace o2 ::framework;
5854using namespace o2 ::framework::expressions;
@@ -86,7 +82,6 @@ struct upcPolarisationJPsiIncorr {
8682 Produces<o2::aod::DiMu> dimuSel;
8783 // defining histograms using histogram registry: different histos for the different process functions
8884 HistogramRegistry registry{" registry" , {}, OutputObjHandlingPolicy::AnalysisObject, true , true };
89-
9085 // CONFIGURABLES
9186 static constexpr double Pi = o2::constants::math::PI;
9287 // pT of muon pairs
@@ -109,7 +104,6 @@ struct upcPolarisationJPsiIncorr {
109104 Configurable<int > nBinsPhi{" nBinsPhi" , 600 , " N bins in phi histo" };
110105 Configurable<float > lowPhi{" lowPhi" , -Pi, " lower limit in phi histo" };
111106 Configurable<float > highPhi{" highPhi" , Pi, " upper limit in phi histo" };
112-
113107 // Analysis cuts
114108 Configurable<float > maxJpsiMass{" maxJpsiMass" , 3.18 , " Maximum of the jpsi peak for peak cut" };
115109 Configurable<float > minJpsiMass{" minJpsiMass" , 3.0 , " Minimum of the jpsi peak for peak cut" };
@@ -127,7 +121,6 @@ struct upcPolarisationJPsiIncorr {
127121 const AxisSpec axisEta{nBinsEta, lowEta, highEta, " #eta" };
128122 const AxisSpec axisRapidity{nBinsRapidity, lowRapidity, highRapidity, " Rapidity" };
129123 const AxisSpec axisPhi{nBinsPhi, lowPhi, highPhi, " #varphi" };
130-
131124 // histos
132125 // data and reco MC
133126 registry.add (" hMass" , " Invariant mass of muon pairs;;#counts" , kTH1D , {axisMass});
@@ -136,8 +129,6 @@ struct upcPolarisationJPsiIncorr {
136129 registry.add (" hRapidity" , " Rapidty of muon pairs;;#counts" , kTH1D , {axisRapidity});
137130 registry.add (" hPhi" , " #varphi of muon pairs;;#counts" , kTH1D , {axisPhi});
138131 }
139-
140-
141132 // retrieve particle mass (GeV/c^2) from TDatabasePDG
142133 float particleMass (int pid)
143134 {
@@ -223,7 +214,7 @@ struct upcPolarisationJPsiIncorr {
223214 p.SetXYZM (fwdTrack.px (), fwdTrack.py (), fwdTrack.pz (), mMu );
224215 float eta = p.Eta ();
225216 float pt = p.Pt ();
226-
217+
227218 if (eta < kEtaMin || eta > kEtaMax )
228219 return false ;
229220 if (pt < kPtMin )
@@ -252,7 +243,6 @@ struct upcPolarisationJPsiIncorr {
252243 return ;
253244 }
254245 }
255-
256246 // MCH-MID match selection
257247 int nMIDs = 0 ;
258248 if (tr1.chi2MatchMCHMID () > 0 )
@@ -261,7 +251,6 @@ struct upcPolarisationJPsiIncorr {
261251 nMIDs++;
262252 if (nMIDs != kReqMatchMIDTracks )
263253 return ;
264-
265254 // MFT-MID match selection (if MFT is requested by the trackType)
266255 if (myTrackType == 0 ) {
267256 // if MFT is requested check that the tracks is inside the MFT acceptance
@@ -305,8 +294,6 @@ struct upcPolarisationJPsiIncorr {
305294 return ;
306295 if (p.Rapidity () > highRapidity)
307296 return ;
308-
309-
310297 // fill the histos without looking at neutron emission
311298 registry.fill (HIST (" hMass" ), p.M ());
312299 registry.fill (HIST (" hPt" ), p.Pt ());
@@ -315,11 +302,7 @@ struct upcPolarisationJPsiIncorr {
315302 registry.fill (HIST (" hPhi" ), p.Phi ());
316303
317304 dimuSel (cand.runNumber (),p.M (), p.Pt (), p.Rapidity (), p.Phi ());
318-
319-
320305 }
321-
322-
323306 // PROCESS FUNCTION
324307 void processData (CandidatesFwd const & eventCandidates,
325308 o2::aod::UDZdcsReduced const & ZDCs,
0 commit comments