CgemSim-01-00-41/CgemSim/BesCgemHit.hh

Go to the documentation of this file.

//-------------------------------------------------------------------------------------//
//                BOOST --- BESIII Object_Oriented Simulation Tool                     //
//-------------------------------------------------------------------------------------//
/*
 * =====================================================================================
 *               
 *     Filename: BesCgemHit.hh 
 *  Description: 
 *  Conventions:
 *               f_    : variable used in function parameter list
 *               gv_   : global variable
 *               lv_   : local variable used in function
 *               lvd_  : local variable double 
 *               m_*   : normal member of class
 *               sm_*  : static member of class
 *               m_M_* : class data member, Material of each layer
 *               m_N_* : class data member, Number of layers (CgemLayer,GemFoil)
 *               m_R_* : class data member, Radius of each (material) layer, and so on
 *               m_L_* : class data member, Length of each layer or hole pitch
 *               m_T_* : class data member, Thickness of each (material) layer 
 *               m_A_* : class data member, Angle of anode VStrip
 *               m_E_* : class data member, Energy
 *               m_P_* : class data member, P(Momentum)
 *      Version: CgemSim-01-00-00
 *      Created: 01/02/2014 16:59:21 PM
 *     Revision: CgemSim-00-00-01(in CMT version CgemBoss-0.0.1 written by xiuql)
 *     Compiler: gcc
 *       Author: [email protected]
 * Organization: DG1,EPC,IHEP
 *      History:
 *               <Num> <Author>  <Time>     <Version>   <remark>
 *                 0    juxd    20140102    01-00-00   created,
 *
 * =====================================================================================
 */
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
#ifndef BesCgemHit_h
#define BesCgemHit_h 1
 
/* Header file: Geant4 */
#include "G4VHit.hh"
#include "G4THitsCollection.hh"
#include "G4Allocator.hh"
#include "G4ThreeVector.hh"
 
/* Header file: C++ */
#include <vector>
#include "TArrayI.h"
 
using namespace std;
 
class BesCgemHit : public G4VHit
{
  public:
    BesCgemHit();
    BesCgemHit(const BesCgemHit&);
   ~BesCgemHit();
 
   /* Operator */
    const BesCgemHit& operator=(const BesCgemHit&);
    G4int operator==(const BesCgemHit&) const;
 
    inline void* operator new(size_t);
    inline void  operator delete(void*);
 
    /* Function inherit from base class */
    virtual void Draw(){};
    virtual void Print();
 
    /* Set methods */
   //Add by sunxh
    void SetRdtEl               ( G4int         f_RdtElectron ) {m_RdtElectron = f_RdtElectron ;};
    void SetHitID               ( G4int         f_ID_hit      ) {m_ID_hit      = f_ID_hit      ;};
   //end
    void SetTrackID             ( G4int         f_ID_track    ) {m_ID_track    = f_ID_track    ;};
    void SetIsSecondary         ( G4int         isSec         ) {m_isSecondary = isSec         ;};
    void SetLayerID             ( G4int         f_ID_layer    ) {m_ID_layer    = f_ID_layer    ;};
    void SetPDGCode             ( G4int         f_pdg_code    ) {m_pdg_code    = f_pdg_code    ;};
    void SetParentID            ( G4int         f_ID_parent   ) {m_ID_parent    = f_ID_parent  ;};
    void SetGlobalTime          ( G4double      f_global_time ) {m_global_time = f_global_time ;};
    void SetTotalEnergyDeposit  ( G4double      f_E_deposit   ) {m_E_deposit   = f_E_deposit   ;};
    void SetStepLength          ( G4double      f_L_step      ) {m_L_step      = f_L_step      ;};
    void SetPositionOfPrePoint  ( G4ThreeVector f_XYZ_pre     ) {m_XYZ_pre     = f_XYZ_pre     ;};
    void SetPositionOfPostPoint ( G4ThreeVector f_XYZ_post    ) {m_XYZ_post    = f_XYZ_post    ;};
    void SetPositionOfPrePointAlign  ( G4ThreeVector f_XYZ_pre  ) {m_XYZ_pre_align = f_XYZ_pre;};
    void SetPositionOfPostPointAlign ( G4ThreeVector f_XYZ_post ) {m_XYZ_post_align= f_XYZ_post;};
    void SetMomentumOfPrePoint  ( G4ThreeVector f_P_pre       ) {m_P_pre       = f_P_pre       ;};
    void SetMomentumOfPostPoint ( G4ThreeVector f_P_post      ) {m_P_post      = f_P_post      ;};
    //void SetPrePositionInCu     ( G4ThreeVector f_Cu_pre      ) {m_Cu_pre      = f_Cu_pre      ;};
    //void SetPostPositionInCu    ( G4ThreeVector f_Cu_post     ) {m_Cu_post     = f_Cu_post     ;};
    //void SetMomentumOfCuPre     ( G4ThreeVector f_P_Cu_pre    ) {m_P_Cu_pre    = f_P_Cu_pre    ;};
    //void SetMomentumOfCuPost    ( G4ThreeVector f_P_Cu_post   ) {m_P_Cu_post   = f_P_Cu_post   ;};
   //Add by sunxh
    void AddIdentifier          ( G4int  f_ID_Identifier[2000] , G4int N_dim  ) {m_ID_Identifier.Set(N_dim, f_ID_Identifier); };
   //end
    void SetCreatorProcess(G4String proName) {m_creatorProcess=proName;};
    //void SetCurrentTrkID(G4int id)           {m_ID_current_track=id;};
    //void SetCurrentPdgCode(G4int code)       {m_current_pdgcode=id;};
    void SetFlightLengthPrePoint(G4double len)  {m_flightLength_pre =len;};
    void SetFlightLengthPostPoint(G4double len) {m_flightLength_post=len;};
    void AddDigiIdxQ(G4int id) { m_vecDigiIdx_Q.push_back(id); };
    void SetDigiIdxT(G4int id) { m_digiIdx_T=id; };
 
