Skip to content

[PWGLF] Add additional track info for systematics #14646

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
skundu692 merged 1 commit into from
Jan 27, 2026
+185 −6
Merged

[PWGLF] Add additional track info for systematics #14646

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
54 changes: 52 additions & 2 deletions PWGLF/DataModel/ReducedF1ProtonTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -77,6 +77,31 @@ DECLARE_SOA_COLUMN(ProtonNsigmaTPC, protonNsigmaTPC, float); //! P
DECLARE_SOA_COLUMN(ProtonTOFHit, protonTOFHit, int); //! Proton TOF Hit
DECLARE_SOA_COLUMN(ProtonNsigmaTOF, protonNsigmaTOF, float); //! Proton TOF nsigma
DECLARE_SOA_COLUMN(F1ProtonIndex, f1ProtonIndex, int64_t); //! F1 proton index
DECLARE_SOA_COLUMN(PionTOF, pionTOF, float); //! Pion TOF
DECLARE_SOA_COLUMN(PionDcaxy, pionDcaxy, float); //! Pion DCAxy
DECLARE_SOA_COLUMN(PionDcaz, pionDcaz, float); //! Pion DCAz
DECLARE_SOA_COLUMN(PionTPCNcls, pionTPCNcls, float); //! Pion TPC Ncls
DECLARE_SOA_COLUMN(PionTPCNcrs, pionTPCNcrs, float); //! Pion TPC Ncrs
DECLARE_SOA_COLUMN(KaonTOF, kaonTOF, float); //! Kaon TOF
DECLARE_SOA_COLUMN(KaonDcaxy, kaonDcaxy, float); //! Kaon DCAxy
DECLARE_SOA_COLUMN(KaonDcaz, kaonDcaz, float); //! Kaon DCAz
DECLARE_SOA_COLUMN(KaonTPCNcls, kaonTPCNcls, float); //! Kaon TPC Ncls
DECLARE_SOA_COLUMN(KaonTPCNcrs, kaonTPCNcrs, float); //! Kaon TPC Ncrs
DECLARE_SOA_COLUMN(K0D1Dcaxy, k0D1Dcaxy, float); //! K0 daughter 1 DCAxy
DECLARE_SOA_COLUMN(K0D1TPCNcls, k0D1TPCNcls, float); //! K0 daughter 1 TPC Ncls
DECLARE_SOA_COLUMN(K0D1TPCNcrs, k0D1TPCNcrs, float); //! K0 daughter 1 TPC Ncrs
DECLARE_SOA_COLUMN(K0D2Dcaxy, k0D2Dcaxy, float); //! K0 daughter 2 DCAxy
DECLARE_SOA_COLUMN(K0D2TPCNcls, k0D2TPCNcls, float); //! K0 daughter 2 TPC Ncls
DECLARE_SOA_COLUMN(K0D2TPCNcrs, k0D2TPCNcrs, float); //! K0 daughter 2 TPC Ncrs
DECLARE_SOA_COLUMN(K0Cpa, k0Cpa, float); //! K0 CPA
DECLARE_SOA_COLUMN(K0Radius, k0Radius, float); //! K0 decay radius
DECLARE_SOA_COLUMN(K0DcaDaughters, k0DcaDaughters, float); //! K0 DCA between daughters
DECLARE_SOA_COLUMN(K0Dca, k0Dca, float); //! K0 DCA to PV
DECLARE_SOA_COLUMN(K0LifeTime, k0LifeTime, float); //! K0 proper lifetime
DECLARE_SOA_COLUMN(ProtonDcaxy, protonDcaxy, float); //! Proton DCAxy
DECLARE_SOA_COLUMN(ProtonDcaz, protonDcaz, float); //! Proton DCAz
DECLARE_SOA_COLUMN(ProtonTPCNcls, protonTPCNcls, float); //! Proton TPC Ncls
DECLARE_SOA_COLUMN(ProtonTPCNcrs, protonTPCNcrs, float); //! Proton TPC Ncrs
} // namespace f1protondaughter
DECLARE_SOA_TABLE(F1Tracks, "AOD", "F1TRACK",
o2::soa::Index<>,
Expand Down Expand Up @@ -104,7 +129,28 @@ DECLARE_SOA_TABLE(F1Tracks, "AOD", "F1TRACK",
f1protondaughter::F1PionIndex,
f1protondaughter::F1KaonIndex,
f1protondaughter::F1KshortPositiveIndex,
f1protondaughter::F1KshortNegativeIndex);
f1protondaughter::F1KshortNegativeIndex,
f1protondaughter::PionTOF,
f1protondaughter::PionDcaxy,
f1protondaughter::PionDcaz,
f1protondaughter::PionTPCNcls,
f1protondaughter::PionTPCNcrs,
f1protondaughter::KaonTOF,
f1protondaughter::KaonDcaxy,
f1protondaughter::KaonDcaz,
f1protondaughter::KaonTPCNcls,
f1protondaughter::KaonTPCNcrs,
f1protondaughter::K0D1Dcaxy,
f1protondaughter::K0D1TPCNcls,
f1protondaughter::K0D1TPCNcrs,
f1protondaughter::K0D2Dcaxy,
f1protondaughter::K0D2TPCNcls,
f1protondaughter::K0D2TPCNcrs,
f1protondaughter::K0Cpa,
f1protondaughter::K0Radius,
f1protondaughter::K0DcaDaughters,
f1protondaughter::K0Dca,
f1protondaughter::K0LifeTime);
using F1Track = F1Tracks::iterator;

