WireSagCorFunc.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
 
#include "MdtCalibData/WireSagCorFunc.h"
 
#include <cmath>
 
#include "AthenaKernel/getMessageSvc.h"
#include "GaudiKernel/IMessageSvc.h"
#include "GaudiKernel/MsgStream.h"
 
 
using namespace MuonCalib;
 
void WireSagCorFunc::init(const CalibFunc::ParVec& params) {
    MsgStream log(Athena::getMessageSvc(), "WireSagCorFunc");
    m_params = params;
    if (!m_params.empty()) { log << MSG::ERROR << "Wrong number of parameters: " << m_params.size() << endmsg; }
}
 
std::string WireSagCorFunc::name() const { return std::string("WireSagCorFunc"); }
 
double WireSagCorFunc::correction(double signedDriftRadius, double effectiveSag) const {
 
    double coeff[6] = {-4.47741E-3, 1.75541E-2, -1.32913E-2, 2.57938E-3, -4.55015E-5, -1.70821E-7};
 
    double scaleFactor = effectiveSag / 100.;
 
    double dT = 0;
    for (int i = 0; i < 5; i++) { dT = dT + coeff[i] * std::pow(signedDriftRadius, i); }
 
    return scaleFactor * dT;
}