Skip to content

[PWGLF] added cosphipsivseta and option to resocorr run by run along with c factor for GE #15533

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
romainschotter merged 3 commits into from
Mar 27, 2026
+155 −21
Merged

[PWGLF] added cosphipsivseta and option to resocorr run by run along with c factor for GE #15533

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
88c5c11
added cosphipsivseta and option to resocorr run by run
Mar 26, 2026
98e3b20
added c factor for gain equalization
Mar 27, 2026
2b87c99
c factor over gain equalized energy
Mar 27, 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
60 changes: 56 additions & 4 deletions PWGLF/TableProducer/Common/zdcvector.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -85,6 +85,7 @@ struct zdcvector {
Configurable<float> cfgCutCentralityMin{"cfgCutCentralityMin", 0.0f, "Centrality cut Min"};
Configurable<bool> additionalEvSel{"additionalEvSel", false, "additionalEvSel"};
Configurable<bool> usemem{"usemem", true, "usemem"};
Configurable<bool> usecfactor{"usecfactor", false, "use c factor"};

struct : ConfigurableGroup {
Configurable<int> vzFineNbins{"vzFineNbins", 20, "Number of bins in Vz fine histograms"};
Expand Down Expand Up @@ -119,6 +120,7 @@ struct zdcvector {
AxisSpec vzfineAxis = {configbins.vzFineNbins, configbins.lfinebinVz, configbins.hfinebinVz, "vzfine"};
AxisSpec centfineAxis = {configbins.centFineNbins, configbins.lfinebinCent, configbins.hfinebinCent, "V0M (%) fine"};
AxisSpec VxyAxis = {2, 0, 2, "Vxy"};
// AxisSpec EAxis = {500, 0, 1000, "Eaxis"};

histos.add("htpcnsigmapi", "htpcnsigmapi", kTH1F, {{50, -10, 10.0}});
histos.add("hEvtSelInfo", "hEvtSelInfo", kTH1F, {{10, 0, 10.0}});
Expand All @@ -128,6 +130,14 @@ struct zdcvector {
histos.add("ZDCAmp", "ZDCAmp", kTProfile2D, {channelZDCAxis, vzfineAxis});
histos.add("ZDCAmpCommon", "ZDCAmpCommon", kTProfile2D, {{2, 0.0, 2.0}, vzfineAxis});
histos.add("AvgVxy", "AvgVxy", kTProfile, {VxyAxis});
/*histos.add("hZNA0", "hZNA0", kTH2D, {{centfineAxis}, {EAxis}}, true);
histos.add("hZNA1", "hZNA1", kTH2D, {{centfineAxis}, {EAxis}}, true);
histos.add("hZNA2", "hZNA2", kTH2D, {{centfineAxis}, {EAxis}}, true);
histos.add("hZNA3", "hZNA3", kTH2D, {{centfineAxis}, {EAxis}}, true);
histos.add("hZNC0", "hZNC0", kTH2D, {{centfineAxis}, {EAxis}}, true);
histos.add("hZNC1", "hZNC1", kTH2D, {{centfineAxis}, {EAxis}}, true);
histos.add("hZNC2", "hZNC2", kTH2D, {{centfineAxis}, {EAxis}}, true);
histos.add("hZNC3", "hZNC3", kTH2D, {{centfineAxis}, {EAxis}}, true);*/

ccdb->setURL(cfgCcdbParam.cfgURL);
ccdbApi.init("http://alice-ccdb.cern.ch");
Expand All @@ -141,6 +151,11 @@ struct zdcvector {
TH2D* gainprofile;
TProfile* gainprofilevxy;

// int lastRunNumberTimeRec = -999;
// for time since start of run
// int runForStartTime = -999;
// uint64_t runStartTime = 0;

using MyCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::FT0sCorrected, aod::CentFT0Cs>;
using AllTrackCandidates = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTPCFullPr, aod::pidTPCFullKa>;

Expand Down Expand Up @@ -180,13 +195,22 @@ struct zdcvector {
return;
}

// -------- define time since run start --------
/*uint64_t ts = bc.timestamp();
if (currentRunNumber != runForStartTime) {
runStartTime = ts;
runForStartTime = currentRunNumber;
}
double tsec = 1.e-3 * static_cast<double>(ts - runStartTime);
*/
histos.fill(HIST("hEvtSelInfo"), 1.5);

auto zdc = bc.zdc();
auto zncEnergy = zdc.energySectorZNC();
auto znaEnergy = zdc.energySectorZNA();
auto zncEnergycommon = zdc.energyCommonZNC();
auto znaEnergycommon = zdc.energyCommonZNA();
auto beamEne = 5.36 * 0.5;

if (znaEnergycommon <= 0.0 || zncEnergycommon <= 0.0) {
triggerevent = false;
Expand Down Expand Up @@ -222,6 +246,14 @@ struct zdcvector {
}

histos.fill(HIST("hEvtSelInfo"), 5.5);
/*histos.fill(HIST("hZNA0"), centrality, znaEnergy[0]);
histos.fill(HIST("hZNA1"), centrality, znaEnergy[1]);
histos.fill(HIST("hZNA2"), centrality, znaEnergy[2]);
histos.fill(HIST("hZNA3"), centrality, znaEnergy[3]);
histos.fill(HIST("hZNC0"), centrality, zncEnergy[0]);
histos.fill(HIST("hZNC1"), centrality, zncEnergy[1]);
histos.fill(HIST("hZNC2"), centrality, zncEnergy[2]);
histos.fill(HIST("hZNC3"), centrality, zncEnergy[3]);*/

if (additionalEvSel && (!collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
triggerevent = false;
Expand All @@ -244,6 +276,8 @@ struct zdcvector {
auto alphaZDC = 0.395;
constexpr double x[4] = {-1.75, 1.75, -1.75, 1.75};
constexpr double y[4] = {-1.75, -1.75, 1.75, 1.75};
double zncEnergycommonsum = 0.0;
double znaEnergycommonsum = 0.0;

histos.fill(HIST("ZDCAmpCommon"), 0.5, vz, znaEnergycommon);
histos.fill(HIST("ZDCAmpCommon"), 1.5, vz, zncEnergycommon);
Expand All @@ -252,6 +286,7 @@ struct zdcvector {
constexpr std::size_t ntow = 8;
for (std::size_t iChA = 0; iChA < ntow; iChA++) {
auto chanelid = iChA;
gainequal = 1.0;
if (useGainCallib && gainprofile) {
gainequal = gainprofile->GetBinContent(gainprofile->FindBin(vz + 0.00001, chanelid + 0.5));
}
Expand All @@ -264,6 +299,7 @@ struct zdcvector {
return;
} else {
double ampl = gainequal * znaEnergy[iChA];
znaEnergycommonsum += ampl;
if (followpub) {
ampl = std::pow(ampl, alphaZDC);
}
Expand All @@ -280,6 +316,7 @@ struct zdcvector {
return;
} else {
double ampl = gainequal * zncEnergy[iChA - 4];
zncEnergycommonsum += ampl;
if (followpub) {
ampl = std::pow(ampl, alphaZDC);
}
Expand All @@ -291,13 +328,26 @@ struct zdcvector {
}
}

auto cZNC = 1.0;
auto cZNA = 1.0;

if (sumA > 0) {
qxZDCA = qxZDCA / sumA;
qyZDCA = qyZDCA / sumA;
float nSpecnA = znaEnergycommonsum / beamEne;
if (usecfactor)
cZNA = 1.89358 - 0.71262 / (nSpecnA + 0.71789);
else
cZNA = 1.0;
qxZDCA = cZNA * (qxZDCA / sumA);
qyZDCA = cZNA * (qyZDCA / sumA);
}
if (sumC > 0) {
qxZDCC = qxZDCC / sumC;
qyZDCC = qyZDCC / sumC;
float nSpecnC = zncEnergycommonsum / beamEne;
if (usecfactor)
cZNC = 1.89358 - 0.71262 / (nSpecnC + 0.71789);
else
cZNC = 1.0;
qxZDCC = cZNC * (qxZDCC / sumC);
qyZDCC = cZNC * (qyZDCC / sumC);
}

if (sumA <= 1e-4 || sumC <= 1e-4) {
Expand All @@ -322,6 +372,8 @@ struct zdcvector {
gainprofilevxy = ccdb->getForTimeStamp<TProfile>(confGainPathVxy.value, bc.timestamp());
}

// LOG(info)<<"*****time stamp is:"<<tsec;

if (useCallibvertex) {
vx = vx - gainprofilevxy->GetBinContent(1);
vy = vy - gainprofilevxy->GetBinContent(2);
Expand Down
Loading
Loading