PndRichGeo.h

Go to the documentation of this file.

#ifndef PNDRICHGEO_H
#define PNDRICHGEO_H

#include "FairGeoSet.h"
#include "TVector3.h"
#include "TGraph.h"

class PndRichGeo : public FairGeoSet {

 private:
  TVector3 fRichOffset; 
  TVector3 fAlBoxSize;  
  double fAlBoxWallThickness;

  TVector3 fAerogelSize;                
  TVector3 fAerogelOffset;              
  std::vector<Double_t> fnOpt;          
  std::vector<Double_t> fAerogelLayers; 

  double fAngleExtansionInner;   
  double fAngleExtansionOuter;   
  double fMirrorCurvature;       
  double fAngleOfMirrorPosition; 
  double fMirrorThickness;       
  double fMirrorLength;          

  double fPhDetLength;    
  double fPhDetThickness; 

  double fBeamPipeHoleX; 
  double fBeamPipeHoleY; 

  size_t fMirrorType;

  double fAerogelEntrancePositionZ;
  double fMirrorThetaMin;
  double fMirrorThetaMax;
  double fMirrorRadius;
  TVector3 fMirrorAxis;
  TVector3 fMirrorAxisGlob;

  std::vector<Double_t> fFlatMirrorZ;     
  std::vector<Double_t> fFlatMirrorY;     
  std::vector<Double_t> fFlatMirrorZGlob; 
  std::vector<Double_t> fFlatMirrorYGlob; 

  std::vector<Double_t> fPhDetZ; 
  std::vector<Double_t> fPhDetY; 
  std::vector<Double_t> fWlPhoton;
  std::vector<Double_t> fPDE;
  TGraph *fPhDetEff;

  UInt_t fSenseLevel;
  UInt_t fSensorsPerDevice;
  Int_t fSensorIndexX;
  Int_t fSensorIndexY;
  Int_t fSensorIndex;
  TVector3 fSensorPosition;

  Double_t fPhDetSizeX;
  Double_t fPhDetSizeY;
  Double_t fPhDetGapX;
  Double_t fPhDetGapY;
  UInt_t fPhDetNumX;
  UInt_t fPhDetNumY;
  UInt_t fPhDetPixelNumX;
  UInt_t fPhDetPixelNumY;

  UInt_t fPhDetDev;
  Double_t fPhDetAngle;
  TVector3 fPhDetP0U, fPhDetNxU, fPhDetNyU, fPhDetNzU;
  TVector3 fPhDetP0D, fPhDetNxD, fPhDetNyD, fPhDetNzD;
  Double_t fdX, fdY, fdZ;
  UInt_t fiXmax, fiYmax;

 protected:
  char modName[20]; // name of module
  char eleName[20]; // substring for elements in module
 public:
  PndRichGeo();
  ~PndRichGeo() {}
  const char *getModuleName(Int_t);
  const char *getEleName(Int_t);
  inline Int_t getModNumInMod(const TString &);

  void init(size_t ver = 0);
  Double_t aerogelEntrancePositionZ() { return fAerogelEntrancePositionZ; }
  Double_t mirrorThetaMin() { return fMirrorThetaMin; }
  Double_t mirrorThetaMax() { return fMirrorThetaMax; }
  Double_t mirrorRadius() { return fMirrorRadius; }
  TVector3 mirrorAxis() { return fMirrorAxis; }
  TVector3 mirrorAxisGlob() { return fMirrorAxisGlob; }
  Double_t phDetAngle() { return fPhDetAngle; };
  TVector3 PhDetPositionLocal(TVector3 pos);
  TVector3 PhDetPositionGlobal(TVector3 pos);
  TVector3 PositionDiscretization(TVector3 pos, bool cell = true);
  TVector3 LocalPositionDiscretization(TVector3 pos, Double_t dX = -1, Double_t dY = -1, Double_t dZ = -1);
  UInt_t IndexX(TVector3 pos);
  UInt_t IndexY(TVector3 pos);
  TVector3 PixelPosition(UInt_t ix, UInt_t iy);
  TVector3 PixelPositionLocal(UInt_t ix, UInt_t iy);
  TVector3 PixelPositionGlobal(UInt_t ix, UInt_t iy);

  TVector3 richOffset() { return fRichOffset; }

  TVector3 alBoxSize() { return fAlBoxSize; }

  Double_t alBoxWallThickness() { return fAlBoxWallThickness; }

  TVector3 aerogelSize() { return fAerogelSize; }

  TVector3 aerogelOffset() { return fAerogelOffset; }

  std::vector<Double_t> nOpt() { return fnOpt; }

  std::vector<Double_t> aerogelLayers() { return fAerogelLayers; }

  Double_t angleExtansionInner() { return fAngleExtansionInner; }

  Double_t angleExtansionOuter() { return fAngleExtansionOuter; }

  Double_t mirrorCurvature() { return fMirrorCurvature; }

  Double_t angleOfMirrorPosition() { return fAngleOfMirrorPosition; }

  Double_t mirrorThickness() { return fMirrorThickness; }

  Double_t mirrorLength() { return fMirrorLength; }

  Double_t phDetLength() { return fPhDetLength; }

  Double_t phDetThickness() { return fPhDetThickness; }

  Double_t beamPipeHoleX() { return fBeamPipeHoleX; }

  Double_t beamPipeHoleY() { return fBeamPipeHoleY; }

  std::vector<Double_t> flatMirrorY() { return fFlatMirrorY; }

  std::vector<Double_t> flatMirrorZ() { return fFlatMirrorZ; }

  std::vector<Double_t> flatMirrorYGlob() { return fFlatMirrorYGlob; }

  std::vector<Double_t> flatMirrorZGlob() { return fFlatMirrorZGlob; }

  std::vector<Double_t> phDetY() { return fPhDetY; }

  std::vector<Double_t> phDetZ() { return fPhDetZ; }

  UInt_t sensorIndex() { return fSensorIndex; }

  UInt_t phDetNPixelMaxX() { return fPhDetPixelNumX; }

  UInt_t phDetNPixelMaxY() { return fPhDetPixelNumY; }

  Double_t phDetQEff(Double_t wl);

  Double_t phDetSizeX() { return fPhDetSizeX; }

  Double_t phDetSizeY() { return fPhDetSizeY; }

  Double_t phDetGapX() { return fPhDetGapX; }

  Double_t phDetGapY() { return fPhDetGapY; }

  UInt_t phDetNumX() { return fPhDetNumX; }

  UInt_t phDetNumY() { return fPhDetNumY; }

  UInt_t sensorsPerDevice() { return fSensorsPerDevice; }

  ClassDef(PndRichGeo, 1)
};

inline Int_t PndRichGeo::getModNumInMod(const TString &name)
{
  return (Int_t)(name[4] - '0') - 1; //
}

#endif