diff --git a/PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx b/PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx index 99c753698d0..8e622f577ca 100644 --- a/PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx +++ b/PWGCF/TwoParticleCorrelations/Tasks/neutronProtonCorrZdc.cxx @@ -43,9 +43,9 @@ struct NeutronProtonCorrZdc { Configurable cfgAlphaZmin{"cfgAlphaZmin", -1, "Minimum value for ZDC asymmetry"}; Configurable cfgAlphaZmax{"cfgAlphaZmax", 1, "Maximum value for ZDC asymmetry"}; Configurable cfgProcessRun2{"cfgProcessRun2", false, "Analyse Run 2 converted data"}; - //Configurable cfgCentralityEstimator{"cfgCentralityEstimator", 0, "Choice of centrality estimator"};//0 for FTOC, 1 for FTOA, 2 for FTOM, 3 for FVOA. //To be included at a later stage + // Configurable cfgCentralityEstimator{"cfgCentralityEstimator", 0, "Choice of centrality estimator"};//0 for FTOC, 1 for FTOA, 2 for FTOM, 3 for FVOA. //To be included at a later stage - ConfigurableAxis cfgAxisCent{"cfgAxisCent", {VARIABLE_WIDTH, 0.0f, 0.1f, 0.2f, 0.3f, 0.4f, 0.5f, 0.6f, 0.7f, 0.8f, 0.9f, 1.0f, 1.1f, 1.2f, 1.3f, 1.4f, 1.5f, 1.6f, 1.7f, 1.8f, 1.9f, 2.0f, 2.2f, 2.4f, 2.6f, 2.8f, 3.0f, 3.2f, 3.4f, 3.6f, 3.8f, 4.0f, 4.2f, 4.4f, 4.6f, 4.8f, 5.0f, 5.5f, 6.0f, 6.5f, 7.0f, 7.5f, 8.0f, 8.5f, 9.0f, 9.5f, 10.0f, 11.0f, 12.0f, 13.0f, 14.0f, 15.0f, 16.0f, 17.0f, 18.0f, 19.0f, 20.0f, 21.0f, 22.0f, 23.0f, 24.0f , 25.0f, 26.0f, 27.0f, 28.0f, 29.0f, 30.0f, 31.0f, 32.0f, 33.0f, 34.0f, 35.0f, 36.0f, 37.0f, 38.0f, 39.0f, 40.0f, 41.0f, 42.0f, 43.0f, 44.0f, 45.0f, 46.0f, 47.0f, 48.0f, 49.0f, 50.0f, 51.0f, 52.0f, 53.0f, 54.0f, 55.0f, 56.0f, 57.0f, 58.0f, 59.0f, 60.0f, 61.0f, 62.0f, 63.0f, 64.0f, 65.0f, 66.0f, 67.0f, 68.0f, 69.0f, 70.0f, 71.0f, 72.0f, 73.0f, 74.0f, 75.0f, 76.0f, 77.0f, 78.0f, 79.0f, 80.0f, 81.0f, 82.0f, 83.0f, 84.0f, 85.0f, 86.0f, 87.0f, 88.0f, 89.0f, 90.0f, 91.0f, 92.0f, 93.0f, 94.0f, 95.0f, 96.0f, 97.0f, 98.0f, 99.0f, 100.0f}, "Centrality [%]"}; + ConfigurableAxis cfgAxisCent{"cfgAxisCent", {VARIABLE_WIDTH, 0.0f, 0.1f, 0.2f, 0.3f, 0.4f, 0.5f, 0.6f, 0.7f, 0.8f, 0.9f, 1.0f, 1.1f, 1.2f, 1.3f, 1.4f, 1.5f, 1.6f, 1.7f, 1.8f, 1.9f, 2.0f, 2.2f, 2.4f, 2.6f, 2.8f, 3.0f, 3.2f, 3.4f, 3.6f, 3.8f, 4.0f, 4.2f, 4.4f, 4.6f, 4.8f, 5.0f, 5.5f, 6.0f, 6.5f, 7.0f, 7.5f, 8.0f, 8.5f, 9.0f, 9.5f, 10.0f, 11.0f, 12.0f, 13.0f, 14.0f, 15.0f, 16.0f, 17.0f, 18.0f, 19.0f, 20.0f, 21.0f, 22.0f, 23.0f, 24.0f, 25.0f, 26.0f, 27.0f, 28.0f, 29.0f, 30.0f, 31.0f, 32.0f, 33.0f, 34.0f, 35.0f, 36.0f, 37.0f, 38.0f, 39.0f, 40.0f, 41.0f, 42.0f, 43.0f, 44.0f, 45.0f, 46.0f, 47.0f, 48.0f, 49.0f, 50.0f, 51.0f, 52.0f, 53.0f, 54.0f, 55.0f, 56.0f, 57.0f, 58.0f, 59.0f, 60.0f, 61.0f, 62.0f, 63.0f, 64.0f, 65.0f, 66.0f, 67.0f, 68.0f, 69.0f, 70.0f, 71.0f, 72.0f, 73.0f, 74.0f, 75.0f, 76.0f, 77.0f, 78.0f, 79.0f, 80.0f, 81.0f, 82.0f, 83.0f, 84.0f, 85.0f, 86.0f, 87.0f, 88.0f, 89.0f, 90.0f, 91.0f, 92.0f, 93.0f, 94.0f, 95.0f, 96.0f, 97.0f, 98.0f, 99.0f, 100.0f}, "Centrality [%]"}; Filter collisionVtxZ = nabs(aod::collision::posZ) < 10.f; @@ -57,77 +57,75 @@ struct NeutronProtonCorrZdc { { // define axes you want to use const AxisSpec axisCounter{4, 0, +2, ""}; - const AxisSpec axisZNSectorSignal{cfgNBinsZN,cfgZNmin,cfgZNmax/3.}; - const AxisSpec axisZPSectorSignal{cfgNBinsZP,cfgZPmin,cfgZPmax/3.}; - const AxisSpec axisZNASignal{cfgNBinsZN,cfgZNmin,cfgZNmax, "ZNA (a.u.)"}; - const AxisSpec axisZNCSignal{cfgNBinsZN,cfgZNmin,cfgZNmax, "ZNC (a.u.)"}; - const AxisSpec axisZPASignal{cfgNBinsZP,cfgZPmin,cfgZPmax, "ZPA (a.u.)"}; - const AxisSpec axisZPCSignal{cfgNBinsZP,cfgZPmin,cfgZPmax, "ZPC (a.u.)"}; - const AxisSpec axisZNSignal{2*cfgNBinsZN, cfgZNmin, 1.5*cfgZNmax,"ZN (a.u.)"}; - const AxisSpec axisZPSignal{2*cfgNBinsZP, cfgZPmin, 1.5*cfgZPmax,"ZP (a.u.)"}; - const AxisSpec axisAlphaZ{cfgNBinsAlpha,cfgAlphaZmin,cfgAlphaZmax, "#alpha_{spec}"}; - const AxisSpec axisZDiffSignal{cfgNBinsZN,cfgDiffZNmin,cfgDiffZNmax,"#Delta E"}; - + const AxisSpec axisZNSectorSignal{cfgNBinsZN, cfgZNmin, cfgZNmax / 3.}; + const AxisSpec axisZPSectorSignal{cfgNBinsZP, cfgZPmin, cfgZPmax / 3.}; + const AxisSpec axisZNASignal{cfgNBinsZN, cfgZNmin, cfgZNmax, "ZNA (a.u.)"}; + const AxisSpec axisZNCSignal{cfgNBinsZN, cfgZNmin, cfgZNmax, "ZNC (a.u.)"}; + const AxisSpec axisZPASignal{cfgNBinsZP, cfgZPmin, cfgZPmax, "ZPA (a.u.)"}; + const AxisSpec axisZPCSignal{cfgNBinsZP, cfgZPmin, cfgZPmax, "ZPC (a.u.)"}; + const AxisSpec axisZNSignal{2 * cfgNBinsZN, cfgZNmin, 1.5 * cfgZNmax, "ZN (a.u.)"}; + const AxisSpec axisZPSignal{2 * cfgNBinsZP, cfgZPmin, 1.5 * cfgZPmax, "ZP (a.u.)"}; + const AxisSpec axisAlphaZ{cfgNBinsAlpha, cfgAlphaZmin, cfgAlphaZmax, "#alpha_{spec}"}; + const AxisSpec axisZDiffSignal{cfgNBinsZN, cfgDiffZNmin, cfgDiffZNmax, "#Delta E"}; // create histograms histos.add("eventCounter", "eventCounter", kTH1F, {axisCounter}); - histos.add("ZNASector0Signal","ZNASector0Signal",kTH1F,{axisZNSectorSignal}); - histos.add("ZNASector1Signal","ZNASector1Signal",kTH1F,{axisZNSectorSignal}); - histos.add("ZNASector2Signal","ZNASector2Signal",kTH1F,{axisZNSectorSignal}); - histos.add("ZNASector3Signal","ZNASector3Signal",kTH1F,{axisZNSectorSignal}); - - histos.add("ZNCSector0Signal","ZNCSector0Signal",kTH1F,{axisZNSectorSignal}); - histos.add("ZNCSector1Signal","ZNCSector1Signal",kTH1F,{axisZNSectorSignal}); - histos.add("ZNCSector2Signal","ZNCSector2Signal",kTH1F,{axisZNSectorSignal}); - histos.add("ZNCSector3Signal","ZNCSector3Signal",kTH1F,{axisZNSectorSignal}); - - histos.add("ZPASector0Signal","ZPASector0Signal",kTH1F,{axisZPSectorSignal}); - histos.add("ZPASector1Signal","ZPASector1Signal",kTH1F,{axisZPSectorSignal}); - histos.add("ZPASector2Signal","ZPASector2Signal",kTH1F,{axisZPSectorSignal}); - histos.add("ZPASector3Signal","ZPASector3Signal",kTH1F,{axisZPSectorSignal}); - - histos.add("ZPCSector0Signal","ZPCSector0Signal",kTH1F,{axisZPSectorSignal}); - histos.add("ZPCSector1Signal","ZPCSector1Signal",kTH1F,{axisZPSectorSignal}); - histos.add("ZPCSector2Signal","ZPCSector2Signal",kTH1F,{axisZPSectorSignal}); - histos.add("ZPCSector3Signal","ZPCSector3Signal",kTH1F,{axisZPSectorSignal}); - - histos.add("ZNASignal","ZNASignal",kTH1F,{axisZNASignal}); - histos.add("ZNCSignal","ZNCSignal",kTH1F,{axisZNCSignal}); - histos.add("ZNSignal", "ZNSignal", kTH1F,{axisZNSignal}); - histos.add("ZPASignal","ZPASignal",kTH1F,{axisZPASignal}); - histos.add("ZPCSignal","ZPCSignal",kTH1F,{axisZPCSignal}); - histos.add("ZPSignal","ZPSignal", kTH1F, {axisZPSignal}); + histos.add("ZNASector0Signal", "ZNASector0Signal", kTH1F, {axisZNSectorSignal}); + histos.add("ZNASector1Signal", "ZNASector1Signal", kTH1F, {axisZNSectorSignal}); + histos.add("ZNASector2Signal", "ZNASector2Signal", kTH1F, {axisZNSectorSignal}); + histos.add("ZNASector3Signal", "ZNASector3Signal", kTH1F, {axisZNSectorSignal}); + + histos.add("ZNCSector0Signal", "ZNCSector0Signal", kTH1F, {axisZNSectorSignal}); + histos.add("ZNCSector1Signal", "ZNCSector1Signal", kTH1F, {axisZNSectorSignal}); + histos.add("ZNCSector2Signal", "ZNCSector2Signal", kTH1F, {axisZNSectorSignal}); + histos.add("ZNCSector3Signal", "ZNCSector3Signal", kTH1F, {axisZNSectorSignal}); + + histos.add("ZPASector0Signal", "ZPASector0Signal", kTH1F, {axisZPSectorSignal}); + histos.add("ZPASector1Signal", "ZPASector1Signal", kTH1F, {axisZPSectorSignal}); + histos.add("ZPASector2Signal", "ZPASector2Signal", kTH1F, {axisZPSectorSignal}); + histos.add("ZPASector3Signal", "ZPASector3Signal", kTH1F, {axisZPSectorSignal}); + + histos.add("ZPCSector0Signal", "ZPCSector0Signal", kTH1F, {axisZPSectorSignal}); + histos.add("ZPCSector1Signal", "ZPCSector1Signal", kTH1F, {axisZPSectorSignal}); + histos.add("ZPCSector2Signal", "ZPCSector2Signal", kTH1F, {axisZPSectorSignal}); + histos.add("ZPCSector3Signal", "ZPCSector3Signal", kTH1F, {axisZPSectorSignal}); + + histos.add("ZNASignal", "ZNASignal", kTH1F, {axisZNASignal}); + histos.add("ZNCSignal", "ZNCSignal", kTH1F, {axisZNCSignal}); + histos.add("ZNSignal", "ZNSignal", kTH1F, {axisZNSignal}); + histos.add("ZPASignal", "ZPASignal", kTH1F, {axisZPASignal}); + histos.add("ZPCSignal", "ZPCSignal", kTH1F, {axisZPCSignal}); + histos.add("ZPSignal", "ZPSignal", kTH1F, {axisZPSignal}); histos.add("alphaZN", "alphaZN", kTH1F, {axisAlphaZ}); histos.add("alphaZP", "alphaZP", kTH1F, {axisAlphaZ}); - histos.add("diffZNASignal","diffZNASignal",kTH1F,{axisZDiffSignal}); - histos.add("diffZNCSignal","diffZNCSignal",kTH1F,{axisZDiffSignal}); - histos.add("diffZPASignal","diffZPASignal",kTH1F,{axisZDiffSignal}); - histos.add("diffZPCSignal","diffZPCSignal",kTH1F,{axisZDiffSignal}); - histos.add("diffZNSignal","diffZNSignal",kTH1F,{axisZDiffSignal}); - histos.add("diffZPSignal","diffZPSignal",kTH1F,{axisZDiffSignal}); - - histos.add("CentralityPercentile","CentralityPercentile",kTH1F,{cfgAxisCent}); - - histos.add("CentvsZNASignal","FT0CvsZNASignal",kTH2F,{cfgAxisCent,axisZNASignal}); - histos.add("CentvsZNCSignal","FT0CvsZNCSignal",kTH2F,{cfgAxisCent,axisZNCSignal}); - histos.add("CentvsZPASignal","FT0CvsZPASignal",kTH2F,{cfgAxisCent,axisZPASignal}); - histos.add("CentvsZPCSignal","FT0CvsZPCSignal",kTH2F,{cfgAxisCent,axisZPCSignal}); - histos.add("CentvsZNSignal","FT0CvsZNSignal",kTH2F,{cfgAxisCent,axisZNSignal}); - histos.add("CentvsZPSignal","FT0CvsZPSignal",kTH2F,{cfgAxisCent,axisZPSignal}); - + histos.add("diffZNASignal", "diffZNASignal", kTH1F, {axisZDiffSignal}); + histos.add("diffZNCSignal", "diffZNCSignal", kTH1F, {axisZDiffSignal}); + histos.add("diffZPASignal", "diffZPASignal", kTH1F, {axisZDiffSignal}); + histos.add("diffZPCSignal", "diffZPCSignal", kTH1F, {axisZDiffSignal}); + histos.add("diffZNSignal", "diffZNSignal", kTH1F, {axisZDiffSignal}); + histos.add("diffZPSignal", "diffZPSignal", kTH1F, {axisZDiffSignal}); + + histos.add("CentralityPercentile", "CentralityPercentile", kTH1F, {cfgAxisCent}); + + histos.add("CentvsZNASignal", "FT0CvsZNASignal", kTH2F, {cfgAxisCent, axisZNASignal}); + histos.add("CentvsZNCSignal", "FT0CvsZNCSignal", kTH2F, {cfgAxisCent, axisZNCSignal}); + histos.add("CentvsZPASignal", "FT0CvsZPASignal", kTH2F, {cfgAxisCent, axisZPASignal}); + histos.add("CentvsZPCSignal", "FT0CvsZPCSignal", kTH2F, {cfgAxisCent, axisZPCSignal}); + histos.add("CentvsZNSignal", "FT0CvsZNSignal", kTH2F, {cfgAxisCent, axisZNSignal}); + histos.add("CentvsZPSignal", "FT0CvsZPSignal", kTH2F, {cfgAxisCent, axisZPSignal}); } - void processRun3(soa::Filtered>::iterator const& collision , BCsRun3 const&, aod::Zdcs const&) + void processRun3(soa::Filtered>::iterator const& collision, BCsRun3 const&, aod::Zdcs const&) { - if(!collision.sel8()){ + if (!collision.sel8()) { histos.fill(HIST("eventCounter"), 0.25); return; } const auto& foundBC = collision.foundBC_as(); - if (foundBC.has_zdc()){ + if (foundBC.has_zdc()) { const auto& zdcread = foundBC.zdc(); const auto cent = collision.centFT0C(); @@ -136,22 +134,22 @@ struct NeutronProtonCorrZdc { histos.fill(HIST("ZNCSignal"), zdcread.energyCommonZNC()); histos.fill(HIST("ZPASignal"), zdcread.energyCommonZPA()); histos.fill(HIST("ZPCSignal"), zdcread.energyCommonZPC()); - histos.fill(HIST("ZNASector0Signal"),(zdcread.energySectorZNA())[0]); - histos.fill(HIST("ZNASector1Signal"),(zdcread.energySectorZNA())[1]); - histos.fill(HIST("ZNASector2Signal"),(zdcread.energySectorZNA())[2]); - histos.fill(HIST("ZNASector3Signal"),(zdcread.energySectorZNA())[3]); - histos.fill(HIST("ZNCSector0Signal"),(zdcread.energySectorZNC())[0]); - histos.fill(HIST("ZNCSector1Signal"),(zdcread.energySectorZNC())[1]); - histos.fill(HIST("ZNCSector2Signal"),(zdcread.energySectorZNC())[2]); - histos.fill(HIST("ZNCSector3Signal"),(zdcread.energySectorZNC())[3]); - histos.fill(HIST("ZPASector0Signal"),(zdcread.energySectorZPA())[0]); - histos.fill(HIST("ZPASector1Signal"),(zdcread.energySectorZPA())[1]); - histos.fill(HIST("ZPASector2Signal"),(zdcread.energySectorZPA())[2]); - histos.fill(HIST("ZPASector3Signal"),(zdcread.energySectorZPA())[3]); - histos.fill(HIST("ZPCSector0Signal"),(zdcread.energySectorZPC())[0]); - histos.fill(HIST("ZPCSector1Signal"),(zdcread.energySectorZPC())[1]); - histos.fill(HIST("ZPCSector2Signal"),(zdcread.energySectorZPC())[2]); - histos.fill(HIST("ZPCSector3Signal"),(zdcread.energySectorZPC())[3]); + histos.fill(HIST("ZNASector0Signal"), (zdcread.energySectorZNA())[0]); + histos.fill(HIST("ZNASector1Signal"), (zdcread.energySectorZNA())[1]); + histos.fill(HIST("ZNASector2Signal"), (zdcread.energySectorZNA())[2]); + histos.fill(HIST("ZNASector3Signal"), (zdcread.energySectorZNA())[3]); + histos.fill(HIST("ZNCSector0Signal"), (zdcread.energySectorZNC())[0]); + histos.fill(HIST("ZNCSector1Signal"), (zdcread.energySectorZNC())[1]); + histos.fill(HIST("ZNCSector2Signal"), (zdcread.energySectorZNC())[2]); + histos.fill(HIST("ZNCSector3Signal"), (zdcread.energySectorZNC())[3]); + histos.fill(HIST("ZPASector0Signal"), (zdcread.energySectorZPA())[0]); + histos.fill(HIST("ZPASector1Signal"), (zdcread.energySectorZPA())[1]); + histos.fill(HIST("ZPASector2Signal"), (zdcread.energySectorZPA())[2]); + histos.fill(HIST("ZPASector3Signal"), (zdcread.energySectorZPA())[3]); + histos.fill(HIST("ZPCSector0Signal"), (zdcread.energySectorZPC())[0]); + histos.fill(HIST("ZPCSector1Signal"), (zdcread.energySectorZPC())[1]); + histos.fill(HIST("ZPCSector2Signal"), (zdcread.energySectorZPC())[2]); + histos.fill(HIST("ZPCSector3Signal"), (zdcread.energySectorZPC())[3]); float sumZNC = (zdcread.energySectorZNC())[0] + (zdcread.energySectorZNC())[1] + (zdcread.energySectorZNC())[2] + (zdcread.energySectorZNC())[3]; float sumZNA = (zdcread.energySectorZNA())[0] + (zdcread.energySectorZNA())[1] + (zdcread.energySectorZNA())[2] + (zdcread.energySectorZNA())[3]; @@ -162,15 +160,15 @@ struct NeutronProtonCorrZdc { histos.fill(HIST("alphaZN"), alphaZN); histos.fill(HIST("alphaZP"), alphaZP); - - histos.fill(HIST("diffZNASignal"),sumZNA - zdcread.energyCommonZNA()); - histos.fill(HIST("diffZNCSignal"),sumZNC - zdcread.energyCommonZNC()); - histos.fill(HIST("diffZPASignal"),sumZPA - zdcread.energyCommonZPA()); - histos.fill(HIST("diffZPCSignal"),sumZPC - zdcread.energyCommonZPC()); + + histos.fill(HIST("diffZNASignal"), sumZNA - zdcread.energyCommonZNA()); + histos.fill(HIST("diffZNCSignal"), sumZNC - zdcread.energyCommonZNC()); + histos.fill(HIST("diffZPASignal"), sumZPA - zdcread.energyCommonZPA()); + histos.fill(HIST("diffZPCSignal"), sumZPC - zdcread.energyCommonZPC()); histos.fill(HIST("ZNSignal"), sumZNA + sumZNC); histos.fill(HIST("ZPSignal"), sumZPA + sumZPC); - histos.fill(HIST("diffZNSignal"), (sumZNA + sumZNC) - (zdcread.energyCommonZNA()+zdcread.energyCommonZNC())); - histos.fill(HIST("diffZPSignal"), (sumZPA + sumZPC) - (zdcread.energyCommonZPA()+zdcread.energyCommonZPC())); + histos.fill(HIST("diffZNSignal"), (sumZNA + sumZNC) - (zdcread.energyCommonZNA() + zdcread.energyCommonZNC())); + histos.fill(HIST("diffZPSignal"), (sumZPA + sumZPC) - (zdcread.energyCommonZPA() + zdcread.energyCommonZPC())); histos.fill(HIST("CentralityPercentile"), cent); histos.fill(HIST("CentvsZNASignal"), cent, sumZNA); histos.fill(HIST("CentvsZNCSignal"), cent, sumZNC); @@ -182,14 +180,14 @@ struct NeutronProtonCorrZdc { } PROCESS_SWITCH(NeutronProtonCorrZdc, processRun3, "Process analysis for Run 3 data", true); - void processRun2(soa::Filtered>::iterator const& collision ,aod::BCsWithTimestamps const&, aod::Zdcs const&) + void processRun2(soa::Filtered>::iterator const& collision, aod::BCsWithTimestamps const&, aod::Zdcs const&) { - if(!collision.alias_bit(kINT7)){ + if (!collision.alias_bit(kINT7)) { histos.fill(HIST("eventCounter"), 0.25); return; } - - if (collision.has_zdc()){ + + if (collision.has_zdc()) { const auto& zdcread = collision.zdc(); const auto cent = collision.centRun2V0M(); histos.fill(HIST("eventCounter"), 1.25); @@ -197,22 +195,22 @@ struct NeutronProtonCorrZdc { histos.fill(HIST("ZNCSignal"), zdcread.energyCommonZNC()); histos.fill(HIST("ZPASignal"), zdcread.energyCommonZPA()); histos.fill(HIST("ZPCSignal"), zdcread.energyCommonZPC()); - histos.fill(HIST("ZNASector0Signal"),(zdcread.energySectorZNA())[0]); - histos.fill(HIST("ZNASector1Signal"),(zdcread.energySectorZNA())[1]); - histos.fill(HIST("ZNASector2Signal"),(zdcread.energySectorZNA())[2]); - histos.fill(HIST("ZNASector3Signal"),(zdcread.energySectorZNA())[3]); - histos.fill(HIST("ZNCSector0Signal"),(zdcread.energySectorZNC())[0]); - histos.fill(HIST("ZNCSector1Signal"),(zdcread.energySectorZNC())[1]); - histos.fill(HIST("ZNCSector2Signal"),(zdcread.energySectorZNC())[2]); - histos.fill(HIST("ZNCSector3Signal"),(zdcread.energySectorZNC())[3]); - histos.fill(HIST("ZPASector0Signal"),(zdcread.energySectorZPA())[0]); - histos.fill(HIST("ZPASector1Signal"),(zdcread.energySectorZPA())[1]); - histos.fill(HIST("ZPASector2Signal"),(zdcread.energySectorZPA())[2]); - histos.fill(HIST("ZPASector3Signal"),(zdcread.energySectorZPA())[3]); - histos.fill(HIST("ZPCSector0Signal"),(zdcread.energySectorZPC())[0]); - histos.fill(HIST("ZPCSector1Signal"),(zdcread.energySectorZPC())[1]); - histos.fill(HIST("ZPCSector2Signal"),(zdcread.energySectorZPC())[2]); - histos.fill(HIST("ZPCSector3Signal"),(zdcread.energySectorZPC())[3]); + histos.fill(HIST("ZNASector0Signal"), (zdcread.energySectorZNA())[0]); + histos.fill(HIST("ZNASector1Signal"), (zdcread.energySectorZNA())[1]); + histos.fill(HIST("ZNASector2Signal"), (zdcread.energySectorZNA())[2]); + histos.fill(HIST("ZNASector3Signal"), (zdcread.energySectorZNA())[3]); + histos.fill(HIST("ZNCSector0Signal"), (zdcread.energySectorZNC())[0]); + histos.fill(HIST("ZNCSector1Signal"), (zdcread.energySectorZNC())[1]); + histos.fill(HIST("ZNCSector2Signal"), (zdcread.energySectorZNC())[2]); + histos.fill(HIST("ZNCSector3Signal"), (zdcread.energySectorZNC())[3]); + histos.fill(HIST("ZPASector0Signal"), (zdcread.energySectorZPA())[0]); + histos.fill(HIST("ZPASector1Signal"), (zdcread.energySectorZPA())[1]); + histos.fill(HIST("ZPASector2Signal"), (zdcread.energySectorZPA())[2]); + histos.fill(HIST("ZPASector3Signal"), (zdcread.energySectorZPA())[3]); + histos.fill(HIST("ZPCSector0Signal"), (zdcread.energySectorZPC())[0]); + histos.fill(HIST("ZPCSector1Signal"), (zdcread.energySectorZPC())[1]); + histos.fill(HIST("ZPCSector2Signal"), (zdcread.energySectorZPC())[2]); + histos.fill(HIST("ZPCSector3Signal"), (zdcread.energySectorZPC())[3]); float sumZNC = (zdcread.energySectorZNC())[0] + (zdcread.energySectorZNC())[1] + (zdcread.energySectorZNC())[2] + (zdcread.energySectorZNC())[3]; float sumZNA = (zdcread.energySectorZNA())[0] + (zdcread.energySectorZNA())[1] + (zdcread.energySectorZNA())[2] + (zdcread.energySectorZNA())[3]; float sumZPC = (zdcread.energySectorZPC())[0] + (zdcread.energySectorZPC())[1] + (zdcread.energySectorZPC())[2] + (zdcread.energySectorZPC())[3]; @@ -221,15 +219,15 @@ struct NeutronProtonCorrZdc { float alphaZP = (sumZPA - sumZPC) / (sumZPA + sumZPC); histos.fill(HIST("alphaZN"), alphaZN); histos.fill(HIST("alphaZP"), alphaZP); - - histos.fill(HIST("diffZNASignal"),sumZNA - zdcread.energyCommonZNA()); - histos.fill(HIST("diffZNCSignal"),sumZNC - zdcread.energyCommonZNC()); - histos.fill(HIST("diffZPASignal"),sumZPA - zdcread.energyCommonZPA()); - histos.fill(HIST("diffZPCSignal"),sumZPC - zdcread.energyCommonZPC()); + + histos.fill(HIST("diffZNASignal"), sumZNA - zdcread.energyCommonZNA()); + histos.fill(HIST("diffZNCSignal"), sumZNC - zdcread.energyCommonZNC()); + histos.fill(HIST("diffZPASignal"), sumZPA - zdcread.energyCommonZPA()); + histos.fill(HIST("diffZPCSignal"), sumZPC - zdcread.energyCommonZPC()); histos.fill(HIST("ZNSignal"), sumZNA + sumZNC); histos.fill(HIST("ZPSignal"), sumZPA + sumZPC); - histos.fill(HIST("diffZNSignal"), (sumZNA + sumZNC) - (zdcread.energyCommonZNA()+zdcread.energyCommonZNC())); - histos.fill(HIST("diffZPSignal"), (sumZPA + sumZPC) - (zdcread.energyCommonZPA()+zdcread.energyCommonZPC())); + histos.fill(HIST("diffZNSignal"), (sumZNA + sumZNC) - (zdcread.energyCommonZNA() + zdcread.energyCommonZNC())); + histos.fill(HIST("diffZPSignal"), (sumZPA + sumZPC) - (zdcread.energyCommonZPA() + zdcread.energyCommonZPC())); histos.fill(HIST("CentralityPercentile"), cent); histos.fill(HIST("CentvsZNASignal"), cent, sumZNA); histos.fill(HIST("CentvsZNCSignal"), cent, sumZNC); @@ -246,4 +244,5 @@ WorkflowSpec defineDataProcessing(ConfigContext const& cfgc) { return WorkflowSpec{ adaptAnalysisTask(cfgc)}; -} \ No newline at end of file +} + \ No newline at end of file