SamplingGTS.h

Go to the documentation of this file.

#ifndef SAMPLINGGTS_H
#define SAMPLINGGTS_H
 
#include "CgemGeomSvc/ICgemGeomSvc.h"
#include "CgemDigitizerSvc/DriftAndAvalanche.h"
#include "CgemDigitizerSvc/DiffusionGTS.h"
 
#include <vector>
#include <string>
 
#include "TFile.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TNtuple.h"
 
#define m_testing_sam false
 
 
class G4Svc;
 
class SamplingGTS: public DriftAndAvalanche {
public:
     SamplingGTS();
     ~SamplingGTS();
 
     void init(ICgemGeomSvc* geomSvc, double magConfig);
 
     void setIonElectrons(int layer, int nElectrons, std::vector<double> x, std::vector<double> y, std::vector<double> z, std::vector<double> t);
 
     /* output info of multiplied eletrons */
     int    getNelectrons() const {return m_nMulElec;}
     Float_t getX(int n) const {return m_eX[n];}
     Float_t getY(int n) const {return m_eY[n];}
     Float_t getZ(int n) const {return m_eZ[n];}
     Float_t getT(int n) const {return m_eT[n];}
 
     const std::vector<Float_t> & getXContainer() const {return m_eX;}
     const std::vector<Float_t> & getYContainer() const {return m_eY;}
     const std::vector<Float_t> & getZContainer() const {return m_eZ;}
     const std::vector<Float_t> & getTContainer() const {return m_eT;}
 
     void setDebugging(bool debugging) {m_debugging = debugging;}
 
     // --------------------------------
     // CHECK 11) per ora li metto hard coded
     bool readGemParameters();
     bool readGasParameters();
     
     bool is_survived();
     double compute_gain();
     void compute_diffusion_on_GEM3(double xi, double yi, double zi, double ti, double &xf, double &yf, double &zf, double &tf);
 
     DiffusionGTS *diffusion;
 
private:
     void clear();
     
     bool m_debugging;
     int    m_hv;
     bool   m_field;
     double m_tuning_factor_gain;
     double m_tuning_factor_diff_perp;
     double m_tuning_factor_diff_paral;
 
     double m_eff_col[3]; // efficiency of collection
     double m_eff_ext[3]; // efficiency of extraction
     TH1F *h_gain_eff;
 
     ICgemGeomSvc* m_geomSvc;
     double m_magConfig;        /* magnetic field */
     // x_local
     int n_shift_x_drift;
     double m_shift_x_drift_par[10];
     double m_shift_x_transf_par;
     double m_shift_x_induct_par;
     int n_sigma_x_drift;
     double m_sigma_x_drift_par[10];
     double m_sigma_x_transf_par;
     double m_sigma_x_induct_par;
     // y_local only diffusion, no field
     int n_sigma_y_drift;
     double m_sigma_y_drift_par[10];
     double m_sigma_y_transf_par;
     double m_sigma_y_induct_par;
     // t
     int n_shift_t_drift;
     double m_shift_t_drift_par[10];
     double m_shift_t_transf_par;
     double m_shift_t_induct_par;
     int n_sigma_t_drift;
     double m_sigma_t_drift_par[10];
     double m_sigma_t_transf_par;
     double m_sigma_t_induct_par;
 
     /* info of multiplied electrons */
     int m_nMulElec;
     std::vector<Float_t> m_eX;
     std::vector<Float_t> m_eY;
     std::vector<Float_t> m_eZ;
     std::vector<Float_t> m_eT;
 
     TFile *output;
     TNtuple *tree;
 
 
};
 
#endif