/data/refman/pasoursint/CMSSW_4_2_9_HLT1_bphpatch4/src/MagneticField/VolumeBasedEngine/interface/MagGeometry.h

Go to the documentation of this file.

#ifndef MagGeometry_H
#define MagGeometry_H

#include "DataFormats/GeometrySurface/interface/BoundPlane.h"
#include "MagneticField/Layers/src/MagBinFinders.h"
#include "DetectorDescription/Core/interface/DDCompactView.h"

#include <vector>

class MagBLayer;
class MagESector;
class MagVolume;
class MagVolume6Faces;
template <class T> class PeriodicBinFinderInPhi;
namespace edm {class ParameterSet;}

class MagGeometry {
public:

  typedef Surface::GlobalVector   GlobalVector;
  typedef Surface::GlobalPoint    GlobalPoint;

  MagGeometry(const edm::ParameterSet& config, std::vector<MagBLayer *> ,
                             std::vector<MagESector *> ,
                             std::vector<MagVolume6Faces*> ,
                             std::vector<MagVolume6Faces*> );

  ~MagGeometry();

  GlobalVector fieldInTesla(const GlobalPoint & gp) const;

  MagVolume * findVolume(const GlobalPoint & gp, double tolerance=0.) const;

  bool isZSymmetric() const;

  // FIXME: only for temporary tests, should be removed.
  const std::vector<MagVolume6Faces*> & barrelVolumes() const {return theBVolumes;}
  const std::vector<MagVolume6Faces*> & endcapVolumes() const {return theEVolumes;}

private:

  friend class MagGeometryExerciser; // for debug purposes

  // Linear search (for debug purposes only)
  MagVolume * findVolume1(const GlobalPoint & gp, double tolerance=0.) const;


  bool inBarrel(const GlobalPoint& gp) const;

  mutable MagVolume * lastVolume; // Cache last volume found

  std::vector<MagBLayer *> theBLayers;
  std::vector<MagESector *> theESectors;

  // FIXME: only for temporary tests, should be removed.
  std::vector<MagVolume6Faces*> theBVolumes;
  std::vector<MagVolume6Faces*> theEVolumes;

  MagBinFinders::GeneralBinFinderInR<double>* theBarrelBinFinder;
  PeriodicBinFinderInPhi<float> * theEndcapBinFinder;

  bool cacheLastVolume;
  bool v_85l;
};
#endif