DECLARE_SOA_TABLE(ProtonTracks, "AOD", "PROTONTRACK",
Expand All @@ -117,7 +163,11 @@ DECLARE_SOA_TABLE(ProtonTracks, "AOD", "PROTONTRACK",
f1protondaughter::ProtonNsigmaTPC,
f1protondaughter::ProtonTOFHit,
f1protondaughter::ProtonNsigmaTOF,
f1protondaughter::F1ProtonIndex);
f1protondaughter::F1ProtonIndex,
f1protondaughter::ProtonDcaxy,
f1protondaughter::ProtonDcaz,
f1protondaughter::ProtonTPCNcls,
f1protondaughter::ProtonTPCNcrs);
using ProtonTrack = ProtonTracks::iterator;
} // namespace o2::aod
#endif // PWGLF_DATAMODEL_REDUCEDF1PROTONTABLES_H_
135 changes: 132 additions & 3 deletions PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -40,13 +40,13 @@

#include <Math/GenVector/Boost.h>
#include <Math/Vector4D.h>
#include <TDatabasePDG.h> // FIXME

Check failure on line 43 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
#include <TMath.h>
#include <TPDGCode.h> // FIXME

#include <fairlogger/Logger.h>

#include <iostream>

Check failure on line 49 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[include-iostream] 

Do not include iostream. Use O2 logging instead.
#include <iterator>
#include <map>
#include <string>
Expand Down Expand Up @@ -351,7 +351,7 @@
template <typename Collision, typename V0>
bool SelectionV0(Collision const& collision, V0 const& candidate)
{
if (fabs(candidate.dcav0topv()) > cMaxV0DCA) {

Check failure on line 354 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix] 

Use std:: prefix for names from the std namespace.
return false;
}

Expand All @@ -361,7 +361,7 @@
const float dcaDaughv0 = candidate.dcaV0daughters();
const float cpav0 = candidate.v0cosPA();

float CtauK0s = candidate.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * TDatabasePDG::Instance()->GetParticle(kK0Short)->Mass(); // FIXME: Get from the common header

Check failure on line 364 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
float lowmasscutks0 = 0.497 - 2.0 * cSigmaMassKs0;
float highmasscutks0 = 0.497 + 2.0 * cSigmaMassKs0;

Expand All @@ -385,7 +385,7 @@
return false;
}
}
if (fabs(CtauK0s) > cMaxV0LifeTime || candidate.mK0Short() < lowmasscutks0 || candidate.mK0Short() > highmasscutks0) {

Check failure on line 388 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix] 

Use std:: prefix for names from the std namespace.
return false;
}
return true;
Expand Down Expand Up @@ -488,10 +488,10 @@

std::vector<double> BBProton, BBAntiproton, BBPion, BBAntipion, BBKaon, BBAntikaon;
ROOT::Math::PtEtaPhiMVector F1Vector, F1VectorDummy, F1d1dummy, F1d2dummy, F1d3dummy, KKs0Vector, ProtonVectorDummy, ProtonVectorDummy2;
double massPi = TDatabasePDG::Instance()->GetParticle(kPiPlus)->Mass(); // FIXME: Get from the common header

Check failure on line 491 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
double massKa = TDatabasePDG::Instance()->GetParticle(kKPlus)->Mass(); // FIXME: Get from the common header

