IonizationGTS Class Reference

#include <IonizationGTS.h>

Inheritance diagram for IonizationGTS:

Public Member Functions
	IonizationGTS ()
	~IonizationGTS ()
void	init (unsigned int random, ICgemGeomSvc *geomSvc, double magConfig)
void	setDebugging (bool debugging)
void	setTrack (int particle, int charge, double p, double trkPosIn[], double trkPosOut[])
int	getNumberIonE ()
double	getEx (int nElec)
double	getEy (int nElec)
double	getEz (int nElec)
double	getEt (int nElec)
int	getNumberIonC ()
double	getCx (int nClus)
double	getCy (int nClus)
double	getCz (int nClus)
double	getCt (int 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)
Public Member Functions inherited from Ionization
	Ionization ()
virtual	~Ionization ()
virtual void	init (unsigned int random, ICgemGeomSvc *geomSvc, double magConfig)=0
virtual void	setDebugging (bool debugging)=0
virtual void	setTrack (int particle, int charge, double p, double trkPosIn[], double trkPosOut[])=0
virtual int	getNumberIonE ()=0
virtual double	getEx (int nElec)=0
virtual double	getEy (int nElec)=0
virtual double	getEz (int nElec)=0
virtual double	getEt (int nElec)=0

Detailed Description

Definition at line 13 of file IonizationGTS.h.

Constructor & Destructor Documentation

IonizationGTS()

IonizationGTS::IonizationGTS

(

)

Definition at line 42 of file IonizationGTS.cxx.

                            {
}

~IonizationGTS()

IonizationGTS::~IonizationGTS

(

)

Definition at line 45 of file IonizationGTS.cxx.

                             {
  if(m_testing == true) {
    output->Write();
    output->Close();
  }
}

Member Function Documentation

compute_pos()

void IonizationGTS::compute_pos	(	double	trkPosIn[],
		double	trkPosOut[],
		double &	x,
		double &	y,
		double &	z
	)

Definition at line 205 of file IonizationGTS.cxx.

                                                                                                {
  
  double t = (m_track_length_limit - m_track_length)/m_track_length_limit;
  x = posin[0] * t + posout[0] * (1 - t);
  y = posin[1] * t + posout[1] * (1 - t); 
  z = posin[2] * t + posout[2] * (1 - t);
 
}

Referenced by setTrack().

from_glo_to_loc()

void IonizationGTS::from_glo_to_loc	(	double	xg,
		double	yg,
		double	zg,
		double &	xl,
		double &	yl,
		double &	zl
	)

generate_primary_ele()

bool IonizationGTS::generate_primary_ele

(

)

Definition at line 169 of file IonizationGTS.cxx.

                                         {
  
  double u = gRandom->Uniform(0, 1); // CHECK 6) TRandom
  // step along the trajectory, that for now is a straight line [x0, z0] -> [x1, z1]
  double dl_extracted = -(1./m_n_ion_mm) * TMath::Log(1 - u);                                                                                
  m_track_length += dl_extracted;
 
  //  cout << "m_track_length " <<   m_track_length << " dl_extracted " << dl_extracted << endl;
 
 
  if(m_track_length > m_track_length_limit) return false;
  
  m_nIonC++;
  m_nIonE++;
  return true;
}

Referenced by setTrack().

generate_secondary_ele()

int IonizationGTS::generate_secondary_ele

(

)

Definition at line 188 of file IonizationGTS.cxx.

                                          {
  int n_elec=0;
  double prob = gRandom->Uniform(0, 1);  // CHECK 6) TRandom
 
  for(int iele = 0; iele < 100; iele++) { // CHECK 4) cut at 100
    n_elec = iele + 1;
    if(prob < int_prob_electron_per_cluster[iele]) break;
  }
 
  if(n_elec>100) { 
    std::cout << "IonizationGTS::generate_secondary_ele: increase the size above 100"<<endl; // CHECK 4)
    n_elec=100;
  }
 
  return n_elec;
}

Referenced by setTrack().

getCt()

double IonizationGTS::getCt ( int  nClus )

inline

Definition at line 39 of file IonizationGTS.h.

{return m_ct[nClus];}

getCx()

double IonizationGTS::getCx ( int  nClus )

inline

Definition at line 36 of file IonizationGTS.h.

{return m_cx[nClus];}

getCy()

double IonizationGTS::getCy ( int  nClus )

inline

Definition at line 37 of file IonizationGTS.h.

{return m_cy[nClus];}

getCz()

double IonizationGTS::getCz ( int  nClus )

inline

Definition at line 38 of file IonizationGTS.h.

{return m_cz[nClus];}

getEt()

double IonizationGTS::getEt ( int  nElec )

inlinevirtual

Implements Ionization.

Definition at line 32 of file IonizationGTS.h.

{return m_et[nElec];}

getEx()

double IonizationGTS::getEx ( int  nElec )

inlinevirtual

Implements Ionization.

Definition at line 29 of file IonizationGTS.h.

{return m_ex[nElec];}

getEy()

double IonizationGTS::getEy ( int  nElec )

inlinevirtual

Implements Ionization.

Definition at line 30 of file IonizationGTS.h.

{return m_ey[nElec];}

getEz()

double IonizationGTS::getEz ( int  nElec )

inlinevirtual

Implements Ionization.

Definition at line 31 of file IonizationGTS.h.

{return m_ez[nElec];}

getNumberIonC()

int IonizationGTS::getNumberIonC ( )

inline

Definition at line 35 of file IonizationGTS.h.

{return m_nIonC;}

getNumberIonE()

int IonizationGTS::getNumberIonE ( )

inlinevirtual

