hcal.h

Go to the documentation of this file.

 
#ifndef HCAL_H
#define HCAL_H
 
 
#include "hcalPoint.h"
#include "hcalStructure.h"
#include "hcalInf.h"
 
#include "FairDetector.h"
 
#include "TClonesArray.h"
#include "TLorentzVector.h"
#include "TVector3.h"
 
#include <list>
 
class hcalPoint; 
class FairVolume;
class TGeoTranslation;
class hcalLightMap;
 
#define kNumberOfHCALSensitiveVolumes 6
 
class hcal : public FairDetector
{
 
public:
 
  hcal();
 
 
  hcal(const char* name, Bool_t active, 
      const char* fileGeo="hcal.geo");
 
 
  virtual ~hcal();
 
 
  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(hcalPoint* 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& section);
  static Bool_t GetCellCoordInf(Int_t fVolumeID, Float_t &x, Float_t &y, Int_t& section);
protected:
  hcalPoint* AddHit(Int_t trackID, Int_t detID, TVector3 pos,
                       TVector3 mom, Double_t time, Double_t length,
                       Double_t eLoss, Int_t pdgcode);
  hcalPoint* 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 SetHcalCuts(Int_t medium);
  hcalPoint* FindHit(Int_t VolId, Int_t TrackId);
 
private:
  hcalInf*  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*  fHcalCollection;   
  TClonesArray*  fLiteCollection;   
  Float_t fHcalSize[3];         
  Int_t fSimpleGeo;         
  Int_t fFastMC;            
  Int_t fXSize;             
  Int_t fYSize;             
  Float_t fDX;              
  Float_t fDY;              
  Float_t fModuleSize;          
  Float_t fZHcal;           
  Float_t fSemiX;
  Float_t fSemiY;
  TString fAbsorber;            
  Float_t fThicknessAbsorber;       
  Float_t fThicknessScin;       
  Float_t fThicknessTyvk;       
  Float_t fThicknessLayer;      
  Float_t fThicknessSteel;      
  Float_t fEdging;          
  Float_t fHoleRad;         
  Float_t fFiberRad;            
  Float_t fXCell;           
  Float_t fYCell;           
  Int_t fNH;                
  Int_t fCF;                
  TString fLightMapName;        
  hcalLightMap* fLightMap;      
  Int_t   fNLayers;         
  Int_t   fNLayers1;            
  Float_t fModuleLength;        
  Int_t   fVolIdMax;            
  Int_t fFirstNumber;           
  Int_t fVolArr[kNumberOfHCALSensitiveVolumes];     
 
  TGeoVolume* ConstructRaw(Int_t number);
  void ConstructModule();
  void ConstructTile(Int_t material);
  void ConstructModuleSimple();
  void ConstructTileSimple(Int_t material);
  TGeoVolume* fModule;          
  TGeoVolume* fScTile;          
  TGeoVolume* fTileEdging;      
  TGeoVolume* fPbTile;          
  TGeoVolume* fTvTile;          
  TGeoVolume* fHoleVol[3];              
  TGeoVolume* fFiberVol[3];             
  TGeoVolume* fSteelTapes[2];               
  TGeoTranslation** fHolePos;       
  Int_t fModules;           
  std::list<std::pair<Int_t, TGeoVolume*> > fRawNumber; 
 
  Int_t fStructureId;           
  Int_t InitMedium(const char* name);
  void InitMedia();
 
  hcal(const hcal&);
  hcal& operator=(const hcal&);
 
  ClassDef(hcal,1)
};
 
inline void hcal::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