Check failure on line 492 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
double massPr = TDatabasePDG::Instance()->GetParticle(kProton)->Mass(); // FIXME: Get from the common header

Check failure on line 493 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
double massK0s = TDatabasePDG::Instance()->GetParticle(kK0Short)->Mass(); // FIXME: Get from the common header

Check failure on line 494 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.

double massF1{0.};
double masskKs0{0.};
Expand Down Expand Up @@ -556,12 +556,16 @@
std::vector<int> PionTOFHitFinal = {};
std::vector<float> PionTPC = {};
std::vector<float> PionTPCFinal = {};
std::vector<float> PionTOFNsigma = {};
std::vector<float> PionTOFNsigmaFinal = {};

// keep TOF Hit of kaon
std::vector<int> KaonTOFHit = {};
std::vector<int> KaonTOFHitFinal = {};
std::vector<float> KaonTPC = {};
std::vector<float> KaonTPCFinal = {};
std::vector<float> KaonTOFNsigma = {};
std::vector<float> KaonTOFNsigmaFinal = {};
std::vector<float> KaonTPCPionHypo = {};
std::vector<float> KaonTPCPionHypoFinal = {};

Expand All @@ -571,6 +575,68 @@
// keep status of F1 signal unlike or wrong sign
std::vector<float> f1signal = {};

// keep track dcaxy
std::vector<float> PionDcaxy = {};
std::vector<float> KaonDcaxy = {};
std::vector<float> ProtonDcaxy = {};
std::vector<float> K0D1Dcaxy = {};
std::vector<float> K0D2Dcaxy = {};

// keep track ncls
std::vector<float> PionTPCNcls = {};
std::vector<float> KaonTPCNcls = {};
std::vector<float> ProtonTPCNcls = {};
std::vector<float> K0D1TPCNcls = {};
std::vector<float> K0D2TPCNcls = {};

// keep track ncrs
std::vector<float> PionTPCNcrs = {};
std::vector<float> KaonTPCNcrs = {};
std::vector<float> ProtonTPCNcrs = {};
std::vector<float> K0D1TPCNcrs = {};
std::vector<float> K0D2TPCNcrs = {};

// keep track dcaz
std::vector<float> PionDcaz = {};
std::vector<float> KaonDcaz = {};
std::vector<float> ProtonDcaz = {};

// keep v0 topology
std::vector<float> K0Cpa = {};
std::vector<float> K0radius = {};
std::vector<float> K0DcaBetweenDaughters = {};
std::vector<float> K0Dca = {};
std::vector<float> K0LifeTime = {};

// keep track dcaxy
std::vector<float> PionDcaxyFinal = {};
std::vector<float> KaonDcaxyFinal = {};
std::vector<float> K0D1DcaxyFinal = {};
std::vector<float> K0D2DcaxyFinal = {};

// keep track ncls
std::vector<float> PionTPCNclsFinal = {};
std::vector<float> KaonTPCNclsFinal = {};
std::vector<float> K0D1TPCNclsFinal = {};
std::vector<float> K0D2TPCNclsFinal = {};

// keep track ncrs
std::vector<float> PionTPCNcrsFinal = {};
std::vector<float> KaonTPCNcrsFinal = {};
std::vector<float> K0D1TPCNcrsFinal = {};
std::vector<float> K0D2TPCNcrsFinal = {};

// keep track dcaz
std::vector<float> PionDcazFinal = {};
std::vector<float> KaonDcazFinal = {};

// keep v0 topology
std::vector<float> K0CpaFinal = {};
std::vector<float> K0radiusFinal = {};
std::vector<float> K0DcaBetweenDaughtersFinal = {};
std::vector<float> K0DcaFinal = {};
std::vector<float> K0LifeTimeFinal = {};

// Prepare vectors for different species
std::vector<ROOT::Math::PtEtaPhiMVector> protons, kaons, pions, kshorts, f1resonance, f1resonanced1, f1resonanced2, f1resonanced3;
// , protonsfinal;
Expand Down Expand Up @@ -649,7 +715,14 @@
pions.push_back(temp);
PionIndex.push_back(track.globalIndex());
PionCharge.push_back(track.sign());

PionDcaxy.push_back(std::abs(track.dcaXY()));
PionDcaz.push_back(std::abs(track.dcaZ()));
PionTPCNcls.push_back(std::abs(track.tpcNClsFound()));
PionTPCNcrs.push_back(std::abs(track.tpcNClsCrossedRows()));

