A class to describe a MARS particle read from an input ASCII file. More...

#include <MarsParticle.h>

Collaboration diagram for MarsParticle:

Public Member Functions
	MarsParticle ()
int	read (std::vector< std::string > *eventAsStringVector)
	A function to read the values from a single row definition and fill the data members.
void	print (bool beamGas)
	A function to print the contents of this particle.
long	eventNumber () const
	A function to return the event number of this particle.
int	particleId () const
	A function to return the particle identify of this particle.
double	kineticEnergy () const
	A function to return the relativistic kinetic energy of this particle.
double	weight () const
	A function to return the weight associated with this particle in the input ASCII file.
HepMC::FourVector	positionAtScoringPlane () const
	A function to the position of this particle with respect to the scoring plane.
HepMC::FourVector	directionalCosines () const
	A function to return the directional cosines of this particle.
double	timeOfFlight () const
	A function to return the time of flight for this particle.
double	primaryProtonZ () const
	A function to return the primary proton Z position for this particle.
const MarsShowerPrimary &	showerPrimary () const
	A function to return the shower primary infomation for this particle.
MarsBeamGasInteraction	beamGasInteraction () const
	A function to return the beam gas interaction information for this particle.
int	pdgId ()
	A function to return the PDG id for this particle.

Private Attributes
long	m_eventNumber
int	m_particleId
double	m_kineticEnergy
double	m_weight
HepMC::FourVector	m_positionAtScoringPlane
HepMC::FourVector	m_directionalCosines
double	m_timeOfFlight
double	m_primaryProtonZ
MarsShowerPrimary	m_showerPrimary
MarsBeamGasInteraction	m_beamGasInteraction

Detailed Description

A class to describe a MARS particle read from an input ASCII file.

Author: W. H. Bell W.Bel.nosp@m.l@ce.nosp@m.rn.ch

Definition at line 55 of file MarsParticle.h.

Constructor & Destructor Documentation

◆ MarsParticle()

MarsParticle::MarsParticle ( )

Definition at line 23 of file MarsParticle.cxx.

                          : m_eventNumber(0),
                  m_particleId(0),
                  m_kineticEnergy(0.),
                  m_weight(0.),
                  m_positionAtScoringPlane(),
                  m_directionalCosines(),
                  m_timeOfFlight(0.),
                  m_primaryProtonZ(0.),
                  m_showerPrimary(),
                  m_beamGasInteraction() {
}

Member Function Documentation

◆ beamGasInteraction()

MarsBeamGasInteraction MarsParticle::beamGasInteraction ( ) const

inline

A function to return the beam gas interaction information for this particle.

Definition at line 103 of file MarsParticle.h.

103{ return m_beamGasInteraction; }

◆ directionalCosines()

HepMC::FourVector MarsParticle::directionalCosines ( ) const

inline

A function to return the directional cosines of this particle.

Definition at line 88 of file MarsParticle.h.

88{ return m_directionalCosines; }

◆ eventNumber()

long MarsParticle::eventNumber ( ) const

inline

A function to return the event number of this particle.

Definition at line 69 of file MarsParticle.h.

69{ return m_eventNumber; }

◆ kineticEnergy()

double MarsParticle::kineticEnergy ( ) const

inline

A function to return the relativistic kinetic energy of this particle.

Definition at line 76 of file MarsParticle.h.

76{ return m_kineticEnergy; }

◆ particleId()

int MarsParticle::particleId ( ) const

inline

A function to return the particle identify of this particle.

Definition at line 72 of file MarsParticle.h.

72{ return m_particleId; }

◆ pdgId()

int MarsParticle::pdgId ( )

A function to return the PDG id for this particle.

