IonizationGar Class Reference

#include <IonizationGar.h>

Inheritance diagram for IonizationGar:

Public Member Functions
	IonizationGar ()
	~IonizationGar ()
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)
Public Member Functions inherited from Ionization
	Ionization ()
virtual	~Ionization ()

Detailed Description

Definition at line 21 of file IonizationGar.h.

Constructor & Destructor Documentation

IonizationGar()

IonizationGar::IonizationGar

(

)

Definition at line 19 of file IonizationGar.cxx.

                            {
}

~IonizationGar()

IonizationGar::~IonizationGar

(

)

Definition at line 22 of file IonizationGar.cxx.

                             {
     delete m_gas;
     delete m_geo;
     delete m_comp;
     delete m_box;
     delete m_sensor;
}

Member Function Documentation

getEt()

double IonizationGar::getEt ( int nElec )

inlinevirtual

Implements Ionization.

Definition at line 40 of file IonizationGar.h.

{return m_et[nElec];}

getEx()

double IonizationGar::getEx ( int nElec )

inlinevirtual

Implements Ionization.

Definition at line 37 of file IonizationGar.h.

{return m_ex[nElec];}

getEy()

double IonizationGar::getEy ( int nElec )

inlinevirtual

Implements Ionization.

Definition at line 38 of file IonizationGar.h.

{return m_ey[nElec];}

getEz()

double IonizationGar::getEz ( int nElec )

inlinevirtual

Implements Ionization.

Definition at line 39 of file IonizationGar.h.

{return m_ez[nElec];}

getNumberIonE()

int IonizationGar::getNumberIonE ( )

inlinevirtual

Implements Ionization.

Definition at line 36 of file IonizationGar.h.

{return m_nIonE;}

init()

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

virtual

Implements Ionization.

Definition at line 30 of file IonizationGar.cxx.

                                                                                    {
     m_random = random;
     m_geomSvc = geomSvc;
     m_magConfig = magConfig;
 
     initGeoGas();
 
     //TrackHeed* 
     m_track = new TrackHeed();
     m_track->SetSensor(m_sensor);
     if(m_debugging) m_track->EnableDebugging();
     else m_track->DisableDebugging();
     //m_track->EnableMagneticField();
 
}

setDebugging()

void IonizationGar::setDebugging ( bool debugging )

inlinevirtual

Implements Ionization.

Definition at line 27 of file IonizationGar.h.

{m_debugging = debugging;}

setTrack()

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

virtual

Implements Ionization.

