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

Go to the documentation of this file.

//-------------------------------------------------------------------------------------//
//                BOOST --- BESIII Object_Oriented Simulation Tool                     //
//-------------------------------------------------------------------------------------//
/*
 * =====================================================================================
 *               
 *     Filename: BesCgemDigitizer.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
 *      Version: CgemSim-01-00-00
 *      Created: 01/05/2014 09:49: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    20140105    01-00-00   created,
 *
 * =====================================================================================
 */
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
#ifndef BesCgemDigitizer_h
#define BesCgemDigitizer_h 1
 
class BesCgemDigitizerMessenger;
class TFile;
class TH1F;
class G4Svc;
 
/* Header file: BOSS */
#include "BesCgemHit.hh"
#include "BesCgemDigi.hh"
#include "CgemGeomSvc/CgemGeomSvc.h"
 
#include "GaudiKernel/NTuple.h"
#include "GaudiKernel/INTupleSvc.h"
#include "CgemDigitizerSvc/CgemDigitizerSvc.h"
 
/* Header file: Geant4 */
#include "globals.hh"
#include "G4VDigitizerModule.hh"
 
/* Header file: C++ */
#include <vector>
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
class BesCgemDigitizer : public G4VDigitizerModule
{
public:
     BesCgemDigitizer(G4String modName);
     ~BesCgemDigitizer();
 
     /* Realize digitize */
     virtual void Digitize();
     void Digitize_v1();
     void Digitize_v2();
     void Digitize_v3();
 
     /* Lorentz transform */
     void DoLorentzDiffusion(const G4int f_ID_layer, G4ThreeVector &f_XYZ_pre, G4ThreeVector &f_XYZ_post, const int f_F_print=0);
 
     /* Diffusion of hit(or step) */
     void Smear(void);
 
     /* Add noise */
     void AddNoise(void);
 
     /* Return ID from XYZ of readout layer */
     void GetIDFromXYZ(const G4int f_ID_layer, const G4ThreeVector f_XYZ,                /* Input */
               G4double &f_x, G4double &f_v, G4int &f_ID_sheet, G4int &f_ID_x, G4int &f_ID_v); /* Output */
 
     /* Return ID from segment in the same sheet of readout layer */
     void GetIDFromSegmentInSameSheet(const G4int f_ID_x_pre , const G4int f_ID_v_pre, 
                      const G4int f_ID_x_post, const G4int f_ID_v_post, 
                      G4int &f_N_hit_strip_x , G4int &f_N_hit_strip_v, 
                      G4int &f_ID_x_start    , G4int &f_ID_v_start    );
    
     /* Return X and V strip ID of the middle point by the segment which crossed sheet */
     void GetMiddleIDFromSegmentCrossedSheet(const G4int f_ID_layer,
                         const G4int f_ID_sheet_pre , const G4double f_x_pre , const G4double f_v_pre , 
                         const G4int f_ID_sheet_post, const G4double f_x_post, const G4double f_v_post,
                         G4int &f_ID_sheet_mid_pre  , G4int &f_ID_x_mid_pre  , G4int &f_ID_v_mid_pre,
                         G4int &f_ID_sheet_mid_post , G4int &f_ID_x_mid_post , G4int &f_ID_v_mid_post);
 
     /* Return XYZ of readout layer from assigned ID */
     void GetXYZFromID();
 
private:
     BesCgemDigisCollection *m_digi_collection;         
     CgemGeomSvc            *m_cgem_geomsvc;
     CgemDigitizerSvc       *m_cgemDigiSvc;
     G4int m_DigitizerVer;/* Flag of Digitizer version: 2: Geo projection, 3: full digitization by CgemDigitizerSvc */
 
     G4double m_E_threshold; /* Energy threshold */
 
     G4int m_F_lorentz;            /* Flag of Lorentzdiffusion(): 0-OFF; 1-ON */
     G4int m_F_smear;              /* Flag of Smear(): 0-OFF; 1-ON */
     G4int m_F_noise;              /* Flag of AddNoise(): 0-OFF; 1-ON */
     G4int m_F_hit[3][2][2][1500]; /* [layer][sheet][x/v-view][strip] Flag of X/V strip hit: -1(1 hit); >-1(>1 hit) */
     BesCgemHit* m_cgemHitTimePnt[3][2][2][1500]; /* [layer][sheet][x/v-view][strip] pointer to the Cgem hit giving the time */
 
     G4int m_F_printStrip;    /* Flag of Print Strip information: 0-OFF; 1-ON */
     G4int m_F_printHitStrip; /* Flag of Print Hit Strip information: 0-OFF; 1-ON */
     G4int m_F_printDigi;     /* Flag of Print Digi information: 0-OFF; 1-ON */
     G4int m_F_ntuple;        /* Flag of fill ntuple: 0-OFF; 1-ON */
 
     NTuple::Tuple* m_nt1;
     NTuple::Item<int> m_evt;
     NTuple::Item<int> m_nhit;
     NTuple::Array<int> m_layer;
     NTuple::Array<int> m_sheet;
     NTuple::Array<double> m_phi;
     NTuple::Array<double> m_v;
};
#endif
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......