99// granted to it by virtue of its status as an Intergovernmental Organization
1010// or submit itself to any jurisdiction.
1111//
12- // \brief tau tau analysis 1e+3pi topology
13- 14- // \since January 2024
12+ // / \file upcTauTau13topo.cxx
13+ // / \brief tau tau analysis 1e+3pi topology
14+ 15+ // / \since January 2024
1516// to run it execute:
1617// copts="--configuration json://tautauConfig.json -b"
1718// o2-analysis-ud-tautau13topo $copts > output.log
@@ -2883,7 +2884,7 @@ struct TauTau13topo {
28832884 registryMC.get <TH1>(HIST (" electronMC/hMCpt" ))->Fill (daughter.pt ());
28842885
28852886 electronFound = !electronFound;
2886- partPt = static_cast <float > daughter.pt ();
2887+ partPt = static_cast <float > ( daughter.pt () );
28872888 // singleElectronGlobalIndex = daughter.globalIndex();
28882889 // LOGF(info,"e pt %f",daughter.pt());
28892890 if (std::abs (daughter.eta ()) > 0.9 )
@@ -2894,7 +2895,7 @@ struct TauTau13topo {
28942895 if (daughter.pdgCode () == 13 )
28952896 flagMuPlusMuMinus = true ;
28962897 muonFound = !muonFound;
2897- partPt = static_cast <float > daughter.pt ();
2898+ partPt = static_cast <float > ( daughter.pt () );
28982899 // LOGF(info,"mu pt %f",daughter.pt());
28992900 if (std::abs (daughter.eta ()) > 0.9 )
29002901 partFromTauInEta = false ;
@@ -2909,7 +2910,7 @@ struct TauTau13topo {
29092910 auto mcPartTmp = mcParticle.daughters_as <aod::McParticles>().begin () + singlePionIndex;
29102911 if (mcPartTmp.pdgCode () == -211 )
29112912 flagPiPlusPiMinus = true ;
2912- partPt = static_cast <float > mcPartTmp.pt ();
2913+ partPt = static_cast <float > ( mcPartTmp.pt () );
29132914 // motherOfSinglePionIndex = mcParticle.index();
29142915 if (std::abs (mcPartTmp.eta ()) > 0.9 )
29152916 partFromTauInEta = false ;
@@ -3481,10 +3482,10 @@ struct TauTau13topo {
34813482 float tmpPhiData = phi (tmptrack.px (), tmptrack.py ());
34823483 registryMC.get <TH2>(HIST (" global1MCrec/hTrackEtaPhiPV" ))->Fill (tmpEtaData, tmpPhiData);
34833484 registryMC.get <TH1>(HIST (" global1MCrec/hTrackPtPV" ))->Fill (tmptrack.pt ());
3484- p1.SetXYZM (v1.X (), v1.Y (), v1.Z (), pdg-> Mass ( 211 ) ); // in case of ghost
3485+ p1.SetXYZM (v1.X (), v1.Y (), v1.Z (), pionMass ); // in case of ghost
34853486
34863487 if (trackMCId[i] >= 0 ) {
3487- p1.SetXYZM (v1.X (), v1.Y (), v1.Z (), (std::abs (tmptrack.udMcParticle ().pdgCode ()) == 211 ? pdg-> Mass ( 211 ) : pdg-> Mass ( 11 ) ));
3488+ p1.SetXYZM (v1.X (), v1.Y (), v1.Z (), (std::abs (tmptrack.udMcParticle ().pdgCode ()) == 211 ? pionMass : electronMass ));
34883489 float tmpPt = pt (tmptrack.udMcParticle ().px (), tmptrack.udMcParticle ().py ());
34893490 float tmpEta = eta (tmptrack.udMcParticle ().px (), tmptrack.udMcParticle ().py (), tmptrack.udMcParticle ().pz ());
34903491 float tmpPhi = phi (tmptrack.udMcParticle ().px (), tmptrack.udMcParticle ().py ());
@@ -3640,9 +3641,9 @@ struct TauTau13topo {
36403641 auto const tmptrack = groupedTracks.begin () + trackId[i];
36413642 // if (tmptrack.hasTOF()) trkHasTof[i] = true;
36423643 v1.SetXYZ (tmptrack.px (), tmptrack.py (), tmptrack.pz ());
3643- p1.SetXYZM (v1.X (), v1.Y (), v1.Z (), pdg-> Mass ( 211 ) ); // in case of ghost
3644+ p1.SetXYZM (v1.X (), v1.Y (), v1.Z (), pionMass ); // in case of ghost
36443645 if (trackMCId[i] >= 0 ) {
3645- p1.SetXYZM (v1.X (), v1.Y (), v1.Z (), (i == matchedElIndexToData ? pdg-> Mass ( 11 ) : pdg-> Mass ( 211 ) ));
3646+ p1.SetXYZM (v1.X (), v1.Y (), v1.Z (), (i == matchedElIndexToData ? electronMass : pionMass ));
36463647 }
36473648
36483649 nSigmaEl[counterTmp] = tmptrack.tpcNSigmaEl ();
0 commit comments