Skip to content

Commit

Permalink
Please consider the following formatting changes
Browse files Browse the repository at this point in the history
  • Loading branch information
alibuild committed Apr 18, 2024
1 parent 13eca33 commit c13786d
Show file tree
Hide file tree
Showing 2 changed files with 5 additions and 6 deletions.
3 changes: 1 addition & 2 deletions Detectors/MUON/MCH/Simulation/src/DEDigitizer.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -20,7 +20,6 @@
#include <TGeoGlobalMagField.h>
#include "Field/MagneticField.h"


/// Convert collision time to ROF time (ROF duration = 4 BC)
std::pair<o2::InteractionRecord, uint8_t> time2ROFtime(const o2::InteractionRecord& time)
{
Expand Down Expand Up @@ -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.};
Expand Down
8 changes: 4 additions & 4 deletions Detectors/MUON/MCH/Simulation/src/Response.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -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 {
Expand All @@ -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);
}
}
Expand All @@ -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
Expand All @@ -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);
}

0 comments on commit c13786d

Please sign in to comment.