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