From c13786df0b4e6e962d4821906737d3c33b1698a1 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 18 Apr 2024 20:08:52 +0000 Subject: [PATCH] Please consider the following formatting changes --- Detectors/MUON/MCH/Simulation/src/DEDigitizer.cxx | 3 +-- Detectors/MUON/MCH/Simulation/src/Response.cxx | 8 ++++---- 2 files changed, 5 insertions(+), 6 deletions(-) diff --git a/Detectors/MUON/MCH/Simulation/src/DEDigitizer.cxx b/Detectors/MUON/MCH/Simulation/src/DEDigitizer.cxx index 49379bbb7ab1c..ff9ebdd127223 100644 --- a/Detectors/MUON/MCH/Simulation/src/DEDigitizer.cxx +++ b/Detectors/MUON/MCH/Simulation/src/DEDigitizer.cxx @@ -20,7 +20,6 @@ #include #include "Field/MagneticField.h" - /// Convert collision time to ROF time (ROF duration = 4 BC) std::pair time2ROFtime(const o2::InteractionRecord& time) { @@ -84,7 +83,7 @@ void DEDigitizer::processHit(const Hit& hit, const InteractionRecord& collisionT auto localY = lpos.Y(); // calculate angle between track and wire assuming wire perpendicular to z-axis - auto thetawire = atan2(lentrance.Y() - lexit.Y(), hitlengthZ);//to be checked! + auto thetawire = atan2(lentrance.Y() - lexit.Y(), hitlengthZ); // to be checked! // local b-field double b[3] = {0., 0., 0.}; diff --git a/Detectors/MUON/MCH/Simulation/src/Response.cxx b/Detectors/MUON/MCH/Simulation/src/Response.cxx index 50519201c4497..dd55810ec9085 100644 --- a/Detectors/MUON/MCH/Simulation/src/Response.cxx +++ b/Detectors/MUON/MCH/Simulation/src/Response.cxx @@ -107,7 +107,7 @@ float Response::inclandbfield(float thetawire, float betagamma, float bx) const sigmaEffectThetadegrees = sigmaEffect10degrees / angleEffectNorma(thetawire * TMath::RadToDeg()); if (o2::mch::Station() == o2::mch::Station::Type1) { sigmaEffectThetadegrees /= 1.09833 + 0.017 * (thetawire * TMath::RadToDeg()); - } + } yAngleEffect = 0.0001 * gRandom->Gaus(0, sigmaEffectThetadegrees); // error due to the angle effect in cm } } else { @@ -118,7 +118,7 @@ float Response::inclandbfield(float thetawire, float betagamma, float bx) const sigmaEffectThetadegrees = sigmaEffect10degrees / magAngleEffectNorma(thetawire * TMath::RadToDeg(), bx / 10.); // check b-field unit in aliroot and O2 if (o2::mch::Station() == o2::mch::Station::Type1) { sigmaEffectThetadegrees /= 1.09833 + 0.017 * (thetawire * TMath::RadToDeg()); - } + } yAngleEffect = 0.0001 * gRandom->Gaus(0, sigmaEffectThetadegrees); } } @@ -131,7 +131,7 @@ float Response::eLossRatio(float logbetagamma) const // Ratio of particle mean eloss with respect MIP's Khalil Boudjemline, sep 2003, PhD.Thesis and Particle Data Book /// copied from aliroot AliMUONv1.cxx float eLossRatioParam[5] = {1.02138, -9.54149e-02, +7.83433e-02, -9.98208e-03, +3.83279e-04}; - return eLossRatioParam[0] + eLossRatioParam[1] * logbetagamma + eLossRatioParam[2] * logbetagamma * logbetagamma + eLossRatioParam[3] * logbetagamma * logbetagamma * logbetagamma + eLossRatioParam[4] * logbetagamma * logbetagamma * logbetagamma * logbetagamma; + return eLossRatioParam[0] + eLossRatioParam[1] * logbetagamma + eLossRatioParam[2] * logbetagamma * logbetagamma + eLossRatioParam[3] * logbetagamma * logbetagamma * logbetagamma + eLossRatioParam[4] * logbetagamma * logbetagamma * logbetagamma * logbetagamma; } //_____________________________________________________________________ float Response::angleEffect10(float elossratio) const @@ -158,5 +158,5 @@ float Response::magAngleEffectNorma(float angle, float bfield) const /// copied from aliroot AliMUONv1.cxx float angleEffectParam[7] = {8.6995, 25.4022, 13.8822, 2.4717, 1.1551, -0.0624, 0.0012}; float aux = std::abs(angle - angleEffectParam[0] * bfield); - return 121.24 / ((angleEffectParam[1] + angleEffectParam[2] * std::abs(bfield)) + angleEffectParam[3] * aux + angleEffectParam[4] * aux * aux + angleEffectParam[5] * aux * aux * aux + angleEffectParam[6] * aux * aux * aux *aux); + return 121.24 / ((angleEffectParam[1] + angleEffectParam[2] * std::abs(bfield)) + angleEffectParam[3] * aux + angleEffectParam[4] * aux * aux + angleEffectParam[5] * aux * aux * aux + angleEffectParam[6] * aux * aux * aux * aux); }