Skip to content

[PWGCF] DptDpt - Tracking the insertion of pT bins of interest #10864

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
victor-gonzalez merged 1 commit into from
Apr 15, 2025
Merged

[PWGCF] DptDpt - Tracking the insertion of pT bins of interest #10864

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
8 changes: 6 additions & 2 deletions PWGCF/TwoParticleCorrelations/Tasks/dptDptEfficiencyAndQc.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -95,7 +95,7 @@
int noOfNSigmaBins = static_cast<int>((maxNSigma - minNSigma) / widthNSigmaBin);

/* the pT bins of interest for the relative separation within TPC sectors data collection */
std::vector<int> ptBinsOfInterest{1, 2, 3};
std::vector<int> ptBinsOfInterest{};

/* the PID selector object to help with the configuration and the id of the selected particles */
o2::analysis::dptdptfilter::PIDSpeciesSelection pidselector;
Expand Down Expand Up @@ -411,7 +411,7 @@

float phiInTpcSector = std::fmod(track.phi(), kTpcPhiSectorWidth);
float phiShiftedPercentInTpcSector = phiInTpcSector * 100 / kTpcPhiSectorWidth;
phiShiftedPercentInTpcSector = (phiShiftedPercentInTpcSector > 50.0f) ? (phiShiftedPercentInTpcSector - 100.0f) : phiShiftedPercentInTpcSector;

Check warning on line 414 in PWGCF/TwoParticleCorrelations/Tasks/dptDptEfficiencyAndQc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
if (track.sign() > 0) {
fhPhiVsPtPosB->Fill(track.pt(), phiInTpcSector);
fhPerColNchVsPhiVsPtPosB->Fill(track.pt(), phiInTpcSector);
Expand Down Expand Up @@ -490,7 +490,7 @@
float genid = findgenid(mcparticle);

bool isprimary = mcparticle.isPhysicalPrimary();
bool issecdecay = !isprimary && (mcparticle.getProcess() == 4);

Check warning on line 493 in PWGCF/TwoParticleCorrelations/Tasks/dptDptEfficiencyAndQc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
bool isfrommaterial = !isprimary && !issecdecay;
fillpurityhistos(fhPtPurityPosPrimA, fhPtPurityNegPrimA, genid, track, isprimary);
fillpurityhistos(fhPtPurityPosSecA, fhPtPurityNegSecA, genid, track, issecdecay);
Expand Down Expand Up @@ -524,7 +524,7 @@
/* pure generator level */
if (track.isPhysicalPrimary()) {
fhPtVsEtaPrimA[track.trackacceptedid()]->Fill(track.eta(), track.pt());
} else if (track.getProcess() == 4) {

Check warning on line 527 in PWGCF/TwoParticleCorrelations/Tasks/dptDptEfficiencyAndQc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
fhPtVsEtaSecA[track.trackacceptedid()]->Fill(track.eta(), track.pt());
} else {
fhPtVsEtaMatA[track.trackacceptedid()]->Fill(track.eta(), track.pt());
Expand Down Expand Up @@ -761,7 +761,7 @@
} else {
ix = 2 * ix;
}
for (uint when = 0; when < 2; ++when) {

Check warning on line 764 in PWGCF/TwoParticleCorrelations/Tasks/dptDptEfficiencyAndQc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
fhTPCnSigmasVsP[when][ix]->Fill(tpcmom, o2::aod::pidutils::tpcNSigma<id>(track));
fhTOFnSigmasVsP[when][ix]->Fill(tofmom, o2::aod::pidutils::tofNSigma<id>(track));
if (track.trackacceptedid() < 0) {
Expand All @@ -779,7 +779,7 @@
} else {
ix = 2 * ix;
}
for (uint when = 0; when < 2; ++when) {

Check warning on line 782 in PWGCF/TwoParticleCorrelations/Tasks/dptDptEfficiencyAndQc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
fhTPCdEdxSignalDiffVsP[when][ix]->Fill(tpcmom, o2::aod::pidutils::tpcExpSignalDiff<id>(track));
fhTOFSignalDiffVsP[when][ix]->Fill(tofmom, o2::aod::pidutils::tofExpSignalDiff<id>(track));
fhTPCTOFSigmaVsP[when][ix]->Fill(tpcmom, o2::aod::pidutils::tpcNSigma<id>(track), o2::aod::pidutils::tofNSigma<id>(track));
Expand All @@ -793,7 +793,7 @@
template <typename TrackObject>
void fillPID(TrackObject const& track, float tpcmom, float tofmom)
{
for (uint when = 0; when < 2; ++when) {

Check warning on line 796 in PWGCF/TwoParticleCorrelations/Tasks/dptDptEfficiencyAndQc.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
if constexpr (framework::has_type_v<o2::aod::mcpidtpc::DeDxTunedMc, typename TrackObject::all_columns>) {
fhTPCdEdxSignalVsP[when]->Fill(tpcmom, track.mcTunedTPCSignal());
} else {
Expand Down Expand Up @@ -1114,7 +1114,11 @@
noOfNSigmaBins = static_cast<int>((maxNSigma - minNSigma) / widthNSigmaBin);

/* configure the pT bins of interest */
ptBinsOfInterest = cfgPtBinsOfInterest.value;
TAxis ptAxis{ptbins, ptlow, ptup};
for (const int& bin : cfgPtBinsOfInterest.value) {
ptBinsOfInterest.push_back(bin);
LOGF(info, "Inserted pT bin %d: %.2f < pT < %.2f GeV/c", bin, ptAxis.GetBinLowEdge(bin), ptAxis.GetBinUpEdge(bin));
}

bool doBasicAnalysis = doprocessDetectorLevelNotStored || doprocessReconstructedNotStored || doprocessGeneratorLevelNotStored;
bool doExtraAnalysis = doprocessExtraDetectorLevelNotStored || doprocessExtraReconstructedNotStored || doprocessExtraGeneratorLevelNotStored;
Expand Down
Loading