Skip to content
Merged
Changes from 2 commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 6 additions & 2 deletions EventFiltering/PWGLF/strangenessFilter.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -146,6 +146,7 @@
Configurable<bool> hastof{"hastof", 1, "Has TOF (OOB condition)"};
Configurable<float> ptthrtof{"ptthrtof", 1.0, "Pt threshold to apply TOF condition"};
Configurable<bool> sel8{"sel8", 0, "Apply sel8 event selection"};
Configurable<bool> isTriggerTVX{"isTriggerTVX", 1, "Require TVX"};
Configurable<int> LowLimitFT0MMult{"LowLimitFT0MMult", 3100, "FT0M selection for omega + high multiplicity trigger"};
Configurable<int> LowLimitFT0MMultNorm{"LowLimitFT0MMultNorm", 70, "FT0M selection for omega + high multiplicity trigger with Normalised FT0M"};
Configurable<bool> useNormalisedMult{"useNormalisedMult", 1, "Use avarage multiplicity for HM omega like in multFilter.cxx"};
Expand Down Expand Up @@ -551,7 +552,7 @@

float getMassWindow(const stfilter::species s, const float pt, const float nsigma = 6)
{
const auto sigma = parSigmaMass->get(0u, s) * exp(parSigmaMass->get(1, s) * pt) + parSigmaMass->get(2, s) * exp(parSigmaMass->get(3, s) * pt);

Check failure on line 555 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
return nsigma * sigma;
}

Expand All @@ -574,12 +575,15 @@

initCCDB(collision.bc().runNumber());

hProcessedEvents->Fill(-0.5);
if (sel8 && !collision.sel8()) {
fillTriggerTable(keepEvent);
return;
}
hProcessedEvents->Fill(-0.5);

if (isTriggerTVX && !collision.selection_bit(aod::evsel::kIsTriggerTVX)){
fillTriggerTable(keepEvent);
return;
}
if (isTimeFrameBorderCut && !collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) {
fillTriggerTable(keepEvent);
return;
Expand Down Expand Up @@ -686,7 +690,7 @@
o2::dataformats::DCA impactParameterTrk;

std::vector<std::tuple<int64_t, int64_t, TVector3, TVector3>> v0sSelTuple;
for (auto& v00 : v0Base) { // loop over v0 for pre selection

Check failure on line 693 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
hCandidate->Fill(0.5); // All candidates

if (v00.v0Type() != 1) {
Expand Down Expand Up @@ -725,7 +729,7 @@
if (std::hypot(mStraHelper.v0.momentum[0], mStraHelper.v0.momentum[1]) < cfgLLCuts.cfgV0PtMin) {
continue;
}
double yLambda = RecoDecay::y(array{mStraHelper.v0.momentum[0], mStraHelper.v0.momentum[1], mStraHelper.v0.momentum[2]}, o2::constants::physics::MassLambda0);

Check failure on line 732 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
if (yLambda < cfgLLCuts.cfgV0RapMin) {
continue;
}
Expand Down Expand Up @@ -781,7 +785,7 @@
}
}

for (auto& casc : cascadesBase) { // loop over cascades

Check failure on line 788 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
hCandidate->Fill(0.5); // All candidates

const auto bachTrack = casc.bachelor_as<TrackCandidates>();
Expand Down Expand Up @@ -818,7 +822,7 @@
// Rapidity
double etaCasc = RecoDecay::eta(std::array{mStraHelper.cascade.cascadeMomentum[0], mStraHelper.cascade.cascadeMomentum[1], mStraHelper.cascade.cascadeMomentum[2]});
// pointing angle
double v0DauCPA = RecoDecay::cpa(pvPos, array{mStraHelper.cascade.v0Position[0], mStraHelper.cascade.v0Position[1], mStraHelper.cascade.v0Position[2]}, array{mStraHelper.cascade.positiveMomentum[0] + mStraHelper.cascade.negativeMomentum[0], mStraHelper.cascade.positiveMomentum[1] + mStraHelper.cascade.negativeMomentum[1], mStraHelper.cascade.positiveMomentum[2] + mStraHelper.cascade.negativeMomentum[2]});

Check failure on line 825 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
double cascCPA = RecoDecay::cpa(
pvPos,
array{mStraHelper.cascade.cascadePosition[0], mStraHelper.cascade.cascadePosition[1], mStraHelper.cascade.cascadePosition[2]},
Expand All @@ -833,14 +837,14 @@
// massLambda
double LambdaMass = 0;
if (mStraHelper.cascade.charge < 0) {
LambdaMass = RecoDecay::m(array{array{mStraHelper.cascade.positiveMomentum[0], mStraHelper.cascade.positiveMomentum[1], mStraHelper.cascade.positiveMomentum[2]}, array{mStraHelper.cascade.negativeMomentum[0], mStraHelper.cascade.negativeMomentum[1], mStraHelper.cascade.negativeMomentum[2]}}, array{o2::constants::physics::MassProton, o2::constants::physics::MassPionCharged});

Check failure on line 840 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
} else {
LambdaMass = RecoDecay::m(array{array{mStraHelper.cascade.positiveMomentum[0], mStraHelper.cascade.positiveMomentum[1], mStraHelper.cascade.positiveMomentum[2]}, array{mStraHelper.cascade.negativeMomentum[0], mStraHelper.cascade.negativeMomentum[1], mStraHelper.cascade.negativeMomentum[2]}}, array{o2::constants::physics::MassPionCharged, o2::constants::physics::MassProton});

Check failure on line 842 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
}

// rapidity
double yXi = RecoDecay::y(array{mStraHelper.cascade.bachelorMomentum[0] + mStraHelper.cascade.positiveMomentum[0] + mStraHelper.cascade.negativeMomentum[0], mStraHelper.cascade.bachelorMomentum[1] + mStraHelper.cascade.positiveMomentum[1] + mStraHelper.cascade.negativeMomentum[1], mStraHelper.cascade.bachelorMomentum[2] + mStraHelper.cascade.positiveMomentum[2] + mStraHelper.cascade.negativeMomentum[2]}, o2::constants::physics::MassXiMinus);

Check failure on line 846 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
double yOmega = RecoDecay::y(array{mStraHelper.cascade.bachelorMomentum[0] + mStraHelper.cascade.positiveMomentum[0] + mStraHelper.cascade.negativeMomentum[0], mStraHelper.cascade.bachelorMomentum[1] + mStraHelper.cascade.positiveMomentum[1] + mStraHelper.cascade.negativeMomentum[1], mStraHelper.cascade.bachelorMomentum[2] + mStraHelper.cascade.positiveMomentum[2] + mStraHelper.cascade.negativeMomentum[2]}, o2::constants::physics::MassOmegaMinus);

Check failure on line 847 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.

if (mStraHelper.cascade.charge > 0) {
if (std::fabs(mStraHelper.cascade.positiveDCAxy) < dcamesontopv) {
Expand Down Expand Up @@ -1069,7 +1073,7 @@

// QA tracks
int triggcounterAllEv = 0;
for (auto track : tracks) { // start loop over tracks

Check failure on line 1076 in EventFiltering/PWGLF/strangenessFilter.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop]

Use constant references for non-modified iterators in range-based for loops.
if (isTrackFilter && !selectTrack(track)) {
continue;
}
Expand Down
Loading