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>
#include <atomic>

class MagBLayer;
class MagESector;
class MagVolume;
class MagVolume6Faces;
template <class T> class PeriodicBinFinderInPhi;

class MagGeometry {
public:

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

  MagGeometry(int geomVersion, const std::vector<MagBLayer *>& ,
                 const std::vector<MagESector *>& ,
                 const std::vector<MagVolume6Faces*>& ,
                 const std::vector<MagVolume6Faces*>& );
  MagGeometry(int geomVersion, const std::vector<MagBLayer const*>& ,
                 const std::vector<MagESector const*>& ,
                 const std::vector<MagVolume6Faces const*>& ,
                 const std::vector<MagVolume6Faces const*>& );

  ~MagGeometry();

  GlobalVector fieldInTesla(const GlobalPoint & gp) const;

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

  // Deprecated, will be removed
  bool isZSymmetric() const {return false;}

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

private:

  friend class MagGeometryExerciser; // for debug purposes

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


  bool inBarrel(const GlobalPoint& gp) const;

  mutable std::atomic<MagVolume const*> lastVolume; // Cache last volume found

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

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

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

  bool cacheLastVolume;
  int geometryVersion;
};
#endif