Skip to content

Commit 981ff92

Browse files
hernasabSabrina Hernandez
hernasab
and
Sabrina Hernandez
authored
[PWGCF] remove obsolete variables (#11614)
Co-authored-by: Sabrina Hernandez <sabrina.gabrielle.hernandez@cern.ch>
1 parent 7b5501c commit 981ff92

File tree

1 file changed

+100
-150
lines changed

1 file changed

+100
-150
lines changed

‎PWGCF/Flow/Tasks/flowZdcTask.cxx‎

Lines changed: 100 additions & 150 deletions
Original file line numberDiff line numberDiff line change
@@ -189,10 +189,6 @@ struct FlowZdcTask {
189189
{
190190
// define axes
191191
const AxisSpec axisCounter{1, 0, +1, ""};
192-
const AxisSpec axisQ{100, -1, 1, "Q"};
193-
const AxisSpec axisQZNA{100, -1, 1, "Q"};
194-
const AxisSpec axisREQ{100, -1, 1, "real Q"};
195-
const AxisSpec axisIMQ{100, -1, 1, "imag Q"};
196192
const AxisSpec axisEvent{18, 0.5, 18.5, ""};
197193
const AxisSpec axisZpos{48, -12., 12., "Vtx_{z} (cm)"};
198194

@@ -203,41 +199,22 @@ struct FlowZdcTask {
203199
// create histograms
204200
histos.add("hEventCounter", "Event counter", kTH1F, {axisEvent});
205201
histos.add("zPos", ";;Entries;", kTH1F, {axisZpos});
206-
histos.add("etaHistogram", "etaHistogram", kTH1F, {axisEta});
207-
histos.add("ptHistogram", "ptHistogram", kTH1F, {axisPt});
208202

209203
histos.add("eventCounter", "eventCounter", kTH1F, {axisCounter});
210204
histos.add("centHistogram", "centHistogram", kTH1F, {axisCent});
211205
histos.add("multHistogram", "multHistogram", kTH1F, {axisMultiplicity});
212-
histos.add("multvsCent", "centrality vs multiplicity", kTH2F, {axisCent, axisMultiplicity});
213206
histos.add("phiHistogram", "phiHistogram", kTH1F, {axisPhi});
214-
histos.add("TPCmultiplicity", "TPCmultiplicity", kTH1F, {axisMultTpc});
215-
216-
histos.add("REqHistogram", "REqHistogram", kTH1F, {axisQ});
217-
histos.add("IMqHistogram", "IMqHistogram", kTH1F, {axisQ});
218-
219-
histos.add("REqHistogramZNA", "REqHistogramZNA", kTH1F, {axisQZNA});
220-
histos.add("IMqHistogramZNA", "IMqHistogramZNA", kTH1F, {axisQZNA});
221-
222-
histos.add("REqHistogramZNC", "REqHistogramZNC", kTH1F, {axisQZNA});
223-
histos.add("IMqHistogramZNC", "IMqHistogramZNC", kTH1F, {axisQZNA});
224-
225-
histos.add("EnergyZNA", "ZNA Sector Energy", kTH1F, {axisEnergy});
226-
histos.add("EnergyZNC", "ZNC Sector Energy", kTH1F, {axisEnergy});
227-
histos.add("hCentFT0C", "FT0C Centrality Distribution", kTH1F, {{100, 0, 105}});
228207
histos.add("hZNvsFT0Ccent",
229208
"ZN Energy vs FT0C Centrality",
230209
kTH2F,
231210
{axisCent, axisZN});
232-
233211
histos.add("hZPvsFT0Ccent",
234212
"ZP Energy vs FT0C Centrality;Centrality [%];ZP Energy",
235213
kTH2F,
236214
{axisCent, axisZP});
237215
histos.add("hNchvsNPV", ";NPVTracks (|#eta|<1);N_{ch} (|#eta|<0.8);",
238216
kTH2F,
239217
{{{nBinsNch, -0.5, maxNch}, {nBinsNch, -0.5, maxNch}}});
240-
histos.add("revsimag", "revsimag", kTH2F, {axisREQ, axisIMQ}); // for q vector recentering
241218
histos.add("hYield", "Nch vs pT", kTH2F, {axisMultiplicity, axisPt});
242219
histos.add("hGlobalTracks", "hGlobalTracks", kTH1F, {axisMultiplicity});
243220
// event selection steps
@@ -268,17 +245,13 @@ struct FlowZdcTask {
268245
if (doprocessZdcCollAssoc) { // Check if the process function for ZDCCollAssoc is enabled
269246
histos.add("ZNAcoll", "ZNAcoll; ZNA amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZn}}});
270247
histos.add("ZNCcoll", "ZNCcoll; ZNC amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZn}}});
248+
histos.add("ZPCcoll", "ZPCcoll; ZPC amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZn}}});
249+
histos.add("ZPAcoll", "ZPAcoll; ZPA amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZn}}});
271250
histos.add("ZEM1coll", "ZEM1coll; ZEM1 amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZem}}});
272251
histos.add("ZEM2coll", "ZEM2coll; ZEM2 amplitude; Entries", {HistType::kTH1F, {{nBinsAmp, -0.5, maxZem}}});
273252
histos.add("ZNvsZEMcoll", "ZNvsZEMcoll; ZEM; ZNA+ZNC", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZem}, {nBinsAmp, -0.5, 2. * maxZn}}}});
274253
histos.add("ZNAvsZNCcoll", "ZNAvsZNCcoll; ZNC; ZNA", {HistType::kTH2F, {{{nBinsAmp, -0.5, maxZn}, {nBinsAmp, -0.5, maxZn}}}});
275254

276-
histos.add("RealQHistogramZNA", "RealQHistogramZNA", kTH1F, {axisQZNA});
277-
histos.add("ImagQHistogramZNA", "ImagQHistogramZNA", kTH1F, {axisQZNA});
278-
histos.add("RealQHistogramZNC", "RealQHistogramZNC", kTH1F, {axisQZNA});
279-
histos.add("ImagQHistogramZNC", "ImagQHistogramZNC", kTH1F, {axisQZNA});
280-
281-
histos.add("Acorrelations", "Acorrelations", kTH2F, {{axisQZNA}, {axisQZNA}});
282255
histos.add("SPAngleZNA", "Spectator Plane Angle ZNA;Angle (radians);Entries", {HistType::kTH1F, {{100, -o2::constants::math::PI, o2::constants::math::PI}}});
283256
histos.add("SPAngleZNC", "Spectator Plane Angle ZNC;Angle (radians);Entries", {HistType::kTH1F, {{100, -o2::constants::math::PI, o2::constants::math::PI}}});
284257

@@ -351,15 +324,6 @@ struct FlowZdcTask {
351324
histos.add("ZNVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNA+ZNC;", kTH2F, {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
352325
histos.add("ZNDifVsNch", ";#it{N}_{ch} (|#eta|<0.8);ZNA-ZNC;", kTH2F, {{{nBinsNch, minNch, maxNch}, {100, -50., 50.}}});
353326
}
354-
LOG(info) << "\tccdbNoLaterThan=" << ccdbNoLaterThan.value;
355-
LOG(info) << "\tapplyEff=" << applyEff.value;
356-
LOG(info) << "\tpaTH=" << paTH.value;
357-
LOG(info) << "\tuseMidRapNchSel=" << useMidRapNchSel.value;
358-
LOG(info) << "\tpaTHmeanNch=" << paTHmeanNch.value;
359-
LOG(info) << "\tpaTHsigmaNch=" << paTHsigmaNch.value;
360-
LOG(info) << "\tminPt=" << minPt.value;
361-
LOG(info) << "\tmaxPt=" << maxPt.value;
362-
LOG(info) << "\tmaxPtSpectra=" << maxPtSpectra.value;
363327

364328
ccdb->setURL("http://alice-ccdb.cern.ch");
365329
// Enabling object caching, otherwise each call goes to the CCDB server
@@ -491,6 +455,17 @@ struct FlowZdcTask {
491455
float tZPC{zdc.timeZPC()};
492456
float tZDCdif{tZNC + tZPC - tZNA - tZPA};
493457
float tZDCsum{tZNC + tZPC + tZNA + tZPA};
458+
const double normT0M{(aT0A + aT0C) / 100.};
459+
float znA = zdc.amplitudeZNA() / cfgCollisionEnergy;
460+
float znC = zdc.amplitudeZNC() / cfgCollisionEnergy;
461+
float zpA = zdc.amplitudeZPA() / cfgCollisionEnergy;
462+
float zpC = zdc.amplitudeZPC() / cfgCollisionEnergy;
463+
float aZEM1{zdc.amplitudeZEM1()};
464+
float aZEM2{zdc.amplitudeZEM2()};
465+
float sumZEMs{aZEM1 + aZEM2};
466+
float tZEM1{zdc.timeZEM1()};
467+
float tZEM2{zdc.timeZEM2()};
468+
float sumZNs{znA + znC};
494469

495470
// TDC cut
496471
if (isTDCcut) {
@@ -500,10 +475,6 @@ struct FlowZdcTask {
500475
histos.fill(HIST("hEventCounter"), EvCutLabel::Tdc);
501476
}
502477

503-
float aZEM1{zdc.amplitudeZEM1()};
504-
float aZEM2{zdc.amplitudeZEM2()};
505-
float sumZEMs{aZEM1 + aZEM2};
506-
507478
// ZEM cut
508479
if (isZEMcut) {
509480
if (sumZEMs < zemCut) {
@@ -512,29 +483,22 @@ struct FlowZdcTask {
512483
histos.fill(HIST("hEventCounter"), EvCutLabel::Zem);
513484
}
514485

515-
const double normT0M{(aT0A + aT0C) / 100.};
516-
float znA = zdc.amplitudeZNA() / cfgCollisionEnergy;
517-
float znC = zdc.amplitudeZNC() / cfgCollisionEnergy;
518-
float zpA = zdc.amplitudeZPA() / cfgCollisionEnergy;
519-
float zpC = zdc.amplitudeZPC() / cfgCollisionEnergy;
520-
521-
float tZEM1{zdc.timeZEM1()};
522-
float tZEM2{zdc.timeZEM2()};
523-
float sumZNs{znA + znC};
524-
525486
float sumZNC = (zdc.energySectorZNC())[0] + (zdc.energySectorZNC())[1] + (zdc.energySectorZNC())[2] + (zdc.energySectorZNC())[3];
526487
float sumZNA = (zdc.energySectorZNA())[0] + (zdc.energySectorZNA())[1] + (zdc.energySectorZNA())[2] + (zdc.energySectorZNA())[3];
527488
float sumZPC = (zdc.energySectorZPC())[0] + (zdc.energySectorZPC())[1] + (zdc.energySectorZPC())[2] + (zdc.energySectorZPC())[3];
528489
float sumZPA = (zdc.energySectorZPA())[0] + (zdc.energySectorZPA())[1] + (zdc.energySectorZPA())[2] + (zdc.energySectorZPA())[3];
529490

530-
int itsTracks = 0, glbTracks = 0;
531491
float et = 0., meanpt = 0.;
492+
int itsTracks = 0, glbTracks = 0;
532493
for (const auto& track : tracks) {
533-
if (track.hasITS() && track.itsNCls() >= minITSnCls) {
494+
if (track.hasITS()) {
534495
itsTracks++;
535496
}
536497
// Track Selection
537-
if (track.isGlobalTrack()) {
498+
if (!track.isGlobalTrack()) {
499+
continue;
500+
}
501+
if ((track.pt() < minPt) || (track.pt() > maxPt)) {
538502
continue;
539503
}
540504
glbTracks++;
@@ -637,42 +601,6 @@ struct FlowZdcTask {
637601
}
638602
}
639603

640-
void processQVector(AodCollisions::iterator const& collision, aod::BCsWithTimestamps const&, AodTracks const& tracks, BCsRun3 const& /*bcs*/, aod::Zdcs const& /*zdcsData*/, aod::ZDCMults const& /*zdcMults*/)
641-
{
642-
if (!isEventSelected(collision)) {
643-
return;
644-
}
645-
histos.fill(HIST("eventCounter"), 0.5);
646-
histos.fill(HIST("centHistogram"), collision.centFT0C());
647-
const auto& tracksGrouped = tracksIUWithTPC->sliceByCached(aod::track::collisionId, collision.globalIndex(), cache);
648-
const int multTPC = tracksGrouped.size();
649-
const auto cent = collision.centFT0C();
650-
int globalTracks = tracks.size();
651-
if (globalTracks < 1)
652-
return;
653-
int nTot{0}; // Tracks are already filtered with GlobalTrack || GlobalTrackSDD
654-
for (const auto& track : tracks) {
655-
nTot++;
656-
histos.fill(HIST("etaHistogram"), track.eta());
657-
histos.fill(HIST("phiHistogram"), track.phi());
658-
histos.fill(HIST("ptHistogram"), track.pt());
659-
} // end track loop 1
660-
double pT{0};
661-
for (const auto& track : tracks) {
662-
if (track.tpcNClsCrossedRows() < minTpcNcrossedRows)
663-
continue;
664-
if (std::fabs(track.dcaXY()) > cfgCutDCAxy)
665-
continue;
666-
pT = track.pt();
667-
pT++;
668-
} // end track loop 2
669-
histos.fill(HIST("multvsCent"), cent, nTot);
670-
histos.fill(HIST("hYield"), nTot, pT);
671-
histos.fill(HIST("multHistogram"), nTot);
672-
673-
histos.fill(HIST("TPCmultiplicity"), multTPC);
674-
histos.fill(HIST("hGlobalTracks"), globalTracks);
675-
}
676604
void processZdcCollAssoc(
677605
AodCollisions::iterator const& collision,
678606
AodTracks const& tracks,
@@ -683,10 +611,10 @@ struct FlowZdcTask {
683611
if (!isEventSelected(collision)) {
684612
return;
685613
}
614+
const auto& foundBC = collision.foundBC_as<BCsRun3>();
686615
int nTot = tracks.size();
687616
double ft0aAmp = 0;
688617
double ft0cAmp = 0;
689-
const auto& foundBC = collision.foundBC_as<BCsRun3>();
690618
if (collision.has_foundFT0()) {
691619
auto ft0 = collision.foundFT0();
692620
for (const auto& amplitude : ft0.amplitudeA()) {
@@ -696,66 +624,89 @@ struct FlowZdcTask {
696624
ft0cAmp += amplitude;
697625
}
698626
}
627+
const double normT0M{(ft0aAmp + ft0aAmp) / 100.};
699628
histos.fill(HIST("hFT0AAmp"), ft0aAmp);
700-
histos.fill(HIST("hFT0CAmp"), ft0cAmp);
701-
double ft0mAmp = ft0aAmp + ft0cAmp;
702-
histos.fill(HIST("hFT0MAmp"), ft0mAmp);
703-
if (foundBC.has_zdc()) {
704-
const auto& zdcread = foundBC.zdc();
705-
const auto cent = collision.centFT0C();
706-
707-
// ZDC data and histogram filling
708-
histos.get<TH1>(HIST("ZNAcoll"))->Fill(zdcread.amplitudeZNA());
709-
histos.get<TH1>(HIST("ZNCcoll"))->Fill(zdcread.amplitudeZNC());
710-
histos.get<TH2>(HIST("ZNvsZEMcoll"))->Fill(zdcread.amplitudeZEM1() + zdcread.amplitudeZEM2(), zdcread.amplitudeZNA() + zdcread.amplitudeZNC());
711-
histos.get<TH2>(HIST("ZNAvsZNCcoll"))->Fill(zdcread.amplitudeZNC(), zdcread.amplitudeZNA());
712-
713-
histos.get<TH1>(HIST("ZEM1coll"))->Fill(zdcread.amplitudeZEM1());
714-
histos.get<TH1>(HIST("ZEM2coll"))->Fill(zdcread.amplitudeZEM2());
715-
716-
float sumZNC = (zdcread.energySectorZNC())[0] + (zdcread.energySectorZNC())[1] + (zdcread.energySectorZNC())[2] + (zdcread.energySectorZNC())[3];
717-
float sumZNA = (zdcread.energySectorZNA())[0] + (zdcread.energySectorZNA())[1] + (zdcread.energySectorZNA())[2] + (zdcread.energySectorZNA())[3];
718-
float sumZPC = (zdcread.energySectorZPC())[0] + (zdcread.energySectorZPC())[1] + (zdcread.energySectorZPC())[2] + (zdcread.energySectorZPC())[3];
719-
float sumZPA = (zdcread.energySectorZPA())[0] + (zdcread.energySectorZPA())[1] + (zdcread.energySectorZPA())[2] + (zdcread.energySectorZPA())[3];
720-
float sumZDC = sumZPA + sumZPC + sumZNA + sumZNC;
721-
float sumZEM = zdcread.amplitudeZEM1() + zdcread.amplitudeZEM2();
722-
// ZEM cut
723-
if (isZEMcut) {
724-
if (sumZEM < zemCut) {
725-
return;
726-
}
629+
histos.fill(HIST("hFT0CAmp"), ft0aAmp);
630+
histos.fill(HIST("hFT0MAmp"), normT0M);
631+
632+
if (!foundBC.has_zdc()) {
633+
return;
634+
}
635+
const auto& zdcread = foundBC.zdc();
636+
const auto cent = collision.centFT0C();
637+
638+
// ZDC data and histogram filling
639+
float znA = zdcread.amplitudeZNA();
640+
float znC = zdcread.amplitudeZNC();
641+
float zpA = zdcread.amplitudeZPA();
642+
float zpC = zdcread.amplitudeZPC();
643+
float tZNA{zdcread.timeZNA()};
644+
float tZNC{zdcread.timeZNC()};
645+
float tZPA{zdcread.timeZPA()};
646+
float tZPC{zdcread.timeZPC()};
647+
float tZDCdif{tZNC + tZPC - tZNA - tZPA};
648+
float tZDCsum{tZNC + tZPC + tZNA + tZPA};
649+
float sumZNC = (zdcread.energySectorZNC())[0] + (zdcread.energySectorZNC())[1] + (zdcread.energySectorZNC())[2] + (zdcread.energySectorZNC())[3];
650+
float sumZNA = (zdcread.energySectorZNA())[0] + (zdcread.energySectorZNA())[1] + (zdcread.energySectorZNA())[2] + (zdcread.energySectorZNA())[3];
651+
float sumZPC = (zdcread.energySectorZPC())[0] + (zdcread.energySectorZPC())[1] + (zdcread.energySectorZPC())[2] + (zdcread.energySectorZPC())[3];
652+
float sumZPA = (zdcread.energySectorZPA())[0] + (zdcread.energySectorZPA())[1] + (zdcread.energySectorZPA())[2] + (zdcread.energySectorZPA())[3];
653+
float sumZDC = sumZPA + sumZPC + sumZNA + sumZNC;
654+
float sumZEM = zdcread.amplitudeZEM1() + zdcread.amplitudeZEM2();
655+
znA /= cfgCollisionEnergy;
656+
znC /= cfgCollisionEnergy;
657+
zpA /= cfgCollisionEnergy;
658+
zpC /= cfgCollisionEnergy;
659+
// TDC cut
660+
if (isTDCcut) {
661+
if (std::sqrt(std::pow(tZDCdif, 2.) + std::pow(tZDCsum, 2.)) > tdcCut) {
662+
return;
663+
}
664+
histos.fill(HIST("hEventCounter"), EvCutLabel::Tdc);
665+
}
666+
// ZEM cut
667+
if (isZEMcut) {
668+
if (sumZEM < zemCut) {
669+
return;
727670
}
728-
// common energies
729-
float commonSumZnc = (zdcread.energyCommonZNC());
730-
float commonSumZna = (zdcread.energyCommonZNA());
731-
float commonSumZpc = (zdcread.energyCommonZPC());
732-
float commonSumZpa = (zdcread.energyCommonZPA());
733-
float sumZN = (sumZNC) + (sumZNA);
734-
float sumZP = (sumZPC) + (sumZPA);
735-
736-
histos.fill(HIST("ZNenergy"), sumZN);
737-
histos.fill(HIST("ZPenergy"), sumZP);
738-
histos.fill(HIST("ZNCenergy"), commonSumZnc);
739-
histos.fill(HIST("ZNAenergy"), commonSumZna);
740-
histos.fill(HIST("ZPAenergy"), commonSumZpa);
741-
histos.fill(HIST("ZPCenergy"), commonSumZpc);
742-
histos.fill(HIST("hZNvsFT0Ccent"), cent, sumZN);
743-
histos.fill(HIST("hZPvsFT0Ccent"), cent, sumZP);
744-
histos.fill(HIST("hZNvsFT0CAmp"), ft0cAmp, sumZN);
745-
histos.fill(HIST("hZPvsFT0CAmp"), ft0cAmp, sumZP);
746-
histos.fill(HIST("hZNvsMult"), nTot, sumZN);
747-
histos.fill(HIST("hZPvsMult"), nTot, sumZP);
748-
histos.fill(HIST("hNchvsNPV"), collision.multNTracksPVeta1(), nTot);
749-
750-
float ratioZN = sumZNC / sumZNA;
751-
float ratioZP = sumZPC / sumZPA;
752-
pZNratiovscent->Fill(cent, ratioZN);
753-
pZPratiovscent->Fill(cent, ratioZP);
754-
pZNvsFT0Ccent->Fill(cent, sumZN);
755-
pZPvsFT0Ccent->Fill(cent, sumZP);
756-
757-
histos.get<TH2>(HIST("ZDC_energy_vs_ZEM"))->Fill(sumZEM, sumZDC);
758671
}
672+
// common energies
673+
float commonSumZnc = (zdcread.energyCommonZNC());
674+
float commonSumZna = (zdcread.energyCommonZNA());
675+
float commonSumZpc = (zdcread.energyCommonZPC());
676+
float commonSumZpa = (zdcread.energyCommonZPA());
677+
float sumZN = (sumZNC) + (sumZNA);
678+
float sumZP = (sumZPC) + (sumZPA);
679+
histos.fill(HIST("zPos"), collision.posZ());
680+
histos.fill(HIST("T0Ccent"), collision.centFT0C());
681+
histos.fill(HIST("ZNAcoll"), znA);
682+
histos.fill(HIST("ZNCcoll"), znC);
683+
histos.fill(HIST("ZPAcoll"), zpA);
684+
histos.fill(HIST("ZPCcoll"), zpC);
685+
histos.get<TH2>(HIST("ZNvsZEMcoll"))->Fill(zdcread.amplitudeZEM1() + zdcread.amplitudeZEM2(), zdcread.amplitudeZNA() + zdcread.amplitudeZNC());
686+
histos.get<TH2>(HIST("ZNAvsZNCcoll"))->Fill(zdcread.amplitudeZNC(), zdcread.amplitudeZNA());
687+
histos.get<TH1>(HIST("ZEM1coll"))->Fill(zdcread.amplitudeZEM1());
688+
histos.get<TH1>(HIST("ZEM2coll"))->Fill(zdcread.amplitudeZEM2());
689+
histos.fill(HIST("ZNenergy"), sumZN);
690+
histos.fill(HIST("ZPenergy"), sumZP);
691+
histos.fill(HIST("ZNCenergy"), commonSumZnc);
692+
histos.fill(HIST("ZNAenergy"), commonSumZna);
693+
histos.fill(HIST("ZPAenergy"), commonSumZpa);
694+
histos.fill(HIST("ZPCenergy"), commonSumZpc);
695+
histos.fill(HIST("hZNvsFT0Ccent"), cent, sumZN);
696+
histos.fill(HIST("hZPvsFT0Ccent"), cent, sumZP);
697+
histos.fill(HIST("hZNvsFT0CAmp"), ft0cAmp, sumZN);
698+
histos.fill(HIST("hZPvsFT0CAmp"), ft0cAmp, sumZP);
699+
histos.fill(HIST("hZNvsMult"), nTot, sumZN);
700+
histos.fill(HIST("hZPvsMult"), nTot, sumZP);
701+
histos.fill(HIST("hNchvsNPV"), collision.multNTracksPVeta1(), nTot);
702+
703+
float ratioZN = sumZNC / sumZNA;
704+
float ratioZP = sumZPC / sumZPA;
705+
pZNratiovscent->Fill(cent, ratioZN);
706+
pZPratiovscent->Fill(cent, ratioZP);
707+
pZNvsFT0Ccent->Fill(cent, sumZN);
708+
pZPvsFT0Ccent->Fill(cent, sumZP);
709+
histos.get<TH2>(HIST("ZDC_energy_vs_ZEM"))->Fill(sumZEM, sumZDC);
759710
}
760711

761712
void processCorrelation(CollisionDataTable::iterator const& collision, FilTrackDataTable const& tracks)
@@ -779,7 +730,6 @@ struct FlowZdcTask {
779730

780731
PROCESS_SWITCH(FlowZdcTask, processZdcCollAssoc, "Processing ZDC w. collision association", true);
781732
PROCESS_SWITCH(FlowZdcTask, processQA, "Process QA", true);
782-
PROCESS_SWITCH(FlowZdcTask, processQVector, "Process before recentering", true);
783733
PROCESS_SWITCH(FlowZdcTask, processCorrelation, "Process correlations", true);
784734

785735
}; // end of struct function

0 commit comments

Comments
 (0)