Skip to content

Commit 84b7275

Browse files
alcalivaalcaliva
committed
fixed function
1 parent 7018c24 commit 84b7275

File tree

1 file changed

+18
-16
lines changed

1 file changed

+18
-16
lines changed

‎PWGLF/Tasks/Strangeness/strangenessInJets.cxx‎

Lines changed: 18 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -508,8 +508,8 @@ struct StrangenessInJets {
508508
}
509509

510510
// Xi Selections
511-
template <typename Xi, typename TrackPos, typename TrackNeg, typename TrackBac, typename Coll>
512-
bool passedXiSelection(const Xi& casc, const TrackPos& ptrack, const TrackNeg& ntrack, const TrackBac& btrack, const Coll& coll)
511+
template <typename Xi, typename TrackPos, typename TrackNeg, typename TrackBac, typename vZero, typename Coll>
512+
bool passedXiSelection(const Xi& casc, const TrackPos& ptrack, const TrackNeg& ntrack, const TrackBac& btrack, const vZero& v0, const Coll& coll)
513513
{
514514
// Single-track selections on cascade daughters
515515
if (!passedSingleTrackSelection(ptrack))
@@ -541,7 +541,7 @@ struct StrangenessInJets {
541541
}
542542

543543
// Require that V0 is compatible with Lambda
544-
if (std::fabs(casc.v0().mAntiLambda() - MassLambda0) > deltaMassLambda)
544+
if (std::fabs(v0.mAntiLambda() - MassLambda0) > deltaMassLambda)
545545
return false;
546546
}
547547

@@ -567,7 +567,7 @@ struct StrangenessInJets {
567567
}
568568

569569
// Require that V0 is compatible with Lambda
570-
if (std::fabs(casc.v0().mLambda() - MassLambda0) > deltaMassLambda)
570+
if (std::fabs(v0.mLambda() - MassLambda0) > deltaMassLambda)
571571
return false;
572572
}
573573

@@ -612,8 +612,8 @@ struct StrangenessInJets {
612612
}
613613

614614
// Omega selections
615-
template <typename Omega, typename TrackPos, typename TrackNeg, typename TrackBac, typename Coll>
616-
bool passedOmegaSelection(const Omega& casc, const TrackPos& ptrack, const TrackNeg& ntrack, const TrackBac& btrack, const Coll& coll)
615+
template <typename Omega, typename TrackPos, typename TrackNeg, typename TrackBac, typename vZero, typename Coll>
616+
bool passedOmegaSelection(const Omega& casc, const TrackPos& ptrack, const TrackNeg& ntrack, const TrackBac& btrack, const vZero& v0, const Coll& coll)
617617
{
618618
// Single-track selections on cascade daughters
619619
if (!passedSingleTrackSelection(ptrack))
@@ -645,7 +645,7 @@ struct StrangenessInJets {
645645
}
646646

647647
// Require that V0 is compatible with Lambda
648-
if (std::fabs(casc.v0().mAntiLambda() - MassLambda0) > deltaMassLambda)
648+
if (std::fabs(v0.mAntiLambda() - MassLambda0) > deltaMassLambda)
649649
return false;
650650
}
651651

@@ -671,7 +671,7 @@ struct StrangenessInJets {
671671
}
672672

673673
// Require that V0 is compatible with Lambda
674-
if (std::fabs(casc.v0().mLambda() - MassLambda0) > deltaMassLambda)
674+
if (std::fabs(v0.mLambda() - MassLambda0) > deltaMassLambda)
675675
return false;
676676
}
677677

