AtlasBComponent.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#include "CLHEP/Units/SystemOfUnits.h"
#include "AtlasBComponent.h"
 
 
using CLHEP::kilogauss;
 
  
namespace Genfun {
  FUNCTION_OBJECT_IMP(AtlasBComponent)
    // This class is an adaptor:
    
  // Constructor:
  AtlasBComponent::AtlasBComponent(unsigned int index, MagField::AtlasFieldCache* fieldCache)
  : AbsFunction()
  , m_index(index)
    , m_fieldCache(fieldCache) {}
  
  // Destructor:
  AtlasBComponent::~AtlasBComponent() {}
  
  // Copy Constructor:
  AtlasBComponent::AtlasBComponent(const AtlasBComponent &right)
   : AbsFunction()
   , m_index(right.m_index)
   , m_fieldCache(right.m_fieldCache) {}
  
  // Dimensionality:
  unsigned int AtlasBComponent::dimensionality() const { return 3;}
  
  double AtlasBComponent::operator () (double) const { throw "Invalid parameter!"; }   
 
  double AtlasBComponent::operator () (const Argument & a) const {
    // Changes from doubles to floats and back
    double XYZ_in_mm[3] , BXYZ_in_kgmm[3];
    for (int i=0;i<3;++i) *(XYZ_in_mm+i) = a[i];
    m_fieldCache->getField( XYZ_in_mm , BXYZ_in_kgmm ); // const-violation, see comment in header file
    return BXYZ_in_kgmm[m_index]*kilogauss;
  }
 
}