Skip to content
Merged
Show file tree
Hide file tree
Changes from all 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
2 changes: 1 addition & 1 deletion PWGCF/Femto/DataModel/PionNucleiTables.h
Original file line number Diff line number Diff line change
Expand Up @@ -15,8 +15,8 @@
/// \author CMY
/// \date 2025-04-10

#include "Framework/AnalysisDataModel.h"
#include "Framework/ASoAHelpers.h"
#include "Framework/AnalysisDataModel.h"

#ifndef PWGCF_FEMTO_DATAMODEL_PIONNUCLEITABLES_H_
#define PWGCF_FEMTO_DATAMODEL_PIONNUCLEITABLES_H_
Expand Down
40 changes: 20 additions & 20 deletions PWGCF/Femto/TableProducer/PiNucleiFemto.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -606,22 +606,22 @@
bool fillCandidateInfoHyper(const aod::DataHypCandsWColl::iterator& V0Hyper, const Ttrack& trackPi, PiNucandidate& piHypercand, bool isMixedEvent)
{
piHypercand.collisionID = V0Hyper.collisionId();
//get hypertriton information
//constexpr double mHe3 = o2::constants::physics::MassHelium3;
//constexpr double mPi = o2::constants::physics::MassPiPlus;
// --- He3
// get hypertriton information
// constexpr double mHe3 = o2::constants::physics::MassHelium3;
// constexpr double mPi = o2::constants::physics::MassPiPlus;
// --- He3
float pxHe3 = V0Hyper.ptHe3() * std::cos(V0Hyper.phiHe3());
float pyHe3 = V0Hyper.ptHe3() * std::sin(V0Hyper.phiHe3());
float pzHe3 = V0Hyper.ptHe3() * std::sinh(V0Hyper.etaHe3());
//float pHe3 = V0Hyper.ptHe3() * std::cosh(V0Hyper.etaHe3());
//float enHe3 = std::sqrt(pHe3 * pHe3 + mHe3 * mHe3);
// --- pi
// float pHe3 = V0Hyper.ptHe3() * std::cosh(V0Hyper.etaHe3());
// float enHe3 = std::sqrt(pHe3 * pHe3 + mHe3 * mHe3);
// --- pi
float pxPi = V0Hyper.ptPi() * std::cos(V0Hyper.phiPi());
float pyPi = V0Hyper.ptPi() * std::sin(V0Hyper.phiPi());
float pzPi = V0Hyper.ptPi() * std::sinh(V0Hyper.etaPi());
//float pPi = V0Hyper.ptPi() * std::cosh(V0Hyper.etaPi());
//float enPi = std::sqrt(pPi * pPi + mPi * mPi);
// --- hypertriton
// float pPi = V0Hyper.ptPi() * std::cosh(V0Hyper.etaPi());
// float enPi = std::sqrt(pPi * pPi + mPi * mPi);
// --- hypertriton
float px = pxHe3 + pxPi;
float py = pyHe3 + pyPi;
float pz = pzHe3 + pzPi;
Expand All @@ -639,12 +639,12 @@
}

piHypercand.signPi = trackPi.sign();
if(V0Hyper.isMatter()){
if (V0Hyper.isMatter()) {
piHypercand.signNu = 1;
}else{
} else {
piHypercand.signNu = -1;
}

piHypercand.dcaxyPi = trackPi.dcaXY();
piHypercand.dcazPi = trackPi.dcaZ();
piHypercand.tpcSignalPi = trackPi.tpcSignal();
Expand Down Expand Up @@ -721,14 +721,14 @@
void pairTracksSameEventHyper(const Ttrack& piTracks, const Thypers& V0Hypers)
{
for (const auto& piTrack : piTracks) {

mQaRegistry.fill(HIST("hTrackSel"), Selections::kNoCuts);

if (!selectTrack(piTrack)) {
continue;
}
mQaRegistry.fill(HIST("hTrackSel"), Selections::kTrackCuts);

if (!selectionPIDPion(piTrack)) {
continue;
}
Expand All @@ -738,7 +738,7 @@

SVCand pair;
pair.tr0Idx = piTrack.globalIndex();
pair.tr1Idx = V0Hyper.globalIndex();
pair.tr1Idx = V0Hyper.globalIndex();
const int collIdx = V0Hyper.collisionId();
CollBracket collBracket{collIdx, collIdx};
pair.collBracket = collBracket;
Expand Down Expand Up @@ -856,7 +856,7 @@
mQaRegistry.fill(HIST("hNuPitInvMass"), piNucand.invMass);
mQaRegistry.fill(HIST("hdcaxyNu"), piNucand.dcaxyNu);
mQaRegistry.fill(HIST("hdcazNu"), piNucand.dcazNu);
mQaRegistry.fill(HIST("hdcazNu_min"), (abs(piNucand.dcazNu) - settingCutDeDCAzMin));

Check failure on line 859 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
mQaRegistry.fill(HIST("hNClsNuITS"), piNucand.nClsItsNu);
mQaRegistry.fill(HIST("hNClsPiITS"), piNucand.nClsItsPi);
mQaRegistry.fill(HIST("hisBkgEM"), piNucand.isBkgEM);
Expand All @@ -864,10 +864,10 @@

double computePrTPCnsig(double InnerParamTPCHad, double SignalTPCHad)
{
double m_BBparamsProton[6] = {-54.42066571222577, 0.2857381250239097, 1.247140602468868, 0.6297483918147729, 2.985438833884555, 0.09};

Check failure on line 867 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.

float TPCinnerParam = InnerParamTPCHad;

Check failure on line 869 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
float expTPCSignal = o2::tpc::BetheBlochAleph((TPCinnerParam / 0.9382721), m_BBparamsProton[0], m_BBparamsProton[1], m_BBparamsProton[2], m_BBparamsProton[3], m_BBparamsProton[4]);

Check failure on line 870 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[pdg/explicit-mass]

Avoid hard-coded particle masses. Use o2::constants::physics::Mass... instead.
double resoTPC{expTPCSignal * m_BBparamsProton[5]};
return ((SignalTPCHad - expTPCSignal) / resoTPC);
}
Expand All @@ -878,7 +878,7 @@
const float kp0 = 1.22204e-02;
const float kp1 = 7.48467e-01;

double fSigmaTOFMassHad = (kp0 * TMath::Exp(kp1 * TMath::Abs(ptHad))) * fExpTOFMassHad;

Check failure on line 881 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity]

Replace ROOT entities with equivalents from standard C++ or from O2.
double fNSigmaTOFHad = (MassTOFHad - fExpTOFMassHad) / fSigmaTOFMassHad;
return fNSigmaTOFHad;
}
Expand All @@ -888,7 +888,7 @@
constexpr double massDe = o2::constants::physics::MassDeuteron;
constexpr double massHad = o2::constants::physics::MassPiPlus;

const ROOT::Math::PtEtaPhiMVector De(std::abs(piNucand.recoPtNu()), piNucand.recoEtaNu(), piNucand.recoPhiNu(), massDe);

Check failure on line 891 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/function-variable]

