IonizationPoisson.h

Go to the documentation of this file.

#ifndef IONIZATIONPOISSON_H
#define IONIZATIONPOISSON_H
 
#include "CgemGeomSvc/ICgemGeomSvc.h"
 
#include "CgemDigitizerSvc/Ionization.h"
 
#include <vector>
 
class IonizationPoisson: public Ionization{
public:
     IonizationPoisson();
     ~IonizationPoisson();
 
     void init(unsigned int random, ICgemGeomSvc* geomSvc);
     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;
 
     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;
     int m_n_ion_cm;
};
 
#endif