TouschekReaderTURTLE Class Reference

Class to read Touschek files and store their content in a MCParticle graph. More...

#include <TouschekReaderTURTLE.h>

Public Member Functions
	BELLE2_DEFINE_EXCEPTION (TouschekCouldNotOpenFileError, "Could not open file %1% !")
	Exception is thrown if the Touschek file could not be opened.
	BELLE2_DEFINE_EXCEPTION (TouschekConvertFieldError, "Line %1%: Could not convert field %2%: %3%")
	Exception is thrown if a field in the Touschek file could not be converted to a number.
	TouschekReaderTURTLE (TGeoHMatrix *transMatrix, int pdg)
	Constructor of the TouschekReader class.
	~TouschekReaderTURTLE ()
	Destructor.
void	open (const std::string &filename)
	Opens an ascii file and prepares it for reading.
int	getParticles (int number, MCParticleGraph &graph)
	Reads the specified number of particles from the file and stores the result in the given MCParticle graph.

Protected Attributes
TGeoHMatrix *	m_transMatrix
	Transformation matrix from local Touschek to global geant4 space.
int	m_pdg
	The pdg value of the type of particle that is read (e.g.
std::ifstream	m_input
	The input stream of the ascii file.
int	m_lineNum
	The line number in the ascii file of the last particle which was read.

Detailed Description

Class to read Touschek files and store their content in a MCParticle graph.

The reader supports retrieving the Touschek information from an ascii text file.

The file contains just a list of particles, defined at plane +-50cm from the IP. The reader reads either a specified number of particles from the file, or all particles.

Definition at line 34 of file TouschekReaderTURTLE.h.

Constructor & Destructor Documentation

TouschekReaderTURTLE()

TouschekReaderTURTLE	(	TGeoHMatrix *	transMatrix,
		int	pdg )

Constructor of the TouschekReader class.

Parameters

transMatrix	Pointer to the matrix which transforms the particles from the local Touschek to the global geant4 coordinate system.
pdg	The pdg value of the type of particle that is read (e.g. 11 for e-, -11 for e+).

Definition at line 23 of file TouschekReaderTURTLE.cc.

                                                                           : m_transMatrix(transMatrix), m_pdg(pdg),
  m_lineNum(0)
{
}

~TouschekReaderTURTLE()

~TouschekReaderTURTLE

(

)

Destructor.

Definition at line 29 of file TouschekReaderTURTLE.cc.

{
  if (m_input) m_input.close();
}

Member Function Documentation

getParticles()

int getParticles	(	int	number,
		MCParticleGraph &	graph )

Reads the specified number of particles from the file and stores the result in the given MCParticle graph.

The number of the last particle which was read is saved. So the next time this method is called the reading continues at the position where the last particle which was read.

Parameters

number	The number of particles that should be read from the file. Set it to -1 to read all particles.
graph	Reference to the graph which should be filled with the information from the Touschek file.

Returns

The number of particles which were read.

Definition at line 43 of file TouschekReaderTURTLE.cc.

{
  if (m_transMatrix == NULL) {
    B2ERROR("The transformation matrix is NULL !");
  }
 
  int numParticles = 0;
  string currLine;
  typedef boost::tokenizer<boost::char_separator<char> > tokenizer;
  boost::char_separator<char> sep(",; \t");
  double fields[7];
  double particlePosTouschek[3] = {0.0, 0.0, 0.0};
  double particlePosGeant4[3] = {0.0, 0.0, 0.0};
  double particleMomTouschek[3] = {0.0, 0.0, 0.0};
  double particleMomGeant4[3] = {0.0, 0.0, 0.0};
 
  do {
    getline(m_input, currLine);
    m_lineNum++;
 
    boost::trim(currLine);
    if (currLine == "") continue;
 
    tokenizer tokens(currLine, sep);
 
    int index = 0;
    for (const string& tok : tokens) {
      try {
        if (index >= 7) {
          B2WARNING("This Touschek file has more than 7 fields ! Only the first 7 fields were read.");
          break;
        }
        fields[index] = boost::lexical_cast<double>(tok);
        index++;
      } catch (boost::bad_lexical_cast&) {
        throw (TouschekConvertFieldError() << m_lineNum << index << tok);
      }
    }
 
    if (index < 7) {
      B2ERROR("TouschekReaderTURTLE: incomplete line at " << m_lineNum);
      continue;
    }
 
    //Convert the position of the particle from local Touschek plane space to global geant4 space.
    //Flip the sign for the y and z component to go from the accelerator to the detector coordinate system
    particlePosTouschek[0] = fields[2] * Unit::m;
    particlePosTouschek[1] = -fields[3] * Unit::m;
    if (m_pdg == -11) particlePosTouschek[2] = 400.0; //forLER
    else if (m_pdg == 11) particlePosTouschek[2] = -400.0; //forHER
    else { B2WARNING("m_pdg is not 11/-11"); break;}
 
    m_transMatrix->LocalToMaster(particlePosTouschek, particlePosGeant4);
 
    //Convert the momentum of the particle from local Touschek plane space to global geant4 space.
    //Flip the sign for the y and z component to go from the accelerator to the detector coordinate system
    particleMomTouschek[0] = fields[4];
    particleMomTouschek[1] = -fields[5];
    if (m_pdg == -11) particleMomTouschek[2] = -fields[6]; //forLER
    else if (m_pdg == 11) particleMomTouschek[2] = fields[6]; //forHER
    else { B2WARNING("m_pdg is not 11/-11"); break;}
 
    m_transMatrix->LocalToMasterVect(particleMomTouschek, particleMomGeant4);
 
    double totalMom2 = particleMomGeant4[0] * particleMomGeant4[0];
    totalMom2 += particleMomGeant4[1] * particleMomGeant4[1];
    totalMom2 += particleMomGeant4[2] * particleMomGeant4[2];
 
    //Add particles to MCParticle collection
    MCParticleGraph::GraphParticle& particle = graph.addParticle();
    particle.setStatus(MCParticle::c_PrimaryParticle);
    particle.setPDG(m_pdg);
    particle.setMassFromPDG();
    particle.setMomentum(ROOT::Math::XYZVector(particleMomGeant4[0], particleMomGeant4[1], particleMomGeant4[2]));
    particle.setProductionVertex(ROOT::Math::XYZVector(particlePosGeant4[0], particlePosGeant4[1], particlePosGeant4[2]));
    particle.setEnergy(sqrt(totalMom2 + particle.getMass()*particle.getMass()));
    particle.setProductionTime(0.0);
    particle.setValidVertex(true);
 
    numParticles++;
  } while ((!m_input.eof()) && ((number > -1) && (numParticles < number)));
 
  return numParticles;
}

open()

void open

(

const std::string &

filename

)

Opens an ascii file and prepares it for reading.

Parameters

filename

The filename of the Touschek ascii file which should be read.

Definition at line 35 of file TouschekReaderTURTLE.cc.

{
  m_lineNum = 0;
  m_input.open(filename.c_str());
  if (!m_input) throw(TouschekCouldNotOpenFileError() << filename);
}

Member Data Documentation

m_input

std::ifstream m_input

protected

The input stream of the ascii file.

Definition at line 78 of file TouschekReaderTURTLE.h.

m_lineNum

int m_lineNum

protected

The line number in the ascii file of the last particle which was read.

Definition at line 79 of file TouschekReaderTURTLE.h.

m_pdg

int m_pdg

protected

The pdg value of the type of particle that is read (e.g.

11 for e-, -11 for e+).

Definition at line 77 of file TouschekReaderTURTLE.h.

m_transMatrix

TGeoHMatrix* m_transMatrix

protected

Transformation matrix from local Touschek to global geant4 space.

Definition at line 76 of file TouschekReaderTURTLE.h.

---

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

• generators/touschek/include/TouschekReaderTURTLE.h
• generators/touschek/src/TouschekReaderTURTLE.cc