@@ -1641,7 +1641,7 @@ struct TauTau13topo {
16411641 {
16421642 return std::sqrt (px * px + py * py + pz * pz + mass * mass);
16431643 }
1644-
1644+
16451645 // float calculateDeltaPhi(TLorentzVector p, TLorentzVector p1)
16461646 float calculateDeltaPhi (ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double >> p, ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double >> p1)
16471647 {
@@ -2395,7 +2395,7 @@ struct TauTau13topo {
23952395 if (trk1.hasTPC ())
23962396 nPiHasTPC[trk.index ()]++;
23972397 // p1.SetXYZM(trk1.px(), trk1.py(), trk1.pz(), MassElectron);
2398- p1.SetXYZT (trk1.px (), trk1.py (), trk1.pz (), energy (trk1.px (), trk1.py (), trk1.pz (), MassElectron));
2398+ p1.SetXYZT (trk1.px (), trk1.py (), trk1.pz (), energy (trk1.px (), trk1.py (), trk1.pz (), MassElectron));
23992399 invMass2El[(counterTmp < 3 ? counterTmp : 5 - counterTmp)][(counterTmp < 3 ? 0 : 1 )] = (p + p1).mag2 ();
24002400 gammaPair[(counterTmp < 3 ? counterTmp : 5 - counterTmp)][(counterTmp < 3 ? 0 : 1 )] = (p + p1);
24012401 registry.get <TH1>(HIST (" control/cut0/hInvMass2ElAll" ))->Fill ((p + p1).mag2 ());
@@ -2492,7 +2492,7 @@ struct TauTau13topo {
24922492 trkTimeRes[counterTmp] = trk.trackTimeRes ();
24932493
24942494 // p1.SetXYZM(trk.px(), trk.py(), trk.pz(), MassPiPlus);
2495- p1.SetXYZT (trk.px (), trk.py (), trk.pz (), energy (trk.px (), trk.py (), trk.pz (),MassPiPlus));
2495+ p1.SetXYZT (trk.px (), trk.py (), trk.pz (), energy (trk.px (), trk.py (), trk.pz (), MassPiPlus));
24962496 tmpMomentum[counterTmp] = p1.P ();
24972497 tmpPt[counterTmp] = p1.Pt ();
24982498 tmpDedx[counterTmp] = trk.tpcSignal ();
@@ -3363,7 +3363,7 @@ struct TauTau13topo {
33633363 if (mcParticle.isPhysicalPrimary ()) {
33643364 countBoth++;
33653365 // if (mcParticle.pdgCode() != 22 && std::abs(mcParticle.pdgCode()) != 12 && std::abs(mcParticle.pdgCode()) != 14 && std::abs(mcParticle.pdgCode()) != 16 && mcParticle.pdgCode() != 130 && mcParticle.pdgCode() != 111) {
3366- if (mcParticle.pdgCode () != kGamma && std::abs (mcParticle.pdgCode ()) != kNuE && std::abs (mcParticle.pdgCode ()) != kNuMu && std::abs (mcParticle.pdgCode ()) != kNuTau && mcParticle.pdgCode () != kK0Long && mcParticle.pdgCode () != kPi0 ) {
3366+ if (mcParticle.pdgCode () != kGamma && std::abs (mcParticle.pdgCode ()) != kNuE && std::abs (mcParticle.pdgCode ()) != kNuMu && std::abs (mcParticle.pdgCode ()) != kNuTau && mcParticle.pdgCode () != kK0Long && mcParticle.pdgCode () != kPi0 ) {
33673367 countCharged++;
33683368
33693369 registryMC.get <TH1>(HIST (" globalMC/hMCetaGen" ))->Fill (mcParticle.eta ());
@@ -3374,7 +3374,7 @@ struct TauTau13topo {
33743374 if (mcParticle.has_mothers ()) {
33753375 auto const & mother = mcParticle.mothers_first_as <aod::McParticles>();
33763376 // if (std::abs(mother.pdgCode()) == 15) {
3377- if (std::abs (mother.pdgCode ()) == kTauMinus ) {
3377+ if (std::abs (mother.pdgCode ()) == kTauMinus ) {
33783378 countChargedFromTau++;
33793379 } // mother is tau
33803380 } // mc particle has mother
@@ -3427,7 +3427,7 @@ struct TauTau13topo {
34273427 } // end of electron check
34283428 // muon from tau
34293429 if (std::abs (daughter.pdgCode ()) == kMuonMinus ) { // 13
3430- if (daughter.pdgCode () == kMuonMinus ) // 13
3430+ if (daughter.pdgCode () == kMuonMinus ) // 13
34313431 flagMuPlusMuMinus = true ;
34323432 muonFound = !muonFound;
34333433 partPt = static_cast <float >(daughter.pt ());
@@ -3599,7 +3599,7 @@ struct TauTau13topo {
35993599 // LOGF(info, "<tautau13topo_MC> mcParticle pdg %d", mcParticle.pdgCode());
36003600 if (mcParticle.isPhysicalPrimary ()) {
36013601 // if (mcParticle.pdgCode() != 22 && std::abs(mcParticle.pdgCode()) != 12 && std::abs(mcParticle.pdgCode()) != 14 && std::abs(mcParticle.pdgCode()) != 16 && mcParticle.pdgCode() != 130 && mcParticle.pdgCode() != 111) {
3602- if (mcParticle.pdgCode () != kGamma && std::abs (mcParticle.pdgCode ()) != kNuE && std::abs (mcParticle.pdgCode ()) != kNuMu && std::abs (mcParticle.pdgCode ()) != kNuTau && mcParticle.pdgCode () != kK0Long && mcParticle.pdgCode () != kPi0 ) {
3602+ if (mcParticle.pdgCode () != kGamma && std::abs (mcParticle.pdgCode ()) != kNuE && std::abs (mcParticle.pdgCode ()) != kNuMu && std::abs (mcParticle.pdgCode ()) != kNuTau && mcParticle.pdgCode () != kK0Long && mcParticle.pdgCode () != kPi0 ) {
36033603 if (mcParticle.has_mothers ()) {
36043604 auto const & mother = mcParticle.mothers_first_as <aod::UDMcParticles>();
36053605 // LOGF(info, "<tautau13topo_MC> mcParticle has mother %d",mother.pdgCode());
@@ -3679,8 +3679,8 @@ struct TauTau13topo {
36793679 // if (!onlyPi) LOGF(info, "ERROR: should be 4 pions, but they are not!");
36803680 } // end of special check for pi + 3pi
36813681
3682- int index3ProngMC[3 ] = {0 , 0 , 0 }; // initialised of request
3683- if (index1ProngMC > 0 ) { // electron or muon case + 3pi
3682+ int index3ProngMC[3 ] = {0 , 0 , 0 }; // initialised of request
3683+ if (index1ProngMC > 0 ) { // electron or muon case + 3pi
36843684 int index3pi = 0 ;
36853685 for (int i = 0 ; i < 4 ; i++) {
36863686 if (index1ProngMC == indexProngMC[i])
@@ -3698,7 +3698,7 @@ struct TauTau13topo {
36983698 auto const & tmpPion2MC = mcParticles.begin () + index3ProngMC[1 ];
36993699 auto const & tmpPion3MC = mcParticles.begin () + index3ProngMC[2 ];
37003700
3701- // if (std::abs(tmpPion1MC.pdgCode()) == 211 && std::abs(tmpPion2MC.pdgCode()) == 211 && std::abs(tmpPion3MC.pdgCode()) == 211)
3701+ // if (std::abs(tmpPion1MC.pdgCode()) == 211 && std::abs(tmpPion2MC.pdgCode()) == 211 && std::abs(tmpPion3MC.pdgCode()) == 211)
37023702 if (std::abs (tmpPion1MC.pdgCode ()) == kPiPlus && std::abs (tmpPion2MC.pdgCode ()) == kPiPlus && std::abs (tmpPion3MC.pdgCode ()) == kPiPlus ) // 211 211 211
37033703 is3prong3PiMC = true ;
37043704
@@ -4030,9 +4030,9 @@ struct TauTau13topo {
40304030 registryMC.get <TH1>(HIST (" global1MCrec/hNTracksPV" ))->Fill (nPVTracks);
40314031
40324032 // TLorentzVector p, p1;
4033- ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double >> p, p1;
4033+ ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double >> p, p1;
40344034 // p.SetXYZM(0., 0., 0., 0.);
4035- p.SetXYZT (0 ., 0 ., 0 ., 0 .);
4035+ p.SetXYZT (0 ., 0 ., 0 ., 0 .);
40364036 TVector3 v1 (0 , 0 , 0 );
40374037 TVector3 v2 (0 , 0 , 0 );
40384038 float scalarPtsum = 0 ;
@@ -4067,12 +4067,12 @@ struct TauTau13topo {
40674067 registryMC.get <TH2>(HIST (" global1MCrec/hTrackEtaPhiPV" ))->Fill (tmpEtaData, tmpPhiData);
40684068 registryMC.get <TH1>(HIST (" global1MCrec/hTrackPtPV" ))->Fill (tmptrack.pt ());
40694069 // p1.SetXYZM(v1.X(), v1.Y(), v1.Z(), MassPiPlus); // in case of ghost
4070- p1.SetXYZT (v1.X (), v1.Y (), v1.Z (), energy (v1.X (), v1.Y (), v1.Z (), MassPiPlus)); // in case of ghost
4070+ p1.SetXYZT (v1.X (), v1.Y (), v1.Z (), energy (v1.X (), v1.Y (), v1.Z (), MassPiPlus)); // in case of ghost
40714071
40724072 if (trackMCId[i] >= 0 ) {
40734073 // p1.SetXYZM(v1.X(), v1.Y(), v1.Z(), (std::abs(tmptrack.udMcParticle().pdgCode()) == 211 ? MassPiPlus : MassElectron));
4074- // p1.SetXYZT(v1.X(), v1.Y(), v1.Z(), energy(v1.X(), v1.Y(), v1.Z(), (std::abs(tmptrack.udMcParticle().pdgCode()) == 211 ? MassPiPlus : MassElectron)));
4075- p1.SetXYZT (v1.X (), v1.Y (), v1.Z (), energy (v1.X (), v1.Y (), v1.Z (), (std::abs (tmptrack.udMcParticle ().pdgCode ()) == kPiPlus ? MassPiPlus : MassElectron))); // 211
4074+ // p1.SetXYZT(v1.X(), v1.Y(), v1.Z(), energy(v1.X(), v1.Y(), v1.Z(), (std::abs(tmptrack.udMcParticle().pdgCode()) == 211 ? MassPiPlus : MassElectron)));
4075+ p1.SetXYZT (v1.X (), v1.Y (), v1.Z (), energy (v1.X (), v1.Y (), v1.Z (), (std::abs (tmptrack.udMcParticle ().pdgCode ()) == kPiPlus ? MassPiPlus : MassElectron))); // 211
40764076 float tmpPt = pt (tmptrack.udMcParticle ().px (), tmptrack.udMcParticle ().py ());
40774077 float tmpEta = eta (tmptrack.udMcParticle ().px (), tmptrack.udMcParticle ().py (), tmptrack.udMcParticle ().pz ());
40784078 float tmpPhi = phi (tmptrack.udMcParticle ().px (), tmptrack.udMcParticle ().py ());
@@ -4231,10 +4231,10 @@ struct TauTau13topo {
42314231 // if (tmptrack.hasTOF()) trkHasTof[i] = true;
42324232 v1.SetXYZ (tmptrack.px (), tmptrack.py (), tmptrack.pz ());
42334233 // p1.SetXYZM(v1.X(), v1.Y(), v1.Z(), MassPiPlus); // in case of ghost
4234- p1.SetXYZT (v1.X (), v1.Y (), v1.Z (), energy (v1.X (), v1.Y (), v1.Z (), MassPiPlus)); // in case of ghost
4234+ p1.SetXYZT (v1.X (), v1.Y (), v1.Z (), energy (v1.X (), v1.Y (), v1.Z (), MassPiPlus)); // in case of ghost
42354235 if (trackMCId[i] >= 0 ) {
42364236 // p1.SetXYZM(v1.X(), v1.Y(), v1.Z(), (i == matchedElIndexToData ? MassElectron : MassPiPlus));
4237- p1.SetXYZT (v1.X (), v1.Y (), v1.Z (), energy (v1.X (), v1.Y (), v1.Z (), (i == matchedElIndexToData ? MassElectron : MassPiPlus)));
4237+ p1.SetXYZT (v1.X (), v1.Y (), v1.Z (), energy (v1.X (), v1.Y (), v1.Z (), (i == matchedElIndexToData ? MassElectron : MassPiPlus)));
42384238 }
42394239
42404240 nSigmaEl[counterTmp] = tmptrack.tpcNSigmaEl ();
0 commit comments