Skip to content

Commit 9680383

Browse files
lgansbartlLaura Gansbartl
andauthored
[PWGEM] add option for at least one photon in Dalitz Skimmer (#17996)
Co-authored-by: Laura Gansbartl <lgansbartl@pcikf5.ikf.uni-frankfurt.de>
1 parent 1d65360 commit 9680383

1 file changed

Lines changed: 9 additions & 0 deletions

File tree

PWGEM/PhotonMeson/TableProducer/skimmerPrimaryElectronFromDalitzEE.h

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -166,6 +166,7 @@ struct skimmerPrimaryElectronFromDalitzEE {
166166
Configurable<float> maxMee{"maxMee", 0.5, "max. mee to store dalitz ee pairs"};
167167
Configurable<float> minMeegamma{"minMeegamma", 0.3, "min. mee to store eegamma candidates"};
168168
Configurable<float> maxMeegamma{"maxMeegamma", 0.8, "max. mee to store eegamma candidates"};
169+
Configurable<bool> photonsNeeded{"photonsNeeded", true, "require a minimum of one photons per collision"};
169170
Configurable<bool> fillLS{"fillLS", true, "flag to fill LS histograms for QA"};
170171
Configurable<int> fillingMode{"fillingMode", 1, "Filling mode| 1: fill tracks without pair selection, 2: fill tracks from selected pairs, 3: fill tracks from selected pairs that can be combined with a photon"};
171172
Configurable<bool> fillWithEtaMassCut{"fillWithEtaMassCut", true, "only valid for fillingmode 3; true: filling tabled based on eta candidate selection with minMeegamma < M < maxMeegamma, false: fill identical to fillingmode 2"};
@@ -745,6 +746,10 @@ struct skimmerPrimaryElectronFromDalitzEE {
745746
auto varFillingMode = static_cast<enumFillingMode>(fillingMode.value);
746747

747748
const auto& v0photons_per_coll = v0photons.sliceBy(perCol_pcm, collision.globalIndex());
749+
if (photonsNeeded && v0photons_per_coll.size() == 0) {
750+
continue;
751+
}
752+
748753
const auto& posTracks_per_coll = posTracks->sliceByCached(o2::aod::track::collisionId, collision.globalIndex(), cache);
749754
const auto& negTracks_per_coll = negTracks->sliceByCached(o2::aod::track::collisionId, collision.globalIndex(), cache);
750755
const auto& slicedTracks = tracks.sliceBy(perTracksCollision, collision.globalIndex());
@@ -927,6 +932,10 @@ struct skimmerPrimaryElectronFromDalitzEE {
927932
auto varFillingMode = static_cast<enumFillingMode>(fillingMode.value);
928933

929934
const auto& v0photons_per_coll = v0photons.sliceBy(perCol_pcm, collision.globalIndex());
935+
if (photonsNeeded && v0photons_per_coll.size() == 0) {
936+
continue;
937+
}
938+
930939
const auto& posTracks_per_coll = posTracksMC->sliceByCached(o2::aod::track::collisionId, collision.globalIndex(), cache);
931940
const auto& negTracks_per_coll = negTracksMC->sliceByCached(o2::aod::track::collisionId, collision.globalIndex(), cache);
932941
const auto& slicedTracks = tracks.sliceBy(perTracksCollision, collision.globalIndex());

0 commit comments

Comments
 (0)