Skip to content

[PWGCF] Net charge fluctuations3 #16314

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
SwatiSaha-1997 merged 16 commits into from
May 20, 2026
+128 −0
Merged

[PWGCF] Net charge fluctuations3 #16314

Changes from all commits
Commits
Show all changes
16 commits
Select commit Hold shift + click to select a range
2f570eb
Efficiency Correction added
pravatp98 Mar 29, 2026
2bd91c1
Remove commented-out track QA functions
pravatp98 Mar 29, 2026
bd2b742
Remove unused track identification variables
pravatp98 Mar 29, 2026
42b9c74
Add configurable parameters for charge histograms
pravatp98 Mar 30, 2026
e0eef10
Merge branch 'AliceO2Group:master' into NetChargeFluctuations3
pravatp98 Apr 5, 2026
7c2ea78
Merge branch 'AliceO2Group:master' into NetChargeFluctuations3
pravatp98 May 13, 2026
bae294b
Rename file and update includes in nchCumulantsId.cxx
pravatp98 May 13, 2026
c5d1c10
Refactor conditional checks for particle types
pravatp98 May 13, 2026
574598e
Add string header to nchCumulantsId.cxx
pravatp98 May 13, 2026
e781be7
Move vector include to a new line
pravatp98 May 13, 2026
660ec7b
Merge branch 'AliceO2Group:master' into NetChargeFluctuations3
pravatp98 May 13, 2026
200a765
Remove unused track identification variables
pravatp98 May 13, 2026
90ebd03
Add string header to nchCumulantsId.cxx
pravatp98 May 13, 2026
1869483
Merge branch 'AliceO2Group:master' into NetChargeFluctuations3
pravatp98 May 14, 2026
eab4dea
Merge branch 'AliceO2Group:master' into NetChargeFluctuations3
pravatp98 May 20, 2026
60776c0
Enhance PID selection with rejection cuts
pravatp98 May 20, 2026
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
128 changes: 128 additions & 0 deletions PWGCF/EbyEFluctuations/Tasks/nchCumulantsId.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -48,6 +48,12 @@ using namespace std;
#define ID_BIT_EL 3
#define ID_BIT_DE 4

#define MC_BIT_PI 0 // MC particle identification bits
#define MC_BIT_KA 1
#define MC_BIT_PR 2
#define MC_BIT_EL 3
#define MC_BIT_DE 4

#define BITSET(mask, ithBit) ((mask) |= (1 << (ithBit))) // avoid name bitset as std::bitset is already there
#define BITCHECK(mask, ithBit) ((mask) & (1 << (ithBit))) // bit check will return int value, not bool, use BITCHECK != 0 in Analysi

