Skip to content

Commit 4664832

Browse files
sawankumawatSawan Sawan
sawankumawat
and
Sawan Sawan
authored
[PWGLF] Applied rapidity cut for rotational mother (#12029)
Co-authored-by: Sawan Sawan <sawan.sawan@cern.ch>
1 parent 91f8904 commit 4664832

File tree

1 file changed

+169
-102
lines changed

1 file changed

+169
-102
lines changed

‎PWGLF/Tasks/Resonances/kstarqa.cxx‎

Lines changed: 169 additions & 102 deletions
Original file line numberDiff line numberDiff line change
@@ -119,8 +119,8 @@ struct Kstarqa {
119119
Configurable<float> nsigmaCutCombined{"nsigmaCutCombined", 3.0, "Combined Nsigma cut"};
120120

121121
// Event selection configurables
122-
Configurable<bool> timFrameEvsel{"timFrameEvsel", true, "TPC Time frame boundary cut"};
123-
Configurable<bool> cTVXEvsel{"cTVXEvsel", true, "Triggger selection"};
122+
// Configurable<bool> timFrameEvsel{"timFrameEvsel", true, "TPC Time frame boundary cut"};
123+
// Configurable<bool> cTVXEvsel{"cTVXEvsel", true, "Triggger selection"};
124124
Configurable<float> cutzvertex{"cutzvertex", 10.0f, "Accepted z-vertex range (cm)"};
125125
// Configurable<bool> cMID{"cMID", false, "Misidentification of tracks"};
126126

@@ -213,6 +213,8 @@ struct Kstarqa {
213213
rEventSelection.add("events_check", "No. of events in the generated MC", kTH1I, {{20, 0, 20}});
214214
rEventSelection.add("events_checkrec", "No. of events in the reconstructed MC", kTH1I, {{20, 0, 20}});
215215
hInvMass.add("h1KSRecsplit", "KS meson Rec split", kTH1F, {{100, 0.0f, 10.0f}});
216+
hInvMass.add("kstargenBeforeEvtSel", "Kstar generated before event selection", kTH1F, {ptAxis});
217+
hInvMass.add("kstargenAfterEvtSel", "Kstar generated after event selection", kTH1F, {ptAxis});
216218

217219
// Multplicity distribution
218220
if (cQAevents) {
@@ -229,6 +231,25 @@ struct Kstarqa {
229231
double massPi = o2::constants::physics::MassPiPlus;
230232
double massKa = o2::constants::physics::MassKPlus;
231233

234+
template <typename Coll>
235+
bool selectionEvent(const Coll& collision)
236+
{
237+
if (std::abs(collision.posZ()) > cutzvertex)
238+
return false;
239+
if (!collision.sel8())
240+
return false;
241+
if (!collision.selection_bit(aod::evsel::kIsTriggerTVX)) {
242+
return false;
243+
}
244+
if (!collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kNoITSROFrameBorder)) {
245+
return false;
246+
}
247+
if (rctCut.requireRCTFlagChecker && !rctChecker(collision)) {
248+
return false;
249+
}
250+
return true;
251+
}
252+
232253
template <typename T>
233254
bool selectionTrack(const T& candidate)
234255
{
@@ -451,7 +472,7 @@ struct Kstarqa {
451472

452473
using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::MultZeqs, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs>>;
453474
using TrackCandidates = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTOFFullPi, aod::pidTPCFullKa, aod::pidTOFFullKa, aod::pidTOFbeta>>;
454-
using EventCandidatesMC = soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::CentFT0Ms>;
475+
using EventCandidatesMC = soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs>;
455476

456477
using TrackCandidatesMC = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullKa, aod::pidTOFFullKa, aod::pidTPCFullPi, aod::pidTOFFullPi, aod::McTrackLabels, aod::pidTOFbeta>>;
457478

@@ -468,114 +489,122 @@ struct Kstarqa {
468489
ROOT::Math::Boost boost{mother.BoostToCM()}; // boost mother to center of mass frame
469490
fourVecDauCM = boost(daughterSelected); // boost the frame of daughter same as mother
470491

471-
if (std::abs(mother.Rapidity()) < 0.5) {
472-
if (activateTHnSparseCosThStarHelicity) {
473-
auto cosThetaStarHelicity = mother.Vect().Dot(fourVecDauCM.Vect()) / (std::sqrt(fourVecDauCM.Vect().Mag2()) * std::sqrt(mother.Vect().Mag2()));
492+
// if (std::abs(mother.Rapidity()) < 0.5) {
493+
if (activateTHnSparseCosThStarHelicity) {
494+
auto cosThetaStarHelicity = mother.Vect().Dot(fourVecDauCM.Vect()) / (std::sqrt(fourVecDauCM.Vect().Mag2()) * std::sqrt(mother.Vect().Mag2()));
474495

475-
if (track1.sign() * track2.sign() < 0) {
476-
if (!isMix) {
496+
if (track1.sign() * track2.sign() < 0) {
497+
if (!isMix) {
498+
if (std::abs(mother.Rapidity()) < 0.5) {
477499
hInvMass.fill(HIST("h3KstarInvMassUnlikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarHelicity);
500+
}
478501

479-
for (int i = 0; i < cRotations; i++) {
480-
theta2 = rn->Uniform(o2::constants::math::PI - o2::constants::math::PI / rotationalCut, o2::constants::math::PI + o2::constants::math::PI / rotationalCut);
502+
for (int i = 0; i < cRotations; i++) {
503+
theta2 = rn->Uniform(o2::constants::math::PI - o2::constants::math::PI / rotationalCut, o2::constants::math::PI + o2::constants::math::PI / rotationalCut);
481504

482-
daughterRot = ROOT::Math::PxPyPzMVector(daughter1.Px() * std::cos(theta2) - daughter1.Py() * std::sin(theta2), daughter1.Px() * std::sin(theta2) + daughter1.Py() * std::cos(theta2), daughter1.Pz(), daughter1.M());
505+
daughterRot = ROOT::Math::PxPyPzMVector(daughter1.Px() * std::cos(theta2) - daughter1.Py() * std::sin(theta2), daughter1.Px() * std::sin(theta2) + daughter1.Py() * std::cos(theta2), daughter1.Pz(), daughter1.M());
483506

484-
motherRot = daughterRot + daughter2;
507+
motherRot = daughterRot + daughter2;
485508

486-
ROOT::Math::Boost boost2{motherRot.BoostToCM()};
487-
daughterRotCM = boost2(daughterRot);
509+
ROOT::Math::Boost boost2{motherRot.BoostToCM()};
510+
daughterRotCM = boost2(daughterRot);
488511

489-
auto cosThetaStarHelicityRot = motherRot.Vect().Dot(daughterRotCM.Vect()) / (std::sqrt(daughterRotCM.Vect().Mag2()) * std::sqrt(motherRot.Vect().Mag2()));
512+
auto cosThetaStarHelicityRot = motherRot.Vect().Dot(daughterRotCM.Vect()) / (std::sqrt(daughterRotCM.Vect().Mag2()) * std::sqrt(motherRot.Vect().Mag2()));
490513

491-
if (calcRotational)
492-
hInvMass.fill(HIST("h3KstarInvMassRotated"), multiplicity, motherRot.Pt(), motherRot.M(), cosThetaStarHelicityRot);
493-
}
494-
} else {
495-
hInvMass.fill(HIST("h3KstarInvMassMixed"), multiplicity, mother.Pt(), mother.M(), cosThetaStarHelicity);
496-
}
497-
} else {
498-
if (!isMix) {
499-
if (calcLikeSign)
500-
hInvMass.fill(HIST("h3KstarInvMasslikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarHelicity);
514+
if (calcRotational && motherRot.Rapidity() < 0.5)
515+
hInvMass.fill(HIST("h3KstarInvMassRotated"), multiplicity, motherRot.Pt(), motherRot.M(), cosThetaStarHelicityRot);
501516
}
517+
} else if (isMix && std::abs(mother.Rapidity()) < 0.5) {
518+
hInvMass.fill(HIST("h3KstarInvMassMixed"), multiplicity, mother.Pt(), mother.M(), cosThetaStarHelicity);
519+
}
520+
} else {
521+
if (!isMix) {
522+
if (calcLikeSign && std::abs(mother.Rapidity()) < 0.5)
523+
hInvMass.fill(HIST("h3KstarInvMasslikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarHelicity);
502524
}
525+
}
503526

504-
} else if (activateTHnSparseCosThStarProduction) {
505-
normalVec = ROOT::Math::XYZVector(mother.Py(), -mother.Px(), 0.f);
506-
auto cosThetaStarProduction = normalVec.Dot(fourVecDauCM.Vect()) / (std::sqrt(fourVecDauCM.Vect().Mag2()) * std::sqrt(normalVec.Mag2()));
527+
} else if (activateTHnSparseCosThStarProduction) {
528+
normalVec = ROOT::Math::XYZVector(mother.Py(), -mother.Px(), 0.f);
529+
auto cosThetaStarProduction = normalVec.Dot(fourVecDauCM.Vect()) / (std::sqrt(fourVecDauCM.Vect().Mag2()) * std::sqrt(normalVec.Mag2()));
507530

508-
if (track1.sign() * track2.sign() < 0) {
509-
if (!isMix) {
531+
if (track1.sign() * track2.sign() < 0) {
532+
if (!isMix) {
533+
if (std::abs(mother.Rapidity()) < 0.5) {
510534
hInvMass.fill(HIST("h3KstarInvMassUnlikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarProduction);
511-
for (int i = 0; i < cRotations; i++) {
512-
theta2 = rn->Uniform(0, o2::constants::math::PI);
513-
daughterRot = ROOT::Math::PxPyPzMVector(daughter1.Px() * std::cos(theta2) - daughter1.Py() * std::sin(theta2), daughter1.Px() * std::sin(theta2) + daughter1.Py() * std::cos(theta2), daughter1.Pz(), daughter1.M());
514-
515-
motherRot = daughterRot + daughter2;
516-
if (calcRotational)
517-
hInvMass.fill(HIST("h3KstarInvMassRotated"), multiplicity, motherRot.Pt(), motherRot.M(), cosThetaStarProduction);
518-
}
519-
} else {
520-
hInvMass.fill(HIST("h3KstarInvMassMixed"), multiplicity, mother.Pt(), mother.M(), cosThetaStarProduction);
521535
}
522-
} else {
523-
if (!isMix) {
524-
if (calcLikeSign)
525-
hInvMass.fill(HIST("h3KstarInvMasslikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarProduction);
536+
for (int i = 0; i < cRotations; i++) {
537+
theta2 = rn->Uniform(0, o2::constants::math::PI);
538+
daughterRot = ROOT::Math::PxPyPzMVector(daughter1.Px() * std::cos(theta2) - daughter1.Py() * std::sin(theta2), daughter1.Px() * std::sin(theta2) + daughter1.Py() * std::cos(theta2), daughter1.Pz(), daughter1.M());
539+
540+
motherRot = daughterRot + daughter2;
541+
if (calcRotational && abs(motherRot.Rapidity()) < 0.5)
542+
hInvMass.fill(HIST("h3KstarInvMassRotated"), multiplicity, motherRot.Pt(), motherRot.M(), cosThetaStarProduction);
526543
}
544+
} else if (isMix && std::abs(mother.Rapidity()) < 0.5) {
545+
hInvMass.fill(HIST("h3KstarInvMassMixed"), multiplicity, mother.Pt(), mother.M(), cosThetaStarProduction);
527546
}
528-
} else if (activateTHnSparseCosThStarBeam) {
529-
beamVec = ROOT::Math::XYZVector(0.f, 0.f, 1.f);
530-
auto cosThetaStarBeam = beamVec.Dot(fourVecDauCM.Vect()) / std::sqrt(fourVecDauCM.Vect().Mag2());
547+
} else {
548+
if (!isMix) {
549+
if (calcLikeSign && std::abs(mother.Rapidity()) < 0.5)
550+
hInvMass.fill(HIST("h3KstarInvMasslikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarProduction);
551+
}
552+
}
553+
} else if (activateTHnSparseCosThStarBeam) {
554+
beamVec = ROOT::Math::XYZVector(0.f, 0.f, 1.f);
555+
auto cosThetaStarBeam = beamVec.Dot(fourVecDauCM.Vect()) / std::sqrt(fourVecDauCM.Vect().Mag2());
531556

532-
if (track1.sign() * track2.sign() < 0) {
533-
if (!isMix) {
557+
if (track1.sign() * track2.sign() < 0) {
558+
if (!isMix) {
559+
if (std::abs(mother.Rapidity()) < 0.5) {
534560
hInvMass.fill(HIST("h3KstarInvMassUnlikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarBeam);
535-
for (int i = 0; i < cRotations; i++) {
536-
theta2 = rn->Uniform(0, o2::constants::math::PI);
537-
daughterRot = ROOT::Math::PxPyPzMVector(daughter1.Px() * std::cos(theta2) - daughter1.Py() * std::sin(theta2), daughter1.Px() * std::sin(theta2) + daughter1.Py() * std::cos(theta2), daughter1.Pz(), daughter1.M());
561+
}
562+
for (int i = 0; i < cRotations; i++) {
563+
theta2 = rn->Uniform(0, o2::constants::math::PI);
564+
daughterRot = ROOT::Math::PxPyPzMVector(daughter1.Px() * std::cos(theta2) - daughter1.Py() * std::sin(theta2), daughter1.Px() * std::sin(theta2) + daughter1.Py() * std::cos(theta2), daughter1.Pz(), daughter1.M());
538565

539-
motherRot = daughterRot + daughter2;
540-
if (calcRotational)
541-
hInvMass.fill(HIST("h3KstarInvMassRotated"), multiplicity, motherRot.Pt(), motherRot.M(), cosThetaStarBeam);
542-
}
543-
} else {
544-
hInvMass.fill(HIST("h3KstarInvMassMixed"), multiplicity, mother.Pt(), mother.M(), cosThetaStarBeam);
566+
motherRot = daughterRot + daughter2;
567+
if (calcRotational && std::abs(motherRot.Rapidity()) < 0.5)
568+
hInvMass.fill(HIST("h3KstarInvMassRotated"), multiplicity, motherRot.Pt(), motherRot.M(), cosThetaStarBeam);
545569
}
546-
} else {
547-
if (calcLikeSign)
548-
hInvMass.fill(HIST("h3KstarInvMasslikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarBeam);
570+
} else if (isMix && std::abs(mother.Rapidity()) < 0.5) {
571+
hInvMass.fill(HIST("h3KstarInvMassMixed"), multiplicity, mother.Pt(), mother.M(), cosThetaStarBeam);
549572
}
550-
} else if (activateTHnSparseCosThStarRandom) {
551-
auto phiRandom = gRandom->Uniform(0.f, constants::math::TwoPI);
552-
auto thetaRandom = gRandom->Uniform(0.f, constants::math::PI);
573+
} else {
574+
if (calcLikeSign && std::abs(mother.Rapidity()) < 0.5)
575+
hInvMass.fill(HIST("h3KstarInvMasslikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarBeam);
576+
}
577+
} else if (activateTHnSparseCosThStarRandom) {
578+
auto phiRandom = gRandom->Uniform(0.f, constants::math::TwoPI);
579+
auto thetaRandom = gRandom->Uniform(0.f, constants::math::PI);
553580

554-
randomVec = ROOT::Math::XYZVector(std::sin(thetaRandom) * std::cos(phiRandom), std::sin(thetaRandom) * std::sin(phiRandom), std::cos(thetaRandom));
555-
auto cosThetaStarRandom = randomVec.Dot(fourVecDauCM.Vect()) / std::sqrt(fourVecDauCM.Vect().Mag2());
581+
randomVec = ROOT::Math::XYZVector(std::sin(thetaRandom) * std::cos(phiRandom), std::sin(thetaRandom) * std::sin(phiRandom), std::cos(thetaRandom));
582+
auto cosThetaStarRandom = randomVec.Dot(fourVecDauCM.Vect()) / std::sqrt(fourVecDauCM.Vect().Mag2());
556583

557-
if (track1.sign() * track2.sign() < 0) {
558-
if (!isMix) {
584+
if (track1.sign() * track2.sign() < 0) {
585+
if (!isMix) {
586+
if (std::abs(mother.Rapidity()) < 0.5) {
559587
hInvMass.fill(HIST("h3KstarInvMassUnlikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarRandom);
560-
for (int i = 0; i < cRotations; i++) {
561-
theta2 = rn->Uniform(0, o2::constants::math::PI);
562-
daughterRot = ROOT::Math::PxPyPzMVector(daughter1.Px() * std::cos(theta2) - daughter1.Py() * std::sin(theta2), daughter1.Px() * std::sin(theta2) + daughter1.Py() * std::cos(theta2), daughter1.Pz(), daughter1.M());
563-
564-
motherRot = daughterRot + daughter2;
565-
if (calcRotational)
566-
hInvMass.fill(HIST("h3KstarInvMassRotated"), multiplicity, motherRot.Pt(), motherRot.M(), cosThetaStarRandom);
567-
}
568-
} else {
569-
hInvMass.fill(HIST("h3KstarInvMassMixed"), multiplicity, mother.Pt(), mother.M(), cosThetaStarRandom);
570588
}
571-
} else {
572-
if (!isMix) {
573-
if (calcLikeSign)
574-
hInvMass.fill(HIST("h3KstarInvMasslikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarRandom);
589+
for (int i = 0; i < cRotations; i++) {
590+
theta2 = rn->Uniform(0, o2::constants::math::PI);
591+
daughterRot = ROOT::Math::PxPyPzMVector(daughter1.Px() * std::cos(theta2) - daughter1.Py() * std::sin(theta2), daughter1.Px() * std::sin(theta2) + daughter1.Py() * std::cos(theta2), daughter1.Pz(), daughter1.M());
592+
593+
motherRot = daughterRot + daughter2;
594+
if (calcRotational && std::abs(motherRot.Rapidity()) < 0.5)
595+
hInvMass.fill(HIST("h3KstarInvMassRotated"), multiplicity, motherRot.Pt(), motherRot.M(), cosThetaStarRandom);
575596
}
597+
} else if (isMix && std::abs(mother.Rapidity()) < 0.5) {
598+
hInvMass.fill(HIST("h3KstarInvMassMixed"), multiplicity, mother.Pt(), mother.M(), cosThetaStarRandom);
599+
}
600+
} else {
601+
if (!isMix) {
602+
if (calcLikeSign && std::abs(mother.Rapidity()) < 0.5)
603+
hInvMass.fill(HIST("h3KstarInvMasslikeSign"), multiplicity, mother.Pt(), mother.M(), cosThetaStarRandom);
576604
}
577605
}
578606
}
607+
// }
579608
}
580609

581610
// int counter = 0;
@@ -585,23 +614,24 @@ struct Kstarqa {
585614
{
586615
rEventSelection.fill(HIST("events_check_data"), 0.5);
587616

588-
if (cTVXEvsel && (!collision.selection_bit(aod::evsel::kIsTriggerTVX))) {
589-
return;
590-
}
591-
rEventSelection.fill(HIST("events_check_data"), 1.5);
617+
// if (cTVXEvsel && (!collision.selection_bit(aod::evsel::kIsTriggerTVX))) {
618+
// return;
619+
// }
620+
// rEventSelection.fill(HIST("events_check_data"), 1.5);
592621

593-
if (timFrameEvsel && (!collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kNoITSROFrameBorder))) {
594-
return;
595-
}
622+
// if (timFrameEvsel && (!collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kNoITSROFrameBorder))) {
623+
// return;
624+
// }
596625

597-
if (!collision.sel8()) {
598-
return;
599-
}
600-
rEventSelection.fill(HIST("events_check_data"), 2.5);
626+
// if (!collision.sel8()) {
627+
// return;
628+
// }
629+
// rEventSelection.fill(HIST("events_check_data"), 2.5);
601630

602-
if (rctCut.requireRCTFlagChecker && !rctChecker(collision)) {
631+
if (!selectionEvent(collision)) {
603632
return;
604633
}
634+
605635
rEventSelection.fill(HIST("events_check_data"), 3.5);
606636

607637
multiplicity = -1;
@@ -613,7 +643,7 @@ struct Kstarqa {
613643
} else if (cSelectMultEstimator == 2) {
614644
multiplicity = collision.centFT0C();
615645
} else {
616-
multiplicity = collision.multFT0M();
646+
multiplicity = collision.centFT0M();
617647
}
618648

619649
// Fill the event counter
@@ -747,11 +777,12 @@ struct Kstarqa {
747777
// Map estimator to pair and multiplicity accessor
748778
auto runMixing = [&](auto& pair, auto multiplicityGetter) {
749779
for (const auto& [c1, tracks1, c2, tracks2] : pair) {
750-
if (!c1.sel8() || !c2.sel8())
751-
continue;
780+
// if (!c1.sel8() || !c2.sel8())
781+
// continue;
752782

753-
if (rctCut.requireRCTFlagChecker && (!rctChecker(c1) || !rctChecker(c2)))
783+
if (!selectionEvent(c1) || !selectionEvent(c2)) {
754784
continue;
785+
}
755786

756787
multiplicity = multiplicityGetter(c1);
757788

@@ -815,10 +846,10 @@ struct Kstarqa {
815846
continue;
816847
}
817848

818-
if (timFrameEvsel && !collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) {
849+
if (!collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) {
819850
continue;
820851
}
821-
if (cTVXEvsel && (!collision.selection_bit(aod::evsel::kIsTriggerTVX))) {
852+
if (!collision.selection_bit(aod::evsel::kIsTriggerTVX)) {
822853
continue;
823854
}
824855
if (!collision.sel8()) {
@@ -882,6 +913,42 @@ struct Kstarqa {
882913
}
883914
PROCESS_SWITCH(Kstarqa, processGen, "Process Generated", false);
884915

916+
void processEvtLossSigLossMC(aod::McCollisions::iterator const&, aod::McParticles const& mcParticles, const soa::SmallGroups<EventCandidatesMC>& recCollisions)
917+
{
918+
919+
bool isSel = false;
920+
// auto multiplicity1 = -999.;
921+
for (const auto& RecCollision : recCollisions) {
922+
if (!selectionEvent(RecCollision))
923+
continue;
924+
925+
// if (cSelectMultEstimator == 0) {
926+
// multiplicity1 = RecCollision.centFT0M();
927+
// } else if (cSelectMultEstimator == 1) {
928+
// multiplicity1 = RecCollision.centFT0A();
929+
// } else if (cSelectMultEstimator == 2) {
930+
// multiplicity1 = RecCollision.centFT0C();
931+
// } else {
932+
// multiplicity1 = RecCollision.centFT0M();
933+
// }
934+
935+
isSel = true;
936+
}
937+
938+
// Generated MC
939+
for (const auto& mcPart : mcParticles) {
940+
if (std::abs(mcPart.y()) >= 0.5 || std::abs(mcPart.pdgCode()) != 313)
941+
continue;
942+
943+
// signal loss estimation
944+
hInvMass.fill(HIST("kstargenBeforeEvtSel"), mcPart.pt());
945+
if (isSel) {
946+
hInvMass.fill(HIST("kstargenAfterEvtSel"), mcPart.pt());
947+
}
948+
} // end loop on gen particles
949+
}
950+
PROCESS_SWITCH(Kstarqa, processEvtLossSigLossMC, "Process Signal Loss, Event Loss", false);
951+
885952
void processRec(EventCandidatesMC::iterator const& collision, TrackCandidatesMC const& tracks, aod::McParticles const&, aod::McCollisions const& /*mcCollisions*/)
886953
{
887954

@@ -898,12 +965,12 @@ struct Kstarqa {
898965
}
899966
rEventSelection.fill(HIST("events_checkrec"), 2.5);
900967

901-
if (timFrameEvsel && !collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) {
968+
if (!collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) {
902969
return;
903970
}
904971
rEventSelection.fill(HIST("events_checkrec"), 3.5);
905972

906-
if (cTVXEvsel && (!collision.selection_bit(aod::evsel::kIsTriggerTVX))) {
973+
if (!collision.selection_bit(aod::evsel::kIsTriggerTVX)) {
907974
return;
908975
}
909976
rEventSelection.fill(HIST("events_checkrec"), 4.5);

0 commit comments

Comments
 (0)