DiffusionGTS Class Reference

#include <DiffusionGTS.h>

Public Member Functions
	DiffusionGTS (ICgemGeomSvc *geomSvc)
	~DiffusionGTS ()
void	readGasPerpParameters (std::string fileName)
void	readGasParalParameters (std::string fileName)
double	shift_x_drift (double z)
double	shift_x_transf ()
double	shift_x_induct ()
double	sigma_x_drift (double z)
double	sigma_x_transf ()
double	sigma_x_induct ()
double	shift_y_drift (double z)
double	shift_y_transf ()
double	shift_y_induct ()
double	sigma_y_drift (double z)
double	sigma_y_transf ()
double	sigma_y_induct ()
double	shift_t_drift (double z)
double	shift_t_transf ()
double	shift_t_induct ()
double	sigma_t_drift (double z)
double	sigma_t_transf ()
double	sigma_t_induct ()
void	diffusion_drift (double xi, double yi, double zi, double ti, double &xf, double &yf, double &zf, double &tf)
void	diffusion_transf (double xi, double yi, double zi, double ti, double &xf, double &yf, double &zf, double &tf)
void	diffusion_induct (double xi, double yi, double zi, double ti, double &xf, double &yf, double &zf, double &tf)

Detailed Description

Definition at line 7 of file DiffusionGTS.h.

Constructor & Destructor Documentation

DiffusionGTS()

DiffusionGTS::DiffusionGTS

(

ICgemGeomSvc *

geomSvc

)

Definition at line 14 of file DiffusionGTS.cxx.

                                                {
  m_geomSvc = geomSvc;
}

~DiffusionGTS()

DiffusionGTS::~DiffusionGTS

(

)

Definition at line 17 of file DiffusionGTS.cxx.

{}

Member Function Documentation

diffusion_drift()

void DiffusionGTS::diffusion_drift	(	double	xi,
		double	yi,
		double	zi,
		double	ti,
		double &	xf,
		double &	yf,
		double &	zf,
		double &	tf
	)

Definition at line 302 of file DiffusionGTS.cxx.

                                                                                                                             {
 
  // x
  double dx = shift_x_drift(zi);
  double sx = sigma_x_drift(zi);
  xf = xi + gRandom->Gaus(dx, sx);
 
  // y
  double dy = shift_y_drift(zi);
  double sy = sigma_y_drift(zi);
  yf = yi + gRandom->Gaus(dy, sy);
 
  // t
  double dt = shift_t_drift(zi);
  double st = sigma_t_drift(zi);
  tf = ti + gRandom->Gaus(dt, st);
}

diffusion_induct()

void DiffusionGTS::diffusion_induct	(	double	xi,
		double	yi,
		double	zi,
		double	ti,
		double &	xf,
		double &	yf,
		double &	zf,
		double &	tf
	)

Definition at line 339 of file DiffusionGTS.cxx.

                                                                                                                              {
 
  // x
  double dx = shift_x_induct();
  double sx = sigma_x_induct();
  xf = xi + gRandom->Gaus(dx, sx);
 
  // y
  double dy = shift_y_induct();
  double sy = sigma_y_induct();
  yf = yi + gRandom->Gaus(dy, sy);
 
  // t                      
  double dt = shift_t_induct();
  double st = sigma_t_induct();
  tf = ti + gRandom->Gaus(dt, st);
 
}

diffusion_transf()

void DiffusionGTS::diffusion_transf	(	double	xi,
		double	yi,
		double	zi,
		double	ti,
		double &	xf,
		double &	yf,
		double &	zf,
		double &	tf
	)

Definition at line 320 of file DiffusionGTS.cxx.

                                                                                                                              {
 
  // x
  double dx = shift_x_transf();
  double sx = sigma_x_transf();
  xf = xi + gRandom->Gaus(dx, sx);
 
  // y
  double dy = shift_y_transf();
  double sy = sigma_y_transf();
  yf = yi + gRandom->Gaus(dy, sy);
 
  // t
  double dt = shift_t_transf();
  double st = sigma_t_transf();
  tf = ti + gRandom->Gaus(dt, st);
}

readGasParalParameters()

void DiffusionGTS::readGasParalParameters

(

std::string

fileName

)

