/home/mweber/Dokumente/CMS/Alignment/kalmanalignment/trunk/simulation/Detector.h

Go to the documentation of this file.

#ifndef _Detector_h
#define _Detector_h

// STL includes
#include <vector>

// C/C++ includes
#include <math.h>

// local includes
#include "Det.h"

class Detector;

// defined in Detector.cc
extern Detector * gDetector;

class Detector
{
 public:
  Detector(const char * fname);
  virtual ~Detector();

  const int    GetNLayer() const { return fNLayer; };
  const int    GetNRod() const { return fNRod; };
  const int    GetNModule() const { return fNModule; };
  const double GetLayerDistance() const { return fLayerDistance; };

  // detector size
  const double GetWidth() const { return fNRod*fModuleWidth; };
  const double GetLength() const { return fNModule*fModuleLength; };
  const double GetHeight() const { return fNLayer*fLayerDistance; };

  const double GetMisalignPosition() const { return fMisalignPosition; };
  const double GetMisalignRotation() const { return fMisalignRotation; };
  const CLHEP::HepVector & GetTrueMisalign() const { return fTrueMisalignment; }

  int MakeID(int layer, int rod, int module) const;

  const Det & GetDet(int layer, int rod, int module) const;
  void ResetDet(int moduleID);

 protected:
  // Detector layout
  int    fNLayer;
  int    fNRod;
  int    fNModule;
  double fLayerDistance;
  bool   fRandomPlacement;

  // Module layout
  double fModuleWidth;
  double fModuleLength;
  double fModuleThickness;
  double fModulePitchX;
  double fModulePitchY;

  // Misalignment
  double fMisalignPosition;
  double fMisalignRotation;

  std::vector<Det *> DetVec;

  CLHEP::HepVector fTrueMisalignment;

 private:
  void ReadConfiguration(const char * fname);
  void CreateModules();
  Det * CreateDetector(int ID, int layer, int rod, int module);
  void MisalignDetector(Det * adet);
};

#endif