    /* Get methods */
   //Add by sunxh
    G4int         GetRdtEl               () const {return m_RdtElectron ;};
    G4int         GetHitID               () const {return m_ID_hit      ;};
   //end
    G4int         GetTrackID             () const {return m_ID_track    ;};
    G4int         GetIsSecondary         () const {return m_isSecondary ;};
    G4int         GetLayerID             () const {return m_ID_layer    ;};
    G4int         GetPDGCode             () const {return m_pdg_code    ;};
    G4int         GetParentID            () const {return m_ID_parent   ;};
    G4double      GetGlobalTime          () const {return m_global_time ;};
    G4double      GetTotalEnergyDeposit  () const {return m_E_deposit   ;};
    G4double      GetStepLength          () const {return m_L_step      ;};
    G4ThreeVector GetPositionOfPrePoint  () const {return m_XYZ_pre     ;};
    G4ThreeVector GetPositionOfPostPoint () const {return m_XYZ_post    ;};
    G4ThreeVector GetPositionOfPrePointAlign  () const {return m_XYZ_pre_align;};
    G4ThreeVector GetPositionOfPostPointAlign () const {return m_XYZ_post_align;};
    G4ThreeVector GetMomentumOfPrePoint  () const {return m_P_pre       ;};
    G4ThreeVector GetMomentumOfPostPoint () const {return m_P_post      ;};
    //G4ThreeVector GetPrePositionInCu     () const {return m_Cu_pre      ;};
    //G4ThreeVector GetPostPositionInCu    () const {return m_Cu_post     ;};
    //G4ThreeVector GetMomentumOfCuPre     () const {return m_P_Cu_pre    ;};
    //G4ThreeVector GetMomentumOfCuPost    () const {return m_P_Cu_post   ;};
   //Add by sunxh
    TArrayI GetIdentifier                () const {return m_ID_Identifier;};
   //end
   G4String GetCreatorProcess() const {return m_creatorProcess;};
   //G4int    GetCurrentTrkID()   const {return m_ID_current_track;};
   //G4int    GetCurrentPdgCode() const {return m_current_pdgcode;};
   G4int    GetFlightLengthPrePoint()  const {return m_flightLength_pre;};
   G4int    GetFlightLengthPostPoint() const {return m_flightLength_post;};
   const    vector<G4int>& GetVecDigiIdxQ() const {return m_vecDigiIdx_Q;};
   G4int    GetDigiIdxT() const {return m_digiIdx_T;};
 
      
  private:
   //Add by sunxh
    G4int         m_RdtElectron;    
    G4int         m_ID_hit;    
   //end
    G4int         m_ID_track;    /* trk index of current track, could be mother or mother's mother or ... */
    //G4int         m_ID_current_track; /* G4 index of current track */
    //G4int         m_current_pdgcode; /*  pdg code of current track */
    G4int         m_isSecondary; /* if the hit from the track index directly */
    G4int         m_ID_layer;    /* ID of CgemLayer Hits belong to */
    G4int         m_pdg_code;    /* PDG code of the track particle */
    G4int         m_ID_parent;   /* parent ID of the track particle */
    G4double      m_global_time; /* Time since the current event began */
    G4double      m_E_deposit;   /* Total energy deposited during the step */
    G4double      m_L_step;      /* Step length */
    G4ThreeVector m_XYZ_pre;     /* Position of Pre  Point */
    G4ThreeVector m_XYZ_post;    /* Position of Post Point */
    G4ThreeVector m_XYZ_pre_align;     /* Position of Pre  Point after alignment conversion */
    G4ThreeVector m_XYZ_post_align;    /* Position of Post Point after alignment conversion */
    G4ThreeVector m_P_pre;       /* Momentum of Pre  Point */
    G4ThreeVector m_P_post;      /* Momentum of post point */
    G4double      m_flightLength_pre;
    G4double      m_flightLength_post;
    G4String      m_creatorProcess; /* Name of the generation process */
    //G4ThreeVector m_Cu_pre;      //position of pre point in Cu//
    //G4ThreeVector m_Cu_post;     //position of post point in Cu//
    //G4ThreeVector m_P_Cu_pre;    //momentum of pre point in Cu//
    //G4ThreeVector m_P_Cu_post;   //momentum of post point in Cu//
   //Add by sunxh
    TArrayI       m_ID_Identifier;
    vector<G4int> m_vecDigiIdx_Q;// cgem digi indexes whose Q are (partly) from this true hit
    G4int         m_digiIdx_T;   // cgem digi index whose T is from this true hit
   //end
    
};
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
typedef G4THitsCollection<BesCgemHit> BesCgemHitsCollection;
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
extern G4Allocator<BesCgemHit> BesCgemHitAllocator;
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
inline void* BesCgemHit::operator new(size_t)
{
  void *hit;
  hit = (void *) BesCgemHitAllocator.MallocSingle();
  return hit;
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
inline void BesCgemHit::operator delete(void *hit)
{
  BesCgemHitAllocator.FreeSingle((BesCgemHit*) hit);
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
#endif