4141#include " ReconstructionDataFormats/Track.h"
4242
4343#include < Math/Vector4D.h>
44- #include < TDatabasePDG.h>
4544#include < TFile.h>
4645#include < TH2F.h>
47- #include < TLorentzVector.h>
4846#include < TPDGCode.h>
4947#include < TProfile.h>
5048
@@ -1011,7 +1009,7 @@ struct sigmaanalysis {
10111009
10121010 // _______________________________________
10131011 // Gamma MC association
1014- if (sigma.photonPDGCode () == 22 ) {
1012+ if (sigma.photonPDGCode () == PDG_t:: kGamma ) {
10151013 if (sigma.photonmcpt () > 0 ) {
10161014 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/Reso/h2dGammaPtResolution" ), sigma.photonmcpt (), (sigma.photonPt ()/sigma.photonmcpt ()) - 1 .f ); // pT resolution
10171015 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/Reso/h2dGammaInvPtResolution" ), 1 .f / sigma.photonmcpt (), 1 .f / sigma.photonPt () - 1 .f / sigma.photonmcpt ()); // pT resolution
@@ -1020,7 +1018,7 @@ struct sigmaanalysis {
10201018
10211019 // _______________________________________
10221020 // Lambda MC association
1023- if (sigma.lambdaPDGCode () == 3122 ) {
1021+ if (sigma.lambdaPDGCode () == PDG_t:: kLambda0 ) {
10241022 if (sigma.lambdamcpt () > 0 ) {
10251023 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/Reso/h2dLambdaPtResolution" ), 1 .f / sigma.lambdamcpt (), 1 .f / sigma.lambdaPt () - 1 .f / sigma.lambdamcpt ()); // 1/pT resolution
10261024 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/Reso/h3dLambdaPtResoVsTPCCR" ), 1 .f / sigma.lambdamcpt (), 1 .f / sigma.lambdaPt () - 1 .f / sigma.lambdamcpt (), -1 * sigma.lambdaNegTPCCrossedRows ()); // 1/pT resolution
@@ -1030,7 +1028,7 @@ struct sigmaanalysis {
10301028
10311029 // _______________________________________
10321030 // AntiLambda MC association
1033- if (sigma.lambdaPDGCode () == - 3122 ) {
1031+ if (sigma.lambdaPDGCode () == PDG_t:: kLambda0Bar ) {
10341032 if (sigma.lambdamcpt () > 0 ) {
10351033 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/Reso/h2dAntiLambdaPtResolution" ), 1 .f / sigma.lambdamcpt (), 1 .f / sigma.lambdaPt () - 1 .f / sigma.lambdamcpt ()); // pT resolution
10361034 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/Reso/h3dAntiLambdaPtResoVsTPCCR" ), 1 .f / sigma.lambdamcpt (), 1 .f / sigma.lambdaPt () - 1 .f / sigma.lambdamcpt (), -1 * sigma.lambdaNegTPCCrossedRows ()); // 1/pT resolution
@@ -1074,24 +1072,24 @@ struct sigmaanalysis {
10741072
10751073 // _______________________________________
10761074 // Real Gamma x Real Lambda - but not from the same sigma0/antisigma0!
1077- if ((PhotonPDGCode == 22 ) && ((LambdaPDGCode == 3122 ) || (LambdaPDGCode == - 3122 )) && (!fIsSigma && !fIsAntiSigma )) {
1075+ if ((PhotonPDGCode == PDG_t:: kGamma ) && ((LambdaPDGCode == PDG_t:: kLambda0 ) || (LambdaPDGCode == PDG_t:: kLambda0Bar )) && (!fIsSigma && !fIsAntiSigma )) {
10781076 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/BkgStudy/h2dPtVsMassSigma_TrueDaughters" ), sigmapT, sigmaMass);
10791077 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/BkgStudy/h2dTrueDaughtersMatrix" ), LambdaPDGCodeMother, PhotonPDGCodeMother);
10801078 }
10811079
10821080 // _______________________________________
10831081 // Real Gamma x fake Lambda
1084- if ((PhotonPDGCode == 22 ) && (LambdaPDGCode != 3122 ) && (LambdaPDGCode != - 3122 ))
1082+ if ((PhotonPDGCode == PDG_t:: kGamma ) && (LambdaPDGCode != PDG_t:: kLambda0 ) && (LambdaPDGCode != PDG_t:: kLambda0Bar ))
10851083 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/BkgStudy/h2dPtVsMassSigma_TrueGammaFakeLambda" ), sigmapT, sigmaMass);
10861084
10871085 // _______________________________________
10881086 // Fake Gamma x Real Lambda
1089- if ((PhotonPDGCode != 22 ) && ((LambdaPDGCode == 3122 ) || (LambdaPDGCode == - 3122 )))
1087+ if ((PhotonPDGCode != PDG_t:: kGamma ) && ((LambdaPDGCode == PDG_t:: kLambda0 ) || (LambdaPDGCode == PDG_t:: kLambda0Bar )))
10901088 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/BkgStudy/h2dPtVsMassSigma_FakeGammaTrueLambda" ), sigmapT, sigmaMass);
10911089
10921090 // _______________________________________
10931091 // Fake Gamma x Fake Lambda
1094- if ((PhotonPDGCode != 22 ) && (LambdaPDGCode != 3122 ) && (LambdaPDGCode != - 3122 ))
1092+ if ((PhotonPDGCode != PDG_t:: kGamma ) && (LambdaPDGCode != PDG_t:: kLambda0 ) && (LambdaPDGCode != PDG_t:: kLambda0Bar ))
10951093 histos.fill (HIST (MainDir[mode]) + HIST (" /MC/BkgStudy/h2dPtVsMassSigma_FakeDaughters" ), sigmapT, sigmaMass);
10961094 }
10971095
0 commit comments