InductionGTS.h

Go to the documentation of this file.

#ifndef INDUCTIONGTS_H
#define INDUCTIONGTS_H
 
#include "CgemGeomSvc/ICgemGeomSvc.h"
#include "CgemDigitizerSvc/Induction.h"
#include "CgemDigitizerSvc/DiffusionGTS.h"
 
#include <vector>
#include <string>
#include "TNtuple.h"
#include "TFile.h"
#include "TF1.h"
#include "TH1F.h"
 
#define m_testing_ind false
#define hnbin 700
#define hxmin 0
#define hxmax 700
#define apv_tstart 0.
 
class G4Svc;
 
class InductionGTS : public Induction {
public:
     InductionGTS();
     ~InductionGTS();
 
     void init(ICgemGeomSvc* geomSvc, double magConfig);
     void setDebugOutput(bool debugOutput){m_debugOutput = debugOutput;}
     void setMultiElectrons(int layer, int nElectrons, const std::vector<Float_t>&  x, const std::vector<Float_t>& y, const std::vector<Float_t> &z, const std::vector<Float_t> &t);
 
     /* output info of fired strips */
     int    getNXstrips() const;
     int    getNVstrips() const;
     int    getXstripSheet(int n) const;
     int    getXstripID(int n) const;
     int    getVstripSheet(int n) const;
     int    getVstripID(int n) const;
     double getXstripQ(int n) const;
     double getVstripQ(int n) const;
     double getXstripT(int n) const;
     double getVstripT(int n) const;
 
     void clear();
     bool drive_to_anode(int, double, double, double, double, double&, double&, double&, double&);
     bool useAPV(int stripid, int view, std::vector< double > stripidvec, 
         std::vector< double > indchargevec, std::vector< double > indtimevec, 
         double &charge, double &time, double &dtime);
 
     void   setVsampleDelay(double delay) {return;}
     void   setStoreFlag(bool flag) {return;}
     void   setLUTFilePath(std::string path) {return;}
     void   setSaturation(bool flag) {return;}
     double getXstripT_Branch(int n) const {return m_XstripT.at(n);}
     double getVstripT_Branch(int n) const {return m_VstripT.at(n);}
     double getXstripQ_Branch(int n) const {return m_XstripQ.at(n);}
     double getVstripQ_Branch(int n) const {return m_VstripQ.at(n);}
     double getXfirstT(int n) const {return m_XstripT.at(n);}
     double getVfirstT(int n) const {return m_VstripT.at(n);}
 
private:
 
     int m_field;
     double m_tuning_factor_diff_perp;
     double m_tuning_factor_diff_paral;
 
     ICgemGeomSvc* m_geomSvc;
     double m_magConfig;        /* magnetic field */
     DiffusionGTS *diffusion;
     bool m_debugOutput;
 
 
     TFile *output;
     TH1F *hcollected_charge;
     TH1F *hintegratore;
     TH1F *hcharge;
     TF1 *f[hnbin];
     TF1 *f_FD;
     TNtuple *tree;
     TNtuple *tree_strip;
 
     // temp
     std::vector< int > stripid_x;
     std::vector< int > sheetid_x;
     std::vector<double> charge_x;
     std::vector<double> time_x;
     std::vector< int > stripid_v;
     std::vector< int > sheetid_v;
     std::vector<double> charge_v;
     std::vector<double> time_v;
 
 
     // final
     int m_NXstrips;
     int m_NVstrips;
 
     std::vector<double> m_XstripSheet; 
     std::vector<double> m_XstripID; 
     std::vector<double> m_VstripSheet; 
     std::vector<double> m_VstripID; 
     std::vector<double> m_XstripQ;
     std::vector<double> m_VstripQ;
     std::vector<double> m_XstripT;
     std::vector<double> m_VstripT;             
 
     //
     int evt;
};
 
#endif