|
17 | 17 | #include "Framework/AnalysisTask.h" |
18 | 18 | #include "Common/DataModel/TrackSelectionTables.h" |
19 | 19 | #include "Framework/ASoAHelpers.h" |
20 | | -#include "TDatabasePDG.h" |
| 20 | +// #include "TDatabasePDG.h" |
21 | 21 | #include "PWGUD/Core/UPCHelpers.h" |
22 | 22 | #include "PWGUD/DataModel/UDTables.h" |
23 | 23 | // #include "TLorentzVector.h" |
@@ -68,18 +68,21 @@ struct ProcessMCDPMJetSGv3 { |
68 | 68 |
|
69 | 69 | void init(InitContext const&) |
70 | 70 | { |
71 | | - TParticlePDG* pionPDG = fPDG->GetParticle(codePion); |
72 | | - if (pionPDG != nullptr) { |
73 | | - massPion = pionPDG->Mass(); |
74 | | - } |
75 | | - TParticlePDG* kaonPDG = fPDG->GetParticle(codeKaon); |
76 | | - if (kaonPDG != nullptr) { |
77 | | - massKaon = kaonPDG->Mass(); |
78 | | - } |
79 | | - TParticlePDG* protonPDG = fPDG->GetParticle(codeProton); |
80 | | - if (protonPDG != nullptr) { |
81 | | - massProton = protonPDG->Mass(); |
82 | | - } |
| 71 | + // TParticlePDG* pionPDG = fPDG->GetParticle(codePion); |
| 72 | + // if (pionPDG != nullptr) { |
| 73 | + // massPion = pionPDG->Mass(); |
| 74 | + // } |
| 75 | + // TParticlePDG* kaonPDG = fPDG->GetParticle(codeKaon); |
| 76 | + // if (kaonPDG != nullptr) { |
| 77 | + // massKaon = kaonPDG->Mass(); |
| 78 | + // } |
| 79 | + // TParticlePDG* protonPDG = fPDG->GetParticle(codeProton); |
| 80 | + // if (protonPDG != nullptr) { |
| 81 | + // massProton = protonPDG->Mass(); |
| 82 | + // } |
| 83 | + massPion = o2::constants::physics::MassPionCharged; |
| 84 | + massKaon = o2::constants::physics::MassKaonCharged; |
| 85 | + massProton = o2::constants::physics::MassProton; |
83 | 86 |
|
84 | 87 | // define axes you want to use |
85 | 88 | const AxisSpec axisCounter{10, 0, 10, ""}; |
@@ -273,7 +276,7 @@ struct ProcessMCDPMJetSGv3 { |
273 | 276 | // continue; |
274 | 277 | // } |
275 | 278 |
|
276 | | - double momentum = TMath::Sqrt(track.px() * track.px() + track.py() * track.py() + track.pz() * track.pz()); |
| 279 | + double momentum = sqrt(track.px() * track.px() + track.py() * track.py() + track.pz() * track.pz()); |
277 | 280 | double dEdx = track.tpcSignal(); |
278 | 281 | histos.fill(HIST("hdEdx"), momentum, dEdx); |
279 | 282 |
|
|
0 commit comments