Definition at line 112 of file DiffusionGTS.cxx.

                                                           {
 
  ifstream fin(fileName.c_str(), ios::in);
  
  string strline;
  string strpar;
  if( ! fin.is_open() ){
    std::cout << "DiffusionGTS::readGasParalParameters ERROR: can not open file " << fileName.c_str() << endl;
  }
  while(fin >> strpar){
    if(strpar[0] == '#'){
      getline(fin, strline);
    }
    else if(strpar == "n_shift_x_drift"){
      fin >> n_shift_y_drift;
    }
    else if(strpar == "shift_x_drift_par"){
      for(int i=0; i<n_shift_y_drift; i++) fin >> m_shift_y_drift_par[i]; 
    }
    else if(strpar == "n_sigma_x_drift"){
      fin >> n_sigma_y_drift;
    }
    else if(strpar == "sigma_x_drift_par"){
      for(int i=0; i<n_sigma_y_drift; i++) fin >> m_sigma_y_drift_par[i];
    }
    else if(strpar == "shift_x_transf_par"){
      fin >> m_shift_y_transf_par;
    }
    else if(strpar == "sigma_x_transf_par"){
      fin >> m_sigma_y_transf_par;
    }
    else if(strpar == "shift_x_induct_par"){
      fin >> m_shift_y_induct_par;
    }
    else if(strpar == "sigma_x_induct_par"){
      fin >> m_sigma_y_induct_par;
    }
    else break;
  }
 
 
  std::cout << "DiffusionGTS::readGasParalParameters: open file " << fileName.c_str() << endl;
 
  cout << "shift y in drift";  
  for(int i=0; i<n_shift_y_drift; i++) cout << " " << m_shift_y_drift_par[i];
  cout << endl;
  cout << "sigma y in drift";
  for(int i=0; i<n_sigma_y_drift; i++) cout << " " << m_sigma_y_drift_par[i];
  cout << endl;
  cout << "shift y in transf " <<  m_shift_y_transf_par << endl;
  cout << "sigma y in transf " << m_sigma_y_transf_par << endl;
  cout << "shift y in induct " << m_shift_y_induct_par << endl;
  cout << "sigma y in induct " << m_sigma_y_induct_par << endl;
 
  cout << "shift t in drift";
  // cout << n_shift_t_drift << endl;
  for(int i=0; i<n_shift_t_drift; i++) cout << " " << m_shift_t_drift_par[i];
  cout << endl;
  
}

Referenced by InductionGTS::init(), and SamplingGTS::readGasParameters().

readGasPerpParameters()

void DiffusionGTS::readGasPerpParameters

(

std::string

fileName

)

Definition at line 19 of file DiffusionGTS.cxx.

                                                          {
 
  ifstream fin(fileName.c_str(), ios::in);
  
  string strline;
  string strpar;
  if( ! fin.is_open() ){
    std::cout << "DiffusionGTS::readGasPerpParameters ERROR: can not open file " << fileName.c_str() << endl;
  }
  while(fin >> strpar){
    if(strpar[0] == '#'){
      getline(fin, strline);
    }
    else if(strpar == "n_shift_x_drift"){
      fin >> n_shift_x_drift;
    }
    else if(strpar == "shift_x_drift_par"){
      for(int i=0; i<n_shift_x_drift; i++) fin >> m_shift_x_drift_par[i]; 
    }
    else if(strpar == "n_sigma_x_drift"){
      fin >> n_sigma_x_drift;
    }
    else if(strpar == "sigma_x_drift_par"){
      for(int i=0; i<n_sigma_x_drift; i++) fin >> m_sigma_x_drift_par[i];
    }
    else if(strpar == "shift_x_transf_par"){
      fin >> m_shift_x_transf_par;
    }
    else if(strpar == "sigma_x_transf_par"){
      fin >> m_sigma_x_transf_par;
    }
    else if(strpar == "shift_x_induct_par"){
      fin >> m_shift_x_induct_par;
    }
    else if(strpar == "sigma_x_induct_par"){
      fin >> m_sigma_x_induct_par;
    }
    else if(strpar == "n_shift_t_drift"){
      fin >> n_shift_t_drift;
    }
    else if(strpar == "shift_t_drift_par"){
      for(int i=0; i<n_shift_t_drift; i++) fin >> m_shift_t_drift_par[i];
    }
    else if(strpar == "n_sigma_t_drift"){
      fin >> n_sigma_t_drift;
    }
    else if(strpar == "sigma_t_drift_par"){
      for(int i=0; i<n_sigma_t_drift; i++) fin >> m_sigma_t_drift_par[i];
    }
    else if(strpar == "shift_t_transf_par"){
      fin >> m_shift_t_transf_par;
    }
    else if(strpar == "sigma_t_transf_par"){
      fin >> m_sigma_t_transf_par;
    }
    else if(strpar == "shift_t_induct_par"){
      fin >> m_shift_t_induct_par;
    }
    else if(strpar == "sigma_t_induct_par"){
      fin >> m_sigma_t_induct_par;
    }
  }
 
  std::cout << "DiffusionGTS::readGasPerpParameters: open file " << fileName.c_str() << endl;
 
  //  cout << n_shift_x_drift  << endl;
  cout << "shift x in drift";  
  for(int i=0; i<n_shift_x_drift; i++) cout << " " << m_shift_x_drift_par[i];
  cout << endl;
  // cout << n_sigma_x_drift << endl;
  cout << "sigma x in drift";
  for(int i=0; i<n_sigma_x_drift; i++) cout << " " << m_sigma_x_drift_par[i];
  cout << endl;
  cout << "shift x in transf " << m_shift_x_transf_par << endl;
  cout << "sigma x in transf " << m_sigma_x_transf_par << endl;
  cout << "shift x in induct " << m_shift_x_induct_par << endl;
  cout << "sigma x in induct " << m_sigma_x_induct_par << endl;
 
  cout << "shift t in drift";
  // cout << n_shift_t_drift << endl;
  for(int i=0; i<n_shift_t_drift; i++) cout << " " << m_shift_t_drift_par[i];
  cout << endl;
  //  cout << n_sigma_t_drift << endl;
  cout << "sigma t in drift";
  for(int i=0; i<n_sigma_t_drift; i++) cout << " " << m_sigma_t_drift_par[i];
  cout << endl;
  cout << "shift t in transf " << m_shift_t_transf_par << endl;
  cout << "sigma t in transf " << m_sigma_t_transf_par << endl;
  cout << "shift t in induct " << m_shift_t_induct_par << endl;
  cout << "sigma t in induct " << m_sigma_t_induct_par << endl;
  
}