Use lowerCamelCase for names of functions and variables.
const ROOT::Math::PtEtaPhiMVector Had(std::abs(piNucand.recoPtPi()), piNucand.recoEtaPi(), piNucand.recoPhiPi(), massHad);
const ROOT::Math::PtEtaPhiMVector trackSum = De + Had;

Expand All @@ -912,12 +912,12 @@
{
double PrTPCnsigma = computePrTPCnsig(piNucand.momPiTPC, piNucand.tpcSignalPi);
double PrTOFnsigma = tofNSigmaCalculation(piNucand.massTOFPi, piNucand.recoPtPi());
if (abs(PrTPCnsigma) < settingCutNsigTPCPrMin)

Check failure on line 915 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
return;
if (abs(PrTOFnsigma) < settingCutNsigTOFPrMin)

Check failure on line 917 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
return;
float DeDCAxyMin = 0.015 + 0.0305 / TMath::Power(piNucand.recoPtNu(), 1.1);

Check failure on line 919 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[root/entity]

Replace ROOT entities with equivalents from standard C++ or from O2.
if (abs(piNucand.dcaxyNu) > DeDCAxyMin || abs(piNucand.dcazNu) > settingCutDeDCAzMin || abs(piNucand.dcaxyPi) > settingCutPiDCAxyMin || abs(piNucand.dcazPi) > settingCutPiDCAzMin)

Check failure on line 920 in PWGCF/Femto/TableProducer/PiNucleiFemto.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[std-prefix]

Use std:: prefix for names from the std namespace.
return;
fillHistograms(piNucand);

Expand Down Expand Up @@ -971,7 +971,7 @@

auto v0hyper = V0Hypers.rawIteratorAt(trackPair.tr1Idx);
auto piTrack = piTracks.rawIteratorAt(trackPair.tr0Idx);
//auto collBracket = trackPair.collBracket;
// auto collBracket = trackPair.collBracket;

PiNucandidate piNucand;
if (!fillCandidateInfoHyper(v0hyper, piTrack, piNucand, isMixedEvent)) {
Expand Down Expand Up @@ -1033,7 +1033,7 @@
mGoodCollisions[collision.globalIndex()] = true;
const uint64_t collIdx = collision.globalIndex();
auto trackTableThisCollision = pitracks.sliceBy(mPerCol, collIdx);
auto hypdTableThisCollision = V0Hypers.sliceBy(hypPerCol, collIdx);
auto hypdTableThisCollision = V0Hypers.sliceBy(hypPerCol, collIdx);
trackTableThisCollision.bindExternalIndices(&pitracks);
hypdTableThisCollision.bindExternalIndices(&V0Hypers);

Expand All @@ -1043,7 +1043,7 @@
continue;
}

fillPairsHyper(collisions, pitracks, V0Hypers,/*isMixedEvent*/ false);
fillPairsHyper(collisions, pitracks, V0Hypers, /*isMixedEvent*/ false);
}
}
PROCESS_SWITCH(PiNucleiFemto, processSameEventHyper, "Process Same event", false);
Expand Down
Loading