OAE85lParametrizedMagneticField.h

Go to the documentation of this file.

#ifndef ParametrizedEngine_OAE85lParametrizedMagneticField_h
#define ParametrizedEngine_OAE85lParametrizedMagneticField_h

#include "MagneticField/Engine/interface/MagneticField.h"

namespace edm { class ParameterSet; }

class OAE85lParametrizedMagneticField : public MagneticField {
 public:

  OAE85lParametrizedMagneticField(float b0_ = 40.681, 
                  float a_ = 4.6430, 
                  float l_ = 15.284);

  OAE85lParametrizedMagneticField(const edm::ParameterSet& parameters);


  virtual ~OAE85lParametrizedMagneticField();
  
  GlobalVector inTesla (const GlobalPoint& gp) const;

  GlobalVector inTeslaUnchecked (const GlobalPoint& gp) const;

  bool isDefined(const GlobalPoint& gp) const;

 private:

  void ffunkti(float u, float* ff) const;

  void init();

  //phenomen. parameters:
  // b0=field at centre, l=solenoid length, a=radius (m) 
  float b0;
  float l;
  float a;
  // Derived parameters
  float ap2;
  float hb0;
  float hlova;
  float ainv;
};
#endif