Skip to content

[PWGLF] updated data models and table producer for reducing deriveddata size #15177

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 2 commits into from
Feb 26, 2026
+135 −46
Merged

[PWGLF] updated data models and table producer for reducing deriveddata size #15177

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
dec23a2
added few QA histograms
Feb 22, 2026
c3db27a
updated datamodel and tables for reducing size
Feb 26, 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
28 changes: 8 additions & 20 deletions PWGLF/DataModel/ZDCCalTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -34,16 +34,10 @@ DECLARE_SOA_COLUMN(Cent, cent, float);
DECLARE_SOA_COLUMN(Vx, vx, float);
DECLARE_SOA_COLUMN(Vy, vy, float);
DECLARE_SOA_COLUMN(Vz, vz, float);
DECLARE_SOA_COLUMN(ZnaC, znaC, float);
DECLARE_SOA_COLUMN(ZncC, zncC, float);
DECLARE_SOA_COLUMN(ZnaE0, znaE0, float);
DECLARE_SOA_COLUMN(ZnaE1, znaE1, float);
DECLARE_SOA_COLUMN(ZnaE2, znaE2, float);
DECLARE_SOA_COLUMN(ZnaE3, znaE3, float);
DECLARE_SOA_COLUMN(ZncE0, zncE0, float);
DECLARE_SOA_COLUMN(ZncE1, zncE1, float);
DECLARE_SOA_COLUMN(ZncE2, zncE2, float);
DECLARE_SOA_COLUMN(ZncE3, zncE3, float);
DECLARE_SOA_COLUMN(QxA, qxA, float);
DECLARE_SOA_COLUMN(QxC, qxC, float);
DECLARE_SOA_COLUMN(QyA, qyA, float);
DECLARE_SOA_COLUMN(QyC, qyC, float);
} // namespace zdccaltable
DECLARE_SOA_TABLE(ZDCCalTables, "AOD", "ZDCCALTABLE",
zdccaltable::TriggerEventZDC,
Expand All @@ -52,16 +46,10 @@ DECLARE_SOA_TABLE(ZDCCalTables, "AOD", "ZDCCALTABLE",
zdccaltable::Vx,
zdccaltable::Vy,
zdccaltable::Vz,
zdccaltable::ZnaC,
zdccaltable::ZncC,
zdccaltable::ZnaE0,
zdccaltable::ZnaE1,
zdccaltable::ZnaE2,
zdccaltable::ZnaE3,
zdccaltable::ZncE0,
zdccaltable::ZncE1,
zdccaltable::ZncE2,
zdccaltable::ZncE3);
zdccaltable::QxA,
zdccaltable::QxC,
zdccaltable::QyA,
zdccaltable::QyC);
using ZDCCalTable = ZDCCalTables::iterator;
} // namespace o2::aod
#endif // PWGLF_DATAMODEL_ZDCCALTABLES_H_
27 changes: 17 additions & 10 deletions PWGLF/TableProducer/Common/zdcvector.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -176,7 +176,7 @@ struct zdcvector {

if (!bc.has_zdc()) {
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
return;
}

Expand All @@ -190,38 +190,43 @@ struct zdcvector {

if (znaEnergycommon <= 0.0 || zncEnergycommon <= 0.0) {
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
// zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
return;
}

histos.fill(HIST("hEvtSelInfo"), 2.5);

if (znaEnergy[0] <= 0.0 || znaEnergy[1] <= 0.0 || znaEnergy[2] <= 0.0 || znaEnergy[3] <= 0.0) {
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
// zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
return;
}
histos.fill(HIST("hEvtSelInfo"), 3.5);

if (zncEnergy[0] <= 0.0 || zncEnergy[1] <= 0.0 || zncEnergy[2] <= 0.0 || zncEnergy[3] <= 0.0) {
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
// zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
return;
}

histos.fill(HIST("hEvtSelInfo"), 4.5);

if (rctCut.requireRCTFlagChecker && !rctChecker(collision)) {
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
// zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
return;
}

histos.fill(HIST("hEvtSelInfo"), 5.5);

if (additionalEvSel && (!collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
// zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
return;
}

Expand Down Expand Up @@ -255,7 +260,7 @@ struct zdcvector {

if (znaEnergy[iChA] <= 0.0) {
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
return;
} else {
double ampl = gainequal * znaEnergy[iChA];
Expand All @@ -270,7 +275,8 @@ struct zdcvector {
} else {
if (zncEnergy[iChA - 4] <= 0.0) {
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
// zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
return;
} else {
double ampl = gainequal * zncEnergy[iChA - 4];
Expand Down Expand Up @@ -300,7 +306,8 @@ struct zdcvector {
qyZDCA = 0.0;
qyZDCC = 0.0;
triggerevent = false;
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, 0.0, 0.0, 0.0, 0.0);
// zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
return;
}

Expand All @@ -322,7 +329,7 @@ struct zdcvector {

lastRunNumber = currentRunNumber;
}
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, znaEnergycommon, zncEnergycommon, znaEnergy[0], znaEnergy[1], znaEnergy[2], znaEnergy[3], zncEnergy[0], zncEnergy[1], zncEnergy[2], zncEnergy[3]);
zdccaltable(triggerevent, currentRunNumber, centrality, vx, vy, vz, qxZDCA, qxZDCC, qyZDCA, qyZDCC);
}
};

Expand Down
96 changes: 96 additions & 0 deletions PWGLF/Tasks/Strangeness/lambdapolsp.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -190,6 +190,16 @@ struct lambdapolsp {
ConfigurableAxis configbinAxis{"configbinAxis", {VARIABLE_WIDTH, -0.8, -0.4, -0.2, 0, 0.2, 0.4, 0.8}, "BA"};
*/
// ConfigurableAxis configphiAxis{"configphiAxis", {VARIABLE_WIDTH, 0.0, 0.2, 0.4, 0.8, 1.0, 2.0, 2.5, 3.0, 4.0, 5.0, 5.5, 6.28}, "PhiAxis"};

struct : ConfigurableGroup {
Configurable<bool> isQA{"isQA", true, "Flag to fill QA"};
ConfigurableAxis centfineAxis{"centfineAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "Cent V0M"};
ConfigurableAxis vxfineAxis{"vxfineAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "vx fine axis"};
ConfigurableAxis vyfineAxis{"vyfineAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "vy fine axis"};
ConfigurableAxis vzfineAxis{"vzfineAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "vz fine axis"};
ConfigurableAxis qxZDCAxis{"qxZDCAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "qx axis"};
ConfigurableAxis psiAxis{"psiAxis", {VARIABLE_WIDTH, 0.0, 10.0, 40.0, 80.0}, "psi axis"};
} QAgrp;
struct : ConfigurableGroup {
Configurable<bool> requireRCTFlagChecker{"requireRCTFlagChecker", true, "Check event quality in run condition table"};
Configurable<std::string> cfgEvtRCTFlagCheckerLabel{"cfgEvtRCTFlagCheckerLabel", "CBT_hadronPID", "Evt sel: RCT flag checker label"};
Expand Down Expand Up @@ -456,6 +466,32 @@ struct lambdapolsp {
histos.add("hSparseRecLambda", "hSparseRecLambda", HistType::kTHnSparseF, runaxes2, true);
histos.add("hSparseRecAntiLambda", "hSparseRecAntiLambda", HistType::kTHnSparseF, runaxes2, true);

if (QAgrp.isQA) {
histos.add("hCentQxZDCA", "hCentQxZDCA", kTH2F, {{QAgrp.centfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hCentQyZDCA", "hCentQyZDCA", kTH2F, {{QAgrp.centfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hCentQxZDCC", "hCentQxZDCC", kTH2F, {{QAgrp.centfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hCentQyZDCC", "hCentQyZDCC", kTH2F, {{QAgrp.centfineAxis}, {QAgrp.qxZDCAxis}});

histos.add("hvxQxZDCA", "hvxQxZDCA", kTH2F, {{QAgrp.vxfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvxQyZDCA", "hvxQyZDCA", kTH2F, {{QAgrp.vxfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvxQxZDCC", "hvxQxZDCC", kTH2F, {{QAgrp.vxfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvxQyZDCC", "hvxQyZDCC", kTH2F, {{QAgrp.vxfineAxis}, {QAgrp.qxZDCAxis}});

histos.add("hvyQxZDCA", "hvyQxZDCA", kTH2F, {{QAgrp.vyfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvyQyZDCA", "hvyQyZDCA", kTH2F, {{QAgrp.vyfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvyQxZDCC", "hvyQxZDCC", kTH2F, {{QAgrp.vyfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvyQyZDCC", "hvyQyZDCC", kTH2F, {{QAgrp.vyfineAxis}, {QAgrp.qxZDCAxis}});

histos.add("hvzQxZDCA", "hvzQxZDCA", kTH2F, {{QAgrp.vzfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvzQyZDCA", "hvzQyZDCA", kTH2F, {{QAgrp.vzfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvzQxZDCC", "hvzQxZDCC", kTH2F, {{QAgrp.vzfineAxis}, {QAgrp.qxZDCAxis}});
histos.add("hvzQyZDCC", "hvzQyZDCC", kTH2F, {{QAgrp.vzfineAxis}, {QAgrp.qxZDCAxis}});

histos.add("PsiZDCC", "PsiZDCC", kTH2F, {QAgrp.centfineAxis, QAgrp.psiAxis});
histos.add("PsiZDCA", "PsiZDCA", kTH2F, {QAgrp.centfineAxis, QAgrp.psiAxis});
histos.add("PsiZDC", "PsiZDC", kTH2F, {QAgrp.centfineAxis, QAgrp.psiAxis});
}

ccdb->setURL(cfgCcdbParam.cfgURL);
ccdbApi.init("http://alice-ccdb.cern.ch");
ccdb->setCaching(true);
Expand Down Expand Up @@ -939,6 +975,10 @@ struct lambdapolsp {
// currentRunNumber = collision.foundBC_as<BCsRun3>().runNumber();
auto bc = collision.foundBC_as<BCsRun3>();

auto vz = collision.vz();
auto vx = collision.vx();
auto vy = collision.vy();

auto qxZDCA = collision.qxZDCA();
auto qxZDCC = collision.qxZDCC();
auto qyZDCA = collision.qyZDCA();
Expand All @@ -963,6 +1003,32 @@ struct lambdapolsp {
psiZDC = psiZDCC - psiZDCA;
}*/

if (QAgrp.isQA) {
histos.fill(HIST("hCentQxZDCA"), centrality, modqxZDCA);
histos.fill(HIST("hCentQyZDCA"), centrality, modqyZDCA);
histos.fill(HIST("hCentQxZDCC"), centrality, modqxZDCC);
histos.fill(HIST("hCentQyZDCC"), centrality, modqyZDCC);

histos.fill(HIST("hvxQxZDCA"), vx, modqxZDCA);
histos.fill(HIST("hvxQyZDCA"), vx, modqyZDCA);
histos.fill(HIST("hvxQxZDCC"), vx, modqxZDCC);
histos.fill(HIST("hvxQyZDCC"), vx, modqyZDCC);

histos.fill(HIST("hvyQxZDCA"), vy, modqxZDCA);
histos.fill(HIST("hvyQyZDCA"), vy, modqyZDCA);
histos.fill(HIST("hvyQxZDCC"), vy, modqxZDCC);
histos.fill(HIST("hvyQyZDCC"), vy, modqyZDCC);

histos.fill(HIST("hvzQxZDCA"), vz, modqxZDCA);
histos.fill(HIST("hvzQyZDCA"), vz, modqyZDCA);
histos.fill(HIST("hvzQxZDCC"), vz, modqxZDCC);
histos.fill(HIST("hvzQyZDCC"), vz, modqyZDCC);

histos.fill(HIST("PsiZDCA"), centrality, psiZDCA);
histos.fill(HIST("PsiZDCC"), centrality, psiZDCC);
histos.fill(HIST("PsiZDC"), centrality, psiZDC);
}

histos.fill(HIST("hCentrality"), centrality);
if (!checkwithpub) {
// histos.fill(HIST("hVtxZ"), collision.posZ());
Expand Down Expand Up @@ -1260,6 +1326,10 @@ struct lambdapolsp {
else if (centestim == 3)
centrality = collision.centFV0A();

auto vz = collision.posZ();
auto vx = collision.posX();
auto vy = collision.posY();

auto runnumber = collision.runNumber();
if (!collision.triggereventsp()) { // provided by StraZDCSP
return;
Expand Down Expand Up @@ -1319,6 +1389,32 @@ struct lambdapolsp {

auto psiZDC = TMath::ATan2((modqyZDCC - modqyZDCA), (modqxZDCC - modqxZDCA)); // full event plane

if (QAgrp.isQA) {
histos.fill(HIST("hCentQxZDCA"), centrality, modqxZDCA);
histos.fill(HIST("hCentQyZDCA"), centrality, modqyZDCA);
histos.fill(HIST("hCentQxZDCC"), centrality, modqxZDCC);
histos.fill(HIST("hCentQyZDCC"), centrality, modqyZDCC);

histos.fill(HIST("hvxQxZDCA"), vx, modqxZDCA);
histos.fill(HIST("hvxQyZDCA"), vx, modqyZDCA);
histos.fill(HIST("hvxQxZDCC"), vx, modqxZDCC);
histos.fill(HIST("hvxQyZDCC"), vx, modqyZDCC);

histos.fill(HIST("hvyQxZDCA"), vy, modqxZDCA);
histos.fill(HIST("hvyQyZDCA"), vy, modqyZDCA);
histos.fill(HIST("hvyQxZDCC"), vy, modqxZDCC);
histos.fill(HIST("hvyQyZDCC"), vy, modqyZDCC);

histos.fill(HIST("hvzQxZDCA"), vz, modqxZDCA);
histos.fill(HIST("hvzQyZDCA"), vz, modqyZDCA);
histos.fill(HIST("hvzQxZDCC"), vz, modqxZDCC);
histos.fill(HIST("hvzQyZDCC"), vz, modqyZDCC);

histos.fill(HIST("PsiZDCA"), centrality, psiZDCA);
histos.fill(HIST("PsiZDCC"), centrality, psiZDCC);
histos.fill(HIST("PsiZDC"), centrality, psiZDC);
}

// fill histograms
histos.fill(HIST("hCentrality"), centrality);
if (!checkwithpub) {
Expand Down
30 changes: 14 additions & 16 deletions PWGLF/Tasks/Strangeness/zdccalderived.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -116,8 +116,8 @@ struct zdccalderived {
Configurable<bool> useShift{"useShift", false, "shift histograms"};
Configurable<bool> ispolarization{"ispolarization", false, "Flag to check polarization"};
Configurable<bool> followpub{"followpub", true, "flag to use alphaZDC"};
Configurable<bool> useGainCallib{"useGainCallib", false, "use gain calibration"};
Configurable<bool> useCallibvertex{"useCallibvertex", false, "use calibration for vxy"};
// Configurable<bool> useGainCallib{"useGainCallib", false, "use gain calibration"};
// Configurable<bool> useCallibvertex{"useCallibvertex", false, "use calibration for vxy"};
Configurable<bool> coarse1{"coarse1", false, "RE1"};
Configurable<bool> fine1{"fine1", false, "REfine1"};
Configurable<bool> coarse2{"coarse2", false, "RE2"};
Expand All @@ -133,8 +133,6 @@ struct zdccalderived {
Configurable<bool> useRecentereSp{"useRecentereSp", false, "use Recentering with Sparse or THn"};
Configurable<bool> useRecenterefineSp{"useRecenterefineSp", false, "use fine Recentering with THn"};

Configurable<std::string> confGainPath{"confGainPath", "Users/p/prottay/My/Object/NewPbPbpass4_10092024/gaincallib", "Path to gain calibration"};
Configurable<std::string> confGainPathVxy{"confGainPathVxy", "Users/p/prottay/My/Object/swapcoords/PbPbpass4_20112024/recentervert", "Path to gain calibration for vxy"};
Configurable<std::string> confRecentereSp{"confRecentereSp", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn path for recentering"};
Configurable<std::string> confRecentereSp2{"confRecentereSp2", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn path for recentering 2"};
Configurable<std::string> confRecentereSp3{"confRecentereSp3", "Users/p/prottay/My/Object/Testingwithsparse/NewPbPbpass4_17092024/recenter", "Sparse or THn path for recentering 3"};
Expand Down Expand Up @@ -354,20 +352,20 @@ struct zdccalderived {
auto vx = collision.vx();
auto vy = collision.vy();

float psiZDCC = -99;
float psiZDCA = -99;
auto qxZDCA = 0.0;
auto qxZDCC = 0.0;
auto qyZDCA = 0.0;
auto qyZDCC = 0.0;
auto sumA = 0.0;
auto sumC = 0.0;
double psiZDCC = -99;
double psiZDCA = -99;
auto qxZDCA = collision.qxA();
auto qxZDCC = collision.qxC();
auto qyZDCA = collision.qyA();
auto qyZDCC = collision.qyC();
// auto sumA = 0.0;
// auto sumC = 0.0;

auto timestamps = ccdb->getRunDuration(currentRunNumber, true); /// fatalise if timestamps are not found
int64_t sorTimestamp = timestamps.first; // timestamp of the SOR/SOX/STF in ms
int64_t eorTimestamp = timestamps.second; // timestamp of the EOR/EOX/ETF in ms
int64_t ts = eorTimestamp / 2 + sorTimestamp / 2; // timestamp of the middle of the run

/*
std::array<float, 4> znaEnergy = {
collision.znaE0(),
collision.znaE1(),
Expand Down Expand Up @@ -442,14 +440,14 @@ struct zdccalderived {
qyZDCC = 0.0;
return;
}

*/
histos.fill(HIST("hEvtSelInfo"), 8.5);
histos.fill(HIST("hCentrality"), centrality);
histos.fill(HIST("Vz"), vz);

histos.fill(HIST("AvgVxy"), 0.5, vx);
histos.fill(HIST("AvgVxy"), 1.5, vy);

/*
if (useCallibvertex && (currentRunNumber != lastRunNumber)) {
gainprofilevxy = ccdb->getForTimeStamp<TProfile>(confGainPathVxy.value, ts);
}
Expand All @@ -458,7 +456,7 @@ struct zdccalderived {
vx = vx - gainprofilevxy->GetBinContent(1);
vy = vy - gainprofilevxy->GetBinContent(2);
}

*/
bool res = 0;
bool resfine = 0;
int check = 1;
Expand Down
Loading