PndHyp.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                     CbmHyp header file                        -----
// -----               created by A. Sanchez                  -----
// -------------------------------------------------------------------------

#ifndef PNDHYP_H
#define PNDHYP_H


#include "TClonesArray.h"
#include "TVector3.h"
#include "TLorentzVector.h"
#include "FairDetector.h"
#include "PndGeoHypPar.h"

#include "TRandom.h"

#include "PndHypGeoHandling.h"
//#include "PndHypDecayer.h"
#include "HypStatDecay.h"
using namespace std;


class TClonesArray;
class PndHypPoint;
class PndHypSecTarPoint;
//class PndHypSTMatBudPoint;
class FairVolume; 

class PndHyp : public FairDetector 
{

 public:

  PndHyp();


  PndHyp(const char* name, Bool_t active);


  virtual ~PndHyp();


  virtual void Initialize();


  virtual Bool_t ProcessHits(FairVolume* vol = 0);

  virtual void SetSpecialPhysicsCuts();
  virtual void EndOfEvent();


  virtual void BeginEvent();
  virtual void Register();


  virtual TClonesArray* GetCollection(Int_t iColl) const;


  virtual void Print() const;    


  virtual void Reset();


  virtual void CopyClones(TClonesArray* cl1, TClonesArray* cl2,
                          Int_t offset);


  virtual void ConstructGeometry();
  //void ConstructASCIIGeometry();
  virtual void FinishRun();

  PndHypPoint* AddHit(Int_t trackID,Int_t evtID, 
                      Int_t detID, TString detName,
                      TVector3 posin, 
                      TVector3 momin,
                      TVector3 posout, 
                      TVector3 momout,
                      //TVector3 posInLocal, 
                      //TVector3 posOutLocal,
                      Double_t tof, 
                      Double_t length, 
                      Double_t eLoss,
                      Double_t charge, 
                      Double_t mass, 
                      Int_t pdgCode,
                      Double_t dist,
                      Double_t PLin,
                      Double_t PLout);

  PndHypPoint* AddSecTarHit(Int_t trackID, Int_t evtID,
                            Int_t detID,TString detName,
                            TVector3 posin, 
                            TVector3 momin,
                            TVector3 posout, 
                            TVector3 momout,
                            //TVector3 posInLocal, 
                            //TVector3 posOutLocal,
                            Double_t tof, 
                            Double_t length,
                            Double_t eLoss,
                            Double_t charge, 
                            Double_t mass, 
                            Int_t pdgCode,
                            Double_t dist,
                            Double_t PLin,
                            Double_t PLout);


 PndHypPoint* AddSTMatBudHit(Int_t trackID, Int_t evtID,
                            Int_t detID,TString detName,
                            TVector3 posin, 
                            TVector3 momin,
                            TVector3 posout, 
                            TVector3 momout,
                            //TVector3 posInLocal, 
                            //TVector3 posOutLocal,
                            Double_t tof, 
                            Double_t length,
                            Double_t eLoss,
                            Double_t charge, 
                            Double_t mass, 
                            Int_t pdgCode,
                            Double_t dist,
                            Double_t PLin,
                            Double_t PLout);

  
  
  void PreTrack();

  void SetHypSDtoFile(bool onf,bool val)
  {
    fUseRAZHOption = onf;// **switch On/Off HypStatDecay**
    fUseFileOption = val;
  };
  
  void SetHypGamFEm(bool gam )
  {
    fUseGamOption = gam;
  }

   void SetMatbud(bool mbud )
  {
    fMatBud = mbud;
  }

  void SetTreeFName(const Char_t* Name){ fFileName    = Name; };

  void SetSensorVol(TString VolSi){ 
    fVolNamSi=VolSi.Data();

  }
  void SetAbsorberVol(TString VolAb){ 
    fVolNamAb=VolAb.Data();

  }

  void SetStartEvID(Int_t EvID=0){ 
    fStartEvID=EvID;
    cout<<" starting value for EvId "<<fStartEvID<<endl;

  }

  void SetGeoVersion(TString vers="standard"){
    fVers=vers;
  }

  void SetAbsMat(TString abmat="HYPdiamond"){
    fAbsMat=abmat;
  }

 void SetBPipeMat(TString bpmat="HYPcarbon"){
    fBPipeMat=bpmat;
  }
   
 void SetListMat(TString mat="carbon"){
   fListMat=kTRUE;
   fListOfMaterials.push_back(mat.Data());
  }
  

  private:
  
  std::vector<std::string> fListOfSensitives;
  std::vector<TString> fListOfMaterials;

  bool CheckIfSensitive(std::string name);

  TString fVolNamSi,fVolNamAb,fVers,fAbsMat,fSiMat,fBPipeMat;
  Bool_t fRootSensVol,fListMat,fStandard,fCurrent;
  Int_t fStartEvID;
  //PndGeoHypPar *par;
  Int_t          fTrackID;           //  track index
  Int_t          fVolumeID;          //  volume id
  Int_t          fEventID;           //  event id
  TLorentzVector fPosIn;               //   position
  TLorentzVector fMomIn;               //  momentum
  TLorentzVector fPosOut;               //  position
  TLorentzVector fMomOut;               //  momentum
  /*  TLorentzVector fPosInLocal;        // entry position in module frame */
  /*   TLorentzVector fPosOutLocal;       //  exit position in module frame */
  Double_t fPLout,fPLin;               //  total momentum
 
  
  TClonesArray* fEvt; 
  TFile* fFile; 
  TTree* ft;
 
  
  Int_t        activeCnt;
  Double_t     weight;
  Double_t     seed;
  Int_t          fcount;
  
  Double32_t     fTime;              //   time
  Double32_t     fLength;            //   length
  Double32_t     fELoss;             //   energy loss
  Int_t fPosIndex;      // 
  Int_t fpdgCode;      //  MC volume ID of MUO
  Int_t  SiId,CId,alId,beId,CpipeId ;     
  
  Double_t fcharge;  
  Double_t fmass, fdist;

  const Char_t* fFileName;   

   PndHypGeoHandling* fGeoH;          
  
  //PndHypDecayer* fread;          //! Gives Access to the Statistical decay products
  //HypStatDecay* fread;          //! Gives Access to the Statistical decay products

  TClonesArray* fHypCollection;        
  TClonesArray* fHypSecTarCollection;        // Hit collection(Absorver)
  TClonesArray* fHypSTMatBudCollection;        // Hit collection(pipehyp)
  TRandom r; 

  // reset all parameters   
  void ResetParameters();
  Bool_t  fTrackStopNxtStep;
   Bool_t  fUseFileOption;
  Bool_t fUseRAZHOption;
  Bool_t fUseGamOption;
  Bool_t fMatBud;
  void SetHypStatDecay(bool cal,bool active);

  ClassDef(PndHyp,15)

}; 


inline void PndHyp::ResetParameters() {
  fTrackID = fVolumeID = fEventID =-1;
  fPosIn.SetXYZM(0.0, 0.0, 0.0, 0.0);
  fPosOut.SetXYZM(0.0, 0.0, 0.0, 0.0);
  fMomIn.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
  fMomOut.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
 /*  fPosInLocal.SetXYZM(0.0, 0.0, 0.0, 0.0); */
/*   fPosOutLocal.SetXYZM(0.0, 0.0, 0.0, 0.0); */
  fTime = fLength = fELoss = 0;

};



#endif