IonizationGTS.h

Go to the documentation of this file.

#ifndef IONIZATIONGTS_H
#define IONIZATIONGTS_H
 
#include "CgemGeomSvc/ICgemGeomSvc.h"
 
#include "CgemDigitizerSvc/Ionization.h"
 
#include <vector>
#include "TH1F.h"
#include "TFile.h"
 
#define m_testing false
class IonizationGTS: public Ionization{
public:
     IonizationGTS();
     ~IonizationGTS();
 
     void init(unsigned int random, ICgemGeomSvc* geomSvc, double magConfig);
     void setDebugging(bool debugging) {m_debugging = debugging;}
 
     /* return the number of electrons from ionization */
     /* particle = 0,1,2,3,4 (e,mu,pi,K,p) */
     /* p: momentum (GeV/c) */
     /* trkPosIn & trkPosOut:  starting & exit positions (x,y,z) of the track segment in drift area */
     /* ex, ey, ez: x, y, z of the electron from ionization */
     void setTrack(int particle, int charge, double p, double trkPosIn[], double trkPosOut[]);
 
     int getNumberIonE() {return m_nIonE;}
     double getEx(int nElec) {return m_ex[nElec];}
     double getEy(int nElec) {return m_ey[nElec];}
     double getEz(int nElec) {return m_ez[nElec];}
     double getEt(int nElec) {return m_et[nElec];}
 
     // ------------------------------------------
     int getNumberIonC() {return m_nIonC;}
     double getCx(int nClus) {return m_cx[nClus];}
     double getCy(int nClus) {return m_cy[nClus];}
     double getCz(int nClus) {return m_cz[nClus];}
     double getCt(int nClus) {return m_ct[nClus];}
 
     bool generate_primary_ele();
     int generate_secondary_ele();
     void compute_pos(double trkPosIn[], double trkPosOut[], 
              double &x, double &y, double &z);
     void from_glo_to_loc(double xg, double yg, double zg,
              double &xl, double &yl, double &zl);
 
 
private:
     void clear();
 
     int m_random;
     ICgemGeomSvc* m_geomSvc;
     bool m_debugging;
     double m_magConfig;
 
     int m_nIonC;
     std::vector<double> m_cx;
     std::vector<double> m_cy;
     std::vector<double> m_cz;
     std::vector<double> m_ct;
 
     int m_nIonE;
     std::vector<double> m_ex;
     std::vector<double> m_ey;
     std::vector<double> m_ez;
     std::vector<double> m_et;
 
     double m_track_length_limit;
     double m_track_length;
     double m_n_ion_mm;
 
     TH1F * h_distance_cluster;
     TH1F *h_nof_cluster;
     TFile *output;
};
 
#endif