Implements Ionization.

Definition at line 28 of file IonizationGTS.h.

{return m_nIonE;}

init()

void IonizationGTS::init ( unsigned int  random, ICgemGeomSvc *  geomSvc, double  magConfig  )

virtual

Implements Ionization.

Definition at line 52 of file IonizationGTS.cxx.

                                                                                    {
  double time_spent = 0.;
  clock_t begin = clock();
  
  //  std::cout << "IonizationGTS::init" << std::endl;
  m_random = random;
  gRandom->SetSeed(m_random); // CHECK 6) TRandom 
  m_geomSvc = geomSvc;
  m_magConfig = magConfig;
  
  cout << "IonizationGTS::init: m_testing " << m_testing << endl;
  
  
  if(m_testing==true) {
    output = new TFile("ionization.root", "RECREATE");
    h_distance_cluster = new TH1F("h_distance_cluster", "distance between two consecutive clusters", 100, 0, 5);
    h_nof_cluster = new TH1F("h_nof_cluster", "number of clusters", 100, 0, 100);
  }
 
  clock_t end = clock();
  time_spent += (double)(end - begin) / CLOCKS_PER_SEC;
 }

setDebugging()

void IonizationGTS::setDebugging ( bool  debugging )

inlinevirtual

Implements Ionization.

Definition at line 19 of file IonizationGTS.h.

{m_debugging = debugging;}

setTrack()

void IonizationGTS::setTrack ( int  particle, int  charge, double  p, double  trkPosIn[], double  trkPosOut[]  )

virtual

cout << "m_n_ion_mm " << m_n_ion_mm << " N_PRI_AVE " << N_PRI_AVE << " m_track_length_limit " << m_track_length_limit << endl; cout << "xIN " << trkPosIn[0] << " yIN " << trkPosIn[1] << " zIN " << trkPosIn[2] << endl;
cout << "xOUT " << trkPosOut[0] << " yOUT " << trkPosOut[1] << " zOUT " << trkPosOut[2] << endl;

Implements Ionization.

Definition at line 77 of file IonizationGTS.cxx.

                                                                                                      {
   
  double time_spent = 0.;
  clock_t begin = clock();
 
  clear();
  // cout << "IonizationGTS setTrack" << endl;
  
  // CHECK 1) 3D
  m_track_length_limit = TMath::Sqrt((trkPosOut[0] - trkPosIn[0])*(trkPosOut[0] - trkPosIn[0]) + (trkPosOut[1] - trkPosIn[1])*(trkPosOut[1] - trkPosIn[1]) + (trkPosOut[2] - trkPosIn[2])*(trkPosOut[2] - trkPosIn[2]));
  m_track_length = 0;
 
  m_n_ion_mm = N_PRI_AVE/m_track_length_limit; // CHECK 2) N_PRI_AVE
 
  /**
  cout << "m_n_ion_mm " << m_n_ion_mm << " N_PRI_AVE " << N_PRI_AVE << " m_track_length_limit " << m_track_length_limit << endl;
  cout << "xIN  " << trkPosIn[0]  << " yIN  " << trkPosIn[1]  << " zIN  " << trkPosIn[2] << endl;    
  cout << "xOUT " << trkPosOut[0] << " yOUT " << trkPosOut[1] << " zOUT " << trkPosOut[2] << endl;    
  **/
 
  // loop over ionization clusters
  int counter = 0;
 
  double last_x;
  double last_y;
  double last_z;
  double last_t;
 
  while(1) {
    
    // primary ionization
    bool is_inside = generate_primary_ele();
    if(!is_inside) break;
    
    // secondary ionization
    int nof_sec_ele = generate_secondary_ele();
    m_nIonE += nof_sec_ele;
    
    // compute position
    double x = 0;
    double y = 0; 
    double z = 0;
    double t = 0; // <-track_time;  // CHECK 3) add the track time 
    compute_pos(trkPosIn, trkPosOut, x, y, z);
    //    cout << "ionization electron pos " << x << " " << y << " " << z << endl;    
 
    // set cluster pos
    m_cx.push_back(x); // CHECK 5) cluster info
    m_cy.push_back(y); // CHECK 5) cluster info
    m_cz.push_back(z); // CHECK 5) cluster info
    m_ct.push_back(t); // CHECK 5) cluster info
 
    //    cout << "x " << x << " y " << y << " z " << z << " t " << t << endl;
    
    if(m_testing==true && counter > 0) {
      double distance = TMath::Sqrt((last_x - x)*(last_x - x) + (last_y - y)*(last_y - y) + (last_z - z)*(last_z - z));
      h_distance_cluster->Fill(distance);
      last_x = x;
      last_y = y;
      last_z = z;
      last_t = t;
    }
 
    // set ele pos
    for(int iele = 0; iele < nof_sec_ele + 1; iele++)
      {
    m_ex.push_back(x);
    m_ey.push_back(y);
    m_ez.push_back(z);
    m_et.push_back(t);
    
      }
    counter++;
    //    cout << "ion position " << x << " " << y << " " << z << endl;
 
  }
  
  if(m_nIonE != m_ex.size()) std::cout << "IonizationGTS::setTrack ERRRRRRRRRRRRRROR" << std::endl; // CHECK delete this
 
  if(m_testing==true) h_nof_cluster->Fill(counter-1);
 
  clock_t end = clock();
  time_spent += (double)(end - begin) / CLOCKS_PER_SEC;
  cout << "Sampling::ionization " << m_nIonE << " time " << time_spent << " seconds" << endl;
 
 
 
 
}

---

The documentation for this class was generated from the following files:

• IonizationGTS.h
• IonizationGTS.cxx