Skip to content

Commit 4a8f353

Browse files
alcalivaalcaliva
committed
removed cut on V0 mass
1 parent 1ed5dec commit 4a8f353

File tree

1 file changed

+22
-20
lines changed

1 file changed

+22
-20
lines changed

‎PWGLF/Tasks/Strangeness/strangenessInJets.cxx‎

Lines changed: 22 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -509,8 +509,8 @@ struct StrangenessInJets {
509509
}
510510

511511
// Xi Selections
512-
template <typename Xi, typename TrackPos, typename TrackNeg, typename TrackBac, typename vZero, typename Coll>
513-
bool passedXiSelection(const Xi& casc, const TrackPos& ptrack, const TrackNeg& ntrack, const TrackBac& btrack, const vZero& v0, const Coll& coll)
512+
template <typename Xi, typename TrackPos, typename TrackNeg, typename TrackBac, typename Coll>
513+
bool passedXiSelection(const Xi& casc, const TrackPos& ptrack, const TrackNeg& ntrack, const TrackBac& btrack, const Coll& coll)
514514
{
515515
// Single-track selections on cascade daughters
516516
if (!passedSingleTrackSelection(ptrack))
@@ -542,8 +542,9 @@ struct StrangenessInJets {
542542
}
543543

544544
// Require that V0 is compatible with Lambda
545+
/*
545546
if (std::fabs(v0.mAntiLambda() - MassLambda0) > deltaMassLambda)
546-
return false;
547+
return false;*/
547548
}
548549

549550
// Xi- selection (Xi- -> L + pi-)
@@ -568,8 +569,9 @@ struct StrangenessInJets {
568569
}
569570

570571
// Require that V0 is compatible with Lambda
572+
/*
571573
if (std::fabs(v0.mLambda() - MassLambda0) > deltaMassLambda)
572-
return false;
574+
return false;*/
573575
}
574576

575577
// V0 selections
@@ -613,8 +615,8 @@ struct StrangenessInJets {
613615
}
614616

615617
// Omega selections
616-
template <typename Omega, typename TrackPos, typename TrackNeg, typename TrackBac, typename vZero, typename Coll>
617-
bool passedOmegaSelection(const Omega& casc, const TrackPos& ptrack, const TrackNeg& ntrack, const TrackBac& btrack, const vZero& v0, const Coll& coll)
618+
template <typename Omega, typename TrackPos, typename TrackNeg, typename TrackBac, typename Coll>
619+
bool passedOmegaSelection(const Omega& casc, const TrackPos& ptrack, const TrackNeg& ntrack, const TrackBac& btrack, const Coll& coll)
618620
{
619621
// Single-track selections on cascade daughters
620622
if (!passedSingleTrackSelection(ptrack))
@@ -646,8 +648,9 @@ struct StrangenessInJets {
646648
}
647649

648650
// Require that V0 is compatible with Lambda
651+
/*
649652
if (std::fabs(v0.mAntiLambda() - MassLambda0) > deltaMassLambda)
650-
return false;
653+
return false;*/
651654
}
652655

653656
// Omega- selection (Omega- -> L + K-)
@@ -672,8 +675,9 @@ struct StrangenessInJets {
672675
}
673676

674677
// Require that V0 is compatible with Lambda
678+
/*
675679
if (std::fabs(v0.mLambda() - MassLambda0) > deltaMassLambda)
676-
return false;
680+
return false;*/
677681
}
678682

