From a9b0e6a77f4beab5a8f5682cba01afa49e76136d Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Nicol=C3=B2=20Jacazio?= Date: Fri, 30 Jul 2021 14:49:45 +0200 Subject: [PATCH] Clean FPE hack in AO2D writing --- Detectors/AOD/src/AODProducerWorkflowSpec.cxx | 17 ----------------- 1 file changed, 17 deletions(-) diff --git a/Detectors/AOD/src/AODProducerWorkflowSpec.cxx b/Detectors/AOD/src/AODProducerWorkflowSpec.cxx index 0185eb387a131..9059c2621f67a 100644 --- a/Detectors/AOD/src/AODProducerWorkflowSpec.cxx +++ b/Detectors/AOD/src/AODProducerWorkflowSpec.cxx @@ -408,23 +408,6 @@ void AODProducerWorkflowDPL::fillMCParticlesTable(o2::steer::MCKinematicsReader& float pY = (float)mcParticles[particle].Py(); float pZ = (float)mcParticles[particle].Pz(); float energy = (float)mcParticles[particle].GetEnergy(); - // HACK to avoid FPE in expression columns. Affect only particles in the Beam Pipe. - // TO BE REMOVED asap - { - const float limit = 1e-4; - const float mom = TMath::Sqrt(pX * pX + pY * pY + pZ * pZ); - const float eta = 0.5f * TMath::Log((mom + pZ) / (mom - pZ)); - if (TMath::Abs(eta) > 0.9) { - if (TMath::Abs((mom - pZ) / pZ) <= limit) { - pX = truncateFloatFraction(TMath::Sqrt((1.f + limit) * (1.f + limit) - 1.f) * pZ * 0.70710678, mMcParticleMom); - pY = truncateFloatFraction(TMath::Sqrt((1.f + limit) * (1.f + limit) - 1.f) * pZ * 0.70710678, mMcParticleMom); - } - if (TMath::Abs(energy - pZ) < limit) { - energy = truncateFloatFraction(pZ + limit, mMcParticleMom); - } - } - } - // End of HACK mcParticlesCursor(0, mccolid,