Cgem/CgemDigitizerSvc/CgemDigitizerSvc-00-00-24/CgemDigitizerSvc/IonizationGar.h

Go to the documentation of this file.

#ifndef IONIZATIONGAR_H
#define IONIZATIONGAR_H
 
#include "ComponentAnsys123.hh"
#include "ComponentAnalyticField.hh"
#include "MediumMagboltz.hh"
#include "GeometrySimple.hh"
#include "SolidBox.hh"
#include "Sensor.hh"
#include "TrackHeed.hh"
#include "RandomEngine.hh"
#include "Random.hh"
 
#include "CgemGeomSvc/ICgemGeomSvc.h"
 
#include "CgemDigitizerSvc/Ionization.h"
 
#include <vector>
#include <string>
 
class IonizationGar: public Ionization{
public:
     IonizationGar();
     ~IonizationGar();
 
     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];}
 
private:
     void clear();
     void initGeoGas();
     std::string  getParticle(int particle, int charge) const; /* charge = 1 or -1 */
 
     unsigned int m_random;
     ICgemGeomSvc* m_geomSvc;
     bool m_debugging;
     double m_magConfig;
 
     /* Garfield::ComponentAnsys123* m_gem; */
     Garfield::ComponentAnalyticField* m_comp;
     Garfield::MediumMagboltz* m_gas;
     Garfield::GeometrySimple* m_geo;
     Garfield::SolidBox* m_box;
     Garfield::Sensor* m_sensor;
     /* Garfield::TrackHeed* m_track; */
 
     int m_nIonE;
     std::vector<double> m_ex;
     std::vector<double> m_ey;
     std::vector<double> m_ez;
     std::vector<double> m_et;
};
 
#endif