ecal.h

Go to the documentation of this file.

 
#ifndef ECALDETAILED_H
#define ECALDETAILED_H
 
 
#include "ecalPoint.h"
#include "ecalStructure.h"
#include "ecalInf.h"
 
#include "FairDetector.h"
 
#include "TClonesArray.h"
#include "TLorentzVector.h"
#include "TVector3.h"
 
#include <list>
 
class ecalPoint; 
class FairVolume;
class TGeoTranslation;
class ecalLightMap;
 
#define kNumberOfECALSensitiveVolumes 6
const Int_t cMaxModuleType=5;
 
class ecal : public FairDetector
{
 
public:
 
  ecal();
 
 
  ecal(const char* name, Bool_t active, 
      const char* fileGeo="ecal_Detailed.geo");
 
 
  virtual ~ecal();
 
 
  virtual Bool_t  ProcessHits(FairVolume* vol = NULL);
 
 
  virtual void ConstructGeometry();
 
  virtual void EndOfEvent();
  virtual void BeginEvent();
  virtual void Reset();
  virtual void Print() const;
  virtual void CopyClones(TClonesArray* cl1, TClonesArray* cl2, Int_t offset);
  virtual void Register();
  virtual void ChangeHit(ecalPoint* oldHit=NULL);
  virtual void FinishPrimary();
 
  virtual void Initialize();
 
  virtual TClonesArray* GetCollection(Int_t iColl) const;
  virtual void SetSpecialPhysicsCuts();
  
  static Bool_t GetCellCoord(Int_t fVolumeID, Float_t &x, Float_t &y, Int_t& tenergy);
  static Bool_t GetCellCoordInf(Int_t fVolumeID, Float_t &x, Float_t &y, Int_t& tenergy);
  // Get cell for python
  static Bool_t GetCellCoordForPy(Int_t fVolID, TVector3 &all);
protected:
  ecalPoint* AddHit(Int_t trackID, Int_t detID, TVector3 pos,
                       TVector3 mom, Double_t time, Double_t length,
                       Double_t eLoss, Int_t pdgcode);
  ecalPoint* AddLiteHit(Int_t trackID, Int_t detID, Double32_t time, Double32_t eLoss);
 
private:
  Bool_t FillLitePoint(Int_t volnum);
  void FillWallPoint();
  void ResetParameters();
  void SetEcalCuts(Int_t medium);
  ecalPoint* FindHit(Int_t VolId, Int_t TrackId);
 
private:
  ecalInf*  fInf;           
  Option_t* fDebug;         
 
  Int_t GetVolType(Int_t volnum);
  Int_t          fTrackID;              
  Int_t          fVolumeID;             
  TLorentzVector fPos;                  
  TLorentzVector fMom;                  
  Double32_t     fTime;                 
  Double32_t     fLength;               
  Double32_t     fELoss;                
  Int_t          fPosIndex;             
 
  TClonesArray*  fEcalCollection;   
  TClonesArray*  fLiteCollection;   
  Float_t fEcalSize[3];         
  Int_t fSimpleGeo;         
  Int_t fXSize;             
  Int_t fYSize;             
  Float_t fDX;              
  Float_t fDY;              
  Float_t fModuleSize;          
  Float_t fZEcal;           
  Float_t fSemiX;
  Float_t fSemiY;
  Float_t fThicknessLead;       
  Float_t fThicknessScin;       
  Float_t fThicknessTyvk;       
  Float_t fThicknessLayer;      
  Float_t fThicknessSteel;      
  Float_t fEdging;          
  Float_t fHoleRad;         
  Float_t fFiberRad;            
  Float_t fXCell[cMaxModuleType];   
  Float_t fYCell[cMaxModuleType];   
  Int_t fNH[cMaxModuleType];        
  Int_t fCF[cMaxModuleType];        
  TString fLightMapNames[cMaxModuleType];       
  ecalLightMap* fLightMaps[cMaxModuleType];     
  Int_t   fNLayers;         
  Float_t fModuleLenght;        
  Float_t fGeoScale;            
  Int_t   fNColumns1;           
  Int_t   fNRows1;          
  Int_t   fNColumns2;           
  Int_t   fNRows2;          
  Int_t   fNColumns;            
  Int_t   fNRows;           
  Int_t   fVolIdMax;            
  Int_t fFirstNumber;           
  Int_t fVolArr[kNumberOfECALSensitiveVolumes];     
 
  TGeoVolume* ConstructRaw(Int_t number);
  void ConstructModule(Int_t type);
  void ConstructCell(Int_t type);
  void ConstructTile(Int_t type, Int_t material);
  void ConstructModuleSimple(Int_t type);
  void ConstructCellSimple(Int_t type);
  void ConstructTileSimple(Int_t type, Int_t material);
  TGeoVolume* fModules[cMaxModuleType];         
  TGeoVolume* fCells[cMaxModuleType];           
  TGeoVolume* fScTiles[cMaxModuleType];         
  TGeoVolume* fTileEdging[cMaxModuleType];      
  TGeoVolume* fPbTiles[cMaxModuleType];         
  TGeoVolume* fTvTiles[cMaxModuleType];         
  TGeoVolume* fHoleVol[3];              
  TGeoVolume* fFiberVol[3];             
  TGeoVolume* fSteelTapes[2];               
  TGeoTranslation** fHolePos[cMaxModuleType];       
  Int_t fModulesWithType[cMaxModuleType];       
  std::list<std::pair<Int_t, TGeoVolume*> > fRawNumber; 
 
  Int_t fStructureId;           
  Int_t InitMedium(const char* name);
  void InitMedia();
 
  ecal(const ecal&);
  ecal& operator=(const ecal&);
 
  ClassDef(ecal,1)
 
};
 
inline void ecal::ResetParameters()
{
  fTrackID = fVolumeID = 0;
  fPos.SetXYZM(0.0, 0.0, 0.0, 0.0);
  fMom.SetXYZM(0.0, 0.0, 0.0, 0.0);
  fTime = fLength = fELoss = 0;
  fPosIndex = 0;
};
 
 
#endif