@@ -870,74 +870,74 @@ struct FlowZdcTask {
870870 BCsRun3 const & /* bcs*/ ,
871871 aod::Zdcs const & /* zdcs*/ )
872872 {
873- for (auto & collision : cols) {
874- const auto & foundBC = collision.foundBC_as <BCsRun3>();
875- if (foundBC.has_zdc ()) {
876- const auto & zdc = foundBC.zdc ();
877- auto znA = zdc.amplitudeZNA () / cfgCollisionEnergy;
878- auto znC = zdc.amplitudeZNC () / cfgCollisionEnergy;
879- auto zpA = zdc.amplitudeZPA () / cfgCollisionEnergy;
880- auto zpC = zdc.amplitudeZPC () / cfgCollisionEnergy;
881- float sumZNC = ((zdc.energySectorZNC ())[0 ] + (zdc.energySectorZNC ())[1 ] + (zdc.energySectorZNC ())[2 ] + (zdc.energySectorZNC ())[3 ]) / cfgCollisionEnergy;
882- float sumZNA = ((zdc.energySectorZNA ())[0 ] + (zdc.energySectorZNA ())[1 ] + (zdc.energySectorZNA ())[2 ] + (zdc.energySectorZNA ())[3 ]) / cfgCollisionEnergy;
883- float sumZPC = ((zdc.energySectorZPC ())[0 ] + (zdc.energySectorZPC ())[1 ] + (zdc.energySectorZPC ())[2 ] + (zdc.energySectorZPC ())[3 ]) / cfgCollisionEnergy;
884- float sumZPA = ((zdc.energySectorZPA ())[0 ] + (zdc.energySectorZPA ())[1 ] + (zdc.energySectorZPA ())[2 ] + (zdc.energySectorZPA ())[3 ]) / cfgCollisionEnergy;
885- float commonSumZnc = zdc.energyCommonZNC () / cfgCollisionEnergy;
886- float commonSumZna = zdc.energyCommonZNA () / cfgCollisionEnergy;
887- float commonSumZpc = zdc.energyCommonZPC () / cfgCollisionEnergy;
888- float commonSumZpa = zdc.energyCommonZPA () / cfgCollisionEnergy;
889- float aZEM1 = zdc.amplitudeZEM1 ();
890- float aZEM2 = zdc.amplitudeZEM2 ();
891- float sumZEMs = aZEM1 + aZEM2;
892- auto tZNA = zdc.timeZNA ();
893- auto tZNC = zdc.timeZNC ();
894- auto tZPA = zdc.timeZPA ();
895- auto tZPC = zdc.timeZPC ();
896- if (isTDCcut) {
897- if ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn))
898- histos.fill (HIST (" ampZna" ), znA);
899- if ((tZNC >= minTdcZn) && (tZNC <= minTdcZn))
900- histos.fill (HIST (" ampZnc" ), znC);
901- if ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp))
902- histos.fill (HIST (" ampZpa" ), zpA);
903- if ((tZPC >= minTdcZp) && (tZPC <= maxTdcZp))
904- histos.fill (HIST (" ampZpc" ), zpC);
905- if (((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) && ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)))
906- histos.fill (HIST (" ZnVsZem" ), sumZEMs, znC + znA);
907- if (((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) && ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)))
908- histos.fill (HIST (" ZnaVsZnc" ), znA, znC);
909- if (((tZPC >= minTdcZp) && (tZPC <= maxTdcZp)) && ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp)))
910- histos.fill (HIST (" ZpaVsZpc" ), zpA, zpC);
911- if ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn))
912- histos.fill (HIST (" ZnaVsZpa" ), znA, zpA);
913- if ((tZNC >= minTdcZn) && (tZNC <= maxTdcZn))
914- histos.fill (HIST (" ZncVsZpc" ), znC, zpC);
915- } else {
916- histos.fill (HIST (" ampZna" ), znA);
917- histos.fill (HIST (" ampZnc" ), znC);
918- histos.fill (HIST (" ampZpa" ), zpA);
919- histos.fill (HIST (" ampZpc" ), zpC);
920- histos.fill (HIST (" ZnVsZem" ), sumZEMs, znC + znA);
921- histos.fill (HIST (" ZnaVsZnc" ), znA, znC);
922- histos.fill (HIST (" ZpaVsZpc" ), zpA, zpC);
923- histos.fill (HIST (" ZnaVsZpa" ), znA, zpA);
924- histos.fill (HIST (" ZncVsZpc" ), znC, zpC);
925- }
926- histos.fill (HIST (" ampZEM1" ), aZEM1);
927- histos.fill (HIST (" ampZEM2" ), aZEM2);
928- histos.fill (HIST (" ZnccVsZncSum" ), sumZNC, commonSumZnc);
929- histos.fill (HIST (" ZnacVsZnaSum" ), sumZNA, commonSumZna);
930- histos.fill (HIST (" ZpccVsZpcSum" ), sumZPC, commonSumZpc);
931- histos.fill (HIST (" ZpacVsZpaSum" ), sumZPA, commonSumZpa);
932- histos.fill (HIST (" ZncVsTdc" ), zdc.timeZNC (), znC);
933- histos.fill (HIST (" ZnaVsTdc" ), zdc.timeZNA (), znA);
934- histos.fill (HIST (" ZpcVsTdc" ), zdc.timeZPC (), zpC);
935- histos.fill (HIST (" ZpaVsTdc" ), zdc.timeZPA (), zpA);
936- histos.fill (HIST (" Zem1VsTdc" ), zdc.timeZEM1 (), aZEM1);
937- histos.fill (HIST (" Zem2VsTdc" ), zdc.timeZEM2 (), aZEM2);
938- }
939- }
940- }
873+ for (const auto & collision : cols) {
874+ const auto & foundBC = collision.foundBC_as <BCsRun3>();
875+ if (foundBC.has_zdc ()) {
876+ const auto & zdc = foundBC.zdc ();
877+ auto znA = zdc.amplitudeZNA () / cfgCollisionEnergy;
878+ auto znC = zdc.amplitudeZNC () / cfgCollisionEnergy;
879+ auto zpA = zdc.amplitudeZPA () / cfgCollisionEnergy;
880+ auto zpC = zdc.amplitudeZPC () / cfgCollisionEnergy;
881+ float sumZNC = ((zdc.energySectorZNC ())[0 ] + (zdc.energySectorZNC ())[1 ] + (zdc.energySectorZNC ())[2 ] + (zdc.energySectorZNC ())[3 ]) / cfgCollisionEnergy;
882+ float sumZNA = ((zdc.energySectorZNA ())[0 ] + (zdc.energySectorZNA ())[1 ] + (zdc.energySectorZNA ())[2 ] + (zdc.energySectorZNA ())[3 ]) / cfgCollisionEnergy;
883+ float sumZPC = ((zdc.energySectorZPC ())[0 ] + (zdc.energySectorZPC ())[1 ] + (zdc.energySectorZPC ())[2 ] + (zdc.energySectorZPC ())[3 ]) / cfgCollisionEnergy;
884+ float sumZPA = ((zdc.energySectorZPA ())[0 ] + (zdc.energySectorZPA ())[1 ] + (zdc.energySectorZPA ())[2 ] + (zdc.energySectorZPA ())[3 ]) / cfgCollisionEnergy;
885+ float commonSumZnc = zdc.energyCommonZNC () / cfgCollisionEnergy;
886+ float commonSumZna = zdc.energyCommonZNA () / cfgCollisionEnergy;
887+ float commonSumZpc = zdc.energyCommonZPC () / cfgCollisionEnergy;
888+ float commonSumZpa = zdc.energyCommonZPA () / cfgCollisionEnergy;
889+ float aZEM1 = zdc.amplitudeZEM1 ();
890+ float aZEM2 = zdc.amplitudeZEM2 ();
891+ float sumZEMs = aZEM1 + aZEM2;
892+ auto tZNA = zdc.timeZNA ();
893+ auto tZNC = zdc.timeZNC ();
894+ auto tZPA = zdc.timeZPA ();
895+ auto tZPC = zdc.timeZPC ();
896+ if (isTDCcut) {
897+ if ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn))
898+ histos.fill (HIST (" ampZna" ), znA);
899+ if ((tZNC >= minTdcZn) && (tZNC <= minTdcZn))
900+ histos.fill (HIST (" ampZnc" ), znC);
901+ if ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp))
902+ histos.fill (HIST (" ampZpa" ), zpA);
903+ if ((tZPC >= minTdcZp) && (tZPC <= maxTdcZp))
904+ histos.fill (HIST (" ampZpc" ), zpC);
905+ if (((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) && ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)))
906+ histos.fill (HIST (" ZnVsZem" ), sumZEMs, znC + znA);
907+ if (((tZNC >= minTdcZn) && (tZNC <= maxTdcZn)) && ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn)))
908+ histos.fill (HIST (" ZnaVsZnc" ), znA, znC);
909+ if (((tZPC >= minTdcZp) && (tZPC <= maxTdcZp)) && ((tZPA >= minTdcZp) && (tZPA <= maxTdcZp)))
910+ histos.fill (HIST (" ZpaVsZpc" ), zpA, zpC);
911+ if ((tZNA >= minTdcZn) && (tZNA <= maxTdcZn))
912+ histos.fill (HIST (" ZnaVsZpa" ), znA, zpA);
913+ if ((tZNC >= minTdcZn) && (tZNC <= maxTdcZn))
914+ histos.fill (HIST (" ZncVsZpc" ), znC, zpC);
915+ } else {
916+ histos.fill (HIST (" ampZna" ), znA);
917+ histos.fill (HIST (" ampZnc" ), znC);
918+ histos.fill (HIST (" ampZpa" ), zpA);
919+ histos.fill (HIST (" ampZpc" ), zpC);
920+ histos.fill (HIST (" ZnVsZem" ), sumZEMs, znC + znA);
921+ histos.fill (HIST (" ZnaVsZnc" ), znA, znC);
922+ histos.fill (HIST (" ZpaVsZpc" ), zpA, zpC);
923+ histos.fill (HIST (" ZnaVsZpa" ), znA, zpA);
924+ histos.fill (HIST (" ZncVsZpc" ), znC, zpC);
925+ }
926+ histos.fill (HIST (" ampZEM1" ), aZEM1);
927+ histos.fill (HIST (" ampZEM2" ), aZEM2);
928+ histos.fill (HIST (" ZnccVsZncSum" ), sumZNC, commonSumZnc);
929+ histos.fill (HIST (" ZnacVsZnaSum" ), sumZNA, commonSumZna);
930+ histos.fill (HIST (" ZpccVsZpcSum" ), sumZPC, commonSumZpc);
931+ histos.fill (HIST (" ZpacVsZpaSum" ), sumZPA, commonSumZpa);
932+ histos.fill (HIST (" ZncVsTdc" ), zdc.timeZNC (), znC);
933+ histos.fill (HIST (" ZnaVsTdc" ), zdc.timeZNA (), znA);
934+ histos.fill (HIST (" ZpcVsTdc" ), zdc.timeZPC (), zpC);
935+ histos.fill (HIST (" ZpaVsTdc" ), zdc.timeZPA (), zpA);
936+ histos.fill (HIST (" Zem1VsTdc" ), zdc.timeZEM1 (), aZEM1);
937+ histos.fill (HIST (" Zem2VsTdc" ), zdc.timeZEM2 (), aZEM2);
938+ }
939+ }
940+ }
941941
942942 PROCESS_SWITCH (FlowZdcTask, processQA, " Process QA" , true );
943943 PROCESS_SWITCH (FlowZdcTask, processZdcCollAssoc, " Processing ZDC w. collision association" , false );
0 commit comments