Definition at line 115 of file MarsParticle.cxx.

                        {
  int pdgId = 0;
 
  switch (m_particleId) {
  case 1 : pdgId = 2212; break;
  case 2 : pdgId = 2112; break;
  case 3 : pdgId = 211; break;
  case 4 : pdgId = -211; break;
  case 5 : pdgId = 321; break;
  case 6 : pdgId = -321; break;
  case 7 : pdgId = -13; break;
  case 8 : pdgId = 13; break;
  case 9 : pdgId = 22; break;
  case 10 : pdgId = 11; break;
  case 11 : pdgId = -11; break;
  case 12 : pdgId = -2212; break;
  case 13 : pdgId = 111; break;
  case 14 : pdgId = 1000010020; break; // Deuterium
  case 15 : pdgId = 1000010030; break; // Tritium
  case 16 : pdgId = 1000020030; break; // He3
  case 17 : pdgId = 1000020040; break; // Alpha-(He4)
  case 18 : pdgId = 14; break;
  case 19 : pdgId = -14; break;
  case 20 : pdgId = 12; break;
  case 21 : pdgId = -12; break;
  case 22 : pdgId = 130; break;
  case 23 : pdgId = 310; break;
  case 24 : pdgId = 311; break;
  case 25 : pdgId = -311; break;
  case 26 : pdgId = 3122; break;
  case 27 : pdgId = -3122; break;
  case 28 : pdgId = 3222; break;
  case 29 : pdgId = 3212; break;
  case 30 : pdgId = 3112; break;
  case 31 : pdgId = -2112; break;
  case 32 : pdgId = 3322; break;
  case 33 : pdgId = 3312; break;
  case 34 : pdgId = 3334; break;
  case 35 : pdgId = -3112; break;
  case 36 : pdgId = -3212; break;
  case 37 : pdgId = -3222; break;
  case 38 : pdgId = -3322; break;
  case 39 : pdgId = -3312; break;
  case 40 : pdgId = -3334; break;
  default : pdgId = 0; break;
  }
 
  return pdgId;
}

◆ positionAtScoringPlane()

HepMC::FourVector MarsParticle::positionAtScoringPlane ( ) const

inline

A function to the position of this particle with respect to the scoring plane.

Definition at line 84 of file MarsParticle.h.

84{ return m_positionAtScoringPlane; }

◆ primaryProtonZ()

double MarsParticle::primaryProtonZ ( ) const

inline

A function to return the primary proton Z position for this particle.

Definition at line 95 of file MarsParticle.h.

95{ return m_primaryProtonZ; }

◆ print()

void MarsParticle::print ( bool beamGas )

A function to print the contents of this particle.

If the beamGas value is true then the beam gas information is printed.

Definition at line 80 of file MarsParticle.cxx.

                                     {
  std::cout.fill(' ');
  std::cout.width(11); std::cout << m_eventNumber << " ";
  std::cout.width(6); std::cout << m_particleId << " ";
  std::cout.width(12); std::cout.precision(5); std::cout << std::scientific << m_kineticEnergy << " "; 
  std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_weight << " "; 
  std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_positionAtScoringPlane.x() << " ";
  std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_positionAtScoringPlane.y() << " ";
  std::cout.width(14); std::cout.precision(7); std::cout << std::scientific << m_directionalCosines.x() << " ";
  std::cout.width(14); std::cout.precision(7); std::cout << std::scientific << m_directionalCosines.y() << " ";
  std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_timeOfFlight << " "; 
  std::cout.width(12); std::cout.precision(5); std::cout << std::scientific << m_primaryProtonZ << " ";
  std::cout.width(12); std::cout.precision(5); std::cout << std::scientific << m_showerPrimary.position.z() << " "; 
  std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_showerPrimary.position.x() << " "; 
  std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_showerPrimary.position.y() << " "; 
  std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_showerPrimary.kineticEnergy << " "; 
  std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_showerPrimary.weight << " ";
  std::cout.fill(' ');
  std::cout.width(3); std::cout << m_showerPrimary.particleId << " ";
  std::cout.width(3); std::cout << m_showerPrimary.processId;
  
  if(beamGas) {
    std::cout << " ";
    std::cout.width(3); std::cout << m_beamGasInteraction.nevve << " ";
    std::cout.width(3); std::cout << m_beamGasInteraction.nucleusType << " ";
    std::cout.width(11); std::cout.precision(4); std::cout << std::scientific << m_beamGasInteraction.timeOfFlight;
  }
 
  std::cout << std::endl;
 
  std::cout.precision(6);
}

◆ read()

int MarsParticle::read ( std::vector< std::string > * eventAsStringVector )

A function to read the values from a single row definition and fill the data members.