@@ -895,6 +895,7 @@ struct StrangenessInJets {
895895
auto bach = casc.bachelor_as<DaughterTracks>();
896896
auto pos = casc.posTrack_as<DaughterTracks>();
897897
auto neg = casc.negTrack_as<DaughterTracks>();
898+
auto const& v0 = casc.v0();
898899
TVector3 cascadeDir(casc.px(), casc.py(), casc.pz());
899900

900901
// Calculate distance from jet and UE axes
@@ -909,7 +910,7 @@ struct StrangenessInJets {
909910
double deltaRue2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
910911

911912
// Xi+
912-
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() > 0) {
913+
if (passedXiSelection(casc, pos, neg, bach, v0, collision) && bach.sign() > 0) {
913914
if (deltaRjet < rJet) {
914915
registryData.fill(HIST("XiPos_in_jet"), multiplicity, casc.pt(), casc.mXi());
915916
}
@@ -918,7 +919,7 @@ struct StrangenessInJets {
918919
}
919920
}
920921
// Xi-
921-
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() < 0) {
922+
if (passedXiSelection(casc, pos, neg, bach, v0, collision) && bach.sign() < 0) {
922923
if (deltaRjet < rJet) {
923924
registryData.fill(HIST("XiNeg_in_jet"), multiplicity, casc.pt(), casc.mXi());
924925
}
@@ -927,7 +928,7 @@ struct StrangenessInJets {
927928
}
928929
}
929930
// Omega+
930-
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() > 0) {
931+
if (passedOmegaSelection(casc, pos, neg, bach, v0, collision) && bach.sign() > 0) {
931932
if (deltaRjet < rJet) {
932933
registryData.fill(HIST("OmegaPos_in_jet"), multiplicity, casc.pt(), casc.mOmega());
933934
}
@@ -936,7 +937,7 @@ struct StrangenessInJets {
936937
}
937938
}
938939
// Omega-
939-
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() < 0) {
940+
if (passedOmegaSelection(casc, pos, neg, bach, v0, collision) && bach.sign() < 0) {
940941
if (deltaRjet < rJet) {
941942
registryData.fill(HIST("OmegaNeg_in_jet"), multiplicity, casc.pt(), casc.mOmega());
942943
}
@@ -1296,6 +1297,7 @@ struct StrangenessInJets {
12961297
auto bach = casc.bachelor_as<DaughterTracksMC>();
12971298
auto pos = casc.posTrack_as<DaughterTracksMC>();
12981299
auto neg = casc.negTrack_as<DaughterTracksMC>();
1300+
auto const& v0 = casc.v0();
12991301

13001302
// Get MC particles
13011303
if (!bach.has_mcParticle() || !pos.has_mcParticle() || !neg.has_mcParticle())
@@ -1339,7 +1341,7 @@ struct StrangenessInJets {
13391341
double deltaRue2 = std::sqrt(deltaEtaUe2 * deltaEtaUe2 + deltaPhiUe2 * deltaPhiUe2);
13401342

13411343
// Xi+
1342-
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() > 0 && pdgParent == kXiPlusBar) {
1344+
if (passedXiSelection(casc, pos, neg, bach, v0, collision) && bach.sign() > 0 && pdgParent == kXiPlusBar) {
13431345
if (deltaRjet < rJet) {
13441346
registryMC.fill(HIST("XiPos_reconstructed_jet"), multiplicity, casc.pt());
13451347
}
@@ -1348,7 +1350,7 @@ struct StrangenessInJets {
13481350
}
13491351
}
13501352
// Xi-
1351-
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() < 0 && pdgParent == kXiMinus) {
1353+
if (passedXiSelection(casc, pos, neg, bach, v0, collision) && bach.sign() < 0 && pdgParent == kXiMinus) {
13521354
if (deltaRjet < rJet) {
13531355
registryMC.fill(HIST("XiNeg_reconstructed_jet"), multiplicity, casc.pt());
13541356
}
@@ -1357,7 +1359,7 @@ struct StrangenessInJets {
13571359
}
13581360
}
13591361
// Omega+
1360-
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() > 0 && pdgParent == kOmegaPlusBar) {
1362+
if (passedOmegaSelection(casc, pos, neg, bach, v0, collision) && bach.sign() > 0 && pdgParent == kOmegaPlusBar) {
13611363
if (deltaRjet < rJet) {
13621364
registryMC.fill(HIST("OmegaPos_reconstructed_jet"), multiplicity, casc.pt());
13631365
}
@@ -1366,7 +1368,7 @@ struct StrangenessInJets {
13661368
}
13671369
}
13681370
// Omega-
1369-
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() < 0 && pdgParent == kOmegaMinus) {
1371+
if (passedOmegaSelection(casc, pos, neg, bach, v0, collision) && bach.sign() < 0 && pdgParent == kOmegaMinus) {
13701372
if (deltaRjet < rJet) {
13711373
registryMC.fill(HIST("OmegaNeg_reconstructed_jet"), multiplicity, casc.pt());
13721374
}

0 commit comments

Comments
 (0)