PndRichGeo.cxx File Reference

#include "PndRichGeo.h"
#include "FairGeoNode.h"
#include "TRandom.h"

Go to the source code of this file.

Functions
	ClassImp (PndRichGeo) PndRichGeo

Function Documentation

ClassImp

(

PndRichGeo

)

Definition at line 5 of file PndRichGeo.cxx.

References Double_t.

  : FairGeoSet()
{
  // Constructor
  // fName has to be the name used in the geometry for all volumes.
  // If there is a mismatch the geometry cannot be build.
  fName="rich";
  maxSectors=0;
  maxModules=10;
   
  fRichOffset = TVector3( 0, 0, 700-35-10/*+50*/ ); //-70
  fAlBoxSize = TVector3( 600, 600, 100 );
  fAlBoxWallThickness = 0.05;
  fAerogelSize = TVector3( 0*590+1*290, 0*590+1*120, 4 );
  fAerogelOffset = TVector3( 0, 0, 1 );
  fnOpt = std::vector<Double_t>(1,1.05);
  fAerogelLayers = std::vector<Double_t>(1,1);
  fAngleExtansionInner = 1.0;
  fAngleExtansionOuter = 1.0;
  fMirrorCurvature = 20;
  fAngleOfMirrorPosition = 50;
  fMirrorThickness = 0.1;
  fMirrorLength = 290;
  fPhDetLength = 290;
  fPhDetThickness = 3;
  fBeamPipeHoleX = 10;
  fBeamPipeHoleY = 10;

   fMirrorType = 0;
   fFlatMirrorZ = std::vector<Double_t>(2);
   fFlatMirrorY = std::vector<Double_t>(2);
   fFlatMirrorZGlob = std::vector<Double_t>(2);
   fFlatMirrorYGlob = std::vector<Double_t>(2);
   fPhDetZ = std::vector<Double_t>(2);
   fPhDetY = std::vector<Double_t>(2);

   fSenseLevel = 4;
   fSensorsPerDevice = 1<<(fSenseLevel-1); // number of sensores per device in one direction
   // Quantum efficiency
   fPhDetDev = 0; // pde_dpc3200_22.dat
   fPhDetDev = 1; // pde_h12700.dat
   std::string workdir(getenv( "VMCWORKDIR" ));
   std::string effFileName = workdir;
   Double_t keff;
   if (fPhDetDev==0)
   {
      effFileName += "/detectors/rich/pde_dpc3200_22.dat";
      fPhDetSizeX = 3.26;
      fPhDetSizeY = 3.26;
      fPhDetGapX = 0.01;
      fPhDetGapY = 0.01;
      keff = 1.7; // measured difference MC-EXP
   }
   if (fPhDetDev==1)
   {
      effFileName += "/detectors/rich/pde_h12700.dat";
      fPhDetSizeX = 5.2;
      fPhDetSizeY = 5.2;
      fPhDetGapX = 0.01;
      fPhDetGapY = 0.01;
      keff = 1;
   }
   std::ifstream from( effFileName.c_str() );
   Double_t wli, pdei;
   from >> wli >> pdei;
   while( !from.eof() ) {
      fWlPhoton.push_back(wli); // nm
      fPDE.push_back(pdei/100.0/keff); // %/100
      from >> wli >> pdei;
   };
   fPhDetEff = new TGraph(fWlPhoton.size(),fWlPhoton.data(),fPDE.data());

}