Expand Down Expand Up @@ -108,6 +114,7 @@ struct NchCumulantsId {

Configurable<bool> checkCollPosZMc{"checkCollPosZMc", false, "checkCollPosZMc"};
Configurable<bool> flagUnusedVariableError{"flagUnusedVariableError", false, "flagUnusedVariableError"};
Configurable<bool> cfgDoRejectionForId{"cfgDoRejectionForId", false, "Apply rejection cut before PID selection (selTrackForId)"};

Configurable<bool> cfgEvSel01doNoSameBunchPileup{"cfgEvSel01doNoSameBunchPileup", true, "apply kNoSameBunchPileup"};
Configurable<bool> cfgEvSel02doIsGoodZvtxFT0vsPV{"cfgEvSel02doIsGoodZvtxFT0vsPV", true, "apply kIsGoodZvtxFT0vsPV"};
Expand Down Expand Up @@ -150,6 +157,12 @@ struct NchCumulantsId {
Configurable<float> cfgIdPr08NSigmaTOFHighP{"cfgIdPr08NSigmaTOFHighP", 2.0, "cfgIdPr08NSigmaTOFHighP"};
Configurable<float> cfgIdPr09NSigmaRadHighP{"cfgIdPr09NSigmaRadHighP", 4.0, "cfgIdPr09NSigmaRadHighP"};

Configurable<float> cfgIdElRejLowNSigma{"cfgIdElRejLowNSigma", -3.0, "cfgIdElRejLowNSigma"};
Configurable<float> cfgIdElRejHighNSigma{"cfgIdElRejHighNSigma", 5.0, "cfgIdElRejHighNSigma"};
Configurable<float> cfgIdPiRejNSigma{"cfgIdPiRejNSigma", 3.0, "cfgIdPiRejNSigma"};
Configurable<float> cfgIdKaRejNSigma{"cfgIdKaRejNSigma", 3.0, "cfgIdKaRejNSigma"};
Configurable<float> cfgIdPrRejNSigma{"cfgIdPrRejNSigma", 3.0, "cfgIdPrRejNSigma"};

struct : ConfigurableGroup {
Configurable<float> cfgVetoId01PiTPC{"cfgVetoId01PiTPC", 3.0, "cfgVetoId01PiTPC"};
Configurable<float> cfgVetoId02PiTOF{"cfgVetoId02PiTOF", 3.0, "cfgVetoId02PiTOF"};
Expand Down Expand Up @@ -244,8 +257,13 @@ struct NchCumulantsId {
const AxisSpec axisPiCh(101, -1, 100, "Pion_Positive");
const AxisSpec axisAPiCh(101, -1, 100, "Pion_Negative");

const AxisSpec axisIdTag = {32, -0.5f, 31.5f, "idTag"};
const AxisSpec axisMcTag = {32, -0.5f, 31.5f, "mcTag"};

HistogramConfigSpec qnHist1({HistType::kTHnSparseD, {axisNch, axisPosCh, axisNegCh, axisPrCh, axisAPrCh, axisKaCh, axisAKaCh, axisNt, axisCent}});
HistogramConfigSpec qnHist2({HistType::kTHnSparseD, {axisNch, axisPosCh, axisNegCh, axisPiCh, axisAPiCh, axisKaCh, axisAKaCh, axisNt, axisCent}});
HistogramConfigSpec histTPCPIDSparse({HistType::kTHnSparseD, {axisP, axisTPCNSigma, axisIdTag, axisMcTag}});
HistogramConfigSpec histTOFPIDSparse({HistType::kTHnSparseD, {axisP, axisTOFNSigma, axisIdTag, axisMcTag}});

HistogramConfigSpec histPPt({HistType::kTH2F, {axisP, axisPt}});
HistogramConfigSpec histPTpcInnerParam({HistType::kTH2F, {axisP, axisTPCInnerParam}});
Expand All @@ -265,6 +283,17 @@ struct NchCumulantsId {

HistogramConfigSpec histPtMc({HistType::kTH1F, {axisPt}});

// Register histograms for PID validation
// Register TPC spares per species
hist.add("PIDValidation/tpcSparse_Pi", "p vs tpcNSigmaPi vs idTag vs mcTag (Pion)", histTPCPIDSparse);
hist.add("PIDValidation/tpcSparse_Ka", "p vs tpcNSigmaKa vs idTag vs mcTag (Kaon)", histTPCPIDSparse);
hist.add("PIDValidation/tpcSparse_Pr", "p vs tpcNSigmaPr vs idTag vs mcTag (Proton)", histTPCPIDSparse);

// Register TOF spares per species
hist.add("PIDValidation/tofSparse_Pi", "p vs tofNSigmaPi vs idTag vs mcTag (Pion)", histTOFPIDSparse);
hist.add("PIDValidation/tofSparse_Ka", "p vs tofNSigmaKa vs idTag vs mcTag (Kaon)", histTOFPIDSparse);
hist.add("PIDValidation/tofSparse_Pr", "p vs tofNSigmaPr vs idTag vs mcTag (Proton)", histTOFPIDSparse);

// QA check histos

hist.add("QA/events/preSel/h_VtxZ", "V_{Z}", kTH1D, {axisVtxZ});
Expand Down Expand Up @@ -762,8 +791,44 @@ struct NchCumulantsId {
}
}

template <typename T>
int getMCTag(const T& track)
{
int mcTag = 0;
if (!track.has_mcParticle())
return mcTag;
auto mcPart = track.mcParticle();
int pdgCode = std::abs(mcPart.pdgCode());

if (pdgCode == kPiPlus || pdgCode == kPiMinus)
BITSET(mcTag, MC_BIT_PI);
else if (pdgCode == kKPlus || pdgCode == kKMinus)
BITSET(mcTag, MC_BIT_KA);
else if (pdgCode == kProton || pdgCode == kProtonBar)
BITSET(mcTag, MC_BIT_PR);
else if (pdgCode == kElectron || pdgCode == kPositron)
BITSET(mcTag, MC_BIT_EL);
else if (pdgCode == kDeuteron || pdgCode == -kDeuteron)
BITSET(mcTag, MC_BIT_DE);

return mcTag;
}
//
//______________________________Identification Functions________________________________________________________________
// Victor slection cuts for electron and other rejections
template <typename T>
bool selTrackForId(const T& track)
{
if (cfgIdElRejLowNSigma < track.tpcNSigmaEl() && track.tpcNSigmaEl() < cfgIdElRejHighNSigma &&
std::fabs(track.tpcNSigmaPi()) > cfgIdPiRejNSigma &&
std::fabs(track.tpcNSigmaKa()) > cfgIdKaRejNSigma &&
std::fabs(track.tpcNSigmaPr()) > cfgIdPrRejNSigma) {
return false;
} else {
return true;
}
}

// Pion
template <typename T>
bool selPion(const T& track, int& IdMethod)
Expand Down Expand Up @@ -1137,6 +1202,11 @@ struct NchCumulantsId {
nM += hPtEtaForEffCorrection[kCh][kNeg]->GetBinContent(ptEtaBin);
}

// Reject electrons first
if (cfgDoRejectionForId && !selTrackForId(track)) {
continue;
}

int idMethod;
// pion
if (selPion(track, idMethod)) {
Expand Down Expand Up @@ -1235,17 +1305,34 @@ struct NchCumulantsId {
idMethodKa = kUnidentified;
idMethodPr = kUnidentified;

// Get MC tag
int mcTag = getMCTag(track);

// Reject electrons first
if (cfgDoRejectionForId && !selTrackForId(track)) {
continue;
}

if (selPion(track, idMethodPi)) {
trackIsPion = true;
BITSET(trackIdTag, ID_BIT_PI);
hist.fill(HIST("PIDValidation/tpcSparse_Pi"), track.p(), track.tpcNSigmaPi(), trackIdTag, mcTag);
if (track.hasTOF())
hist.fill(HIST("PIDValidation/tofSparse_Pi"), track.p(), track.tofNSigmaPi(), trackIdTag, mcTag);
}
if (selKaon(track, idMethodKa)) {
trackIsKaon = true;
BITSET(trackIdTag, ID_BIT_KA);
hist.fill(HIST("PIDValidation/tpcSparse_Ka"), track.p(), track.tpcNSigmaKa(), trackIdTag, mcTag);
if (track.hasTOF())
hist.fill(HIST("PIDValidation/tofSparse_Ka"), track.p(), track.tofNSigmaKa(), trackIdTag, mcTag);
}
if (selProton(track, idMethodPr)) {
trackIsProton = true;
BITSET(trackIdTag, ID_BIT_PR);
hist.fill(HIST("PIDValidation/tpcSparse_Pr"), track.p(), track.tpcNSigmaPr(), trackIdTag, mcTag);
if (track.hasTOF())
hist.fill(HIST("PIDValidation/tofSparse_Pr"), track.p(), track.tofNSigmaPr(), trackIdTag, mcTag);
}

if constexpr (analysisType == doPurityProcessing) {
Expand Down Expand Up @@ -1555,6 +1642,8 @@ struct NchCumulantsId {
continue;
int pdg = mcPart.pdgCode();

int mcTag = getMCTag(track);

fillTrackQA<qaTracksPostSel>(track);

trackIsPion = false;
Expand All @@ -1565,17 +1654,41 @@ struct NchCumulantsId {
idMethodKa = kUnidentified;
idMethodPr = kUnidentified;

// Reject electrons first
if (cfgDoRejectionForId && !selTrackForId(track)) {
continue;
}

// Fill separate spares for each species if it passes the cut
if (selPion(track, idMethodPi)) {
trackIsPion = true;
BITSET(trackIdTag, ID_BIT_PI);
// Fill TPC sparse for pion
hist.fill(HIST("PIDValidation/tpcSparse_Pi"), track.p(), track.tpcNSigmaPi(), trackIdTag, mcTag);
// Fill TOF sparse for pion if has TOF
if (track.hasTOF()) {
hist.fill(HIST("PIDValidation/tofSparse_Pi"), track.p(), track.tofNSigmaPi(), trackIdTag, mcTag);
}
}
if (selKaon(track, idMethodKa)) {
trackIsKaon = true;
BITSET(trackIdTag, ID_BIT_KA);
// Fill TPC sparse for kaon
hist.fill(HIST("PIDValidation/tpcSparse_Ka"), track.p(), track.tpcNSigmaKa(), trackIdTag, mcTag);
// Fill TOF sparse for kaon if has TOF
if (track.hasTOF()) {
hist.fill(HIST("PIDValidation/tofSparse_Ka"), track.p(), track.tofNSigmaKa(), trackIdTag, mcTag);
}
}
if (selProton(track, idMethodPr)) {
trackIsProton = true;
BITSET(trackIdTag, ID_BIT_PR);
// Fill TPC sparse for proton
hist.fill(HIST("PIDValidation/tpcSparse_Pr"), track.p(), track.tpcNSigmaPr(), trackIdTag, mcTag);
// Fill TOF sparse for proton if has TOF
if (track.hasTOF()) {
hist.fill(HIST("PIDValidation/tofSparse_Pr"), track.p(), track.tofNSigmaPr(), trackIdTag, mcTag);
}
}

// bool isKnownCharged = (std::abs(pdg) == kPiPlus ||
Expand Down Expand Up @@ -1614,6 +1727,11 @@ struct NchCumulantsId {
nAPiRec += hPtEtaForEffCorrection[kPi][kNeg]->GetBinContent(ptEtaBin);
fillRecoTrackQA<recoAnalysisDir, kPi, kNeg>(recoAnalysis, track);
}
// PID band QA for pions
if (idMethodPi == kTPCidentified)
fillIdentificationQA<qaTracksIdfd, kPi, tpcId>(hist, track);
if (idMethodPi == kTPCTOFidentified)
fillIdentificationQA<qaTracksIdfd, kPi, tpctofId>(hist, track);
} else if (trackIsKaon) {
if (track.sign() > 0) {
nKaRec += hPtEtaForEffCorrection[kKa][kPos]->GetBinContent(ptEtaBin);
Expand All @@ -1622,6 +1740,11 @@ struct NchCumulantsId {
nAKaRec += hPtEtaForEffCorrection[kKa][kNeg]->GetBinContent(ptEtaBin);
fillRecoTrackQA<recoAnalysisDir, kKa, kNeg>(recoAnalysis, track);
}
// PID band QA for kaons
if (idMethodKa == kTPCidentified)
fillIdentificationQA<qaTracksIdfd, kKa, tpcId>(hist, track);
if (idMethodKa == kTPCTOFidentified)
fillIdentificationQA<qaTracksIdfd, kKa, tpctofId>(hist, track);
} else if (trackIsProton) {
if (track.sign() > 0) {
nPrRec += hPtEtaForEffCorrection[kPr][kPos]->GetBinContent(ptEtaBin);
Expand All @@ -1630,6 +1753,11 @@ struct NchCumulantsId {
nAPrRec += hPtEtaForEffCorrection[kPr][kNeg]->GetBinContent(ptEtaBin);
fillRecoTrackQA<recoAnalysisDir, kPr, kNeg>(recoAnalysis, track);
}
// PID band QA for protons
if (idMethodPr == kTPCidentified)
fillIdentificationQA<qaTracksIdfd, kPr, tpcId>(hist, track);
if (idMethodPr == kTPCTOFidentified)
fillIdentificationQA<qaTracksIdfd, kPr, tpctofId>(hist, track);
}
// purity check - check pdg aginst sign
bool purityPion = false;
Expand Down
Loading