diff --git a/PWGEM/PhotonMeson/TableProducer/createEMEventPhoton.cxx b/PWGEM/PhotonMeson/TableProducer/createEMEventPhoton.cxx index e3cb64523e2..b83fcf67fbd 100644 --- a/PWGEM/PhotonMeson/TableProducer/createEMEventPhoton.cxx +++ b/PWGEM/PhotonMeson/TableProducer/createEMEventPhoton.cxx @@ -191,10 +191,10 @@ struct CreateEMEvent { const size_t qvecSize = collision.qvecFT0CReVec().size(); if (qvecSize >= 2) { // harmonics 2,3 event_qvec(collision.qvecFT0MReVec()[0], collision.qvecFT0MImVec()[0], collision.qvecFT0AReVec()[0], collision.qvecFT0AImVec()[0], collision.qvecFT0CReVec()[0], collision.qvecFT0CImVec()[0], collision.qvecBPosReVec()[0], collision.qvecBPosImVec()[0], collision.qvecBNegReVec()[0], collision.qvecBNegImVec()[0], collision.qvecBTotReVec()[0], collision.qvecBTotImVec()[0], - collision.qvecFT0MReVec()[1], collision.qvecFT0MImVec()[1], collision.qvecFT0AReVec()[1], collision.qvecFT0AImVec()[1], collision.qvecFT0CReVec()[1], collision.qvecFT0CImVec()[1], collision.qvecBPosReVec()[1], collision.qvecBPosImVec()[1], collision.qvecBNegReVec()[1], collision.qvecBNegImVec()[1], collision.qvecBTotReVec()[1], collision.qvecBTotImVec()[1]); + collision.qvecFT0MReVec()[1], collision.qvecFT0MImVec()[1], collision.qvecFT0AReVec()[1], collision.qvecFT0AImVec()[1], collision.qvecFT0CReVec()[1], collision.qvecFT0CImVec()[1], collision.qvecBPosReVec()[1], collision.qvecBPosImVec()[1], collision.qvecBNegReVec()[1], collision.qvecBNegImVec()[1], collision.qvecBTotReVec()[1], collision.qvecBTotImVec()[1]); } else if (qvecSize >= 1) { // harmonics 2 event_qvec(collision.qvecFT0MReVec()[0], collision.qvecFT0MImVec()[0], collision.qvecFT0AReVec()[0], collision.qvecFT0AImVec()[0], collision.qvecFT0CReVec()[0], collision.qvecFT0CImVec()[0], collision.qvecBPosReVec()[0], collision.qvecBPosImVec()[0], collision.qvecBNegReVec()[0], collision.qvecBNegImVec()[0], collision.qvecBTotReVec()[0], collision.qvecBTotImVec()[0], - qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault); + qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault, qDefault); } } else { event_cent(105.f, 105.f, 105.f); diff --git a/PWGEM/PhotonMeson/Tasks/taskPi0FlowEMC.cxx b/PWGEM/PhotonMeson/Tasks/taskPi0FlowEMC.cxx index 6df505d8b5b..1e2eeee94b1 100644 --- a/PWGEM/PhotonMeson/Tasks/taskPi0FlowEMC.cxx +++ b/PWGEM/PhotonMeson/Tasks/taskPi0FlowEMC.cxx @@ -708,10 +708,10 @@ struct TaskPi0FlowEMC { float dTheta = photon1.Theta() - photon3.Theta(); float dPhi = photon1.Phi() - photon3.Phi(); if (mesonConfig.minTanThetadPhi > std::fabs(getAngleDegree(std::atan(dTheta / dPhi)))) { - registry.fill(HIST("hSparseCalibBack"), mother1.M(), mother1.E()/2., cent); + registry.fill(HIST("hSparseCalibBack"), mother1.M(), mother1.E() / 2., cent); } } else { - registry.fill(HIST("hSparseCalibBack"), mother1.M(), mother1.E()/2., cent); + registry.fill(HIST("hSparseCalibBack"), mother1.M(), mother1.E() / 2., cent); } } } @@ -726,10 +726,10 @@ struct TaskPi0FlowEMC { float dTheta = photon2.Theta() - photon3.Theta(); float dPhi = photon2.Phi() - photon3.Phi(); if (mesonConfig.minTanThetadPhi > std::fabs(getAngleDegree(std::atan(dTheta / dPhi)))) { - registry.fill(HIST("hSparseCalibBack"), mother2.M(), mother2.E()/2., cent); + registry.fill(HIST("hSparseCalibBack"), mother2.M(), mother2.E() / 2., cent); } } else { - registry.fill(HIST("hSparseCalibBack"), mother2.M(), mother2.E()/2., cent); + registry.fill(HIST("hSparseCalibBack"), mother2.M(), mother2.E() / 2., cent); } } } @@ -1096,8 +1096,9 @@ struct TaskPi0FlowEMC { PROCESS_SWITCH(TaskPi0FlowEMC, processResolution, "Process resolution", false); // EMCal calibration - void processEMCalCalib(CollsWithQvecs const& collisions, EMCalPhotons const& clusters){ - if(!correctionConfig.doEMCalCalib) { + void processEMCalCalib(CollsWithQvecs const& collisions, EMCalPhotons const& clusters) + { + if (!correctionConfig.doEMCalCalib) { return; } int nColl = 1; @@ -1197,7 +1198,7 @@ struct TaskPi0FlowEMC { } if ((v1.E() - v2.E()) / (v1.E() + v2.E()) > 0.9) { // only use symmetric decays registry.fill(HIST("hClusterCuts"), 6); - registry.fill(HIST("hSparseCalibSE"), vMeson.M(), vMeson.E()/2., getCentrality(collision)); + registry.fill(HIST("hSparseCalibSE"), vMeson.M(), vMeson.E() / 2., getCentrality(collision)); } } if (cfgDoRotation) {