@@ -543,24 +543,26 @@ struct TableMakerMC {
543543 if (mcflags == 0 ) {
544544 continue ;
545545 }
546-
546+
547547 auto mcCollision = mctrack.template mcCollision_as <MyEventsMcWithMults>();
548-
548+
549549 bool isjpsi = false ;
550- if (mctrack.pdgCode () == 443 ) {
550+ if (mctrack.pdgCode () == 443 ) {
551551 isjpsi = true ;
552552 }
553553
554- // add rapidity and pT cuts for mcparticles
555- if (isjpsi && (fConfigJpsiMinPtMC .value > 0 && std::abs (mctrack.pt ()) < fConfigJpsiMinPtMC .value )) continue ;
556- if (isjpsi && (fConfigJpsiMaxPtMC .value > 0 && std::abs (mctrack.pt ()) > fConfigJpsiMaxPtMC .value )) continue ;
557- if (isjpsi && (fConfigTrackMaxRapidityMC .value > 0 && std::abs (mctrack.y ()) > fConfigTrackMaxRapidityMC .value )) continue ;
554+ // add rapidity and pT cuts for mcparticles
555+ if (isjpsi && (fConfigJpsiMinPtMC .value > 0 && std::abs (mctrack.pt ()) < fConfigJpsiMinPtMC .value ))
556+ continue ;
557+ if (isjpsi && (fConfigJpsiMaxPtMC .value > 0 && std::abs (mctrack.pt ()) > fConfigJpsiMaxPtMC .value ))
558+ continue ;
559+ if (isjpsi && (fConfigTrackMaxRapidityMC .value > 0 && std::abs (mctrack.y ()) > fConfigTrackMaxRapidityMC .value ))
560+ continue ;
558561
559- // Save jpsi events that satisfied the cut conditions
560- if (fConfigApplyselctionMC && isjpsi && ((std::abs (mctrack.pt ()) > fConfigJpsiMinPtMC .value ) && (std::abs (mctrack.pt ()) < fConfigJpsiMaxPtMC .value ) && (std::abs (mctrack.y ()) < fConfigTrackMaxRapidityMC .value )))
561- {
562+ // Save jpsi events that satisfied the cut conditions
563+ if (fConfigApplyselctionMC && isjpsi && ((std::abs (mctrack.pt ()) > fConfigJpsiMinPtMC .value ) && (std::abs (mctrack.pt ()) < fConfigJpsiMaxPtMC .value ) && (std::abs (mctrack.y ()) < fConfigTrackMaxRapidityMC .value ))) {
562564 labelsToReserve.insert (mcCollision.globalIndex ());
563- }
565+ }
564566
565567 // If this MC track was not already added to the map, add it now
566568 if (fLabelsMap .find (mctrack.globalIndex ()) == fLabelsMap .end ()) {
@@ -785,11 +787,12 @@ struct TableMakerMC {
785787 // However, in data analysis one should loop over associations, so this one should not be used.
786788 // In the case of Run2-like analysis, there will be no associations, so this ID will be the one originally assigned in the AO2Ds (updated for the skims)
787789 uint32_t reducedEventIdx = fCollIndexMap [track.collisionId ()];
788-
790+
789791 // apply the event selection for reconstrcted track
790792 if (fConfigApplyselctionMC && track.has_mcParticle ()) {
791- auto mcCollision = collision.template mcCollision_as <MyEventsMcWithMults>();
792- if ((labelsToReserve.count (mcCollision.globalIndex ()) == 0 )) continue ;
793+ auto mcCollision = collision.template mcCollision_as <MyEventsMcWithMults>();
794+ if ((labelsToReserve.count (mcCollision.globalIndex ()) == 0 ))
795+ continue ;
793796 }
794797
795798 // NOTE: trackBarrelInfo stores the index of the collision as in AO2D (for use in some cases where the analysis on skims is done
0 commit comments