diff --git a/PWGLF/DataModel/LFStrangenessPIDTables.h b/PWGLF/DataModel/LFStrangenessPIDTables.h index 5fbc99ae5ca..36a345b87ad 100644 --- a/PWGLF/DataModel/LFStrangenessPIDTables.h +++ b/PWGLF/DataModel/LFStrangenessPIDTables.h @@ -25,7 +25,7 @@ namespace o2::aod namespace dautrack { // ==== TPC INFORMATION === -DECLARE_SOA_COLUMN(TPCSignal, tpcSignal, float); //! track TPC signal +DECLARE_SOA_COLUMN(TPCSignal, tpcSignal, float); //! track TPC signal DECLARE_SOA_COLUMN(TPCNSigmaEl, tpcNSigmaEl, float); //! Nsigma proton DECLARE_SOA_COLUMN(TPCNSigmaPi, tpcNSigmaPi, float); //! Nsigma proton DECLARE_SOA_COLUMN(TPCNSigmaKa, tpcNSigmaKa, float); //! Nsigma proton @@ -35,7 +35,7 @@ DECLARE_SOA_COLUMN(TPCNSigmaHe, tpcNSigmaHe, float); //! Nsigma proton // ==== TOF INFORMATION === DECLARE_SOA_COLUMN(TOFSignal, tofSignal, float); //! track TOF signal DECLARE_SOA_COLUMN(TOFEvTime, tofEvTime, float); //! track TOF signal -DECLARE_SOA_COLUMN(Length, length, float); //! track TOF signal +DECLARE_SOA_COLUMN(Length, length, float); //! track TOF signal } // namespace dautrack DECLARE_SOA_TABLE(DauTrackTPCPIDs, "AOD", "DAUTRACKTPCPID", // nsigma table (for analysis) diff --git a/PWGLF/TableProducer/cascadepid.cxx b/PWGLF/TableProducer/cascadepid.cxx index f2f7e7c3697..dc890a738c5 100644 --- a/PWGLF/TableProducer/cascadepid.cxx +++ b/PWGLF/TableProducer/cascadepid.cxx @@ -59,7 +59,7 @@ using namespace o2::framework::expressions; using std::array; // For original data loops -using CascOriginalDatas = soa::Join; +using CascOriginalDatas = soa::Join; using TracksWithAllExtras = soa::Join; // Cores with references and TOF pid @@ -74,7 +74,8 @@ struct cascadepid { Service ccdb; // For manual sliceBy - Preslice perCollisionOriginal = o2::aod::cascdata::collisionId;; + Preslice perCollisionOriginal = o2::aod::cascdata::collisionId; + ; Preslice perCollisionDerived = o2::aod::cascdata::straCollisionId; HistogramRegistry histos{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject}; @@ -383,7 +384,8 @@ struct cascadepid { } template - void processCascadeCandidate(TCollision const& collision, TCascade const& cascade, TTrack const& pTra, TTrack const &nTra, TTrack const &bTra) { + void processCascadeCandidate(TCollision const& collision, TCascade const& cascade, TTrack const& pTra, TTrack const& nTra, TTrack const& bTra) + { // initialize from positions and momenta as needed o2::track::TrackPar posTrack = o2::track::TrackPar({cascade.xlambda(), cascade.ylambda(), cascade.zlambda()}, {cascade.pxpos(), cascade.pypos(), cascade.pzpos()}, +1); o2::track::TrackPar negTrack = o2::track::TrackPar({cascade.xlambda(), cascade.ylambda(), cascade.zlambda()}, {cascade.pxneg(), cascade.pyneg(), cascade.pzneg()}, -1); @@ -580,8 +582,8 @@ struct cascadepid { auto pTra = cascade.posTrack_as(); auto nTra = cascade.negTrack_as(); auto bTra = cascade.bachelor_as(); - - processCascadeCandidate(collision, cascade, pTra, nTra, bTra); + + processCascadeCandidate(collision, cascade, pTra, nTra, bTra); } } } @@ -600,8 +602,8 @@ struct cascadepid { auto pTra = cascade.posTrackExtra_as(); auto nTra = cascade.negTrackExtra_as(); auto bTra = cascade.bachTrackExtra_as(); - - processCascadeCandidate(collision, cascade, pTra, nTra, bTra); + + processCascadeCandidate(collision, cascade, pTra, nTra, bTra); } } } diff --git a/PWGLF/TableProducer/converters/strangederivedconverter.cxx b/PWGLF/TableProducer/converters/strangederivedconverter.cxx index 9603bf6d01f..4e5a95189bf 100644 --- a/PWGLF/TableProducer/converters/strangederivedconverter.cxx +++ b/PWGLF/TableProducer/converters/strangederivedconverter.cxx @@ -67,8 +67,8 @@ struct strangeDerivedConverter { auto nTra = casc.negTrackExtra_as>(); auto bTra = casc.bachTrackExtra_as>(); casctofs(pTra.length(), nTra.length(), bTra.length(), - pTra.tofSignal(), nTra.tofSignal(), bTra.tofSignal(), - pTra.tofEvTime(), nTra.tofEvTime(), bTra.tofEvTime()); + pTra.tofSignal(), nTra.tofSignal(), bTra.tofSignal(), + pTra.tofEvTime(), nTra.tofEvTime(), bTra.tofEvTime()); } PROCESS_SWITCH(strangeDerivedConverter, processStraRawCents000to001, "from StraRawCents 000 to 001", false); diff --git a/PWGLF/TableProducer/lambdakzeropid.cxx b/PWGLF/TableProducer/lambdakzeropid.cxx index 41cd3c421a0..56f8f308102 100644 --- a/PWGLF/TableProducer/lambdakzeropid.cxx +++ b/PWGLF/TableProducer/lambdakzeropid.cxx @@ -64,7 +64,7 @@ using TracksWithAllExtras = soa::Join; -using V0DerivedDatas = soa::Join; +using V0DerivedDatas = soa::Join; struct lambdakzeropid { // TOF pid for strangeness (recalculated with topology) @@ -76,7 +76,8 @@ struct lambdakzeropid { Service ccdb; // For manual sliceBy - Preslice perCollisionOriginal = o2::aod::v0data::collisionId;; + Preslice perCollisionOriginal = o2::aod::v0data::collisionId; + ; Preslice perCollisionDerived = o2::aod::v0data::straCollisionId; HistogramRegistry histos{"Histos", {}, OutputObjHandlingPolicy::AnalysisObject}; @@ -403,9 +404,10 @@ struct lambdakzeropid { return 0.0299792458 * TMath::Sqrt(lA / (1 + lA)); } - // templatized process function for symmetric operation in derived and original AO2D + // templatized process function for symmetric operation in derived and original AO2D template - void processV0Candidate(TCollision const& collision, TV0 const& v0, TTrack const& pTra, TTrack const &nTra) { + void processV0Candidate(TCollision const& collision, TV0 const& v0, TTrack const& pTra, TTrack const& nTra) + { // time of V0 segment float lengthV0 = std::hypot(v0.x() - collision.posX(), v0.y() - collision.posY(), v0.z() - collision.posZ()); float velocityK0Short = velocity(v0.p(), o2::constants::physics::MassKaonNeutral); @@ -497,10 +499,10 @@ struct lambdakzeropid { betaK0Short = (lengthV0 / decayTimeK0Short) / 0.0299792458; } - v0tofpid( deltaTimePositiveLambdaPi, deltaTimePositiveLambdaPr, - deltaTimeNegativeLambdaPi, deltaTimeNegativeLambdaPr, - deltaTimePositiveK0ShortPi, deltaTimeNegativeK0ShortPi, - deltaDecayTimeLambda, deltaDecayTimeAntiLambda, deltaDecayTimeK0Short); + v0tofpid(deltaTimePositiveLambdaPi, deltaTimePositiveLambdaPr, + deltaTimeNegativeLambdaPi, deltaTimeNegativeLambdaPr, + deltaTimePositiveK0ShortPi, deltaTimeNegativeK0ShortPi, + deltaDecayTimeLambda, deltaDecayTimeAntiLambda, deltaDecayTimeK0Short); v0tofbeta(betaLambda, betaAntiLambda, betaK0Short); v0tofdebugs(timeLambda, timeK0Short, timePositivePr, timePositivePi, timeNegativePr, timeNegativePi); @@ -594,8 +596,8 @@ struct lambdakzeropid { // de-reference interlinks by hand for derived data auto pTra = v0.posTrack_as(); auto nTra = v0.negTrack_as(); - - processV0Candidate(collision, v0, pTra, nTra); + + processV0Candidate(collision, v0, pTra, nTra); } } } @@ -614,8 +616,8 @@ struct lambdakzeropid { // de-reference interlinks by hand for derived data auto pTra = v0.posTrackExtra_as(); auto nTra = v0.negTrackExtra_as(); - - processV0Candidate(collision, v0, pTra, nTra); + + processV0Candidate(collision, v0, pTra, nTra); } } } diff --git a/PWGLF/Tasks/Strangeness/cascadeanalysis.cxx b/PWGLF/Tasks/Strangeness/cascadeanalysis.cxx index 7118dd2464d..baf37b72717 100644 --- a/PWGLF/Tasks/Strangeness/cascadeanalysis.cxx +++ b/PWGLF/Tasks/Strangeness/cascadeanalysis.cxx @@ -275,12 +275,12 @@ struct cascadeAnalysis { bool lCompatiblePID_Xi = (lPIDvalue >> 0 & 1); bool lCompatiblePID_Om = (lPIDvalue >> 1 & 1); - if (casc.v0radius() > v0setting_radius && - casc.cascradius() > cascadesetting_cascradius && - casc.v0cosPA(pvx, pvy, pvz) > v0setting_cospa && - casc.casccosPA(pvx, pvy, pvz) > cascadesetting_cospa && - casc.dcav0topv(pvx, pvy, pvz) > cascadesetting_mindcav0topv && - TMath::Abs(casc.mLambda() - 1.115683) < cascadesetting_v0masswindow) { + if (casc.v0radius() > v0setting_radius && + casc.cascradius() > cascadesetting_cascradius && + casc.v0cosPA(pvx, pvy, pvz) > v0setting_cospa && + casc.casccosPA(pvx, pvy, pvz) > cascadesetting_cospa && + casc.dcav0topv(pvx, pvy, pvz) > cascadesetting_mindcav0topv && + TMath::Abs(casc.mLambda() - 1.115683) < cascadesetting_v0masswindow) { registry.fill(HIST("hCandidateCounter"), 3.5); // pass cascade selections if (casc.sign() < 0) { // FIXME: could be done better... if (TMath::Abs(casc.yXi()) < 0.5 && lCompatiblePID_Xi) {