Skip to content

PWGLF: change in derivedcascadeanalysis.cxx #5842

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
njacazio merged 4 commits into from
Apr 25, 2024
Merged

PWGLF: change in derivedcascadeanalysis.cxx #5842

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
e00f2dd
momentum thr for TOF usage
Apr 25, 2024
aeaf05f
MegaLinter fix
Apr 25, 2024
d19bebc
MegaLinter fix_2
Apr 25, 2024
9ccf0d8
MegaLinter fix_3
Apr 25, 2024
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
95 changes: 58 additions & 37 deletions PWGLF/Tasks/Strangeness/derivedcascadeanalysis.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,13 +32,8 @@
#include <TH2F.h>
#include <TProfile.h>
#include <TLorentzVector.h>
#include <Math/Vector4D.h>
#include <TPDGCode.h>
#include <TDatabasePDG.h>
#include <cmath>
#include <array>
#include <cstdlib>
#include "Framework/ASoAHelpers.h"

// constants
const float ctauxiPDG = 4.91; // from PDG
Expand Down Expand Up @@ -421,7 +416,7 @@ struct derivedCascadeAnalysis {
}

template <typename TCascade>
bool IsCascadeCandidateAccepted(TCascade casc, int counter, float /*centrality*/)
bool IsCascadeCandidateAccepted(TCascade casc, int counter, float centrality)
{
float cut = masswin;
histos.fill(HIST("hCutValue"), 2, cut);
Expand Down Expand Up @@ -462,17 +457,19 @@ struct derivedCascadeAnalysis {
return false;
}
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (doDCAV0DauCut) {
cut = dcav0dau;
histos.fill(HIST("hCutValue"), 5, cut);
if (casc.dcaV0daughters() > dcav0dau)
return false;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (doCascadeRadiusCut) {
if (doPtDepCascRadiusCut) {
Expand All @@ -491,8 +488,9 @@ struct derivedCascadeAnalysis {
if (casc.cascradius() > maxRadius)
return false;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
counter += 2;
}

if (doV0RadiusCut) {
if (doPtDepV0RadiusCut) {
Expand All @@ -510,8 +508,9 @@ struct derivedCascadeAnalysis {
if (casc.v0radius() > maxV0Radius)
return false;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
counter += 2;
}

cut = lambdaMassWin;
histos.fill(HIST("hCutValue"), 10, cut);
Expand All @@ -526,15 +525,17 @@ struct derivedCascadeAnalysis {
return false;
}
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (doV0CosPaCut) {
if (!IsCosPAAccepted(casc, casc.x(), casc.y(), casc.z(), doPtDepV0CosPaCut, false))
return false;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

cut = rejcomp;
histos.fill(HIST("hCutValue"), 13, cut);
Expand Down Expand Up @@ -562,8 +563,9 @@ struct derivedCascadeAnalysis {
if (casc.sign() < 0 && (TMath::Abs(casc.dcapostopv()) < dcaBaryonToPV || TMath::Abs(casc.dcanegtopv()) < dcaMesonToPV))
return false;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

return true;
}
Expand Down Expand Up @@ -613,17 +615,19 @@ struct derivedCascadeAnalysis {
if (!IsCosPAAccepted(casc, coll.posX(), coll.posY(), coll.posZ(), doPtDepCosPaCut, true))
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

cut = dcaV0ToPV;
histos.fill(HIST("hCutValue"), 17, cut);
if (doDCAV0ToPVCut) {
if (TMath::Abs(casc.dcav0topv(coll.posX(), coll.posY(), coll.posZ())) < dcaV0ToPV)
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (doNTPCSigmaCut) {
if (casc.sign() < 0) {
Expand All @@ -639,8 +643,9 @@ struct derivedCascadeAnalysis {
continue;
histos.fill(HIST("hCandidate"), ++counter);
}
} else
} else {
++counter;
}

if (posExtra.tpcCrossedRows() < mintpccrrows || negExtra.tpcCrossedRows() < mintpccrrows || bachExtra.tpcCrossedRows() < mintpccrrows)
continue;
Expand All @@ -667,27 +672,31 @@ struct derivedCascadeAnalysis {
cut = ctau;
histos.fill(HIST("hCutValue"), 22, cut);

double fullmomentum;

if (posExtra.hasTOF()) {
fullmomentum = TMath::Sqrt(TMath::Power(casc.pxpos(), 2) + TMath::Power(casc.pypos(), 2) + TMath::Power(casc.pzpos(), 2));
if (doNTOFSigmaProtonCut && casc.sign() < 0) {
histos.fill(HIST("hNsigmaTOFProton"), casc.tofNSigmaXiLaPr(), TMath::Sqrt(TMath::Power(casc.pxpos(), 2) + TMath::Power(casc.pypos(), 2) + TMath::Power(casc.pzpos(), 2)), coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPr()) > nsigmatofPr)
histos.fill(HIST("hNsigmaTOFProton"), casc.tofNSigmaXiLaPr(), fullmomentum, coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPr()) > nsigmatofPr && fullmomentum > 0.6)
continue;
}
if (doNTOFSigmaV0PionCut && casc.sign() > 0) {
histos.fill(HIST("hNsigmaTOFV0Pion"), casc.tofNSigmaXiLaPi(), TMath::Sqrt(TMath::Power(casc.pxpos(), 2) + TMath::Power(casc.pypos(), 2) + TMath::Power(casc.pzpos(), 2)), coll.centFT0C());
histos.fill(HIST("hNsigmaTOFV0Pion"), casc.tofNSigmaXiLaPi(), fullmomentum, coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPi()) > nsigmatofPion)
continue;
}
}

if (negExtra.hasTOF()) {
fullmomentum = TMath::Sqrt(TMath::Power(casc.pxneg(), 2) + TMath::Power(casc.pyneg(), 2) + TMath::Power(casc.pzneg(), 2));
if (doNTOFSigmaProtonCut && casc.sign() > 0) {
histos.fill(HIST("hNsigmaTOFProton"), casc.tofNSigmaXiLaPr(), TMath::Sqrt(TMath::Power(casc.pxneg(), 2) + TMath::Power(casc.pyneg(), 2) + TMath::Power(casc.pzneg(), 2)), coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPr()) > nsigmatofPr)
histos.fill(HIST("hNsigmaTOFProton"), casc.tofNSigmaXiLaPr(), fullmomentum, coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPr()) > nsigmatofPr && fullmomentum > 0.6)
continue;
}
if (doNTOFSigmaV0PionCut && casc.sign() < 0) {
histos.fill(HIST("hNsigmaTOFV0Pion"), casc.tofNSigmaXiLaPi(), TMath::Sqrt(TMath::Power(casc.pxneg(), 2) + TMath::Power(casc.pyneg(), 2) + TMath::Power(casc.pzneg(), 2)), coll.centFT0C());
histos.fill(HIST("hNsigmaTOFV0Pion"), casc.tofNSigmaXiLaPi(), fullmomentum, coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPi()) > nsigmatofPion)
continue;
}
Expand All @@ -702,8 +711,9 @@ struct derivedCascadeAnalysis {
if (TMath::Abs(bachExtra.tpcNSigmaPi()) > nsigmatpcPi)
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (bachExtra.hasTOF() && doNTOFSigmaBachelorCut) {
histos.fill(HIST("hNsigmaTOFBachelorPion"), casc.tofNSigmaXiPi(), TMath::Sqrt(TMath::Power(casc.pxbach(), 2) + TMath::Power(casc.pybach(), 2) + TMath::Power(casc.pzbach(), 2)), coll.centFT0C());
Expand All @@ -725,8 +735,9 @@ struct derivedCascadeAnalysis {
if (TMath::Abs(bachExtra.tpcNSigmaKa()) > nsigmatpcKa)
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (bachExtra.hasTOF() && doNTOFSigmaBachelorCut) {
histos.fill(HIST("hNsigmaTOFBachelorKaon"), casc.tofNSigmaOmKa(), TMath::Sqrt(TMath::Power(casc.pxbach(), 2) + TMath::Power(casc.pybach(), 2) + TMath::Power(casc.pzbach(), 2)), coll.centFT0C());
Expand All @@ -739,8 +750,9 @@ struct derivedCascadeAnalysis {
if (ctau > proplifetime)
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}
}

if (casc.sign() < 0) {
Expand Down Expand Up @@ -839,15 +851,17 @@ struct derivedCascadeAnalysis {
if (!IsCosPAAccepted(casc, coll.posX(), coll.posY(), coll.posZ(), doPtDepCosPaCut, true))
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (doDCAV0ToPVCut) {
if (TMath::Abs(casc.dcav0topv(coll.posX(), coll.posY(), coll.posZ())) < dcaV0ToPV)
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (doNTPCSigmaCut) {
if (casc.sign() < 0) {
Expand All @@ -863,8 +877,9 @@ struct derivedCascadeAnalysis {
continue;
histos.fill(HIST("hCandidate"), ++counter);
}
} else
} else {
++counter;
}

if (posExtra.tpcCrossedRows() < mintpccrrows || negExtra.tpcCrossedRows() < mintpccrrows || bachExtra.tpcCrossedRows() < mintpccrrows)
continue;
Expand All @@ -888,28 +903,31 @@ struct derivedCascadeAnalysis {
float cascpos = std::hypot(casc.x() - coll.posX(), casc.y() - coll.posY(), casc.z() - coll.posZ());
float cascptotmom = std::hypot(casc.px(), casc.py(), casc.pz());
float ctau = -10;
double fullmomentum;

if (posExtra.hasTOF()) {
fullmomentum = TMath::Sqrt(TMath::Power(casc.pxpos(), 2) + TMath::Power(casc.pypos(), 2) + TMath::Power(casc.pzpos(), 2));
if (doNTOFSigmaProtonCut && casc.sign() < 0) {
histos.fill(HIST("hNsigmaTOFProton"), casc.tofNSigmaXiLaPr(), TMath::Sqrt(TMath::Power(casc.pxpos(), 2) + TMath::Power(casc.pypos(), 2) + TMath::Power(casc.pzpos(), 2)), coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPr()) > nsigmatofPr)
histos.fill(HIST("hNsigmaTOFProton"), casc.tofNSigmaXiLaPr(), fullmomentum, coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPr()) > nsigmatofPr && fullmomentum > 0.6)
continue;
}
if (doNTOFSigmaV0PionCut && casc.sign() > 0) {
histos.fill(HIST("hNsigmaTOFV0Pion"), casc.tofNSigmaXiLaPi(), TMath::Sqrt(TMath::Power(casc.pxpos(), 2) + TMath::Power(casc.pypos(), 2) + TMath::Power(casc.pzpos(), 2)), coll.centFT0C());
histos.fill(HIST("hNsigmaTOFV0Pion"), casc.tofNSigmaXiLaPi(), fullmomentum, coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPi()) > nsigmatofPion)
continue;
}
}

if (negExtra.hasTOF()) {
fullmomentum = TMath::Sqrt(TMath::Power(casc.pxneg(), 2) + TMath::Power(casc.pyneg(), 2) + TMath::Power(casc.pzneg(), 2));
if (doNTOFSigmaProtonCut && casc.sign() > 0) {
histos.fill(HIST("hNsigmaTOFProton"), casc.tofNSigmaXiLaPr(), TMath::Sqrt(TMath::Power(casc.pxneg(), 2) + TMath::Power(casc.pyneg(), 2) + TMath::Power(casc.pzneg(), 2)), coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPr()) > nsigmatofPr)
histos.fill(HIST("hNsigmaTOFProton"), casc.tofNSigmaXiLaPr(), fullmomentum, coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPr()) > nsigmatofPr && fullmomentum > 0.6)
continue;
}
if (doNTOFSigmaV0PionCut && casc.sign() < 0) {
histos.fill(HIST("hNsigmaTOFV0Pion"), casc.tofNSigmaXiLaPi(), TMath::Sqrt(TMath::Power(casc.pxneg(), 2) + TMath::Power(casc.pyneg(), 2) + TMath::Power(casc.pzneg(), 2)), coll.centFT0C());
histos.fill(HIST("hNsigmaTOFV0Pion"), casc.tofNSigmaXiLaPi(), fullmomentum, coll.centFT0C());
if (TMath::Abs(casc.tofNSigmaXiLaPi()) > nsigmatofPion)
continue;
}
Expand All @@ -924,8 +942,9 @@ struct derivedCascadeAnalysis {
if (TMath::Abs(bachExtra.tpcNSigmaPi()) > nsigmatpcPi)
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (bachExtra.hasTOF() && doNTOFSigmaBachelorCut) {
histos.fill(HIST("hNsigmaTOFBachelorPion"), casc.tofNSigmaXiPi(), TMath::Sqrt(TMath::Power(casc.pxbach(), 2) + TMath::Power(casc.pybach(), 2) + TMath::Power(casc.pzbach(), 2)), coll.centFT0C());
Expand All @@ -949,8 +968,9 @@ struct derivedCascadeAnalysis {
if (TMath::Abs(bachExtra.tpcNSigmaKa()) > nsigmatpcKa)
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}

if (bachExtra.hasTOF() && doNTOFSigmaBachelorCut) {
histos.fill(HIST("hNsigmaTOFBachelorKaon"), casc.tofNSigmaOmKa(), TMath::Sqrt(TMath::Power(casc.pxbach(), 2) + TMath::Power(casc.pybach(), 2) + TMath::Power(casc.pzbach(), 2)), coll.centFT0C());
Expand All @@ -963,8 +983,9 @@ struct derivedCascadeAnalysis {
if (ctau > proplifetime)
continue;
histos.fill(HIST("hCandidate"), ++counter);
} else
} else {
++counter;
}
invmass = casc.mOmega();
}

Expand Down