Definition at line 46 of file IonizationGar.cxx.

                                                                                                     {
     clear();
     int nE = 0;
 
     double x0 = 0., y0 = 0., z0 = 0., t0 = 0.;
 
     // //the length of track segment.
     // double d0 = sqrt(pow((trkPosOut[0]-trkPosIn[0]),2)+pow((trkPosOut[1]-trkPosIn[1]),2)+pow((trkPosOut[2]-trkPosIn[2]),2));
     // double d1 = sqrt(pow((trkPosOut[0]-trkPosIn[0]),2)+pow((trkPosOut[1]-trkPosIn[1]),2));
     // double trkL = d0;
 
     // double lambda = acos(d1/d0);
     // double theta = acos((trkPosOut[1]-trkPosIn[1])/d1);
       
     // double dx0 = cos(lambda)*sin(theta)/10.0; // mm->cm
     // double dy0 = cos(lambda)*cos(theta)/10.0; // mm->cm
     // double dz0 = sin(lambda)/10.0;         // mm->cm
 
     double dx0 = (trkPosOut[0] - trkPosIn[0])/10.0; // cm
     double dy0 = (trkPosOut[1] - trkPosIn[1])/10.0; // cm
     double dz0 = (trkPosOut[2] - trkPosIn[2])/10.0; // cm
     double trkL = sqrt(dx0*dx0 + dy0*dy0 + dz0*dz0);
 
     // cout << "track segment: " << setw(15) << trkPosIn[0] << setw(15) << trkPosIn[1] << setw(15) << trkPosIn[2]
     //       << setw(15) << trkPosOut[0] << setw(15) << trkPosOut[1] << setw(15) << trkPosOut[2] << endl;
     // cout << setw(15) << trkL << setw(15) << dx0 << setw(15) << dy0 << setw(15) << dz0 << endl;
 
     double minX = 0. < dx0 ? 0. : dx0;
     double maxX = 0. < dx0 ? dx0 : 0.;
     double minY = 0. < dy0 ? 0. : dy0;
     double maxY = 0. < dy0 ? dy0 : 0.;
     double minZ = 0. < dz0 ? 0. : dz0;
     double maxZ = 0. < dz0 ? dz0 : 0.;
     if(maxX<=minX){
      minX += -0.1;
      maxX += 0.1;
     }
     if(maxY<=minY){
      minY += -0.1;
      maxY += 0.1;
     }
     if(maxZ<=minZ){
      minZ += -0.1;
      maxZ += 0.1;
     }
 
     string particleType = getParticle(particle, charge);
 
     m_track->SetParticle(particleType);
     m_track->SetMomentum(p*1.e9);
     // m_sensor->SetArea(-0.5, 0, -0.5, 0.5, 0.5, 0.5);
     m_sensor->SetArea(minX, minY, minZ, maxX, maxY, maxZ);
     // cout << setw(15) << minX << setw(15) << maxX << setw(15) << minY << setw(15) << maxY
     //       << setw(15) << minZ << setw(15) << maxZ << endl;
 
     // cout << "trkpos" << endl;
     // cout << setw(15) << trkPosIn[0] << setw(15) << trkPosIn[1] << setw(15) << trkPosIn[2]
     //       << setw(15) << trkPosOut[0] << setw(15) << trkPosOut[1] << setw(15) << trkPosOut[2] << endl;
 
     //cout<<"new Track"<<endl;
     m_track->NewTrack(x0, y0, z0, t0, dx0, dy0, dz0);
     //cout<<"after new Track"<<endl;
     // m_track->NewTrack(0, 0, 0, t0, 0, 1, 0);
     // cout << "particle: " << particleType << ", P = " << p << endl;
     // cout << setw(15) << x0 << setw(15) << y0 << setw(15) << z0 << setw(15) << t0
     //       << setw(15) << dx0 << setw(15) << dy0 << setw(15) << dz0 << endl;
 
     double xc = 0., yc = 0., zc = 0., tc = 0.; // cluster position
     int nc = 0;              // number of electrons
     double ec = 0.;              // energy
     double extra = 0.;               // not used
     int nCluster = 0;
 
     while(m_track->GetCluster(xc, yc, zc, tc, nc, ec, extra)){
      //   cout<<"nCluster="<<nCluster<<endl;
          double length = sqrt(pow(xc,2)+pow(yc,2)+pow(zc,2));
      // cout << "length = " << length << ",  trkL = " << trkL << endl;
          if(length>trkL){
           cout << "overflow -------------------------" << endl;
           continue;
      }
      // cout << setw(5) << nCluster << setw(15) << xc << setw(15) << yc << setw(15) << zc << setw(15) << nc << endl;
          nCluster++;
          for(int j = 0; j<nc; ++j) {
               double xe, ye, ze, te;
               double ee, dxe, dye, dze;
           //cout <<" get e "<<j<<endl;
               m_track->GetElectron(j, xe, ye, ze, te, ee, dxe, dye, dze); // unit: cm
 
               // cout << setw(5) << nE << setw(15) << xe << setw(15) << ye << setw(15) << ze << setw(15) << te << endl;
               m_ex.push_back(trkPosIn[0] + xe*10.0); // mm
               m_ey.push_back(trkPosIn[1] + ye*10.0); // mm
               m_ez.push_back(trkPosIn[2] + ze*10.0); // mm
               m_et.push_back(te);
               nE++;
          }
     }
     m_nIonE = nE;
     //delete m_track;
}

---

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

• IonizationGar.h
• IonizationGar.cxx