@@ -533,7 +533,7 @@ struct UpcTauRl {
533533 histos.add (" EventTwoTracks/ElectronMuPi/PID/hTOFnSigmaVsPPofO" , " ;Not-electron #it{p} (GeV/c);n#sigma^{#pi}_{TOF} (arb. units)" , HistType::kTH2D , {confAxis.zzAxisMom , confAxis.zzAxisNsigma });
534534 histos.add (" EventTwoTracks/ElectronMuPi/PID/hTOFnSigmaEvsnSigmaPofO" , " ;Not-electron n#sigma^{e}_{TOF} (arb. units);Not-electron n#sigma^{#pi}_{TOF} (arb. units)" , HistType::kTH2D , {confAxis.zzAxisNsigma , confAxis.zzAxisNsigma });
535535
536- if (doMixedEventsHistos){
536+ if (doMixedEventsHistos) {
537537 histos.add (" EventTwoTracks/MixedEvents/TwoElectrons/hInvariantMass" , " Mixed events;Invariant mass (GeV/c^{2});Number of events (-)" , HistType::kTH1D , {confAxis.zzAxisInvMass });
538538 histos.add (" EventTwoTracks/MixedEvents/TwoElectrons/hInvariantMassWide" , " Mixed events;Invariant mass (GeV/c^{2});Number of events (-)" , HistType::kTH1D , {confAxis.zzAxisInvMassWide });
539539 histos.add (" EventTwoTracks/MixedEvents/TwoElectrons/hAcoplanarity" , " Mixed events;#Delta#phi (rad);Number of events (-)" , HistType::kTH1D , {confAxis.zzAxisAcoplanarity });
@@ -592,7 +592,6 @@ struct UpcTauRl {
592592 histos.add (" EventTwoTracks/MixedEvents/ElectronMuPi/hElectronPhiVsOtherPhi" , " Mixed events;Electron #phi (rad); #mu/#pi #phi (rad)" , HistType::kTH2D , {confAxis.zzAxisPhi , confAxis.zzAxisPhi });
593593 histos.add (" EventTwoTracks/MixedEvents/ElectronMuPi/hElectronRapVsOtherRap" , " Mixed events;Electron #it{y} (-); #mu/#pi #it{y} (-)" , HistType::kTH2D , {confAxis.zzAxisRap , confAxis.zzAxisRap });
594594 }
595-
596595 }
597596
598597 if (doTruthHistos) {
@@ -1191,15 +1190,15 @@ struct UpcTauRl {
11911190 }
11921191 }
11931192 if (isElMuPion) {
1194- if (cutTauEvent.useThresholdsPID ){
1195- if (isElectronCandidate (trkDaug1)){
1196- vecMixElMupion.push_back (std::make_pair (vecPVnewPIDidx[0 ], vecPVnewPIDidx[1 ]));// storing electron first
1193+ if (cutTauEvent.useThresholdsPID ) {
1194+ if (isElectronCandidate (trkDaug1)) {
1195+ vecMixElMupion.push_back (std::make_pair (vecPVnewPIDidx[0 ], vecPVnewPIDidx[1 ])); // storing electron first
11971196 } else {
11981197 vecMixElMupion.push_back (std::make_pair (vecPVnewPIDidx[1 ], vecPVnewPIDidx[0 ]));
11991198 }
12001199 } else {
1201- if (enumMyParticle (trackPDG (trkDaug1, cutPID.cutSiTPC , cutPID.cutSiTOF , cutPID.usePIDwTOF , cutPID.useScutTOFinTPC )) == P_ELECTRON){
1202- vecMixElMupion.push_back (std::make_pair (vecPVidx[0 ], vecPVidx[1 ]));// storing electron first
1200+ if (enumMyParticle (trackPDG (trkDaug1, cutPID.cutSiTPC , cutPID.cutSiTOF , cutPID.usePIDwTOF , cutPID.useScutTOFinTPC )) == P_ELECTRON) {
1201+ vecMixElMupion.push_back (std::make_pair (vecPVidx[0 ], vecPVidx[1 ])); // storing electron first
12031202 } else {
12041203 vecMixElMupion.push_back (std::make_pair (vecPVidx[1 ], vecPVidx[0 ]));
12051204 }
@@ -2013,8 +2012,8 @@ struct UpcTauRl {
20132012 void fillMixedEventHistograms (Ts const & reconstructedBarrelTracks)
20142013 {
20152014 TLorentzVector mother, daug[2 ];
2016- for (int idx = 0 ; idx < static_cast <int >(vecMixElEl.size ()); idx++){
2017- for (int cnt = 0 ; cnt < 5 ; cnt++){
2015+ for (int idx = 0 ; idx < static_cast <int >(vecMixElEl.size ()); idx++) {
2016+ for (int cnt = 0 ; cnt < 5 ; cnt++) {
20182017 std::random_device rand_dev;
20192018 std::mt19937 generator (rand_dev ());
20202019 std::uniform_int_distribution<int > distr (0 , static_cast <int >(vecMixElEl.size ()));
@@ -2055,11 +2054,11 @@ struct UpcTauRl {
20552054 histos.get <TH2>(HIST (" EventTwoTracks/MixedEvents/TwoElectrons/hLeadingPtVsOtherPt" ))->Fill (((daug[0 ].P () > daug[1 ].P ()) ? daug[0 ].Pt () : daug[1 ].Pt ()), ((daug[0 ].P () > daug[1 ].P ()) ? daug[1 ].Pt () : daug[0 ].Pt ()));
20562055 histos.get <TH2>(HIST (" EventTwoTracks/MixedEvents/TwoElectrons/hLeadingPhiVsOtherPhi" ))->Fill (((daug[0 ].P () > daug[1 ].P ()) ? daug[0 ].Phi () : daug[1 ].Phi ()), ((daug[0 ].P () > daug[1 ].P ()) ? daug[1 ].Phi () : daug[0 ].Phi ()));
20572056 histos.get <TH2>(HIST (" EventTwoTracks/MixedEvents/TwoElectrons/hLeadingRapVsOtherRap" ))->Fill (((daug[0 ].P () > daug[1 ].P ()) ? daug[0 ].Rapidity () : daug[1 ].Rapidity ()), ((daug[0 ].P () > daug[1 ].P ()) ? daug[1 ].Rapidity () : daug[0 ].Rapidity ()));
2058- }// cnt
2059- }// idx
2057+ } // cnt
2058+ } // idx
20602059
2061- for (int idx = 0 ; idx < static_cast <int >(vecMixElMupion.size ()); idx++){
2062- for (int cnt = 0 ; cnt < 5 ; cnt++){
2060+ for (int idx = 0 ; idx < static_cast <int >(vecMixElMupion.size ()); idx++) {
2061+ for (int cnt = 0 ; cnt < 5 ; cnt++) {
20632062 std::random_device rand_dev;
20642063 std::mt19937 generator (rand_dev ());
20652064 std::uniform_int_distribution<int > distr (0 , static_cast <int >(vecMixElMupion.size ()));
@@ -2076,7 +2075,7 @@ struct UpcTauRl {
20762075 const auto & electronPt = daug[0 ].Pt ();
20772076 const auto & electronP = daug[0 ].P ();
20782077 const auto & electronE = daug[0 ].E ();
2079- const auto & mupionPt = daug[1 ].Pt ();
2078+ const auto & mupionPt = daug[1 ].Pt ();
20802079 const auto & mupionP = daug[1 ].P ();
20812080 const auto & mupionE = daug[1 ].E ();
20822081 histos.get <TH1>(HIST (" EventTwoTracks/MixedEvents/ElectronMuPi/hInvariantMass" ))->Fill (mother.M ());
@@ -2109,15 +2108,15 @@ struct UpcTauRl {
21092108 histos.get <TH2>(HIST (" EventTwoTracks/MixedEvents/ElectronMuPi/hElectronPtVsOtherPt" ))->Fill (electronPt, mupionPt);
21102109 histos.get <TH2>(HIST (" EventTwoTracks/MixedEvents/ElectronMuPi/hElectronPhiVsOtherPhi" ))->Fill (isElectronCandidate (trkDaug1) ? daug[0 ].Phi () : daug[1 ].Phi (), isElectronCandidate (trkDaug1) ? daug[1 ].Phi () : daug[0 ].Phi ());
21112110 histos.get <TH2>(HIST (" EventTwoTracks/MixedEvents/ElectronMuPi/hElectronRapVsOtherRap" ))->Fill (isElectronCandidate (trkDaug1) ? daug[0 ].Rapidity () : daug[1 ].Rapidity (), isElectronCandidate (trkDaug1) ? daug[1 ].Rapidity () : daug[0 ].Rapidity ());
2112- }// cnt
2113- }// idx
2111+ } // cnt
2112+ } // idx
21142113 } // end fillMixedEventHistograms
21152114
21162115 void processDataDG (FullUDCollisions const & reconstructedCollisions,
21172116 FullUDTracks const & reconstructedBarrelTracks)
21182117 {
21192118
2120- for (const auto & reconstructedCollision : reconstructedCollisions){
2119+ for (const auto & reconstructedCollision : reconstructedCollisions) {
21212120 if (!isGoodROFtime (reconstructedCollision))
21222121 return ;
21232122
@@ -2146,7 +2145,7 @@ struct UpcTauRl {
21462145 void processDataSG (FullSGUDCollisions const & reconstructedCollisions,
21472146 FullUDTracks const & reconstructedBarrelTracks)
21482147 {
2149- for (const auto & reconstructedCollision : reconstructedCollisions){
2148+ for (const auto & reconstructedCollision : reconstructedCollisions) {
21502149
21512150 int gapSide = reconstructedCollision.gapSide ();
21522151 int trueGapSide = sgSelector.trueGap (reconstructedCollision, cutSample.cutTrueGapSideFV0 , cutSample.cutTrueGapSideFT0A , cutSample.cutTrueGapSideFT0C , cutSample.cutTrueGapSideZDC );
@@ -2188,7 +2187,7 @@ struct UpcTauRl {
21882187 aod::UDMcParticles const &)
21892188 {
21902189 isMC = true ;
2191- for (const auto & reconstructedCollision : reconstructedCollisions){
2190+ for (const auto & reconstructedCollision : reconstructedCollisions) {
21922191
21932192 if (!isGoodROFtime (reconstructedCollision))
21942193 return ;
@@ -2231,7 +2230,7 @@ struct UpcTauRl {
22312230 aod::UDMcParticles const &)
22322231 {
22332232 isMC = true ;
2234- for (const auto & reconstructedCollision : reconstructedCollisions){
2233+ for (const auto & reconstructedCollision : reconstructedCollisions) {
22352234
22362235 int gapSide = reconstructedCollision.gapSide ();
22372236 histos.fill (HIST (" Events/UDtableGapSide" ), gapSide);
@@ -2271,7 +2270,6 @@ struct UpcTauRl {
22712270 fillPIDhistograms (reconstructedCollision, reconstructedBarrelTracks);
22722271 fillMCPIDhistograms (reconstructedBarrelTracks);
22732272 }
2274-
22752273 }
22762274
22772275 } // end processMCrecDG
0 commit comments