auto PionTOF = 0;
auto nSigmaPionTOF = -999;
if (track.sign() > 0) {
qaRegistry.fill(HIST("hNsigmaPtpionTPC"), nTPCSigmaP[0], track.pt());
PionTPC.push_back(nTPCSigmaP[0]);
Expand All @@ -660,17 +733,26 @@
}
if (track.hasTOF()) {
qaRegistry.fill(HIST("hNsigmaPtpionTOF"), track.tofNSigmaPi(), track.pt());
nSigmaPionTOF = track.tofNSigmaPi();
PionTOF = 1;
}
PionTOFHit.push_back(PionTOF);
PionTOFNsigma.push_back(nSigmaPionTOF);
}

if ((track.pt() > cMinKaonPt && track.sign() > 0 && SelectionPID(track, strategyPIDKaon, 1, nTPCSigmaP[1]) && (itsResponse.nSigmaITS<o2::track::PID::Kaon>(track) > -3.0 && itsResponse.nSigmaITS<o2::track::PID::Kaon>(track) < 3.0)) || (track.pt() > cMinKaonPt && track.sign() < 0 && SelectionPID(track, strategyPIDKaon, 1, nTPCSigmaN[1]) && (itsResponse.nSigmaITS<o2::track::PID::Kaon>(track) > -3.0 && itsResponse.nSigmaITS<o2::track::PID::Kaon>(track) < 3.0))) {
ROOT::Math::PtEtaPhiMVector temp(track.pt(), track.eta(), track.phi(), massKa);
kaons.push_back(temp);
KaonIndex.push_back(track.globalIndex());
KaonCharge.push_back(track.sign());

KaonDcaxy.push_back(std::abs(track.dcaXY()));
KaonDcaz.push_back(std::abs(track.dcaZ()));
KaonTPCNcls.push_back(std::abs(track.tpcNClsFound()));
KaonTPCNcrs.push_back(std::abs(track.tpcNClsCrossedRows()));

auto KaonTOF = 0;
auto nSigmaKaonTOF = -999.0;
if (track.sign() > 0) {
qaRegistry.fill(HIST("hNsigmaPtkaonTPC"), nTPCSigmaP[1], nTPCSigmaP[0], track.pt());
KaonTPC.push_back(nTPCSigmaP[1]);
Expand All @@ -684,15 +766,23 @@
if (track.hasTOF()) {
qaRegistry.fill(HIST("hNsigmaPtkaonTOF"), track.tofNSigmaKa(), track.pt());
KaonTOF = 1;
nSigmaKaonTOF = track.tofNSigmaKa();
}
KaonTOFHit.push_back(KaonTOF);
KaonTOFNsigma.push_back(nSigmaKaonTOF);
}

if ((track.pt() < cMaxProtonPt && track.sign() > 0 && SelectionPID(track, strategyPIDProton, 2, nTPCSigmaP[2]) && (itsResponse.nSigmaITS<o2::track::PID::Proton>(track) > -3.0 && itsResponse.nSigmaITS<o2::track::PID::Proton>(track) < 3.0)) || (track.pt() < cMaxProtonPt && track.sign() < 0 && SelectionPID(track, strategyPIDProton, 2, nTPCSigmaN[2]) && (itsResponse.nSigmaITS<o2::track::PID::Proton>(track) > -3.0 && itsResponse.nSigmaITS<o2::track::PID::Proton>(track) < 3.0))) {
ROOT::Math::PtEtaPhiMVector temp(track.pt(), track.eta(), track.phi(), massPr);
protons.push_back(temp);
ProtonIndex.push_back(track.globalIndex());
ProtonCharge.push_back(track.sign());

ProtonDcaxy.push_back(std::abs(track.dcaXY()));
ProtonDcaz.push_back(std::abs(track.dcaZ()));
ProtonTPCNcls.push_back(std::abs(track.tpcNClsFound()));
ProtonTPCNcrs.push_back(std::abs(track.tpcNClsCrossedRows()));

if (track.sign() > 0) {
qaRegistry.fill(HIST("hNsigmaPtprotonTPC"), nTPCSigmaP[2], track.pt());
ProtonTPCNsigma.push_back(nTPCSigmaP[2]);
Expand All @@ -707,7 +797,7 @@
ProtonTOFHit.push_back(1);
}
if (!track.hasTOF()) {
ProtonTOFNsigma.push_back(999.0);
ProtonTOFNsigma.push_back(-999.0);
ProtonTOFHit.push_back(0);
}
}
Expand Down Expand Up @@ -739,6 +829,20 @@
kshorts.push_back(temp);
KshortPosDaughIndex.push_back(postrack.globalIndex());
KshortNegDaughIndex.push_back(negtrack.globalIndex());

K0D1Dcaxy.push_back(std::abs(postrack.dcaXY()));
K0D1TPCNcls.push_back(std::abs(postrack.tpcNClsFound()));
K0D1TPCNcrs.push_back(std::abs(postrack.tpcNClsCrossedRows()));

K0D2Dcaxy.push_back(std::abs(negtrack.dcaXY()));
K0D2TPCNcls.push_back(std::abs(negtrack.tpcNClsFound()));
K0D2TPCNcrs.push_back(std::abs(negtrack.tpcNClsCrossedRows()));

K0Cpa.push_back(std::abs(v0.v0cosPA()));
K0radius.push_back(std::abs(v0.v0radius()));
K0DcaBetweenDaughters.push_back(std::abs(v0.dcaV0daughters()));
K0LifeTime.push_back(std::abs(v0.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * TDatabasePDG::Instance()->GetParticle(kK0Short)->Mass()));

Check failure on line 844 in PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/database] 