679683
// V0 selections
@@ -739,7 +743,7 @@ struct StrangenessInJets {
739743

740744
// Process data
741745
void processData(SelCollisions::iterator const& collision, aod::V0Datas const& fullV0s,
742-
CascAndV0 const& Cascades, DaughterTracks const& tracks,
746+
aod::CascDataExt const& Cascades, DaughterTracks const& tracks,
743747
aod::BCsWithTimestamps const&)
744748
{
745749
// Fill event counter before event selection
@@ -896,7 +900,6 @@ struct StrangenessInJets {
896900
auto bach = casc.bachelor_as<DaughterTracks>();
897901
auto pos = casc.posTrack_as<DaughterTracks>();
898902
auto neg = casc.negTrack_as<DaughterTracks>();
899-
auto const& v0 = casc.v0();
900903
TVector3 cascadeDir(casc.px(), casc.py(), casc.pz());
901904

902905
// Calculate distance from jet and UE axes
@@ -911,7 +914,7 @@ struct StrangenessInJets {
911914
double deltaRue2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
912915

913916
// Xi+
914-
if (passedXiSelection(casc, pos, neg, bach, v0, collision) && bach.sign() > 0) {
917+
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() > 0) {
915918
if (deltaRjet < rJet) {
916919
registryData.fill(HIST("XiPos_in_jet"), multiplicity, casc.pt(), casc.mXi());
917920
}
@@ -920,7 +923,7 @@ struct StrangenessInJets {
920923
}
921924
}
922925
// Xi-
923-
if (passedXiSelection(casc, pos, neg, bach, v0, collision) && bach.sign() < 0) {
926+
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() < 0) {
924927
if (deltaRjet < rJet) {
925928
registryData.fill(HIST("XiNeg_in_jet"), multiplicity, casc.pt(), casc.mXi());
926929
}
@@ -929,7 +932,7 @@ struct StrangenessInJets {
929932
}
930933
}
931934
// Omega+
932-
if (passedOmegaSelection(casc, pos, neg, bach, v0, collision) && bach.sign() > 0) {
935+
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() > 0) {
933936
if (deltaRjet < rJet) {
934937
registryData.fill(HIST("OmegaPos_in_jet"), multiplicity, casc.pt(), casc.mOmega());
935938
}
@@ -938,7 +941,7 @@ struct StrangenessInJets {
938941
}
939942
}
940943
// Omega-
941-
if (passedOmegaSelection(casc, pos, neg, bach, v0, collision) && bach.sign() < 0) {
944+
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() < 0) {
942945
if (deltaRjet < rJet) {
943946
registryData.fill(HIST("OmegaNeg_in_jet"), multiplicity, casc.pt(), casc.mOmega());
944947
}
@@ -1109,7 +1112,7 @@ struct StrangenessInJets {
11091112

11101113
// Reconstructed MC events
11111114
void processMCreconstructed(SimCollisions const& collisions, DaughterTracksMC const& mcTracks,
1112-
aod::V0Datas const& fullV0s, CascAndV0 const& Cascades,
1115+
aod::V0Datas const& fullV0s, aod::CascDataExt const& Cascades,
11131116
const aod::McParticles&)
11141117
{
11151118
for (const auto& collision : collisions) {
@@ -1298,7 +1301,6 @@ struct StrangenessInJets {
12981301
auto bach = casc.bachelor_as<DaughterTracksMC>();
12991302
auto pos = casc.posTrack_as<DaughterTracksMC>();
13001303
auto neg = casc.negTrack_as<DaughterTracksMC>();
1301-
auto const& v0 = casc.v0();
13021304

13031305
// Get MC particles
13041306
if (!bach.has_mcParticle() || !pos.has_mcParticle() || !neg.has_mcParticle())
@@ -1342,7 +1344,7 @@ struct StrangenessInJets {
13421344
double deltaRue2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
13431345

13441346
// Xi+
1345-
if (passedXiSelection(casc, pos, neg, bach, v0, collision) && bach.sign() > 0 && pdgParent == kXiPlusBar) {
1347+
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() > 0 && pdgParent == kXiPlusBar) {
13461348
if (deltaRjet < rJet) {
13471349
registryMC.fill(HIST("XiPos_reconstructed_jet"), multiplicity, casc.pt());
13481350
}
@@ -1351,7 +1353,7 @@ struct StrangenessInJets {
13511353
}
13521354
}
13531355
// Xi-
1354-
if (passedXiSelection(casc, pos, neg, bach, v0, collision) && bach.sign() < 0 && pdgParent == kXiMinus) {
1356+
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() < 0 && pdgParent == kXiMinus) {
13551357
if (deltaRjet < rJet) {
13561358
registryMC.fill(HIST("XiNeg_reconstructed_jet"), multiplicity, casc.pt());
13571359
}
@@ -1360,7 +1362,7 @@ struct StrangenessInJets {
13601362
}
13611363
}
13621364
// Omega+
1363-
if (passedOmegaSelection(casc, pos, neg, bach, v0, collision) && bach.sign() > 0 && pdgParent == kOmegaPlusBar) {
1365+
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() > 0 && pdgParent == kOmegaPlusBar) {
13641366
if (deltaRjet < rJet) {
13651367
registryMC.fill(HIST("OmegaPos_reconstructed_jet"), multiplicity, casc.pt());
13661368
}
@@ -1369,7 +1371,7 @@ struct StrangenessInJets {
13691371
}
13701372
}
13711373
// Omega-
1372-
if (passedOmegaSelection(casc, pos, neg, bach, v0, collision) && bach.sign() < 0 && pdgParent == kOmegaMinus) {
1374+
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() < 0 && pdgParent == kOmegaMinus) {
13731375
if (deltaRjet < rJet) {
13741376
registryMC.fill(HIST("OmegaNeg_reconstructed_jet"), multiplicity, casc.pt());
13751377
}

0 commit comments

Comments
 (0)