Skip to content

PWGLF: TOF PID checking hasTOF #5403

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 3 commits into from
Mar 28, 2024
Merged

PWGLF: TOF PID checking hasTOF #5403

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
5c69a1c
PWGLF: TOF PID checking hasTOF
ddobrigk Mar 28, 2024
2f7957a
Please consider the following formatting changes (#272)
alibuild Mar 28, 2024
bdb419d
Change to hasTOF also for cascadePID
ddobrigk Mar 28, 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
14 changes: 7 additions & 7 deletions PWGLF/TableProducer/cascadepid.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -383,19 +383,23 @@ struct cascadepid {
float negDeltaTimeAsOmPi = -1e+6, negDeltaTimeAsOmPr = -1e+6;
float bachDeltaTimeAsOmKa = -1e+6;

if (cascade.posTOFSignal() > 0 && cascade.posTOFEventTime() > 0) {
auto pTra = cascade.posTrackExtra_as<dauTracks>();
auto nTra = cascade.negTrackExtra_as<dauTracks>();
auto bTra = cascade.negTrackExtra_as<dauTracks>();

if (pTra.hasTOF()) {
posDeltaTimeAsXiPi = (cascade.posTOFSignal() - cascade.posTOFEventTime()) - (xiFlight + lambdaFlight + posFlightPi);
posDeltaTimeAsXiPr = (cascade.posTOFSignal() - cascade.posTOFEventTime()) - (xiFlight + lambdaFlight + posFlightPr);
posDeltaTimeAsOmPi = (cascade.posTOFSignal() - cascade.posTOFEventTime()) - (omFlight + lambdaFlight + posFlightPi);
posDeltaTimeAsOmPr = (cascade.posTOFSignal() - cascade.posTOFEventTime()) - (omFlight + lambdaFlight + posFlightPr);
}
if (cascade.negTOFSignal() > 0 && cascade.negTOFEventTime() > 0) {
if (nTra.hasTOF()) {
negDeltaTimeAsXiPi = (cascade.negTOFSignal() - cascade.negTOFEventTime()) - (xiFlight + lambdaFlight + negFlightPi);
negDeltaTimeAsXiPr = (cascade.negTOFSignal() - cascade.negTOFEventTime()) - (xiFlight + lambdaFlight + negFlightPr);
negDeltaTimeAsOmPi = (cascade.negTOFSignal() - cascade.negTOFEventTime()) - (omFlight + lambdaFlight + negFlightPi);
negDeltaTimeAsOmPr = (cascade.negTOFSignal() - cascade.negTOFEventTime()) - (omFlight + lambdaFlight + negFlightPr);
}
if (cascade.bachTOFSignal() > 0 && cascade.bachTOFEventTime() > 0) {
if (bTra.hasTOF()) {
bachDeltaTimeAsXiPi = (cascade.bachTOFSignal() - cascade.bachTOFEventTime()) - (xiFlight + bachFlightPi);
bachDeltaTimeAsOmKa = (cascade.bachTOFSignal() - cascade.bachTOFEventTime()) - (omFlight + bachFlightKa);
}
Expand All @@ -407,10 +411,6 @@ struct cascadepid {
);

if (doQA) {
auto pTra = cascade.posTrackExtra_as<dauTracks>();
auto nTra = cascade.negTrackExtra_as<dauTracks>();
auto bTra = cascade.negTrackExtra_as<dauTracks>();

// fill QA histograms for cross-checking
histos.fill(HIST("hArcDebug"), cascade.pt(), lengthCascade - d3d); // for debugging purposes

Expand Down
44 changes: 35 additions & 9 deletions PWGLF/TableProducer/lambdakzeropid.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -116,6 +116,11 @@ struct lambdakzeropid {
float maxSnp; // max sine phi for propagation
float maxStep; // max step size (cm) for propagation

// enum to keep track of the TOF-related properties for V0s
enum tofEnum { kLength = 0,
kHasTOF,
kNEnums };

/// function to calculate track length of this track up to a certain segment of a detector
/// to be used internally in another funcrtion that calculates length until it finds the proper one
/// warning: this could be optimised further for speed
Expand Down Expand Up @@ -289,6 +294,9 @@ struct lambdakzeropid {
histos.add("h2dDeltaTimePositiveK0ShortPi", "h2dDeltaTimePositiveK0ShortPi", {HistType::kTH3F, {axisPt, axisEta, axisDeltaTime}});
histos.add("h2dDeltaTimeNegativeK0ShortPi", "h2dDeltaTimeNegativeK0ShortPi", {HistType::kTH3F, {axisPt, axisEta, axisDeltaTime}});

histos.add("h2dPositiveTOFProperties", "h2dPositiveTOFProperties", {HistType::kTH2F, {axisPt, {4, -0.5, 3.5f}}});
histos.add("h2dNegativeTOFProperties", "h2dNegativeTOFProperties", {HistType::kTH2F, {axisPt, {4, -0.5, 3.5f}}});

if (doQANSigma) {
// standard NSigma values
histos.add("h2dNSigmaPositiveLambdaPi", "h2dNSigmaPositiveLambdaPi", {HistType::kTH2F, {axisPt, axisNSigma}});
Expand Down Expand Up @@ -411,6 +419,9 @@ struct lambdakzeropid {
o2::track::TrackPar posTrack = o2::track::TrackPar({v0.x(), v0.y(), v0.z()}, {v0.pxpos(), v0.pypos(), v0.pzpos()}, +1);
o2::track::TrackPar negTrack = o2::track::TrackPar({v0.x(), v0.y(), v0.z()}, {v0.pxneg(), v0.pyneg(), v0.pzneg()}, -1);

auto pTra = v0.posTrackExtra_as<dauTracks>();
auto nTra = v0.negTrackExtra_as<dauTracks>();

float deltaTimePositiveLambdaPi = -1e+6;
float deltaTimeNegativeLambdaPi = -1e+6;
float deltaTimePositiveLambdaPr = -1e+6;
Expand All @@ -437,20 +448,38 @@ struct lambdakzeropid {
float timeNegativePr = lengthNegative / velocityNegativePr;
float timeNegativePi = lengthNegative / velocityNegativePi;

if (v0.posTOFSignal() > 0 && v0.posTOFEventTime() > 0 && lengthPositive > 0) {
if (pTra.hasTOF() && lengthPositive > 0) {
deltaTimePositiveLambdaPr = (v0.posTOFSignal() - v0.posTOFEventTime()) - (timeLambda + timePositivePr);
deltaTimePositiveLambdaPi = (v0.posTOFSignal() - v0.posTOFEventTime()) - (timeLambda + timePositivePi);
deltaTimePositiveK0ShortPi = (v0.posTOFSignal() - v0.posTOFEventTime()) - (timeK0Short + timePositivePi);
}
if (v0.negTOFSignal() > 0 && v0.negTOFEventTime() > 0 && lengthNegative > 0) {
if (nTra.hasTOF() && lengthNegative > 0) {
deltaTimeNegativeLambdaPr = (v0.negTOFSignal() - v0.negTOFEventTime()) - (timeLambda + timeNegativePr);
deltaTimeNegativeLambdaPi = (v0.negTOFSignal() - v0.negTOFEventTime()) - (timeLambda + timeNegativePi);
deltaTimeNegativeK0ShortPi = (v0.negTOFSignal() - v0.negTOFEventTime()) - (timeK0Short + timeNegativePi);
}

if (doQA) {
// calculate and pack properties for QA purposes
int posProperties = 0;
if (lengthPositive > 0)
posProperties = posProperties | (int(1) << kLength);
if (pTra.hasTOF())
posProperties = posProperties | (int(1) << kHasTOF);
int negProperties = 0;
if (lengthNegative > 0)
negProperties = negProperties | (int(1) << kLength);
if (nTra.hasTOF())
negProperties = negProperties | (int(1) << kHasTOF);

histos.fill(HIST("h2dPositiveTOFProperties"), v0.pt(), posProperties);
histos.fill(HIST("h2dNegativeTOFProperties"), v0.pt(), negProperties);
}

float deltaDecayTimeLambda = -10e+4;
float deltaDecayTimeAntiLambda = -10e+4;
float deltaDecayTimeK0Short = -10e+4;
if (v0.posTOFSignal() > 0 && v0.negTOFSignal() > 0 && lengthPositive > 0 && lengthNegative > 0) { // does not depend on event time
if (nTra.hasTOF() && pTra.hasTOF() > 0 && lengthPositive > 0 && lengthNegative > 0) { // does not depend on event time
deltaDecayTimeLambda = (v0.posTOFSignal() - timePositivePr) - (v0.negTOFSignal() - timeNegativePi);
deltaDecayTimeAntiLambda = (v0.posTOFSignal() - timePositivePi) - (v0.negTOFSignal() - timeNegativePr);
deltaDecayTimeK0Short = (v0.posTOFSignal() - timePositivePi) - (v0.negTOFSignal() - timeNegativePi);
Expand All @@ -467,7 +496,7 @@ struct lambdakzeropid {
float betaAntiLambda = -1e+6;
float betaK0Short = -1e+6;

if (v0.posTOFSignal() > 0 && v0.negTOFSignal() > 0 && v0.posTOFEventTime() > 0 && v0.negTOFEventTime() > 0) {
if (nTra.hasTOF() && pTra.hasTOF()) {
betaLambda = (lengthV0 / decayTimeLambda) / 0.0299792458;
betaAntiLambda = (lengthV0 / decayTimeAntiLambda) / 0.0299792458;
betaK0Short = (lengthV0 / decayTimeK0Short) / 0.0299792458;
Expand Down Expand Up @@ -505,10 +534,7 @@ struct lambdakzeropid {
}

if (doQA) {
auto pTra = v0.posTrackExtra_as<dauTracks>();
auto nTra = v0.negTrackExtra_as<dauTracks>();

if (v0.posTOFSignal() > 0 && v0.posTOFEventTime() > 0) {
if (pTra.hasTOF()) {
histos.fill(HIST("h2dProtonMeasuredVsExpected"),
(timeLambda + timePositivePr),
(v0.posTOFSignal() - v0.posTOFEventTime()));
Expand All @@ -531,7 +557,7 @@ struct lambdakzeropid {
}
}

if (v0.negTOFSignal() > 0 && v0.negTOFEventTime() > 0) {
if (nTra.hasTOF()) {
histos.fill(HIST("h2dPionMeasuredVsExpected"),
(timeLambda + timeNegativePi),
(v0.negTOFSignal() - v0.negTOFEventTime()));
Expand Down