2020#include " TDatabasePDG.h"
2121#include " PWGUD/Core/UPCHelpers.h"
2222#include " PWGUD/DataModel/UDTables.h"
23- #include " TLorentzVector.h"
23+ // #include "TLorentzVector.h"
2424// #include "TVector3.h"
2525#include " Math/LorentzVector.h" // ROOT::Math::LorentzVector
2626#include " Math/PxPyPzM4D.h" // ROOT::Math::PxPyPzM4D
@@ -180,7 +180,9 @@ struct ProcessMCDPMJetSGv3 {
180180 mcParticle.pz (),
181181 massPion
182182 );
183- if (std::fabs (protoMC.Eta ()) < 0.8 && protoMC.Pt () > 0.1 ) {
183+ double etaMax = 0.8 ;
184+ double ptMin = 0.1 ;
185+ if (std::fabs (protoMC.Eta ()) < etaMax && protoMC.Pt () > ptMin) {
184186 counter += 1 ;
185187 }
186188 if (!mcParticle.isPhysicalPrimary ())
@@ -205,7 +207,8 @@ struct ProcessMCDPMJetSGv3 {
205207 pMC.SetM (massProton);
206208 histos.fill (HIST (" ptGeneratedProton" ), pMC.Pt ());
207209 }
208- if (std::abs (pMC.Rapidity ()) < 0.8 ) {
210+ double yMax = 0.8 ;
211+ if (std::abs (pMC.Rapidity ()) < yMax) {
209212 if (std::abs (mcParticle.pdgCode ()) == codePion)
210213 histos.fill (HIST (" ptGeneratedPionAxE" ), pMC.Pt ());
211214 if (std::abs (mcParticle.pdgCode ()) == codeKaon)
@@ -243,20 +246,23 @@ struct ProcessMCDPMJetSGv3 {
243246 // massProton = protonPDG.Mass();
244247
245248 histos.fill (HIST (" numberOfTracksReco" ), tracks.size ());
249+ double etaMax = 0.8 ;
250+ double ptMin = 0.1 ;
251+ int NFindableMin = 70 ;
246252
247253 int counter = 0 ;
248254 for (auto & track : tracks) {
249255 if (track.isPVContributor ()) {
250256 int NFindable = track.tpcNClsFindable ();
251- if (NFindable < 70 ) {
257+ if (NFindable < ptMin ) {
252258 continue ;
253259 }
254260 // int NMinusFound = track.tpcNClsFindableMinusFound();
255261 // int NCluster = NFindable - NMinusFound;
256262 // if (NCluster < 70) {
257263 // continue;
258264 // }
259- if (track.pt () < 0.1 ) {
265+ if (track.pt () < ptMin ) {
260266 continue ;
261267 }
262268 // if (!(std::abs(track.dcaZ()) < 2.)) {
0 commit comments