Referenced by InductionGTS::init(), and SamplingGTS::readGasParameters().

shift_t_drift()

double DiffusionGTS::shift_t_drift

(

double

z

)

if(z<0||z>5.) { // CHECK 16) gap from geometry! cout<<"DiffusionGTS::shift_t_drift: error in z = " << z << endl; return dt; }

Definition at line 258 of file DiffusionGTS.cxx.

                                          {
  double dt = 0;
  /**
     if(z<0||z>5.) {  // CHECK 16) gap from geometry!
     cout<<"DiffusionGTS::shift_t_drift: error in z = " << z << endl;
     return dt;
     }
  **/
  for(int i=0; i<n_shift_t_drift; i++) dt += pow(z, i) * m_shift_t_drift_par[i];
  return dt;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_drift().

shift_t_induct()

double DiffusionGTS::shift_t_induct

(

)

Definition at line 274 of file DiffusionGTS.cxx.

                                    {
  double dx = m_shift_t_induct_par;
  return dx;
}

Referenced by diffusion_induct(), and InductionGTS::drive_to_anode().

shift_t_transf()

double DiffusionGTS::shift_t_transf

(

)

Definition at line 270 of file DiffusionGTS.cxx.

                                    {
  double dx = m_shift_t_transf_par;
  return dx;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_transf().

shift_x_drift()

double DiffusionGTS::shift_x_drift

(

double

z

)

if(z<0||z>5.) { // CHECK 16) gap from geometry! cout<<"DiffusionGTS::shift_x_drift: error in z = " << z << endl; return dx; }

Definition at line 173 of file DiffusionGTS.cxx.

                                          {
 
  double dx = 0;  
  /**
     if(z<0||z>5.) {  // CHECK 16) gap from geometry! 
     cout<<"DiffusionGTS::shift_x_drift: error in z = " << z << endl;
     return dx;
     } 
  **/
  for(int i=0; i<n_shift_x_drift; i++) dx += pow(z, i) * m_shift_x_drift_par[i];
  return dx;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_drift().

shift_x_induct()

double DiffusionGTS::shift_x_induct

(

)

Definition at line 189 of file DiffusionGTS.cxx.

                                    {
  double dx = m_shift_x_induct_par;
  return dx;
}

Referenced by diffusion_induct(), and InductionGTS::drive_to_anode().

shift_x_transf()

double DiffusionGTS::shift_x_transf

(

)

Definition at line 185 of file DiffusionGTS.cxx.

                                    {
  double dx = m_shift_x_transf_par;
  return dx;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_transf().

shift_y_drift()

double DiffusionGTS::shift_y_drift

(

double

z

)

if(z<0||z>5.) { // CHECK 16) gap from geometry! cout<<"DiffusionGTS::shift_y_drift: error in z = " << z << endl; return dy; }

Definition at line 216 of file DiffusionGTS.cxx.

                                          {
 
  double dy = 0;  
  /**
     if(z<0||z>5.) {  // CHECK 16) gap from geometry! 
     cout<<"DiffusionGTS::shift_y_drift: error in z = " << z << endl;
     return dy;
     } 
  **/
  for(int i=0; i<n_shift_y_drift; i++) dy += pow(z, i) * m_shift_y_drift_par[i];
  return dy;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_drift().

shift_y_induct()

double DiffusionGTS::shift_y_induct

(

)

Definition at line 232 of file DiffusionGTS.cxx.

                                    {
  double dy = m_shift_y_induct_par;
  return dy;
}

Referenced by diffusion_induct(), and InductionGTS::drive_to_anode().

shift_y_transf()

double DiffusionGTS::shift_y_transf

(

)

Definition at line 228 of file DiffusionGTS.cxx.

                                    {
  double dy = m_shift_y_transf_par;
  return dy;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_transf().

sigma_t_drift()

double DiffusionGTS::sigma_t_drift

(

double

z

)

if(z<0||z>5.) { // CHECK 16) gap from geometry! cout<<"DiffusionGTS::sigma_t_drift: error in z = " << z << endl; return dt; }

Definition at line 279 of file DiffusionGTS.cxx.

                                          {
  double dt = 0;
  /**
     if(z<0||z>5.) {  // CHECK 16) gap from geometry! 
     cout<<"DiffusionGTS::sigma_t_drift: error in z = " << z << endl;
     return dt;
     }
  **/
  for(int i=0; i<n_sigma_t_drift; i++) dt += pow(z, i) * m_sigma_t_drift_par[i];
  return dt;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_drift().

sigma_t_induct()

double DiffusionGTS::sigma_t_induct

(

)

Definition at line 294 of file DiffusionGTS.cxx.

                                    {
  double dx = m_sigma_t_induct_par;
  return dx;
}

Referenced by diffusion_induct(), and InductionGTS::drive_to_anode().

sigma_t_transf()

double DiffusionGTS::sigma_t_transf

(

)

Definition at line 290 of file DiffusionGTS.cxx.

                                    {
  double dx = m_sigma_t_transf_par;
  return dx;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_transf().

sigma_x_drift()

double DiffusionGTS::sigma_x_drift

(

double

z

)

if(z<0||z>5.) { // CHECK 16) gap from geometry! cout<<"DiffusionGTS::sigma_x_drift: error in z = " << z << endl; return dx; }

Definition at line 195 of file DiffusionGTS.cxx.

                                          {
  double dx = 0;  
  /**
     if(z<0||z>5.) {  // CHECK 16) gap from geometry! 
     cout<<"DiffusionGTS::sigma_x_drift: error in z = " << z << endl;
     return dx;
     } 
  **/
  for(int i=0; i<n_sigma_x_drift; i++) dx += pow(z, i) * m_sigma_x_drift_par[i];
  return dx;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_drift().

sigma_x_induct()

double DiffusionGTS::sigma_x_induct

(

)

Definition at line 210 of file DiffusionGTS.cxx.

                                    {
  double dx = m_sigma_x_induct_par;
  return dx;
}

Referenced by diffusion_induct(), and InductionGTS::drive_to_anode().

sigma_x_transf()

double DiffusionGTS::sigma_x_transf

(

)

Definition at line 206 of file DiffusionGTS.cxx.

                                    {
  double dx = m_sigma_x_transf_par;
  return dx;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_transf().

sigma_y_drift()

double DiffusionGTS::sigma_y_drift

(

double

z

)

if(z<0||z>5.) { // CHECK 16) gap from geometry! cout<<"DiffusionGTS::sigma_y_drift: error in z = " << z << endl; return dy; }

Definition at line 237 of file DiffusionGTS.cxx.

                                          {
  double dy = 0;  
  /**
     if(z<0||z>5.) {  // CHECK 16) gap from geometry! 
     cout<<"DiffusionGTS::sigma_y_drift: error in z = " << z << endl;
     return dy;
     } 
  **/
  for(int i=0; i<n_sigma_y_drift; i++) dy += pow(z, i) * m_sigma_y_drift_par[i];
  return dy;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_drift().

sigma_y_induct()

double DiffusionGTS::sigma_y_induct

(

)

Definition at line 252 of file DiffusionGTS.cxx.

                                    {
  double dy = m_sigma_y_induct_par;
  return dy;
}

Referenced by diffusion_induct(), and InductionGTS::drive_to_anode().

sigma_y_transf()

double DiffusionGTS::sigma_y_transf

(

)

Definition at line 248 of file DiffusionGTS.cxx.

                                    {
  double dy = m_sigma_y_transf_par;
  return dy;
}

Referenced by SamplingGTS::compute_diffusion_on_GEM3(), and diffusion_transf().

---

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

• DiffusionGTS.h
• DiffusionGTS.cxx