2222#include " PWGUD/DataModel/UDTables.h"
2323// #include "TLorentzVector.h"
2424// #include "TVector3.h"
25- #include " Math/LorentzVector.h" // ROOT::Math::LorentzVector
26- #include " Math/PxPyPzM4D.h" // ROOT::Math::PxPyPzM4D
25+ #include " Math/LorentzVector.h" // ROOT::Math::LorentzVector
26+ #include " Math/PxPyPzM4D.h" // ROOT::Math::PxPyPzM4D
2727#include " TMath.h"
2828
2929using namespace o2 ;
@@ -56,9 +56,6 @@ struct ProcessMCDPMJetSGv3 {
5656 using TC = TCs::iterator;
5757 using LorentzVectorM = ROOT::Math::LorentzVector<ROOT::Math::PxPyPzM4D<double >>;
5858
59-
60-
61-
6259 double massPion = 0 .;
6360 double massKaon = 0 .;
6461 double massProton = 0 .;
@@ -144,6 +141,7 @@ struct ProcessMCDPMJetSGv3 {
144141 histos.add (" allreconstructedPFProtonTOF" , " allreconstructedPFProtonTOF" , kTH1F , {axisPt});
145142
146143 histos.add (" numberOfRecoCollisions" , " numberOfRecoCollisions" , kTH1F , {{100 , -0 .5f , 99 .5f }});
144+ histos.add (" numberOfRecoCollisions2" , " numberOfRecoCollisions2" , kTH1F , {{100 , -0 .5f , 99 .5f }});
147145 histos.add (" numberOfTracksMC" , " numberOfTracksMC" , kTH1F , {{100 , -0 .5f , 99 .5f }});
148146 histos.add (" numberOfTracksReco" , " numberOfTracksReco" , kTH1F , {{100 , -0 .5f , 99 .5f }});
149147
@@ -156,7 +154,7 @@ struct ProcessMCDPMJetSGv3 {
156154 // -----------------------------------------------------------------------------------------------------------------------
157155 void processSim (aod::UDMcCollision const & mcCollision, aod::UDMcParticles const & mcParticles)
158156 {
159- histos.fill (HIST (" eventCounter" ), 0.5 );
157+ // histos.fill(HIST("eventCounter"), 0.5);
160158
161159 // auto massPion = 0.;
162160 // TParticlePDG pionPDG = fPDG->GetParticle(codePion);
@@ -168,6 +166,7 @@ struct ProcessMCDPMJetSGv3 {
168166 // TParticlePDG protonPDG = fPDG->GetParticle(codeProton);
169167 // massProton = protonPDG.Mass();
170168 histos.fill (HIST (" numberOfTracksMC" ), mcParticles.size ());
169+ histos.fill (HIST (" eventCounter" ), mcCollision.size ());
171170 // LOGF(info, "New event! mcParticles.size() = %d", mcParticles.size());
172171
173172 int counterMC = 0 ;
@@ -178,11 +177,10 @@ struct ProcessMCDPMJetSGv3 {
178177 counterMC += 1 ;
179178 // if(mcParticle.isPhysicalPrimary()) counterMC += 1;
180179 LorentzVectorM protoMC (
181- mcParticle.px (),
182- mcParticle.py (),
183- mcParticle.pz (),
184- massPion
185- );
180+ mcParticle.px (),
181+ mcParticle.py (),
182+ mcParticle.pz (),
183+ massPion);
186184 double etaMax = 0.8 ;
187185 double ptMin = 0.1 ;
188186 if (std::fabs (protoMC.Eta ()) < etaMax && protoMC.Pt () > ptMin) {
@@ -234,7 +232,9 @@ struct ProcessMCDPMJetSGv3 {
234232 // aod::UDMcCollisions const& /*mccollisions*/,
235233 aod::UDMcParticles const & mcParticles)
236234 {
237- histos.fill (HIST (" numberOfRecoCollisions" ), 2 .); // number of times coll was reco-ed
235+ histos.fill (HIST (" numberOfRecoCollisions" ), 88 .); // number of times coll was reco-ed
236+ histos.fill (HIST (" numberOfRecoCollisions" ), collision.size ()); // number of times coll was reco-ed
237+ histos.fill (HIST (" numberOfRecoCollisions2" ), mcParticles.size ());
238238 Partition<TCs> pvContributors = aod::udtrack::isPVContributor == true ;
239239 pvContributors.bindTable (tracks);
240240
@@ -249,7 +249,7 @@ struct ProcessMCDPMJetSGv3 {
249249 // massProton = protonPDG.Mass();
250250
251251 histos.fill (HIST (" numberOfTracksReco" ), tracks.size ());
252- double etaMax = 0.8 ;
252+ // double etaMax = 0.8;
253253 double yMax = 0.8 ;
254254 double sigmaMax = 3 .;
255255 double ptMin = 0.1 ;
@@ -270,13 +270,13 @@ struct ProcessMCDPMJetSGv3 {
270270 if (track.pt () < ptMin) {
271271 continue ;
272272 }
273- // if (!(std::abs(track.dcaZ()) < 2.)) {
274- // continue;
275- // }
273+ if (!(std::abs (track.dcaZ ()) < 2 .)) {
274+ continue ;
275+ }
276276 double dcaLimit = 0.0105 + 0.035 / std::pow (track.pt (), 1.1 );
277- // if (!(std::abs(track.dcaXY()) < dcaLimit)) {
278- // continue;
279- // }
277+ if (!(std::abs (track.dcaXY ()) < dcaLimit)) {
278+ continue ;
279+ }
280280
281281 double momentum = std::sqrt (track.px () * track.px () + track.py () * track.py () + track.pz () * track.pz ());
282282 double dEdx = track.tpcSignal ();
0 commit comments