@@ -220,7 +220,7 @@ struct taggingHFE {
220220
221221 struct : ConfigurableGroup {
222222 std::string prefix = " eventCut" ;
223- // Configurable<int> cfgEventGeneratorType{"cfgEventGeneratorType ", -1 , "if positive, select event generator type. i.e. gap or signal "};
223+ Configurable<int > cfgRejectEventGenerator{ " cfgRejectEventGenerator " , 999 , " reject event generator. e.g. reject tracks from gap events " };
224224 Configurable<int > cfgCentEstimator{" cfgCentEstimator" , 2 , " FT0M:0, FT0A:1, FT0C:2" };
225225 Configurable<float > cfgCentMin{" cfgCentMin" , -1 .f , " min. centrality" };
226226 Configurable<float > cfgCentMax{" cfgCentMax" , 999 .f , " max. centrality" };
@@ -1047,10 +1047,6 @@ struct taggingHFE {
10471047 continue ;
10481048 }
10491049 const auto & mcParticle = track.template mcParticle_as <aod::McParticles>();
1050- // const auto& mcCollision = mcParticle.template mcCollision_as<aod::McCollisions>();
1051- // if (eventCut.cfgEventGeneratorType >= 0 && mcCollision.getSubGeneratorId() != eventCut.cfgEventGeneratorType) {
1052- // continue;
1053- // }
10541050 if (!mcParticle.has_mothers ()) {
10551051 continue ;
10561052 }
@@ -1066,6 +1062,11 @@ struct taggingHFE {
10661062 float dcaZ = mDcaInfoCov .getZ ();
10671063
10681064 if (isSelectedTrack (track, trackParCov, dcaXY, dcaZ) && std::abs (mcParticle.pdgCode ()) == cfgPdgLepton && isPrimary) { // keep truely primary electrons
1065+ auto mcCollision_ele = mcParticle.template mcCollision_as <aod::McCollisions>();
1066+ if (mcCollision_ele.getSubGeneratorId () == eventCut.cfgRejectEventGenerator ) {
1067+ continue ;
1068+ }
1069+
10691070 fRegistry .fill (HIST (" Track/Electron/hTPCdEdx" ), track.tpcInnerParam (), track.mcTunedTPCSignal ());
10701071 fRegistry .fill (HIST (" Track/Electron/hTOFbeta" ), track.p (), track.beta ());
10711072 if (track.sign () > 0 ) { // positron
0 commit comments