Skip to content

Commit 9bd1d13

Browse files
committed
Add possibility to recenter without gain eq
1 parent 871f816 commit 9bd1d13

File tree

1 file changed

+17
-12
lines changed

1 file changed

+17
-12
lines changed

PWGCF/Flow/TableProducer/zdcQVectors.cxx

Lines changed: 17 additions & 12 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

@@ -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)