MagVolume.h

Go to the documentation of this file.

#ifndef MagVolume_H
#define MagVolume_H

#include "DataFormats/GeometrySurface/interface/GloballyPositioned.h"
#include "MagneticField/VolumeGeometry/interface/VolumeSide.h"
#include "MagneticField/Engine/interface/MagneticField.h"

#include <vector>

template <class T>
class MagneticFieldProvider;

class MagVolume : public GloballyPositioned<float>, public MagneticField {
public:

  typedef GloballyPositioned<float>    Base;
  typedef GloballyPositioned<float>::LocalPoint     LocalPoint;
  typedef GloballyPositioned<float>::LocalVector    LocalVector;
  typedef GloballyPositioned<float>::GlobalPoint    GlobalPoint;
  typedef GloballyPositioned<float>::GlobalVector   GlobalVector;

  MagVolume( const PositionType& pos, const RotationType& rot, 
         const MagneticFieldProvider<float> * mfp,
         double sf=1.) :
    Base(pos,rot), MagneticField(), theProvider(mfp), 
    theProviderOwned(false), theScalingFactor(sf), isIronFlag(false) {}

  ~MagVolume() override;

  LocalVector fieldInTesla( const LocalPoint& lp) const;
  GlobalVector fieldInTesla( const GlobalPoint& lp) const;

  virtual bool inside( const GlobalPoint& gp, double tolerance=0.) const = 0;
  virtual bool inside( const LocalPoint& lp, double tolerance=0.) const {
    return inside( toGlobal(lp), tolerance);
  }

  const MagneticFieldProvider<float>* provider() const {return theProvider;}

  virtual const std::vector<VolumeSide>& faces() const = 0;

  ::GlobalVector inTesla ( const ::GlobalPoint& gp) const override {
    return fieldInTesla( gp);
  }

  bool isIron() const {return isIronFlag;}
  void setIsIron(bool iron) {isIronFlag = iron;}
  void ownsFieldProvider(bool o) {theProviderOwned=o;}

private:

  const MagneticFieldProvider<float> * theProvider;
  bool theProviderOwned;
  double theScalingFactor;
  // Temporary hack to keep information on material. Will eventually be replaced!
  bool isIronFlag;

};

#endif