@@ -79,12 +79,12 @@ DECLARE_SOA_COLUMN(TrackKaPID, trackKaPID, float[2]);
7979DECLARE_SOA_COLUMN (TrackDcaXY, trackDcaXY, float [2 ]);
8080DECLARE_SOA_COLUMN (TrackDcaZ, trackDcaZ, float [2 ]);
8181DECLARE_SOA_COLUMN (TrackTpcSignal, trackTpcSignal, float [2 ]);
82- } // namespace recoTree
82+ } // namespace reco_tree
8383DECLARE_SOA_TABLE (RecoTree, " AOD" , " RECOTREE" ,
8484 reco_tree::RunNumber, reco_tree::LocalBC, reco_tree::NumContrib, reco_tree::PosX, reco_tree::PosY, reco_tree::PosZ,
8585 reco_tree::TotalFT0AmplitudeA, reco_tree::TotalFT0AmplitudeC, reco_tree::TotalFV0AmplitudeA, reco_tree::TotalFDDAmplitudeA, reco_tree::TotalFDDAmplitudeC,
8686 reco_tree::TimeFT0A, reco_tree::TimeFT0C, reco_tree::TimeFV0A, reco_tree::TimeFDDA, reco_tree::TimeFDDC,
87- reco_tree::EnergyCommonZNA, reco_tree::EnergyCommonZNC, reco_tree::TimeZNA, reco_tree::TimeZNC,
87+ reco_tree::EnergyCommonZNA, reco_tree::EnergyCommonZNC, reco_tree::TimeZNA, reco_tree::TimeZNC,
8888 reco_tree::PhiRandom, reco_tree::PhiCharge, reco_tree::TrackSign, reco_tree::TrackPt, reco_tree::TrackEta, reco_tree::TrackPhi, reco_tree::TrackPiPID, reco_tree::TrackElPID, reco_tree::TrackKaPID, reco_tree::TrackDcaXY, reco_tree::TrackDcaZ, reco_tree::TrackTpcSignal);
8989
9090namespace mc_tree
@@ -102,10 +102,10 @@ DECLARE_SOA_COLUMN(TrackSign, trackSign, int[2]);
102102DECLARE_SOA_COLUMN (TrackPt, trackPt, float [2 ]);
103103DECLARE_SOA_COLUMN (TrackEta, trackEta, float [2 ]);
104104DECLARE_SOA_COLUMN (TrackPhi, trackPhi, float [2 ]);
105- } // namespace mcTree
105+ } // namespace mc_tree
106106DECLARE_SOA_TABLE (McTree, " AOD" , " MCTREE" ,
107107 mc_tree::LocalBc,
108- mc_tree::PosX, mc_tree::PosY, mc_tree::PosZ,
108+ mc_tree::PosX, mc_tree::PosY, mc_tree::PosZ,
109109 mc_tree::PhiRandom, mc_tree::PhiCharge, mc_tree::TrackSign, mc_tree::TrackPt, mc_tree::TrackEta, mc_tree::TrackPhi);
110110} // namespace o2::aod
111111
@@ -157,7 +157,7 @@ struct UpcRhoAnalysis {
157157 ConfigurableAxis momentumFromPhiAxis{" momentumFromPhiAxis" , {400 , -0.1 , 0.1 }, " p (GeV/#it{c})" };
158158 ConfigurableAxis znCommonEnergyAxis{" znCommonEnergyAxis" , {250 , -5.0 , 20.0 }, " ZN common energy (TeV)" };
159159 ConfigurableAxis znTimeAxis{" znTimeAxis" , {200 , -10.0 , 10.0 }, " ZN time (ns)" };
160-
160+
161161 HistogramRegistry rQC{" rQC" , {}, OutputObjHandlingPolicy::AnalysisObject};
162162 HistogramRegistry rTracks{" rTracks" , {}, OutputObjHandlingPolicy::AnalysisObject};
163163 HistogramRegistry rSystem{" rSystem" , {}, OutputObjHandlingPolicy::AnalysisObject};
@@ -244,7 +244,7 @@ struct UpcRhoAnalysis {
244244 rSystem.addClone (" system/selected/no-selection/" , " system/selected/Xn0n/" );
245245 rSystem.addClone (" system/selected/no-selection/" , " system/selected/0nXn/" );
246246 rSystem.addClone (" system/selected/no-selection/" , " system/selected/XnXn/" );
247-
247+
248248 // MC
249249 // collisions
250250 rMC.add (" MC/collisions/hPosXY" , " ;x (cm);y (cm);counts" , kTH2D , {{2000 , -0.1 , 0.1 }, {2000 , -0.1 , 0.1 }});
@@ -481,7 +481,7 @@ struct UpcRhoAnalysis {
481481 return system;
482482 }
483483
484- float getPhiRandom (const std::vector<TLorentzVector>& cutTracksLVs) // decay phi anisotropy
484+ float getPhiRandom (const std::vector<TLorentzVector>& cutTracksLVs) // decay phi anisotropy
485485 { // two possible definitions of phi: randomize the tracks
486486 int indices[2 ] = {0 , 1 };
487487 unsigned seed = std::chrono::system_clock::now ().time_since_epoch ().count (); // get time-based seed
@@ -533,7 +533,7 @@ struct UpcRhoAnalysis {
533533 if (collision.energyCommonZNA () > znCommonEnergyCut && std::abs (collision.timeZNA ()) < znTimeCut &&
534534 collision.energyCommonZNC () > znCommonEnergyCut && std::abs (collision.timeZNC ()) < znTimeCut)
535535 xnxn = true ;
536-
536+
537537 std::vector<decltype (tracks.begin ())> cutTracks; // store selected tracks
538538 for (const auto & track : tracks) {
539539 rQC.fill (HIST (" QC/tracks/hSelectionCounter" ), 0 );
@@ -546,7 +546,7 @@ struct UpcRhoAnalysis {
546546 cutTracks.push_back (track);
547547 }
548548 rQC.fill (HIST (" QC/tracks/selected/hRemainingTracks" ), cutTracks.size ());
549-
549+
550550 if (cutTracks.size () != 2 ) // further consider only two pion systems
551551 return ;
552552 for (int i = 0 ; i < static_cast <int >(cutTracks.size ()); i++)
@@ -562,12 +562,12 @@ struct UpcRhoAnalysis {
562562 trackLV.SetXYZM (track.px (), track.py (), track.pz (), o2::constants::physics::MassPionCharged); // apriori assume pion mass
563563 cutTracksLVs.push_back (trackLV);
564564 }
565-
565+
566566 // differentiate leading- and subleading-momentum tracks
567567 auto leadingMomentumTrack = momentum (cutTracks[0 ].px (), cutTracks[0 ].py (), cutTracks[0 ].pz ()) > momentum (cutTracks[1 ].px (), cutTracks[1 ].py (), cutTracks[1 ].pz ()) ? cutTracks[0 ] : cutTracks[1 ];
568568 auto subleadingMomentumTrack = (leadingMomentumTrack == cutTracks[0 ]) ? cutTracks[1 ] : cutTracks[0 ];
569569 rQC.fill (HIST (" QC/tracks/hTofHitCheck" ), leadingMomentumTrack.hasTOF (), subleadingMomentumTrack.hasTOF ());
570-
570+
571571 float leadingPt = leadingMomentumTrack.pt ();
572572 float subleadingPt = subleadingMomentumTrack.pt ();
573573 float leadingEta = eta (leadingMomentumTrack.px (), leadingMomentumTrack.py (), leadingMomentumTrack.pz ());
@@ -576,7 +576,7 @@ struct UpcRhoAnalysis {
576576 float subleadingPhi = phi (subleadingMomentumTrack.px (), subleadingMomentumTrack.py ());
577577 float phiRandom = getPhiRandom (cutTracksLVs);
578578 float phiCharge = getPhiCharge (cutTracks, cutTracksLVs);
579-
579+
580580 // fill recoTree
581581 int localBc = collision.globalBC () % o2::constants::lhc::LHCMaxBunches;
582582 int trackSigns[2 ] = {leadingMomentumTrack.sign (), subleadingMomentumTrack.sign ()};
@@ -593,9 +593,9 @@ struct UpcRhoAnalysis {
593593 recoTree (collision.runNumber (), localBc, collision.numContrib (), collision.posX (), collision.posY (), collision.posZ (),
594594 collision.totalFT0AmplitudeA (), collision.totalFT0AmplitudeC (), collision.totalFV0AmplitudeA (), collision.totalFDDAmplitudeA (), collision.totalFDDAmplitudeC (),
595595 collision.timeFT0A (), collision.timeFT0C (), collision.timeFV0A (), collision.timeFDDA (), collision.timeFDDC (),
596- collision.energyCommonZNA (), collision.energyCommonZNC (), collision.timeZNA (), collision.timeZNC (),
596+ collision.energyCommonZNA (), collision.energyCommonZNC (), collision.timeZNA (), collision.timeZNC (),
597597 phiRandom, phiCharge, trackSigns, trackPts, trackEtas, trackPhis, trackPiPIDs, trackElPIDs, trackKaPIDs, trackDcaXYs, trackDcaZs, trackTpcSignals);
598-
598+
599599 if (!tracksPassPiPID (cutTracks)) // apply PID cut
600600 return ;
601601 TLorentzVector system = reconstructSystem (cutTracksLVs);
@@ -604,7 +604,7 @@ struct UpcRhoAnalysis {
604604 float pT = system.Pt ();
605605 float rapidity = system.Rapidity ();
606606 float systemPhi = system.Phi () + o2::constants::math::PI;
607-
607+
608608 // fill raw histograms according to total charge
609609 switch (totalCharge) {
610610 case 0 :
@@ -701,7 +701,7 @@ struct UpcRhoAnalysis {
701701 mcParticlesLVs.push_back (pionLV);
702702 }
703703 rMC.fill (HIST (" MC/collisions/hNPions" ), cutMcParticles.size ());
704-
704+
705705 if (cutMcParticles.size () != 2 )
706706 return ;
707707 if (mcParticlesLVs.size () != cutMcParticles.size ())
0 commit comments