@@ -89,6 +89,8 @@ struct Kstarqa {
8989 Configurable<int > cSelectMultEstimator{" cSelectMultEstimator" , 0 , " Select multiplicity estimator: 0 - FT0M, 1 - FT0A, 2 - FT0C" };
9090 Configurable<bool > applyRecMotherRapidity{" applyRecMotherRapidity" , true , " Apply rapidity cut on reconstructed mother track" };
9191 Configurable<bool > applypTdepPID{" applypTdepPID" , false , " Apply pT dependent PID" };
92+ Configurable<bool > ispileupGoodvtxCut{" ispileupGoodvtxCut" , true , " kNoSameBunchPileup, kIsGoodZvtxFT0vsPV cuts" };
93+ Configurable<bool > allLayersGoodITS{" allLayersGoodITS" , true , " Require all ITS layers to be good" };
9294
9395 // Configurables for track selections
9496 Configurable<int > rotationalCut{" rotationalCut" , 10 , " Cut value (Rotation angle pi - pi/cut and pi + pi/cut)" };
@@ -247,6 +249,12 @@ struct Kstarqa {
247249 if (rctCut.requireRCTFlagChecker && !rctChecker (collision)) {
248250 return false ;
249251 }
252+ if (ispileupGoodvtxCut && (!collision.selection_bit (aod::evsel::kNoSameBunchPileup ) || !collision.selection_bit (aod::evsel::kIsGoodZvtxFT0vsPV ))) {
253+ return false ;
254+ }
255+ if (allLayersGoodITS && !collision.selection_bit (o2::aod::evsel::kIsGoodITSLayersAll )) {
256+ return false ;
257+ }
250258 return true ;
251259 }
252260
@@ -855,6 +863,9 @@ struct Kstarqa {
855863 if (!collision.sel8 ()) {
856864 continue ;
857865 }
866+ if (ispileupGoodvtxCut && (!collision.selection_bit (aod::evsel::kNoSameBunchPileup ) || !collision.selection_bit (aod::evsel::kIsGoodZvtxFT0vsPV ))) {
867+ continue ;
868+ }
858869 multiplicity = collision.centFT0M ();
859870 hInvMass.fill (HIST (" h1GenMult" ), multiplicity);
860871 selectedEvents[nevts++] = collision.mcCollision_as <aod::McCollisions>().globalIndex ();
@@ -978,6 +989,9 @@ struct Kstarqa {
978989 if (!collision.sel8 ()) {
979990 return ;
980991 }
992+ if (ispileupGoodvtxCut && (!collision.selection_bit (aod::evsel::kNoSameBunchPileup ) || !collision.selection_bit (aod::evsel::kIsGoodZvtxFT0vsPV ))) {
993+ return ;
994+ }
981995 multiplicity = collision.centFT0M ();
982996 hInvMass.fill (HIST (" h1RecMult" ), multiplicity);
983997
0 commit comments