BesTofDigitizerEcV4.hh

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//      BOOST --- BESIII Object_Oriented Simulation Tool                     //
//---------------------------------------------------------------------------//
//Description:
//Author: Matthias Ullrich
//Created: June, 2011
//Modified:
//Comment: Be carefull the 4 static function will be used within the reconstruction
 
//---------------------------------------------------------------------------//
// $Id: BesTofDigitizerEcV4.hh
 
#ifndef BesTofDigitizerEcV4_h
#define BesTofDigitizerEcV4_h 1
 
#include "BesTofDigitizerV.hh"
#include "BesTofDigi.hh"
#include "globals.hh"
#include "TFile.h"
#include "TTree.h"
#include "ScintSingle.hh"
 
class G4Svc;
 
class BesTofDigitizerEcV4 : public BesTofDigitizerV
{
 
  public:
     BesTofDigitizerEcV4();
    ~BesTofDigitizerEcV4();
    
  virtual void Digitize(ScintSingle*, BesTofDigitsCollection* );
 
 
 
  G4double Smear_gaussian(G4double,G4double,G4double);
  static G4int Calculate_Readoutstrip_number(G4double x,G4double y, G4int partId_f, G4int module_mrpc_f);
  static G4int Calculate_Readoutstrip_number_continuum(G4double x_mm,G4double y_mm, G4int partId_f, G4int module_mrpc_f);
 
  static G4int Produce_unique_identifier(G4int module_mrpc_f, G4int readoutstripnumber_f);
  static G4int Get_module_mrpc_from_unique_identifier(G4int unique_identifier_f);
  static G4int Get_stripnumber_from_unique_identifier(G4int unique_identifier_f);
  static G4double Calculate_strip_transition_time(G4double x_mm,G4double y_mm, G4int partId_f, G4int module_mrpc_f);
  static G4double Calculate_strip_transition_time_cont(G4double x_mm,G4double y_mm, G4int partId_f, G4int module_mrpc_f);
  static G4double Average_transition_time(G4int my_module);
 
  
  G4int Get_gapnumber(G4int partID, G4double z_position);
  G4double Distance_for_avalanche(G4int gap, G4double z_position, G4int partID);
 
  struct adc_info{G4int gap; G4int readoutstrip; G4int ions; G4double zdistance; G4double time; G4int trackindex; G4int pdg_code;G4double transtime;G4double x; G4double y;};
  void simulate_avalanche(G4int actual_strip, vector<adc_info> * avalanche_info, double * avalanche_threshold_time, double *induced_charge);
  double adcsignal_get_sigma(double alpha,double eta, double stepwidth);
  long int adcsignal_get_n_electron(double s, double stepwidth, double alpha, double eta);
  long int adcsignal_simulate_step(int n_elec,double stepwidth,double alpha,double eta, bool saturated);
 
  double chargetotime(double);
  double Get_NINO_pulselength_resolution(double induced_charge_fc);
  double Get_NINO_leadingedge_resolution(double induced_charge_fc);
 
 
  private:
 
  G4double time;   
  G4int module_mrpc; 
  G4int partId;
  G4int trackIndex;
      
  BesTofHitsCollection* m_THC;
 
  G4Svc* m_G4Svc;
 
 
  //In order to store the data into the root-files
  bool m_storedata;
  std::string FileName;
  TFile* m_TreeFile;
  TTree * m_ecv4tree;
  double m_partID;
  double m_stripidentifier;
  double m_trackIndex;
  double m_signal_pc;
  double m_time_threshold;
  double m_time_1sthit;
  double m_time_1;
  double m_time_2;
  double m_firedstrip;
  double m_numberions;
  double m_numberofgaps_with_avalanches;
  double m_numberofstripsfired;
  double m_multihit;
  double m_inefficiency;
  double m_particle_true;
  double m_particle_has_signal;
  double m_trans_time;
  double m_x;
  double m_y;
 
};
 
 
#endif