2121#include " Common/DataModel/TrackSelectionTables.h"
2222
2323#include " CommonConstants/MathConstants.h"
24- #include " Framework/ASoAHelpers.h"
2524#include " Framework/AnalysisTask.h"
2625#include " Framework/HistogramRegistry.h"
2726#include " Framework/RunningWorkflowInfo.h"
2827#include " Framework/runDataProcessing.h"
2928#include < CCDB/BasicCCDBManager.h>
3029
31- #include " TList.h"
32- #include < TComplex.h>
3330#include < TF1.h>
3431#include < TMath.h>
3532#include < TProfile.h>
36- #include < TProfile2D.h>
37- #include < TRandom3.h>
3833
3934#include < cmath>
40- #include < string>
4135#include < vector>
4236
4337using namespace o2 ;
@@ -54,12 +48,12 @@ struct FlowZdcTask {
5448
5549 O2_DEFINE_CONFIGURABLE (cfgCutVertex, float , 10 .0f , " Accepted z-vertex range" )
5650 Configurable<int > eventSelection{" eventSelection" , 1 , " event selection" };
57- Configurable<float > maxZp{" maxZp" , 125.5 , " Max ZP signal" };
5851 Configurable<float > maxZem{" maxZem" , 3099.5 , " Max ZEM signal" };
5952 // for ZDC info and analysis
6053 Configurable<int > nBinsAmp{" nBinsAmp" , 1025 , " nbinsAmp" };
6154 Configurable<int > nBinsADC{" nBinsADC" , 1000 , " nbinsADC" };
6255 Configurable<float > maxZn{" maxZn" , 125.5 , " Max ZN signal" };
56+ Configurable<float > maxZp{" maxZp" , 125.5 , " Max ZP signal" };
6357 // configs for process QA
6458 Configurable<int > nBinsNch{" nBinsNch" , 2501 , " N bins Nch (|eta|<0.8)" };
6559 Configurable<int > nBinsAmpFT0{" nBinsAmpFT0" , 100 , " N bins FT0 amp" };
@@ -71,14 +65,20 @@ struct FlowZdcTask {
7165 Configurable<float > minNch{" minNch" , 0 , " Min Nch (|eta|<0.8)" };
7266 Configurable<float > maxNch{" maxNch" , 2500 , " Max Nch (|eta|<0.8)" };
7367 Configurable<int > nBinsTDC{" nBinsTDC" , 150 , " nbinsTDC" };
74- Configurable<float > minTdcZn{" minTdcZn" , - 4.0 , " minimum TDC for ZN" };
75- Configurable<float > maxTdcZn{" maxTdcZn" , - 4.0 , " maximum TDC for ZN" };
68+ Configurable<float > minTdcZn{" minTdcZn" , 4.0 , " minimum TDC for ZN" };
69+ Configurable<float > maxTdcZn{" maxTdcZn" , 4.0 , " maximum TDC for ZN" };
7670 Configurable<float > minTdcZp{" minTdcZp" , -4.0 , " minimum TDC for ZP" };
7771 Configurable<float > maxTdcZp{" maxTdcZp" , -4.0 , " maximum TDC for ZP" };
78- Configurable<float > cfgCollisionEnergy{" cfgCollisionEnergy" , 2.68 , " cfgCollisionEnergy" };
7972 Configurable<bool > applyZdcCorrection{" applyZdcCorrection" , false , " Apply ZP correction?" };
8073 Configurable<float > zpaCoeff{" zpaCoeff" , 0 .021f , " Coefficient a in zpa correction" };
8174 Configurable<float > zpcCoeff{" zpcCoeff" , 0 .021f , " Coefficient b in zpc correction" };
75+ Configurable<float > posZcut{" posZcut" , +10.0 , " z-vertex position cut" };
76+ Configurable<float > posYcut{" posYcut" , +10.0 , " y-vertex position cut" };
77+ Configurable<float > posXcut{" posXcut" , +10.0 , " x-vertex position cut" };
78+ Configurable<float > minEta{" minEta" , -0.8 , " minimum eta" };
79+ Configurable<float > maxEta{" maxEta" , +0.8 , " maximum eta" };
80+ Configurable<float > minT0CcentCut{" minT0CcentCut" , 0.0 , " Min T0C Cent. cut" };
81+ Configurable<float > maxT0CcentCut{" maxT0CcentCut" , 90.0 , " Max T0C Cent. cut" };
8282 // event selection
8383 Configurable<bool > isNoCollInTimeRangeStrict{" isNoCollInTimeRangeStrict" , true , " isNoCollInTimeRangeStrict?" };
8484 Configurable<bool > isNoCollInTimeRangeStandard{" isNoCollInTimeRangeStandard" , false , " isNoCollInTimeRangeStandard?" };
@@ -102,17 +102,8 @@ struct FlowZdcTask {
102102 Configurable<float > maxPtSpectra{" maxPtSpectra" , 50 ., " maximum pt of the tracks" };
103103 Configurable<float > zemCut{" zemCut" , 100 ., " ZEM cut" };
104104 // axis configs
105- ConfigurableAxis axisZN{" axisZN" , {5000 , 0 , 500 }, " axisZN" };
106- ConfigurableAxis axisZP{" axisZP" , {5000 , 0 , 500 }, " axisZP" };
107105 ConfigurableAxis axisCent{" axisCent" , {10 , 0 , 100 }, " axisCent" };
108106 ConfigurableAxis binsPt{" binsPt" , {VARIABLE_WIDTH, 0.0 , 0.1 , 0.12 }, " pT binning" };
109- Configurable<float > posZcut{" posZcut" , +10.0 , " z-vertex position cut" };
110- Configurable<float > posYcut{" posYcut" , +10.0 , " y-vertex position cut" };
111- Configurable<float > posXcut{" posXcut" , +10.0 , " x-vertex position cut" };
112- Configurable<float > minEta{" minEta" , -0.8 , " minimum eta" };
113- Configurable<float > maxEta{" maxEta" , +0.8 , " maximum eta" };
114- Configurable<float > minT0CcentCut{" minT0CcentCut" , 0.0 , " Min T0C Cent. cut" };
115- Configurable<float > maxT0CcentCut{" maxT0CcentCut" , 90.0 , " Max T0C Cent. cut" };
116107
117108 using ColEvents = soa::Join<aod::Collisions, aod::EvSels>;
118109 Filter collisionFilter = nabs(aod::collision::posZ) < cfgCutVertex;
@@ -237,10 +228,10 @@ struct FlowZdcTask {
237228 histos.add (" ZPAVsZPC" , " ;ZPC;ZPA;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZp}}});
238229 histos.add (" ZNAVsZPA" , " ;ZPA;ZNA;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZn}}});
239230 histos.add (" ZNCVsZPC" , " ;ZPC;ZNC;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZn}}});
240- histos.add (" SectorZNAVsZNC " , " ;Sector ZNC;Sector ZNA" , kTH2F , {{{nBinsZDC, -0.5 , maxZn}, {nBinsZDC, -0.5 , maxZn}}});
241- histos.add (" SectorZPAVsZPC " , " ;Sector ZPC;Sector ZPA;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZp}}});
242- histos.add (" SectorZNAVsZPA " , " ;Sector ZPA;Sector ZNA;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZn}}});
243- histos.add (" SectorZNCVsZPC " , " ;Sector ZPC;Sector ZNC;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZn}}});
231+ histos.add (" CommonZNAVsZNC " , " ;Common ZNC;Common ZNA" , kTH2F , {{{nBinsZDC, -0.5 , maxZn}, {nBinsZDC, -0.5 , maxZn}}});
232+ histos.add (" CommonZPAVsZPC " , " ;Common ZPC;Common ZPA;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZp}}});
233+ histos.add (" CommonZNAVsZPA " , " ;Common ZPA;Common ZNA;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZn}}});
234+ histos.add (" CommonZNCVsZPC " , " ;Common ZPC;Common ZNC;" , kTH2F , {{{nBinsZDC, -0.5 , maxZp}, {nBinsZDC, -0.5 , maxZn}}});
244235 histos.add (" ZNASector" , " ;ZNA;Entries;" , kTH1F , {{nBinsZDC, -0.5 , maxZn}});
245236 histos.add (" ZPASector" , " ;ZPA;Entries;" , kTH1F , {{nBinsZDC, -0.5 , maxZp}});
246237 histos.add (" ZNCSector" , " ;ZNC;Entries;" , kTH1F , {{nBinsZDC, -0.5 , maxZn}});
@@ -276,14 +267,6 @@ struct FlowZdcTask {
276267 histos.add (" ZNAVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);ZNA;" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
277268 histos.add (" ZNVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);ZNA+ZNC;" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
278269 histos.add (" ZNDifVsNch" , " ;#it{N}_{ch} (|#eta|<0.8);ZNA-ZNC;" , kTH2F , {{{nBinsNch, minNch, maxNch}, {100 , -50 ., 50 .}}});
279- histos.add (" ZPAvsCent" , " ;centFT0C;ZPA" , kTH2F , {{{axisCent}, {nBinsZDC, -0.5 , maxZp}}});
280- histos.add (" ZPCvsCent" , " ;centFT0C;ZPC" , kTH2F , {{{axisCent}, {nBinsZDC, -0.5 , maxZp}}});
281- histos.add (" hZPASectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZPA Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZp}}});
282- histos.add (" hZPCSectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZPC Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZp}}});
283- histos.add (" hZNASectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZNA Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
284- histos.add (" hZNCSectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);ZNC Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
285- histos.add (" hZPSectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);(ZPA + ZPC) Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZp}}});
286- histos.add (" hZNSectorvsGlbTrack" , " ;Global Tracks (ITS + TPC);(ZNA + ZNC) Sector Energy" , kTH2F , {{{nBinsNch, minNch, maxNch}, {nBinsZDC, minNch, maxZn}}});
287270 }
288271 if (doprocessZdc) {
289272 histos.add (" ampZna" , " ;ZNA Amplitude;Entries;" , kTH1F , {{nBinsZDC, -0.5 , maxZn}});
@@ -417,7 +400,6 @@ struct FlowZdcTask {
417400 void processQA (ColEvSels::iterator const & collision, BCsRun3 const & /* bcs*/ , aod::Zdcs const & /* zdcsData*/ , aod::FV0As const & /* fv0as*/ , aod::FT0s const & /* ft0s*/ , TheFilteredTracks const & tracks)
418401 {
419402 const auto & foundBC = collision.foundBC_as <BCsRun3>();
420- const auto cent = collision.centFT0C ();
421403 if (!isEventSelected (collision)) {
422404 return ;
423405 }
@@ -447,18 +429,18 @@ struct FlowZdcTask {
447429 float tZPA{zdc.timeZPA ()};
448430 float tZPC{zdc.timeZPC ()};
449431 const double normT0M{(aT0A + aT0C) / 100 .};
450- float znA = zdc.amplitudeZNA () / cfgCollisionEnergy ;
451- float znC = zdc.amplitudeZNC () / cfgCollisionEnergy ;
452- float zpA = zdc.amplitudeZPA () / cfgCollisionEnergy ;
453- float zpC = zdc.amplitudeZPC () / cfgCollisionEnergy ;
454- float commonSumZnc = zdc.energyCommonZNC () / cfgCollisionEnergy ;
455- float commonSumZna = zdc.energyCommonZNA () / cfgCollisionEnergy ;
456- float commonSumZpc = zdc.energyCommonZPC () / cfgCollisionEnergy ;
457- float commonSumZpa = zdc.energyCommonZPA () / cfgCollisionEnergy ;
458- float sumZNC = ((zdc.energySectorZNC ())[0 ] + (zdc.energySectorZNC ())[1 ] + (zdc.energySectorZNC ())[2 ] + (zdc.energySectorZNC ())[3 ]) / cfgCollisionEnergy ;
459- float sumZNA = ((zdc.energySectorZNA ())[0 ] + (zdc.energySectorZNA ())[1 ] + (zdc.energySectorZNA ())[2 ] + (zdc.energySectorZNA ())[3 ]) / cfgCollisionEnergy ;
460- float sumZPC = ((zdc.energySectorZPC ())[0 ] + (zdc.energySectorZPC ())[1 ] + (zdc.energySectorZPC ())[2 ] + (zdc.energySectorZPC ())[3 ]) / cfgCollisionEnergy ;
461- float sumZPA = ((zdc.energySectorZPA ())[0 ] + (zdc.energySectorZPA ())[1 ] + (zdc.energySectorZPA ())[2 ] + (zdc.energySectorZPA ())[3 ]) / cfgCollisionEnergy ;
432+ float znA = zdc.amplitudeZNA ();
433+ float znC = zdc.amplitudeZNC ();
434+ float zpA = zdc.amplitudeZPA ();
435+ float zpC = zdc.amplitudeZPC ();
436+ float commonSumZnc = zdc.energyCommonZNC ();
437+ float commonSumZna = zdc.energyCommonZNA ();
438+ float commonSumZpc = zdc.energyCommonZPC ();
439+ float commonSumZpa = zdc.energyCommonZPA ();
440+ float sumZNC = ((zdc.energySectorZNC ())[0 ] + (zdc.energySectorZNC ())[1 ] + (zdc.energySectorZNC ())[2 ] + (zdc.energySectorZNC ())[3 ]);
441+ float sumZNA = ((zdc.energySectorZNA ())[0 ] + (zdc.energySectorZNA ())[1 ] + (zdc.energySectorZNA ())[2 ] + (zdc.energySectorZNA ())[3 ]);
442+ float sumZPC = ((zdc.energySectorZPC ())[0 ] + (zdc.energySectorZPC ())[1 ] + (zdc.energySectorZPC ())[2 ] + (zdc.energySectorZPC ())[3 ]);
443+ float sumZPA = ((zdc.energySectorZPA ())[0 ] + (zdc.energySectorZPA ())[1 ] + (zdc.energySectorZPA ())[2 ] + (zdc.energySectorZPA ())[3 ]);
462444 if (applyZdcCorrection) {
463445 const float a = zpaCoeff;
464446 const float b = zpcCoeff;
@@ -475,8 +457,6 @@ struct FlowZdcTask {
475457 float tZEM1{zdc.timeZEM1 ()};
476458 float tZEM2{zdc.timeZEM2 ()};
477459 float sumZNs{znA + znC};
478- float sumSectZN = (sumZNC + sumZNA);
479- float sumSectZP = (sumZPC + sumZPA);
480460
481461 if (sumZEMs > zemCut) {
482462 if (isTDCcut) {
@@ -516,7 +496,7 @@ struct FlowZdcTask {
516496 histos.fill (HIST (" ZNVsZEM" ), sumZEMs, sumZNs);
517497 if (((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)) && ((tZNC >= minTdcZn) && (tZNC <= maxTdcZn))) {
518498 histos.fill (HIST (" ZNAVsZNC" ), znC, znA);
519- histos.fill (HIST (" SectorZNAVsZNC " ), sumZNC, sumZNA );
499+ histos.fill (HIST (" CommonZNAVsZNC " ), commonSumZnc, commonSumZna );
520500 histos.fill (HIST (" ZN" ), znA + znC);
521501 histos.fill (HIST (" ZNVsFT0C" ), aT0C / 100 ., znA + znC);
522502 histos.fill (HIST (" ZNVsFT0M" ), (aT0A + aT0C) / 100 ., znA + znC);
@@ -525,7 +505,7 @@ struct FlowZdcTask {
525505 }
526506 if (((tZPA >= minTdcZp) && (tZPA <= maxTdcZp)) && ((tZPC >= minTdcZp) && (tZPC <= maxTdcZp))) {
527507 histos.fill (HIST (" ZPAVsZPC" ), zpC, zpA);
528- histos.fill (HIST (" SectorZPAVsZPC " ), sumZPC, sumZPA );
508+ histos.fill (HIST (" CommonZPAVsZPC " ), commonSumZpc, commonSumZpa );
529509 histos.fill (HIST (" ZPVsFT0A" ), aT0A / 100 ., zpA + zpC);
530510 histos.fill (HIST (" ZPVsFT0C" ), aT0C / 100 ., zpA + zpC);
531511 histos.fill (HIST (" ZPVsFT0M" ), (aT0A + aT0C) / 100 ., zpA + zpC);
@@ -534,11 +514,11 @@ struct FlowZdcTask {
534514 }
535515 if (((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)) && ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp))) {
536516 histos.fill (HIST (" ZNAVsZPA" ), zpA, znA);
537- histos.fill (HIST (" SectorZNAVsZPA " ), sumZPA, sumZNA );
517+ histos.fill (HIST (" CommonZNAVsZPA " ), commonSumZpa, commonSumZna );
538518 }
539519 if (((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) && ((tZPC >= minTdcZp) && (tZPC <= maxTdcZp))) {
540520 histos.fill (HIST (" ZNCVsZPC" ), zpC, znC);
541- histos.fill (HIST (" SectorZNCVsZPC " ), sumZPC, sumZNC );
521+ histos.fill (HIST (" CommonZNCVsZPC " ), commonSumZpc, commonSumZnc );
542522 }
543523 } else {
544524 histos.fill (HIST (" ZNA" ), znA);
@@ -550,10 +530,10 @@ struct FlowZdcTask {
550530 histos.fill (HIST (" ZNAVsZPA" ), zpA, znA);
551531 histos.fill (HIST (" ZNCVsZPC" ), zpC, znC);
552532 histos.fill (HIST (" ZPAVsZPC" ), zpC, zpA);
553- histos.fill (HIST (" SectorZPAVsZPC " ), sumZPC, sumZPA );
554- histos.fill (HIST (" SectorZNAVsZNC " ), sumZNC, sumZNA );
555- histos.fill (HIST (" SectorZNAVsZPA " ), sumZPA, sumZNA );
556- histos.fill (HIST (" SectorZNCVsZPC " ), sumZPC, sumZNC );
533+ histos.fill (HIST (" CommonZPAVsZPC " ), commonSumZpc, commonSumZpa );
534+ histos.fill (HIST (" CommonZNAVsZNC " ), commonSumZnc, commonSumZna );
535+ histos.fill (HIST (" CommonZNAVsZPA " ), commonSumZpa, commonSumZna );
536+ histos.fill (HIST (" CommonZNCVsZPC " ), commonSumZpc, commonSumZnc );
557537 histos.fill (HIST (" ZNACommon" ), commonSumZna);
558538 histos.fill (HIST (" ZNASector" ), sumZNA);
559539 histos.fill (HIST (" ZNCCommon" ), commonSumZnc);
@@ -666,18 +646,6 @@ struct FlowZdcTask {
666646 histos.fill (HIST (" GlbTracks" ), glbTracks);
667647
668648 if (sumZEMs > zemCut) {
669- histos.fill (HIST (" ZPAvsCent" ), cent, zpA);
670- histos.fill (HIST (" ZPCvsCent" ), cent, zpC);
671- if (std::isfinite (zpA) && !std::isnan (zpA) && cent >= minT0CcentCut && cent < maxT0CcentCut && glbTracks >= minNch && glbTracks < maxNch) {
672- histos.fill (HIST (" hZPASectorvsGlbTrack" ), glbTracks, sumZPA);
673- }
674- if (std::isfinite (zpC) && !std::isnan (zpC) && cent >= minT0CcentCut && cent < maxT0CcentCut && glbTracks >= minNch && glbTracks < maxNch) {
675- histos.fill (HIST (" hZPCSectorvsGlbTrack" ), glbTracks, sumZPC);
676- }
677- histos.fill (HIST (" hZNASectorvsGlbTrack" ), glbTracks, sumZNA);
678- histos.fill (HIST (" hZNCSectorvsGlbTrack" ), glbTracks, sumZNC);
679- histos.fill (HIST (" hZPSectorvsGlbTrack" ), glbTracks, sumSectZP);
680- histos.fill (HIST (" hZNSectorvsGlbTrack" ), glbTracks, sumSectZN);
681649 // ZDC Correlations
682650 histos.fill (HIST (" ZNAVsNch" ), glbTracks, znA);
683651 histos.fill (HIST (" ZNCVsNch" ), glbTracks, znC);
@@ -715,18 +683,18 @@ struct FlowZdcTask {
715683 const auto & foundBC = collision.foundBC_as <BCsRun3>();
716684 if (foundBC.has_zdc ()) {
717685 const auto & zdc = foundBC.zdc ();
718- auto znA = zdc.amplitudeZNA () / cfgCollisionEnergy ;
719- auto znC = zdc.amplitudeZNC () / cfgCollisionEnergy ;
720- auto zpA = zdc.amplitudeZPA () / cfgCollisionEnergy ;
721- auto zpC = zdc.amplitudeZPC () / cfgCollisionEnergy ;
722- float sectorSumZNC = ((zdc.energySectorZNC ())[0 ] + (zdc.energySectorZNC ())[1 ] + (zdc.energySectorZNC ())[2 ] + (zdc.energySectorZNC ())[3 ]) / cfgCollisionEnergy ;
723- float sectorSumZNA = ((zdc.energySectorZNA ())[0 ] + (zdc.energySectorZNA ())[1 ] + (zdc.energySectorZNA ())[2 ] + (zdc.energySectorZNA ())[3 ]) / cfgCollisionEnergy ;
724- float sectorSumZPC = ((zdc.energySectorZPC ())[0 ] + (zdc.energySectorZPC ())[1 ] + (zdc.energySectorZPC ())[2 ] + (zdc.energySectorZPC ())[3 ]) / cfgCollisionEnergy ;
725- float sectorSumZPA = ((zdc.energySectorZPA ())[0 ] + (zdc.energySectorZPA ())[1 ] + (zdc.energySectorZPA ())[2 ] + (zdc.energySectorZPA ())[3 ]) / cfgCollisionEnergy ;
726- float commonSumZnc = zdc.energyCommonZNC () / cfgCollisionEnergy ;
727- float commonSumZna = zdc.energyCommonZNA () / cfgCollisionEnergy ;
728- float commonSumZpc = zdc.energyCommonZPC () / cfgCollisionEnergy ;
729- float commonSumZpa = zdc.energyCommonZPA () / cfgCollisionEnergy ;
686+ auto znA = zdc.amplitudeZNA ();
687+ auto znC = zdc.amplitudeZNC ();
688+ auto zpA = zdc.amplitudeZPA ();
689+ auto zpC = zdc.amplitudeZPC ();
690+ float sectorSumZNC = ((zdc.energySectorZNC ())[0 ] + (zdc.energySectorZNC ())[1 ] + (zdc.energySectorZNC ())[2 ] + (zdc.energySectorZNC ())[3 ]);
691+ float sectorSumZNA = ((zdc.energySectorZNA ())[0 ] + (zdc.energySectorZNA ())[1 ] + (zdc.energySectorZNA ())[2 ] + (zdc.energySectorZNA ())[3 ]);
692+ float sectorSumZPC = ((zdc.energySectorZPC ())[0 ] + (zdc.energySectorZPC ())[1 ] + (zdc.energySectorZPC ())[2 ] + (zdc.energySectorZPC ())[3 ]);
693+ float sectorSumZPA = ((zdc.energySectorZPA ())[0 ] + (zdc.energySectorZPA ())[1 ] + (zdc.energySectorZPA ())[2 ] + (zdc.energySectorZPA ())[3 ]);
694+ float commonSumZnc = zdc.energyCommonZNC ();
695+ float commonSumZna = zdc.energyCommonZNA ();
696+ float commonSumZpc = zdc.energyCommonZPC ();
697+ float commonSumZpa = zdc.energyCommonZPA ();
730698 float aZEM1 = zdc.amplitudeZEM1 ();
731699 float aZEM2 = zdc.amplitudeZEM2 ();
732700 float sumZEMs = aZEM1 + aZEM2;
0 commit comments