Definition at line 37 of file MarsParticle.cxx.

                                                                {
  if(eventAsStringVector->size() != 17 &&
     eventAsStringVector->size() != 20) {
    std::cerr << "There are only " << eventAsStringVector->size() << " data words.  This event will be null." << std::endl;
    return 1;
  }
  
  double doubleValue;
  std::vector<std::string>::iterator itr = eventAsStringVector->begin();
  std::vector<std::string>::iterator itr_end = eventAsStringVector->end();
  int wordNumber = 0;
  for(;itr!=itr_end;++itr,++wordNumber) {
    std::istringstream inStr((*itr));
    switch (wordNumber) {
    case 0 : inStr >> m_eventNumber; break;
    case 1 : inStr >> m_particleId; break;
    case 2 : inStr >> m_kineticEnergy; break;
    case 3 : inStr >> m_weight; break;
    case 4 : inStr >> doubleValue; m_positionAtScoringPlane.setX(doubleValue); break;
    case 5 : inStr >> doubleValue; m_positionAtScoringPlane.setY(doubleValue); break;
    case 6 : inStr >> doubleValue; m_directionalCosines.setX(doubleValue); break;
    case 7 : inStr >> doubleValue; m_directionalCosines.setY(doubleValue); break;
    case 8 : inStr >> m_timeOfFlight; break;
    case 9 : inStr >> m_primaryProtonZ; break;
    case 10 : inStr >> doubleValue; m_showerPrimary.position.setZ(doubleValue); break;
    case 11 : inStr >> doubleValue; m_showerPrimary.position.setX(doubleValue); break;
    case 12 : inStr >> doubleValue; m_showerPrimary.position.setY(doubleValue); break;
    case 13 : inStr >> m_showerPrimary.kineticEnergy; break;
    case 14 : inStr >> m_showerPrimary.weight; break;
    case 15 : inStr >> m_showerPrimary.particleId; break;
    case 16 : inStr >> m_showerPrimary.processId; break;
    case 17 : inStr >> m_beamGasInteraction.nevve; break;
    case 18 : inStr >> m_beamGasInteraction.nucleusType; break;
    case 19 : inStr >> m_beamGasInteraction.timeOfFlight; break;
    default : break;
    }
  }
 
  return 0;
}

◆ showerPrimary()

const MarsShowerPrimary & MarsParticle::showerPrimary ( ) const

inline

A function to return the shower primary infomation for this particle.

Definition at line 99 of file MarsParticle.h.

99{ return m_showerPrimary; }

◆ timeOfFlight()

double MarsParticle::timeOfFlight ( ) const

inline

A function to return the time of flight for this particle.

Definition at line 91 of file MarsParticle.h.

91{ return m_timeOfFlight; }

◆ weight()

double MarsParticle::weight ( ) const

inline

A function to return the weight associated with this particle in the input ASCII file.

Definition at line 80 of file MarsParticle.h.

80{ return m_weight; }

Member Data Documentation

◆ m_beamGasInteraction

MarsBeamGasInteraction MarsParticle::m_beamGasInteraction

private

Definition at line 118 of file MarsParticle.h.

◆ m_directionalCosines

HepMC::FourVector MarsParticle::m_directionalCosines

private

Definition at line 114 of file MarsParticle.h.

◆ m_eventNumber

long MarsParticle::m_eventNumber

private

Definition at line 109 of file MarsParticle.h.

◆ m_kineticEnergy

double MarsParticle::m_kineticEnergy

private

Definition at line 111 of file MarsParticle.h.

◆ m_particleId

int MarsParticle::m_particleId

private

Definition at line 110 of file MarsParticle.h.

◆ m_positionAtScoringPlane

HepMC::FourVector MarsParticle::m_positionAtScoringPlane

private

Definition at line 113 of file MarsParticle.h.

◆ m_primaryProtonZ

double MarsParticle::m_primaryProtonZ

private

Definition at line 116 of file MarsParticle.h.

◆ m_showerPrimary

MarsShowerPrimary MarsParticle::m_showerPrimary

private

Definition at line 117 of file MarsParticle.h.

◆ m_timeOfFlight

double MarsParticle::m_timeOfFlight

private

Definition at line 115 of file MarsParticle.h.

◆ m_weight

double MarsParticle::m_weight

private

Definition at line 112 of file MarsParticle.h.

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

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MarsParticle()

Member Function Documentation

◆ beamGasInteraction()

◆ directionalCosines()

◆ eventNumber()

◆ kineticEnergy()

◆ particleId()

◆ pdgId()

◆ positionAtScoringPlane()

◆ primaryProtonZ()

◆ print()

◆ read()

◆ showerPrimary()

◆ timeOfFlight()

◆ weight()

Member Data Documentation

◆ m_beamGasInteraction

◆ m_directionalCosines

◆ m_eventNumber

◆ m_kineticEnergy

◆ m_particleId

◆ m_positionAtScoringPlane

◆ m_primaryProtonZ

◆ m_showerPrimary

◆ m_timeOfFlight

◆ m_weight