DevelopSignal.h

Go to the documentation of this file.

#ifndef DEVELOPSIGNAL_H
#define DEVELOPSIGNAL_H
 
#include "CgemGeomSvc/ICgemGeomSvc.h"
#include "CgemDigitizerSvc/DriftAndAvalanche.h"
 
#include "TF1.h"
 
#include <vector>
#include <string>
 
class G4Svc;
 
class DevelopSignal: public DriftAndAvalanche {
public:
     DevelopSignal();
     ~DevelopSignal();
 
     void init(ICgemGeomSvc* geomSvc);
 
 
     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;}
     double getX(int n) const {return m_eX[n];}
     double getY(int n) const {return m_eY[n];}
     double getZ(int n) const {return m_eZ[n];}
     double getT(int n) const {return m_eT[n];}
 
     // --------------------------------
     // CHECK 11) per ora li metto hard coded
     bool readGemParameters();
     bool readGasParameters();
     // void readFldParameters();
     
     int get_voltage_index(int voltage);
     bool is_survived();
     double shift_x(double z);
     double sigma_x(double z);
 
     double sigma_y(double z);
 
     double shift_t(double z);
     double sigma_t(double z);
 
     void set_eff_polya();
     double compute_gain();
     void compute_drifting(double xi, double yi, double zi, double ti, double &xf, double &yf, double &zf, double &tf);
 
private:
     void clear();
 
     //   GasProp *gas; CHECK 10) maybe obj
     string m_gas; // which gas
 
     // for each GEM CHECK 10) --> GemProp
     int m_hv; // high voltage
     double m_eff_col[3]; // efficiency of collection
     double m_eff_ext[3]; // efficiency of extraction
 
 
     bool m_field; // CHECK 9) might need a full class for the field
 
     TF1 *m_polya;   // CHECK 15)
 
     ICgemGeomSvc* m_geomSvc;
 
     // CHECK set variable number of parameters
     double m_shift_x_parameters[2];
     double m_shift_t_parameters[2];
 
     double m_smear_x_parameters[4];
     double m_smear_y_parameters[4];
     double m_smear_t_parameters[4];
 
     /* info of multiplied electrons */
     int m_nMulElec;
     std::vector<double> m_eX;
     std::vector<double> m_eY;
     std::vector<double> m_eZ;
     std::vector<double> m_eT;
 
     double m_polya_mean[3];
     double m_polya_theta[3];
     
 
 
 
 
};
 
#endif