@@ -126,7 +126,7 @@ struct StudyPnch {
126126 x->SetBinLabel (6 , " INEL > 0" );
127127 x->SetBinLabel (7 , " |vz| < 10" );
128128
129- if (doprocessData || doprocessCorrelation || doprocessMonteCarlo || doprocessTreatedMonteCarlo ) {
129+ if (doprocessData || doprocessCorrelation || doprocessMonteCarlo) {
130130 histos.add (" PhiVsEtaHist" , " PhiVsEtaHist" , kTH2F , {axisPhi, axisEta}, false );
131131 }
132132 if (doprocessData) {
@@ -144,11 +144,6 @@ struct StudyPnch {
144144 histos.add (" hMultiplicityMCgen" , " hMultiplicityMCgen" , kTH1F , {axisMult}, true );
145145 histos.add (" hResponseMatrix" , " hResponseMatrix" , kTH2F , {axisMult, axisMult}, true );
146146 }
147- if (doprocessTreatedMonteCarlo) {
148- histos.add (" hMultiplicityTreatMCrec" , " hMultiplicityTreatMCrec" , kTH1F , {axisMult}, true );
149- histos.add (" hMultiplicityTreatMCgen" , " hMultiplicityTreatMCgen" , kTH1F , {axisMult}, true );
150- histos.add (" hResponseMatrixTreat" , " hResponseMatrixTreat" , kTH2F , {axisMult, axisMult}, true );
151- }
152147 if (doprocessEvtLossSigLossMC) {
153148 histos.add (" MCEventHist" , " MCEventHist" , kTH1F , {axisEvent}, false );
154149 auto hstat = histos.get <TH1>(HIST (" MCEventHist" ));
@@ -338,26 +333,6 @@ struct StudyPnch {
338333 }
339334 }
340335
341- void processTreatedMonteCarlo (ColMCTrueTable::iterator const & mcCollision, ColMCRecTable const & RecCols, TrackMCTrueTable const & GenParticles, FilTrackMCRecTable const & RecTracks)
342- {
343- // Count generated tracks at each iterator
344- auto multgen = countGenTracks (GenParticles, mcCollision);
345- histos.fill (HIST (" hMultiplicityTreatMCgen" ), multgen);
346- for (const auto & RecCol : RecCols) {
347- if (!isEventSelected (RecCol)) {
348- continue ;
349- }
350- // Verify that the reconstructed collision corresponds to the given MC collision
351- if (RecCol.mcCollisionId () != mcCollision.globalIndex ()) {
352- continue ;
353- }
354- auto recTracksPart = RecTracks.sliceBy (perCollision, RecCol.globalIndex ());
355- auto multrec = countNTracksMcCol (recTracksPart, RecCol);
356- histos.fill (HIST (" hMultiplicityTreatMCrec" ), multrec);
357- histos.fill (HIST (" hResponseMatrixTreat" ), multrec, multgen);
358- }
359- }
360-
361336 void processEvtLossSigLossMC (soa::Join<ColMCTrueTable, aod::MultMCExtras>::iterator const & mcCollision, ColMCRecTable const & RecCols, TrackMCTrueTable const & GenParticles)
362337 {
363338 if (isApplyInelgt0 && !mcCollision.isInelGt0 ()) {
@@ -402,7 +377,6 @@ struct StudyPnch {
402377 PROCESS_SWITCH (StudyPnch, processData, " process data CentFT0C" , false );
403378 PROCESS_SWITCH (StudyPnch, processCorrelation, " do correlation study in data" , false );
404379 PROCESS_SWITCH (StudyPnch, processMonteCarlo, " process MC CentFT0C" , false );
405- PROCESS_SWITCH (StudyPnch, processTreatedMonteCarlo, " process Treated MC CentFT0C" , false );
406380 PROCESS_SWITCH (StudyPnch, processEvtLossSigLossMC, " process Signal Loss, Event Loss" , false );
407381};
408382
0 commit comments