Do not use TDatabasePDG directly. Use o2::constants::physics::Mass... or Service<o2::framework::O2DatabasePDG> instead.
K0Dca.push_back(std::abs(v0.dcav0topv()));
}

if (pions.size() != 0 && kaons.size() != 0 && kshorts.size() != 0) {
Expand Down Expand Up @@ -792,6 +896,31 @@
PionTPCFinal.push_back(PionTPC.at(i1)); // Pion TPC
KaonTPCFinal.push_back(KaonTPC.at(i2)); // Kaon TPC
KaonTPCPionHypoFinal.push_back(KaonTPCPionHypo.at(i2)); // Kaon TPC

// Fill variable for systematic

PionTOFNsigmaFinal.push_back(PionTOFNsigma.at(i1));
PionDcaxyFinal.push_back(PionDcaxy.at(i1));
PionDcazFinal.push_back(PionDcaz.at(i1));
PionTPCNcrsFinal.push_back(PionTPCNcrs.at(i1));
PionTPCNclsFinal.push_back(PionTPCNcls.at(i1));
KaonTOFNsigmaFinal.push_back(KaonTOFNsigma.at(i2));
KaonDcaxyFinal.push_back(KaonDcaxy.at(i2));
KaonDcazFinal.push_back(KaonDcaz.at(i2));
KaonTPCNcrsFinal.push_back(KaonTPCNcrs.at(i2));
KaonTPCNclsFinal.push_back(KaonTPCNcls.at(i2));
K0D1DcaxyFinal.push_back(K0D1Dcaxy.at(i3));
K0D1TPCNcrsFinal.push_back(K0D1TPCNcrs.at(i3));
K0D1TPCNclsFinal.push_back(K0D1TPCNcls.at(i3));
K0D2DcaxyFinal.push_back(K0D2Dcaxy.at(i3));
K0D2TPCNcrsFinal.push_back(K0D2TPCNcrs.at(i3));
K0D2TPCNclsFinal.push_back(K0D2TPCNcls.at(i3));
K0CpaFinal.push_back(K0Cpa.at(i3));
K0radiusFinal.push_back(K0radius.at(i3));
K0DcaBetweenDaughtersFinal.push_back(K0DcaBetweenDaughters.at(i3));
K0DcaFinal.push_back(K0Dca.at(i3));
K0LifeTimeFinal.push_back(K0LifeTime.at(i3));

if (pairsign > 0) {
qaRegistry.fill(HIST("hInvMassf1"), F1Vector.M(), F1Vector.Pt());
numberF1 = numberF1 + 1;
Expand Down Expand Up @@ -848,13 +977,13 @@
F1d1dummy = f1resonanced1.at(i5);
F1d2dummy = f1resonanced2.at(i5);
F1d3dummy = f1resonanced3.at(i5);
f1track(indexEvent, f1signal.at(i5), F1VectorDummy.Px(), F1VectorDummy.Py(), F1VectorDummy.Pz(), F1d1dummy.Px(), F1d1dummy.Py(), F1d1dummy.Pz(), F1d2dummy.Px(), F1d2dummy.Py(), F1d2dummy.Pz(), F1d3dummy.Px(), F1d3dummy.Py(), F1d3dummy.Pz(), PionTOFHitFinal.at(i5), KaonTOFHitFinal.at(i5), PionTPCFinal.at(i5), KaonTPCFinal.at(i5), KaonTPCPionHypoFinal.at(i5), F1VectorDummy.M(), f1kaonkshortmass.at(i5), F1PionIndex.at(i5), F1KaonIndex.at(i5), F1KshortDaughterPositiveIndex.at(i5), F1KshortDaughterNegativeIndex.at(i5));
f1track(indexEvent, f1signal.at(i5), F1VectorDummy.Px(), F1VectorDummy.Py(), F1VectorDummy.Pz(), F1d1dummy.Px(), F1d1dummy.Py(), F1d1dummy.Pz(), F1d2dummy.Px(), F1d2dummy.Py(), F1d2dummy.Pz(), F1d3dummy.Px(), F1d3dummy.Py(), F1d3dummy.Pz(), PionTOFHitFinal.at(i5), KaonTOFHitFinal.at(i5), PionTPCFinal.at(i5), KaonTPCFinal.at(i5), KaonTPCPionHypoFinal.at(i5), F1VectorDummy.M(), f1kaonkshortmass.at(i5), F1PionIndex.at(i5), F1KaonIndex.at(i5), F1KshortDaughterPositiveIndex.at(i5), F1KshortDaughterNegativeIndex.at(i5), PionTOFNsigmaFinal.at(i5), PionDcaxyFinal.at(i5), PionDcazFinal.at(i5), PionTPCNclsFinal.at(i5), PionTPCNcrsFinal.at(i5), KaonTOFNsigmaFinal.at(i5), KaonDcaxyFinal.at(i5), KaonDcazFinal.at(i5), KaonTPCNclsFinal.at(i5), KaonTPCNcrsFinal.at(i5), K0D1DcaxyFinal.at(i5), K0D1TPCNclsFinal.at(i5), K0D1TPCNcrsFinal.at(i5), K0D2DcaxyFinal.at(i5), K0D2TPCNclsFinal.at(i5), K0D2TPCNcrsFinal.at(i5), K0CpaFinal.at(i5), K0radiusFinal.at(i5), K0DcaBetweenDaughtersFinal.at(i5), K0DcaFinal.at(i5), K0LifeTimeFinal.at(i5));
}
//// Fill track table for proton//////////////////
for (auto iproton = protons.begin(); iproton != protons.end(); ++iproton) {
auto i6 = std::distance(protons.begin(), iproton);
ProtonVectorDummy2 = protons.at(i6);
protontrack(indexEvent, ProtonCharge.at(i6), ProtonVectorDummy2.Px(), ProtonVectorDummy2.Py(), ProtonVectorDummy2.Pz(), ProtonTPCNsigma.at(i6), ProtonTOFHit.at(i6), ProtonTOFNsigma.at(i6), ProtonIndex.at(i6));
protontrack(indexEvent, ProtonCharge.at(i6), ProtonVectorDummy2.Px(), ProtonVectorDummy2.Py(), ProtonVectorDummy2.Pz(), ProtonTPCNsigma.at(i6), ProtonTOFHit.at(i6), ProtonTOFNsigma.at(i6), ProtonIndex.at(i6), ProtonDcaxy.at(i6), ProtonDcaz.at(i6), ProtonTPCNcls.at(i6), ProtonTPCNcrs.at(i6));
}
}
}
Expand Down
2 changes: 1 addition & 1 deletion PWGLF/Tasks/Resonances/f1protoncorrelation.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -84,7 +84,7 @@ struct f1protoncorrelation {
// Event Mixing
Configurable<int> nEvtMixing{"nEvtMixing", 10, "Number of events to mix"};
Configurable<int> nEvtMixingBkg{"nEvtMixingBkg", 5, "Number of events to mix for background reconstruction"};
ConfigurableAxis CfgVtxBins{"CfgVtxBins", {10, -10, 10}, "Mixing bins - z-vertex"};
ConfigurableAxis CfgVtxBins{"CfgVtxBins", {VARIABLE_WIDTH, 10, -10, 10}, "Mixing bins - z-vertex"};
ConfigurableAxis CfgMultBins{"CfgMultBins", {VARIABLE_WIDTH, 0.0, 40.0, 80.0, 500.0}, "Mixing bins - number of contributor"};

// THnsparse bining
Expand Down
Loading