Skip to content

Commit 3cfe609

Browse files
cnkosteralibuild
andauthored
[PWGCF] ZDCQvectors: add possibility to recenter without gain eq (#14552)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent 0a5042c commit 3cfe609

File tree

1 file changed

+18
-13
lines changed

1 file changed

+18
-13
lines changed

PWGCF/Flow/TableProducer/zdcQVectors.cxx

Lines changed: 18 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -156,6 +156,7 @@ struct ZdcQVectors {
156156
O2_DEFINE_CONFIGURABLE(cfgFillHistRegistry, bool, true, "Fill common registry with histograms");
157157
O2_DEFINE_CONFIGURABLE(cfgFillCutAnalysis, bool, true, "Fill cut analysis with histograms");
158158
O2_DEFINE_CONFIGURABLE(cfgFillNothing, bool, false, "Disable ALL Histograms -> ONLY use to reduce memory");
159+
O2_DEFINE_CONFIGURABLE(cfgNoGain, bool, false, "Do not apply gain correction to ZDC energy calibration");
159160

160161
O2_DEFINE_CONFIGURABLE(cfgCCDBdir_Shift, std::string, "Users/c/ckoster/ZDC/LHC23_PbPb_pass5/Shift", "CCDB directory for Shift ZDC");
161162

@@ -724,7 +725,6 @@ struct ZdcQVectors {
724725
std::vector<float> centralities;
725726

726727
auto cent = collision.centFT0C();
727-
centrality = cent;
728728

729729
centralities.push_back(collision.centFT0C());
730730

@@ -751,6 +751,7 @@ struct ZdcQVectors {
751751

752752
v = {collision.posX(), collision.posY(), collision.posZ()};
753753
cents = centralities;
754+
centrality = cent;
754755

755756
const auto& foundBC = collision.foundBC_as<BCsRun3>();
756757
runnumber = foundBC.runNumber();
@@ -992,6 +993,22 @@ struct ZdcQVectors {
992993
}
993994
}
994995

996+
if (cfgFillHistRegistry && !cfgFillNothing) {
997+
registry.get<TProfile>(HIST("QA/before/ZNA_Qx"))->Fill(Form("%d", runnumber), q[0]);
998+
registry.get<TProfile>(HIST("QA/before/ZNA_Qy"))->Fill(Form("%d", runnumber), q[1]);
999+
registry.get<TProfile>(HIST("QA/before/ZNC_Qx"))->Fill(Form("%d", runnumber), q[2]);
1000+
registry.get<TProfile>(HIST("QA/before/ZNC_Qy"))->Fill(Form("%d", runnumber), q[3]);
1001+
1002+
registry.get<TProfile>(HIST("QA/before/ZNA_Qx_noEq"))->Fill(Form("%d", runnumber), qNoEq[0]);
1003+
registry.get<TProfile>(HIST("QA/before/ZNA_Qy_noEq"))->Fill(Form("%d", runnumber), qNoEq[1]);
1004+
registry.get<TProfile>(HIST("QA/before/ZNC_Qx_noEq"))->Fill(Form("%d", runnumber), qNoEq[2]);
1005+
registry.get<TProfile>(HIST("QA/before/ZNC_Qy_noEq"))->Fill(Form("%d", runnumber), qNoEq[3]);
1006+
}
1007+
1008+
if (cfgNoGain) {
1009+
q = qNoEq;
1010+
}
1011+
9951012
if (cal.calibfilesLoaded[1]) {
9961013
v[0] = v[0] - getCorrection<TProfile, kMeanv>(vnames[0].Data());
9971014
v[1] = v[1] - getCorrection<TProfile, kMeanv>(vnames[1].Data());
@@ -1008,18 +1025,6 @@ struct ZdcQVectors {
10081025

10091026
std::vector<double> qRec(q);
10101027

1011-
if (cfgFillHistRegistry && !cfgFillNothing) {
1012-
registry.get<TProfile>(HIST("QA/before/ZNA_Qx"))->Fill(Form("%d", runnumber), q[0]);
1013-
registry.get<TProfile>(HIST("QA/before/ZNA_Qy"))->Fill(Form("%d", runnumber), q[1]);
1014-
registry.get<TProfile>(HIST("QA/before/ZNC_Qx"))->Fill(Form("%d", runnumber), q[2]);
1015-
registry.get<TProfile>(HIST("QA/before/ZNC_Qy"))->Fill(Form("%d", runnumber), q[3]);
1016-
1017-
registry.get<TProfile>(HIST("QA/before/ZNA_Qx_noEq"))->Fill(Form("%d", runnumber), qNoEq[0]);
1018-
registry.get<TProfile>(HIST("QA/before/ZNA_Qy_noEq"))->Fill(Form("%d", runnumber), qNoEq[1]);
1019-
registry.get<TProfile>(HIST("QA/before/ZNC_Qx_noEq"))->Fill(Form("%d", runnumber), qNoEq[2]);
1020-
registry.get<TProfile>(HIST("QA/before/ZNC_Qy_noEq"))->Fill(Form("%d", runnumber), qNoEq[3]);
1021-
}
1022-
10231028
if (cal.atIteration == 0) {
10241029
if (isSelected && cfgFillHistRegistry)
10251030
fillCommonRegistry<kBefore>(q[0], q[1], q[2], q[3], v, centrality);

0 commit comments

Comments
 (0)