#include <MaterialBudgetData.h>

Public Member Functions
bool	allStepsON ()
void	dataEndTrack (const G4Track *aTrack)
void	dataPerStep (const G4Step *aStep)
void	dataStartTrack (const G4Track *aTrack)
float	getAirDil (int is) const
float	getAirDmb (int is) const
float	getAirFractionIL () const
float	getAirFractionMB () const
float	getAirIL () const
float	getAirMB () const
float	getCablesDil (int is) const
float	getCablesDmb (int is) const
float	getCablesFractionIL () const
float	getCablesFractionMB () const
float	getCablesIL () const
float	getCablesMB () const
float	getCoolingDil (int is) const
float	getCoolingDmb (int is) const
float	getCoolingFractionIL () const
float	getCoolingFractionMB () const
float	getCoolingIL () const
float	getCoolingMB () const
float	getElectronicsDil (int is) const
float	getElectronicsDmb (int is) const
float	getElectronicsFractionIL () const
float	getElectronicsFractionMB () const
float	getElectronicsIL () const
float	getElectronicsMB () const
float	getEnergy () const
float	getEta () const
int	getID () const
float	getIntLen () const
float	getMass () const
int	getNumberOfSteps () const
float	getOtherDil (int is) const
float	getOtherDmb (int is) const
float	getOtherFractionIL () const
float	getOtherFractionMB () const
float	getOtherIL () const
float	getOtherMB () const
float	getPhi () const
float	getPt () const
int	getPVcopyNo () const
std::string	getPVname () const
float	getRadLen () const
float	getSensitiveDil (int is) const
float	getSensitiveDmb (int is) const
float	getSensitiveFractionIL () const
float	getSensitiveFractionMB () const
float	getSensitiveIL () const
float	getSensitiveMB () const
float	getStepDil (int is)
float	getStepDmb (int is)
float	getStepFinalBeta (int is)
float	getStepFinalEnergy (int is)
float	getStepFinalEta (int is)
float	getStepFinalGamma (int is)
float	getStepFinalMass (int is)
float	getStepFinalPhi (int is)
float	getStepFinalPt (int is)
float	getStepFinalPx (int is)
float	getStepFinalPy (int is)
float	getStepFinalPz (int is)
double	getStepFinalX (int is)
double	getStepFinalY (int is)
double	getStepFinalZ (int is)
int	getStepID (int is)
float	getStepInitialBeta (int is)
float	getStepInitialEnergy (int is)
float	getStepInitialEta (int is)
float	getStepInitialGamma (int is)
float	getStepInitialMass (int is)
float	getStepInitialPhi (int is)
float	getStepInitialPt (int is)
float	getStepInitialPx (int is)
float	getStepInitialPy (int is)
float	getStepInitialPz (int is)
double	getStepInitialX (int is)
double	getStepInitialY (int is)
double	getStepInitialZ (int is)
float	getStepMaterialDensity (int is)
int	getStepMaterialID (int is)
float	getStepMaterialLambda0 (int is)
std::string	getStepMaterialName (int is)
float	getStepMaterialX0 (int is)
int	getStepPostProcess (int is)
int	getStepPreProcess (int is)
int	getStepVolumeCopy (int is)
int	getStepVolumeID (int is)
std::string	getStepVolumeName (int is)
float	getStepVolumeX (int is)
CLHEP::HepLorentzVector	getStepVolumeXaxis (int is)
float	getStepVolumeY (int is)
CLHEP::HepLorentzVector	getStepVolumeYaxis (int is)
float	getStepVolumeZ (int is)
CLHEP::HepLorentzVector	getStepVolumeZaxis (int is)
float	getSupportDil (int is) const
float	getSupportDmb (int is) const
float	getSupportFractionIL () const
float	getSupportFractionMB () const
float	getSupportIL () const
float	getSupportMB () const
float	getTotalIL () const
float	getTotalMB () const
float	getTrkLen () const
	MaterialBudgetData ()
void	SetAllStepsToTree ()
	~MaterialBudgetData ()
Private Attributes
bool	allStepsToTree
double	densityConvertionFactor
int	MAXNUMBERSTEPS
MaterialBudgetCategorizer *	myMaterialBudgetCategorizer
int	stepN
float *	theAirDil
float *	theAirDmb
float	theAirFractionIL
float	theAirFractionMB
float	theAirIL
float	theAirMB
float *	theCablesDil
float *	theCablesDmb
float	theCablesFractionIL
float	theCablesFractionMB
float	theCablesIL
float	theCablesMB
float *	theCoolingDil
float *	theCoolingDmb
float	theCoolingFractionIL
float	theCoolingFractionMB
float	theCoolingIL
float	theCoolingMB
float *	theDil
float *	theDmb
float *	theElectronicsDil
float *	theElectronicsDmb
float	theElectronicsFractionIL
float	theElectronicsFractionMB
float	theElectronicsIL
float	theElectronicsMB
float	theEnergy
float	theEta
double *	theFinalX
double *	theFinalY
double *	theFinalZ
int	theID
double *	theInitialX
double *	theInitialY
double *	theInitialZ
float	theIntLen
float	theMass
float *	theMaterialDensity
int *	theMaterialID
float *	theMaterialLambda0
std::string *	theMaterialName
float *	theMaterialX0
float *	theOtherDil
float *	theOtherDmb
float	theOtherFractionIL
float	theOtherFractionMB
float	theOtherIL
float	theOtherMB
float	thePhi
float	thePt
int	thePVcopyNo
std::string	thePVname
float	theRadLen
float *	theSensitiveDil
float *	theSensitiveDmb
float	theSensitiveFractionIL
float	theSensitiveFractionMB
float	theSensitiveIL
float	theSensitiveMB
float *	theStepFinalBeta
float *	theStepFinalEnergy
float *	theStepFinalEta
float *	theStepFinalGamma
float *	theStepFinalMass
float *	theStepFinalPhi
float *	theStepFinalPt
float *	theStepFinalPx
float *	theStepFinalPy
float *	theStepFinalPz
int *	theStepID
float *	theStepInitialBeta
float *	theStepInitialEnergy
float *	theStepInitialEta
float *	theStepInitialGamma
float *	theStepInitialMass
float *	theStepInitialPhi
float *	theStepInitialPt
float *	theStepInitialPx
float *	theStepInitialPy
float *	theStepInitialPz
int	theStepN
int *	theStepPostProcess
int *	theStepPreProcess
float *	theSupportDil
float *	theSupportDmb
float	theSupportFractionIL
float	theSupportFractionMB
float	theSupportIL
float	theSupportMB
float	theTotalIL
float	theTotalMB
float	theTrkLen
int *	theVolumeCopy
int *	theVolumeID
std::string *	theVolumeName
float *	theVolumeX
float *	theVolumeXaxis1
float *	theVolumeXaxis2
float *	theVolumeXaxis3
float *	theVolumeY
float *	theVolumeYaxis1
float *	theVolumeYaxis2
float *	theVolumeYaxis3
float *	theVolumeZ
float *	theVolumeZaxis1
float *	theVolumeZaxis2
float *	theVolumeZaxis3

Detailed Description

Definition at line 18 of file MaterialBudgetData.h.

Constructor & Destructor Documentation

MaterialBudgetData::MaterialBudgetData ( )

Definition at line 12 of file MaterialBudgetData.cc.

References allStepsToTree, densityConvertionFactor, and myMaterialBudgetCategorizer.

{
  //instantiate categorizer to assign an ID to volumes and materials
  myMaterialBudgetCategorizer = 0;
  allStepsToTree = false;
  densityConvertionFactor = 6.24E18;
}

MaterialBudgetData::~MaterialBudgetData ( )

Definition at line 20 of file MaterialBudgetData.cc.

{
}

Member Function Documentation

bool MaterialBudgetData::allStepsON ( ) [inline]

Definition at line 289 of file MaterialBudgetData.h.

References allStepsToTree.

Referenced by MaterialBudgetTree::book(), and MaterialBudgetTree::fillEndTrack().

                    {
    return allStepsToTree;
  }

void MaterialBudgetData::dataEndTrack ( const G4Track * aTrack )

Definition at line 167 of file MaterialBudgetData.cc.

References gather_cfg::cout, theAirIL, theAirMB, theCablesIL, theCablesMB, theCoolingIL, theCoolingMB, theElectronicsIL, theElectronicsMB, theEta, theOtherIL, theOtherMB, thePhi, theSensitiveIL, theSensitiveMB, theStepN, theSupportIL, theSupportMB, theTotalIL, and theTotalMB.

Referenced by MaterialBudgetAction::update().

{
  //-  std::cout << "[OVAL] MaterialBudget " << G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID() << " " << theEta << " " << thePhi << " " << theTotalMB << std::endl;
  // rr
  std::cout << "Recorded steps " << theStepN << std::endl;
  std::cout << " Material Budget: Radiation Length   " << G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID() << " eta " << theEta << " phi " << thePhi << " total X " << theTotalMB << " SUP " << theSupportMB << " SEN " << theSensitiveMB << " CAB " << theCablesMB << " COL " << theCoolingMB << " ELE " << theElectronicsMB << " other " << theOtherMB << " Air " << theAirMB << std::endl;
  std::cout << " Material Budget: Interaction Length " << G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID() << " eta " << theEta << " phi " << thePhi << " total L " << theTotalIL << " SUP " << theSupportIL << " SEN " << theSensitiveIL << " CAB " << theCablesIL << " COL " << theCoolingIL << " ELE " << theElectronicsIL << " other " << theOtherIL << " Air " << theAirIL << std::endl;
  // rr
}

void MaterialBudgetData::dataPerStep ( const G4Step * aStep )

Definition at line 178 of file MaterialBudgetData.cc.

References allStepsToTree, gather_cfg::cout, densityConvertionFactor, MaterialBudgetCategorizer::l0fraction(), MaterialBudgetCategorizer::material(), MAXNUMBERSTEPS, myMaterialBudgetCategorizer, mergeVDriftHistosByStation::name, mathSSE::sqrt(), stepN, matplotRender::t, theAirDil, theAirDmb, theAirFractionIL, theAirFractionMB, theAirIL, theAirMB, theCablesDil, theCablesDmb, theCablesFractionIL, theCablesFractionMB, theCablesIL, theCablesMB, theCoolingDil, theCoolingDmb, theCoolingFractionIL, theCoolingFractionMB, theCoolingIL, theCoolingMB, theDil, theDmb, theElectronicsDil, theElectronicsDmb, theElectronicsFractionIL, theElectronicsFractionMB, theElectronicsIL, theElectronicsMB, theEnergy, theEta, theFinalX, theFinalY, theFinalZ, theID, theInitialX, theInitialY, theInitialZ, theIntLen, theMass, theMaterialDensity, theMaterialID, theMaterialLambda0, theMaterialName, theMaterialX0, theOtherDil, theOtherDmb, theOtherFractionIL, theOtherFractionMB, theOtherIL, theOtherMB, thePhi, thePt, thePVcopyNo, thePVname, theRadLen, theSensitiveDil, theSensitiveDmb, theSensitiveFractionIL, theSensitiveFractionMB, theSensitiveIL, theSensitiveMB, theStepFinalBeta, theStepFinalEnergy, theStepFinalEta, theStepFinalGamma, theStepFinalMass, theStepFinalPhi, theStepFinalPt, theStepFinalPx, theStepFinalPy, theStepFinalPz, theStepID, theStepInitialBeta, theStepInitialEnergy, theStepInitialEta, theStepInitialGamma, theStepInitialMass, theStepInitialPhi, theStepInitialPt, theStepInitialPx, theStepInitialPy, theStepInitialPz, theStepN, theStepPostProcess, theStepPreProcess, theSupportDil, theSupportDmb, theSupportFractionIL, theSupportFractionMB, theSupportIL, theSupportMB, theTotalIL, theTotalMB, theTrkLen, theVolumeCopy, theVolumeID, theVolumeName, theVolumeX, theVolumeXaxis1, theVolumeXaxis2, theVolumeXaxis3, theVolumeY, theVolumeYaxis1, theVolumeYaxis2, theVolumeYaxis3, theVolumeZ, theVolumeZaxis1, theVolumeZaxis2, theVolumeZaxis3, MaterialBudgetCategorizer::volume(), and MaterialBudgetCategorizer::x0fraction().

Referenced by MaterialBudgetAction::update().

{
  G4Material * theMaterialPre = aStep->GetPreStepPoint()->GetMaterial();
  //  G4Material * theMaterialPost = aStep->GetPostStepPoint()->GetMaterial();

  G4StepPoint* prePoint  = aStep->GetPreStepPoint();
  G4StepPoint* postPoint = aStep->GetPostStepPoint();
  
  CLHEP::Hep3Vector prePos  = prePoint->GetPosition();
  CLHEP::Hep3Vector postPos = postPoint->GetPosition();

  G4double steplen = aStep->GetStepLength();

  G4double radlen;
  G4double intlen;
  G4double density;

  radlen  = theMaterialPre->GetRadlen();
  intlen  = theMaterialPre->GetNuclearInterLength();
  density = theMaterialPre->GetDensity() / densityConvertionFactor; // always g/cm3
  
  G4String name = theMaterialPre->GetName();
  //  std::cout << " steplen " << steplen << " radlen " << radlen << " mb " << steplen/radlen << " mate " << theMaterialPre->GetName() << std::endl;
     
  G4LogicalVolume* lv = aStep->GetTrack()->GetVolume()->GetLogicalVolume();
  
  // instantiate the categorizer
  int volumeID   = myMaterialBudgetCategorizer->volume( lv->GetName() );
  int materialID = myMaterialBudgetCategorizer->material( lv->GetMaterial()->GetName() );
  // rr
  std::string volumeName   = lv->GetName();
  std::string materialName = lv->GetMaterial()->GetName();
  // rr
  
  // rr
  /*
    std::cout << " Volume "   << lv->GetName()                << std::endl;
    std::cout << " Material " << lv->GetMaterial()->GetName() << std::endl;
  */    
  theSupportFractionMB     = myMaterialBudgetCategorizer->x0fraction(lv->GetMaterial()->GetName())[0];
  theSensitiveFractionMB   = myMaterialBudgetCategorizer->x0fraction(lv->GetMaterial()->GetName())[1];
  theCablesFractionMB      = myMaterialBudgetCategorizer->x0fraction(lv->GetMaterial()->GetName())[2];
  theCoolingFractionMB     = myMaterialBudgetCategorizer->x0fraction(lv->GetMaterial()->GetName())[3];
  theElectronicsFractionMB = myMaterialBudgetCategorizer->x0fraction(lv->GetMaterial()->GetName())[4];
  theOtherFractionMB       = myMaterialBudgetCategorizer->x0fraction(lv->GetMaterial()->GetName())[5];
  theAirFractionMB         = myMaterialBudgetCategorizer->x0fraction(lv->GetMaterial()->GetName())[6];
  if(theOtherFractionMB!=0) std::cout << " material found with no category " << lv->GetMaterial()->GetName() 
                                      << " in volume " << lv->GetName() << std::endl;
  theSupportFractionIL     = myMaterialBudgetCategorizer->l0fraction(lv->GetMaterial()->GetName())[0];
  theSensitiveFractionIL   = myMaterialBudgetCategorizer->l0fraction(lv->GetMaterial()->GetName())[1];
  theCablesFractionIL      = myMaterialBudgetCategorizer->l0fraction(lv->GetMaterial()->GetName())[2];
  theCoolingFractionIL     = myMaterialBudgetCategorizer->l0fraction(lv->GetMaterial()->GetName())[3];
  theElectronicsFractionIL = myMaterialBudgetCategorizer->l0fraction(lv->GetMaterial()->GetName())[4];
  theOtherFractionIL       = myMaterialBudgetCategorizer->l0fraction(lv->GetMaterial()->GetName())[5];
  theAirFractionIL         = myMaterialBudgetCategorizer->l0fraction(lv->GetMaterial()->GetName())[6];
  if(theOtherFractionIL!=0) std::cout << " material found with no category " << lv->GetMaterial()->GetName() 
                                      << " in volume " << lv->GetName() << std::endl;
  //  if(theOtherFractionMB!=0) LogDebug("MaterialBudgetData") << " material found with no category " << lv->GetMaterial()->GetName() 
  //                             << " in volume " << lv->GetName();
  // rr  
  
  float dmb = steplen/radlen;
  float dil = steplen/intlen;
  
  G4VPhysicalVolume*       pv                = aStep->GetPreStepPoint()->GetPhysicalVolume();
  const G4VTouchable*      t                 = aStep->GetPreStepPoint()->GetTouchable();
  G4ThreeVector            objectTranslation = t->GetTranslation();
  const G4RotationMatrix*  objectRotation    = t->GetRotation();
  const G4VProcess*        interactionPre    = prePoint->GetProcessDefinedStep();
  const G4VProcess*        interactionPost   = postPoint->GetProcessDefinedStep();
  
  G4Track* track = aStep->GetTrack();
  if(theStepN==0) std::cout << " Simulated Particle " << theID << "\tMass " << theMass << " MeV/c2"
                            << "\tPt = " << thePt  << " MeV/c" << "\tEta = " << theEta << "\tPhi = " << thePhi 
                            << "\tEnergy = " << theEnergy << " MeV"
                    //                      << std::endl
                    //                      << "\tMagnetic Field at (0,0,0): (" << B000[0] << "," < B000[1] << "," << B000[2] << ")" 
                            << std::endl;
  
  //fill data per step
  if( allStepsToTree ){
    if( stepN > MAXNUMBERSTEPS ) stepN = MAXNUMBERSTEPS;
    theDmb[theStepN] = dmb; 
    theDil[theStepN] = dil; 
    theSupportDmb[theStepN]     = (dmb * theSupportFractionMB);
    theSensitiveDmb[theStepN]   = (dmb * theSensitiveFractionMB);
    theCablesDmb[theStepN]      = (dmb * theCablesFractionMB);
    theCoolingDmb[theStepN]     = (dmb * theCoolingFractionMB);
    theElectronicsDmb[theStepN] = (dmb * theElectronicsFractionMB);
    theOtherDmb[theStepN]       = (dmb * theOtherFractionMB);
    theAirDmb[theStepN]         = (dmb * theAirFractionMB);
    theSupportDil[theStepN]     = (dil * theSupportFractionIL);
    theSensitiveDil[theStepN]   = (dil * theSensitiveFractionIL);
    theCablesDil[theStepN]      = (dil * theCablesFractionIL);
    theCoolingDil[theStepN]     = (dil * theCoolingFractionIL);
    theElectronicsDil[theStepN] = (dil * theElectronicsFractionIL);
    theOtherDil[theStepN]       = (dil * theOtherFractionIL);
    theAirDil[theStepN]         = (dil * theAirFractionIL);
    theInitialX[theStepN] = prePos.x();
    theInitialY[theStepN] = prePos.y();
    theInitialZ[theStepN] = prePos.z();
    theFinalX[theStepN]   = postPos.x();
    theFinalY[theStepN]   = postPos.y();
    theFinalZ[theStepN]   = postPos.z();
    theVolumeID[theStepN]   = volumeID;
    theVolumeName[theStepN] = volumeName;
    theVolumeCopy[theStepN] = pv->GetCopyNo();
    theVolumeX[theStepN]    = objectTranslation.x();
    theVolumeY[theStepN]    = objectTranslation.y();
    theVolumeZ[theStepN]    = objectTranslation.z();
    theVolumeXaxis1[theStepN] = objectRotation->xx();
    theVolumeXaxis2[theStepN] = objectRotation->xy();
    theVolumeXaxis3[theStepN] = objectRotation->xz();
    theVolumeYaxis1[theStepN] = objectRotation->yx();
    theVolumeYaxis2[theStepN] = objectRotation->yy();
    theVolumeYaxis3[theStepN] = objectRotation->yz();
    theVolumeZaxis1[theStepN] = objectRotation->zx();
    theVolumeZaxis2[theStepN] = objectRotation->zy();
    theVolumeZaxis3[theStepN] = objectRotation->zz();
    theMaterialID[theStepN]      = materialID;
    theMaterialName[theStepN]    = materialName;
    theMaterialX0[theStepN]      = radlen;
    theMaterialLambda0[theStepN] = intlen;
    theMaterialDensity[theStepN] = density;
    theStepID[theStepN]             = track->GetDefinition()->GetPDGEncoding();
    theStepInitialPt[theStepN]      = prePoint->GetMomentum().perp();
    theStepInitialEta[theStepN]     = prePoint->GetMomentum().eta();
    theStepInitialPhi[theStepN]     = prePoint->GetMomentum().phi();
    theStepInitialEnergy[theStepN]  = prePoint->GetTotalEnergy();
    theStepInitialPx[theStepN]      = prePoint->GetMomentum().x();
    theStepInitialPy[theStepN]      = prePoint->GetMomentum().y();
    theStepInitialPz[theStepN]      = prePoint->GetMomentum().z();
    theStepInitialBeta[theStepN]    = prePoint->GetBeta();
    theStepInitialGamma[theStepN]   = prePoint->GetGamma();
    theStepInitialMass[theStepN]    = prePoint->GetMass();
    theStepFinalPt[theStepN]        = postPoint->GetMomentum().perp();
    theStepFinalEta[theStepN]       = postPoint->GetMomentum().eta();
    theStepFinalPhi[theStepN]       = postPoint->GetMomentum().phi();
    theStepFinalEnergy[theStepN]    = postPoint->GetTotalEnergy();
    theStepFinalPx[theStepN]        = postPoint->GetMomentum().x();
    theStepFinalPy[theStepN]        = postPoint->GetMomentum().y();
    theStepFinalPz[theStepN]        = postPoint->GetMomentum().z();
    theStepFinalBeta[theStepN]      = postPoint->GetBeta();
    theStepFinalGamma[theStepN]     = postPoint->GetGamma();
    theStepFinalMass[theStepN]      = postPoint->GetMass();
    int preProcType  = -99;
    int postProcType = -99;
    if (interactionPre) preProcType = interactionPre->GetProcessType();
    theStepPreProcess[theStepN]     = preProcType;
    if (interactionPost) postProcType = interactionPost->GetProcessType();
    theStepPostProcess[theStepN]    = postProcType;
#ifdef TREE_DEBUG
    std::cout << " step " << theStepN
              << "\tDelta MB = " << theDmb[theStepN]
              << std::endl
              << "\t\tSupport "  << theSupportDmb[theStepN]
              << " Sensitive "   << theSensitiveDmb[theStepN]
              << " Cables "      << theCablesDmb[theStepN]
              << " Cooling "     << theCoolingDmb[theStepN]
              << " Electronics " << theElectronicsDmb[theStepN]
              << " Other "       << theOtherDmb[theStepN]
              << " Air "         << theAirDmb[theStepN]
              << std::endl
              << "\tDelta IL = " << theDil[theStepN]
              << std::endl
              << "\t\tSupport "  << theSupportDil[theStepN]
              << " Sensitive "   << theSensitiveDil[theStepN]
              << " Cables "      << theCablesDil[theStepN]
              << " Cooling "     << theCoolingDil[theStepN]
              << " Electronics " << theElectronicsDil[theStepN]
              << " Other "       << theOtherDil[theStepN]
              << " Air "         << theAirDil[theStepN]
              << std::endl;
    if (interactionPre)
      std::cout << "\tProcess Pre " << interactionPre->GetProcessName()
                << " type " << theStepPreProcess[theStepN] << " " << interactionPre->GetProcessTypeName(G4ProcessType(theStepPreProcess[theStepN]))
                << std::endl;
    if (interactionPost)
      std::cout << "\tProcess Post " << interactionPost->GetProcessName()
                << " type " << theStepPostProcess[theStepN] << " "
                << interactionPost->GetProcessTypeName(G4ProcessType(theStepPostProcess[theStepN]))
                << std::endl;
    std::cout << "\tPre x = " << theInitialX[theStepN]
              << "\ty = "     << theInitialY[theStepN]
              << "\tz = "     << theInitialZ[theStepN] 
              << "\tPolar Radius = " << sqrt(prePos.x()*prePos.x()+prePos.y()*prePos.y())
              << "\tPt = "     << theStepInitialPt[theStepN]
              << "\tEnergy = " << theStepInitialEnergy[theStepN]
      //              << std::endl
      //              << "B-field(T) at Pre (cm): " << field->fieldInTesla(GlobalPoint(pos.x()/10.,pos.y()/10.,pos.z()/10.))
              << std::endl;
    std::cout  << "\tPost x = " << theFinalX[theStepN]
               << "\ty = "      << theFinalY[theStepN]
               << "\tz = "      << theFinalZ[theStepN] 
               << "\tPolar Radius = " << sqrt(postPos.x()*postPos.x()+postPos.y()*postPos.y())
              << "\tPt = "     << theStepFinalPt[theStepN]
               << "\tEnergy = " << theStepFinalEnergy[theStepN]
               << std::endl;
    std::cout << "\tvolume " << volumeID << " " << theVolumeName[theStepN] 
              << " copy number " << theVolumeCopy[theStepN]
              << "\tmaterial " << theMaterialID[theStepN] << " " << theMaterialName[theStepN]
              << "\tDensity = " << theMaterialDensity[theStepN] << " g/cm3"
              << "\tX0 = " << theMaterialX0[theStepN] << " mm"
              << "\tLambda0 = " << theMaterialLambda0[theStepN] << " mm"
              << std::endl;
    std::cout << "\t\tParticle "  << theStepID[theStepN] 
              << " Initial Pt = " << theStepInitialPt[theStepN]     << " MeV/c"
              << " eta = "        << theStepInitialEta[theStepN]
              << " phi = "        << theStepInitialPhi[theStepN]
              << " Energy = "     << theStepInitialEnergy[theStepN] << " MeV"
              << " Mass = "       << theStepInitialMass[theStepN]   << " MeV/c2"
              << " Beta = "       << theStepInitialBeta[theStepN]
              << " Gamma = "      << theStepInitialGamma[theStepN]
              << std::endl
              << "\t\tParticle "  << theStepID[theStepN]
              << " Final Pt = "   << theStepFinalPt[theStepN]       << " MeV/c"
              << " eta = "        << theStepFinalEta[theStepN]
              << " phi = "        << theStepFinalPhi[theStepN]
              << " Energy = "     << theStepFinalEnergy[theStepN]   << " MeV"
              << " Mass = "       << theStepFinalMass[theStepN]     << " MeV/c2"
              << " Beta = "       << theStepFinalBeta[theStepN]
              << " Gamma = "      << theStepFinalGamma[theStepN]
              << std::endl;
    std:: cout << "\tVolume Centre x = " << theVolumeX[theStepN]
               << "\ty = "               << theVolumeY[theStepN]
               << "\tz = "               << theVolumeZ[theStepN]
               << "\tPolar Radius = "    << sqrt( theVolumeX[theStepN]*theVolumeX[theStepN] +
                                                  theVolumeY[theStepN]*theVolumeY[theStepN] )
               << std::endl;
    std::cout << "\tx axis = (" 
              << theVolumeXaxis1[theStepN] << "," 
              << theVolumeXaxis2[theStepN] << "," 
              << theVolumeXaxis3[theStepN] << ")"
              << std::endl;
    std::cout << "\ty axis = (" 
              << theVolumeYaxis1[theStepN] << "," 
              << theVolumeYaxis2[theStepN] << "," 
              << theVolumeYaxis3[theStepN] << ")"
              << std::endl;
    std::cout << "\tz axis = (" 
              << theVolumeZaxis1[theStepN] << "," 
              << theVolumeZaxis2[theStepN] << "," 
              << theVolumeZaxis3[theStepN] << ")"
              << std::endl;
    std::cout << "\tSecondaries"
              << std::endl;
    for(G4TrackVector::iterator iSec = aStep->GetSecondary()->begin(); iSec!=aStep->GetSecondary()->end(); iSec++) {
      std::cout << "\t\tid " << (*iSec)->GetDefinition()->GetPDGEncoding()
                << " created through process "
                << " type " << (*iSec)->GetCreatorProcess()->GetProcessType()
                << " " << (*iSec)->GetCreatorProcess()->GetProcessTypeName(G4ProcessType((*iSec)->GetCreatorProcess()->GetProcessType()))
                << std::endl;
    }
#endif
  }
  
  theTrkLen = aStep->GetTrack()->GetTrackLength();
  //-  std::cout << " theTrkLen " << theTrkLen << " theTrkLen2 " << theTrkLen2 << " postPos " << postPos.mag() << postPos << std::endl;
  thePVname = pv->GetName();
  thePVcopyNo = pv->GetCopyNo();
  theRadLen = radlen;
  theIntLen = intlen;
  theTotalMB += dmb;
  theTotalIL += dil;
  
  // rr
  theSupportMB     += (dmb * theSupportFractionMB);
  theSensitiveMB   += (dmb * theSensitiveFractionMB);
  theCablesMB      += (dmb * theCablesFractionMB);
  theCoolingMB     += (dmb * theCoolingFractionMB);
  theElectronicsMB += (dmb * theElectronicsFractionMB);
  theOtherMB       += (dmb * theOtherFractionMB);
  theAirMB         += (dmb * theAirFractionMB);
  theSupportIL     += (dil * theSupportFractionIL);
  theSensitiveIL   += (dil * theSensitiveFractionIL);
  theCablesIL      += (dil * theCablesFractionIL);
  theCoolingIL     += (dil * theCoolingFractionIL);
  theElectronicsIL += (dil * theElectronicsFractionIL);
  theOtherIL       += (dil * theOtherFractionIL);
  theAirIL         += (dil * theAirFractionIL);
  // rr
  
  theStepN++;
  
}

void MaterialBudgetData::dataStartTrack ( const G4Track * aTrack )

Definition at line 100 of file MaterialBudgetData.cc.

References gather_cfg::cout, dir, myMaterialBudgetCategorizer, theAirFractionIL, theAirFractionMB, theAirIL, theAirMB, theCablesFractionIL, theCablesFractionMB, theCablesIL, theCablesMB, theCoolingFractionIL, theCoolingFractionMB, theCoolingIL, theCoolingMB, theElectronicsFractionIL, theElectronicsFractionMB, theElectronicsIL, theElectronicsMB, theEnergy, theEta, theID, theMass, theOtherFractionIL, theOtherFractionMB, theOtherIL, theOtherMB, thePhi, thePt, theSensitiveFractionIL, theSensitiveFractionMB, theSensitiveIL, theSensitiveMB, theStepN, theSupportFractionIL, theSupportFractionMB, theSupportIL, theSupportMB, theTotalIL, and theTotalMB.

Referenced by MaterialBudgetAction::update().

{
  // rr
  std::cout << "MaterialBudget Analysis of Event #" << G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID() << std::endl;
  // rr
  
  const G4ThreeVector& dir = aTrack->GetMomentum() ;
  
  if( myMaterialBudgetCategorizer == 0) myMaterialBudgetCategorizer = new MaterialBudgetCategorizer;
  
  theStepN=0;
  theTotalMB=0;
  theTotalIL=0;
  theEta=0;
  thePhi=0;
  
  // rr
  theID=0;
  thePt=0;
  theEnergy=0;
  theMass=0;
  
  theSupportMB     = 0.;
  theSensitiveMB   = 0.;
  theCablesMB      = 0.;
  theCoolingMB     = 0.;
  theElectronicsMB = 0.;
  theOtherMB       = 0.;
  theAirMB         = 0.;
  theSupportIL     = 0.;
  theSensitiveIL   = 0.;
  theCablesIL      = 0.;
  theCoolingIL     = 0.;
  theElectronicsIL = 0.;
  theOtherIL       = 0.;
  theAirIL         = 0.;
  theSupportFractionMB     = 0.;
  theSensitiveFractionMB   = 0.;
  theCablesFractionMB      = 0.;
  theCoolingFractionMB     = 0.;
  theElectronicsFractionMB = 0.;
  theOtherFractionMB       = 0.;
  theAirFractionMB         = 0.;
  theSupportFractionIL     = 0.;
  theSensitiveFractionIL   = 0.;
  theCablesFractionIL      = 0.;
  theCoolingFractionIL     = 0.;
  theElectronicsFractionIL = 0.;
  theOtherFractionIL       = 0.;
  theAirFractionIL         = 0.;
  // rr
  
  theID = (int)(aTrack->GetDefinition()->GetPDGEncoding());
  thePt = dir.perp();
  if( dir.theta() != 0 ) {
    theEta = dir.eta(); 
  } else {
    theEta = -99;
  }
  //  thePhi = dir.phi()/deg; // better not to store in deg
  thePhi = dir.phi();
  theEnergy = aTrack->GetTotalEnergy();
  theMass = aTrack->GetDefinition()->GetPDGMass();
  
}

float MaterialBudgetData::getAirDil ( int is ) const [inline]

Definition at line 154 of file MaterialBudgetData.h.

References theAirDil.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                  {
    return theAirDil[is]; }

float MaterialBudgetData::getAirDmb ( int is ) const [inline]

Definition at line 136 of file MaterialBudgetData.h.

References theAirDmb.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                  {
    return theAirDmb[is]; }

float MaterialBudgetData::getAirFractionIL ( ) const [inline]

Definition at line 73 of file MaterialBudgetData.h.

References theAirFractionIL.

                                 {
    return theAirFractionIL; }

float MaterialBudgetData::getAirFractionMB ( ) const [inline]

Definition at line 45 of file MaterialBudgetData.h.

References theAirFractionMB.

                                 {
    return theAirFractionMB; }

float MaterialBudgetData::getAirIL ( ) const [inline]

Definition at line 89 of file MaterialBudgetData.h.

References theAirIL.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                         {
    return theAirIL; }

float MaterialBudgetData::getAirMB ( ) const [inline]

Definition at line 59 of file MaterialBudgetData.h.

References theAirMB.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                         {
    return theAirMB; }

float MaterialBudgetData::getCablesDil ( int is ) const [inline]

Definition at line 146 of file MaterialBudgetData.h.

References theCablesDil.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                     {
    return theCablesDil[is]; }

float MaterialBudgetData::getCablesDmb ( int is ) const [inline]

Definition at line 128 of file MaterialBudgetData.h.

References theCablesDmb.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                     {
    return theCablesDmb[is]; }

float MaterialBudgetData::getCablesFractionIL ( ) const [inline]

Definition at line 65 of file MaterialBudgetData.h.

References theCablesFractionIL.

                                    {
    return theCablesFractionIL; }

float MaterialBudgetData::getCablesFractionMB ( ) const [inline]

Definition at line 37 of file MaterialBudgetData.h.

References theCablesFractionMB.

                                    {
    return theCablesFractionMB; }

float MaterialBudgetData::getCablesIL ( ) const [inline]

Definition at line 81 of file MaterialBudgetData.h.

References theCablesIL.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                            {
    return theCablesIL; }

float MaterialBudgetData::getCablesMB ( ) const [inline]

Definition at line 51 of file MaterialBudgetData.h.

References theCablesMB.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                            {
    return theCablesMB; }

float MaterialBudgetData::getCoolingDil ( int is ) const [inline]

Definition at line 148 of file MaterialBudgetData.h.

References theCoolingDil.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                      {
    return theCoolingDil[is]; }

float MaterialBudgetData::getCoolingDmb ( int is ) const [inline]

Definition at line 130 of file MaterialBudgetData.h.

References theCoolingDmb.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                      {
    return theCoolingDmb[is]; }

float MaterialBudgetData::getCoolingFractionIL ( ) const [inline]

Definition at line 67 of file MaterialBudgetData.h.

References theCoolingFractionIL.

                                     {
    return theCoolingFractionIL; }

float MaterialBudgetData::getCoolingFractionMB ( ) const [inline]

Definition at line 39 of file MaterialBudgetData.h.

References theCoolingFractionMB.

                                     {
    return theCoolingFractionMB; }

float MaterialBudgetData::getCoolingIL ( ) const [inline]

Definition at line 83 of file MaterialBudgetData.h.

References theCoolingIL.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                             {
    return theCoolingIL; }

float MaterialBudgetData::getCoolingMB ( ) const [inline]

Definition at line 53 of file MaterialBudgetData.h.

References theCoolingMB.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                             {
    return theCoolingMB; }

float MaterialBudgetData::getElectronicsDil ( int is ) const [inline]

Definition at line 150 of file MaterialBudgetData.h.

References theElectronicsDil.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                          {
    return theElectronicsDil[is]; }

float MaterialBudgetData::getElectronicsDmb ( int is ) const [inline]

Definition at line 132 of file MaterialBudgetData.h.

References theElectronicsDmb.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                          {
    return theElectronicsDmb[is]; }

float MaterialBudgetData::getElectronicsFractionIL ( ) const [inline]

Definition at line 69 of file MaterialBudgetData.h.

References theElectronicsFractionIL.

                                         {
    return theElectronicsFractionIL; }

float MaterialBudgetData::getElectronicsFractionMB ( ) const [inline]

Definition at line 41 of file MaterialBudgetData.h.

References theElectronicsFractionMB.

                                         {
    return theElectronicsFractionMB; }

float MaterialBudgetData::getElectronicsIL ( ) const [inline]

Definition at line 85 of file MaterialBudgetData.h.

References theElectronicsIL.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                                 {
    return theElectronicsIL; }

float MaterialBudgetData::getElectronicsMB ( ) const [inline]

Definition at line 55 of file MaterialBudgetData.h.

References theElectronicsMB.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                                 {
    return theElectronicsMB; }

float MaterialBudgetData::getEnergy ( ) const [inline]

Definition at line 101 of file MaterialBudgetData.h.

References theEnergy.

Referenced by MaterialBudgetTree::fillEndTrack().

                          {
    return theEnergy; }

float MaterialBudgetData::getEta ( ) const [inline]

Definition at line 92 of file MaterialBudgetData.h.

References theEta.

Referenced by MaterialBudgetEcalHistos::fillEndTrack(), MaterialBudgetHistos::fillEndTrack(), MaterialBudgetTree::fillEndTrack(), MaterialBudgetTrackerHistos::fillEndTrack(), and MaterialBudgetTxt::fillStartTrack().

                       {
    return theEta; }

int MaterialBudgetData::getID ( ) const [inline]

Definition at line 97 of file MaterialBudgetData.h.

References theID.

Referenced by MaterialBudgetTree::fillEndTrack().

                    {
    return theID; }

float MaterialBudgetData::getIntLen ( ) const [inline]

Definition at line 118 of file MaterialBudgetData.h.

References theIntLen.

                          {
    return theIntLen; }

float MaterialBudgetData::getMass ( ) const [inline]

Definition at line 103 of file MaterialBudgetData.h.

References theMass.

Referenced by MaterialBudgetTree::fillEndTrack().

                        {
    return theMass; }

int MaterialBudgetData::getNumberOfSteps ( ) const [inline]

Definition at line 107 of file MaterialBudgetData.h.

References theStepN.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                               {
    return theStepN; }

float MaterialBudgetData::getOtherDil ( int is ) const [inline]

Definition at line 152 of file MaterialBudgetData.h.

References theOtherDil.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                    {
    return theOtherDil[is]; }

float MaterialBudgetData::getOtherDmb ( int is ) const [inline]

Definition at line 134 of file MaterialBudgetData.h.

References theOtherDmb.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                    {
    return theOtherDmb[is]; }

float MaterialBudgetData::getOtherFractionIL ( ) const [inline]

Definition at line 71 of file MaterialBudgetData.h.

References theOtherFractionIL.

                                   {
    return theOtherFractionIL; }

float MaterialBudgetData::getOtherFractionMB ( ) const [inline]

Definition at line 43 of file MaterialBudgetData.h.

References theOtherFractionMB.

                                   {
    return theOtherFractionMB; }

float MaterialBudgetData::getOtherIL ( ) const [inline]

Definition at line 87 of file MaterialBudgetData.h.

References theOtherIL.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                           {
    return theOtherIL; }

float MaterialBudgetData::getOtherMB ( ) const [inline]

Definition at line 57 of file MaterialBudgetData.h.

References theOtherMB.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                           {
    return theOtherMB; }

float MaterialBudgetData::getPhi ( ) const [inline]

Definition at line 94 of file MaterialBudgetData.h.

References thePhi.

Referenced by MaterialBudgetEcalHistos::fillEndTrack(), MaterialBudgetHistos::fillEndTrack(), MaterialBudgetTree::fillEndTrack(), MaterialBudgetTrackerHistos::fillEndTrack(), and MaterialBudgetTxt::fillStartTrack().

                       {
    return thePhi; }

float MaterialBudgetData::getPt ( ) const [inline]

Definition at line 99 of file MaterialBudgetData.h.

References thePt.

Referenced by MaterialBudgetTree::fillEndTrack().

                      {
    return thePt; }

int MaterialBudgetData::getPVcopyNo ( ) const [inline]

Definition at line 114 of file MaterialBudgetData.h.

References thePVcopyNo.

Referenced by MaterialBudgetTxt::fillPerStep().

                          {
    return thePVcopyNo; }

std::string MaterialBudgetData::getPVname ( ) const [inline]

Definition at line 112 of file MaterialBudgetData.h.

References thePVname.

Referenced by MaterialBudgetTxt::fillPerStep().

                              {
    return thePVname; }

float MaterialBudgetData::getRadLen ( ) const [inline]

Definition at line 116 of file MaterialBudgetData.h.

References theRadLen.

Referenced by MaterialBudgetTxt::fillPerStep().

                          {
    return theRadLen; }

float MaterialBudgetData::getSensitiveDil ( int is ) const [inline]

Definition at line 144 of file MaterialBudgetData.h.

References theSensitiveDil.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                        {
    return theSensitiveDil[is]; }

float MaterialBudgetData::getSensitiveDmb ( int is ) const [inline]

Definition at line 126 of file MaterialBudgetData.h.

References theSensitiveDmb.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                        {
    return theSensitiveDmb[is]; }

float MaterialBudgetData::getSensitiveFractionIL ( ) const [inline]

Definition at line 63 of file MaterialBudgetData.h.

References theSensitiveFractionIL.

                                       {
    return theSensitiveFractionIL; }

float MaterialBudgetData::getSensitiveFractionMB ( ) const [inline]

Definition at line 35 of file MaterialBudgetData.h.

References theSensitiveFractionMB.

                                       {
    return theSensitiveFractionMB; }

float MaterialBudgetData::getSensitiveIL ( ) const [inline]

Definition at line 79 of file MaterialBudgetData.h.

References theSensitiveIL.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                               {
    return theSensitiveIL; }

float MaterialBudgetData::getSensitiveMB ( ) const [inline]

Definition at line 49 of file MaterialBudgetData.h.

References theSensitiveMB.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                               {
    return theSensitiveMB; }

float MaterialBudgetData::getStepDil ( int is ) [inline]

Definition at line 139 of file MaterialBudgetData.h.

References theDil.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                             {
    return theDil[is];
  }

float MaterialBudgetData::getStepDmb ( int is ) [inline]

Definition at line 121 of file MaterialBudgetData.h.

References theDmb.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                             {
    return theDmb[is];
  }

float MaterialBudgetData::getStepFinalBeta ( int is ) [inline]

Definition at line 229 of file MaterialBudgetData.h.

References theStepFinalBeta.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theStepFinalBeta[is];
  }

float MaterialBudgetData::getStepFinalEnergy ( int is ) [inline]

Definition at line 217 of file MaterialBudgetData.h.

References theStepFinalEnergy.

Referenced by MaterialBudgetTree::fillEndTrack().

                                    {
    return theStepFinalEnergy[is];
  }

float MaterialBudgetData::getStepFinalEta ( int is ) [inline]

Definition at line 211 of file MaterialBudgetData.h.

References theStepFinalEta.

Referenced by MaterialBudgetTree::fillEndTrack().

                                 {
    return theStepFinalEta[is];
  }

float MaterialBudgetData::getStepFinalGamma ( int is ) [inline]

Definition at line 232 of file MaterialBudgetData.h.

References theStepFinalGamma.

Referenced by MaterialBudgetTree::fillEndTrack().

                                   {
    return theStepFinalGamma[is];
  }

float MaterialBudgetData::getStepFinalMass ( int is ) [inline]

Definition at line 235 of file MaterialBudgetData.h.

References theStepFinalMass.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theStepFinalMass[is];
  }

float MaterialBudgetData::getStepFinalPhi ( int is ) [inline]

Definition at line 214 of file MaterialBudgetData.h.

References theStepFinalPhi.

Referenced by MaterialBudgetTree::fillEndTrack().

                                 {
    return theStepFinalPhi[is];
  }

float MaterialBudgetData::getStepFinalPt ( int is ) [inline]

Definition at line 208 of file MaterialBudgetData.h.

References theStepFinalPt.

Referenced by MaterialBudgetTree::fillEndTrack().

                                {
    return theStepFinalPt[is];
  }

float MaterialBudgetData::getStepFinalPx ( int is ) [inline]

Definition at line 220 of file MaterialBudgetData.h.

References theStepFinalPx.

Referenced by MaterialBudgetTree::fillEndTrack().

                                {
    return theStepFinalPx[is];
  }

float MaterialBudgetData::getStepFinalPy ( int is ) [inline]

Definition at line 223 of file MaterialBudgetData.h.

References theStepFinalPy.

Referenced by MaterialBudgetTree::fillEndTrack().

                                {
    return theStepFinalPy[is];
  }

float MaterialBudgetData::getStepFinalPz ( int is ) [inline]

Definition at line 226 of file MaterialBudgetData.h.

References theStepFinalPz.

Referenced by MaterialBudgetTree::fillEndTrack().

                                {
    return theStepFinalPz[is];
  }

double MaterialBudgetData::getStepFinalX ( int is ) [inline]

Definition at line 166 of file MaterialBudgetData.h.

References theFinalX.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                 {
    return theFinalX[is];
  }

double MaterialBudgetData::getStepFinalY ( int is ) [inline]

Definition at line 169 of file MaterialBudgetData.h.

References theFinalY.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                 {
    return theFinalY[is];
  }

double MaterialBudgetData::getStepFinalZ ( int is ) [inline]

Definition at line 172 of file MaterialBudgetData.h.

References theFinalZ.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                 {
    return theFinalZ[is];
  }

int MaterialBudgetData::getStepID ( int is ) [inline]

Definition at line 175 of file MaterialBudgetData.h.

References theStepID.

Referenced by MaterialBudgetTree::fillEndTrack().

                         {
    return theStepID[is];
  }

float MaterialBudgetData::getStepInitialBeta ( int is ) [inline]

Definition at line 199 of file MaterialBudgetData.h.

References theStepInitialBeta.

Referenced by MaterialBudgetTree::fillEndTrack().

                                    {
    return theStepInitialBeta[is];
  }

float MaterialBudgetData::getStepInitialEnergy ( int is ) [inline]

Definition at line 187 of file MaterialBudgetData.h.

References theStepInitialEnergy.

Referenced by MaterialBudgetTree::fillEndTrack().

                                      {
    return theStepInitialEnergy[is];
  }

float MaterialBudgetData::getStepInitialEta ( int is ) [inline]

Definition at line 181 of file MaterialBudgetData.h.

References theStepInitialEta.

Referenced by MaterialBudgetTree::fillEndTrack().

                                   {
    return theStepInitialEta[is];
  }

float MaterialBudgetData::getStepInitialGamma ( int is ) [inline]

Definition at line 202 of file MaterialBudgetData.h.

References theStepInitialGamma.

Referenced by MaterialBudgetTree::fillEndTrack().

                                     {
    return theStepInitialGamma[is];
  }

float MaterialBudgetData::getStepInitialMass ( int is ) [inline]

Definition at line 205 of file MaterialBudgetData.h.

References theStepInitialMass.

Referenced by MaterialBudgetTree::fillEndTrack().

                                    {
    return theStepInitialMass[is];
  }

float MaterialBudgetData::getStepInitialPhi ( int is ) [inline]

Definition at line 184 of file MaterialBudgetData.h.

References theStepInitialPhi.

Referenced by MaterialBudgetTree::fillEndTrack().

                                   {
    return theStepInitialPhi[is];
  }

float MaterialBudgetData::getStepInitialPt ( int is ) [inline]

Definition at line 178 of file MaterialBudgetData.h.

References theStepInitialPt.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theStepInitialPt[is];
  }

float MaterialBudgetData::getStepInitialPx ( int is ) [inline]

Definition at line 190 of file MaterialBudgetData.h.

References theStepInitialPx.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theStepInitialPx[is];
  }

float MaterialBudgetData::getStepInitialPy ( int is ) [inline]

Definition at line 193 of file MaterialBudgetData.h.

References theStepInitialPy.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theStepInitialPy[is];
  }

float MaterialBudgetData::getStepInitialPz ( int is ) [inline]

Definition at line 196 of file MaterialBudgetData.h.

References theStepInitialPz.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theStepInitialPz[is];
  }

double MaterialBudgetData::getStepInitialX ( int is ) [inline]

Definition at line 157 of file MaterialBudgetData.h.

References theInitialX.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                   {
    return theInitialX[is];
  }

double MaterialBudgetData::getStepInitialY ( int is ) [inline]

Definition at line 160 of file MaterialBudgetData.h.

References theInitialY.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                   {
    return theInitialY[is];
  }

double MaterialBudgetData::getStepInitialZ ( int is ) [inline]

Definition at line 163 of file MaterialBudgetData.h.

References theInitialZ.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                   {
    return theInitialZ[is];
  }

float MaterialBudgetData::getStepMaterialDensity ( int is ) [inline]

Definition at line 284 of file MaterialBudgetData.h.

References theMaterialDensity.

Referenced by MaterialBudgetTree::fillEndTrack().

                                         {
    return theMaterialDensity[is];
  }

int MaterialBudgetData::getStepMaterialID ( int is ) [inline]

Definition at line 272 of file MaterialBudgetData.h.

References theMaterialID.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theMaterialID[is];
  }

float MaterialBudgetData::getStepMaterialLambda0 ( int is ) [inline]

Definition at line 281 of file MaterialBudgetData.h.

References theMaterialLambda0.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                         {
    return theMaterialLambda0[is];
  }

std::string MaterialBudgetData::getStepMaterialName ( int is ) [inline]

Definition at line 275 of file MaterialBudgetData.h.

References theMaterialName.

Referenced by MaterialBudgetTree::fillEndTrack().

                                          {
    return theMaterialName[is];
  }

float MaterialBudgetData::getStepMaterialX0 ( int is ) [inline]

Definition at line 278 of file MaterialBudgetData.h.

References theMaterialX0.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                    {
    return theMaterialX0[is];
  }

int MaterialBudgetData::getStepPostProcess ( int is ) [inline]

Definition at line 241 of file MaterialBudgetData.h.

References theStepPostProcess.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theStepPostProcess[is];
  }

int MaterialBudgetData::getStepPreProcess ( int is ) [inline]

Definition at line 238 of file MaterialBudgetData.h.

References theStepPreProcess.

Referenced by MaterialBudgetTree::fillEndTrack().

                                 {
    return theStepPreProcess[is];
  }

int MaterialBudgetData::getStepVolumeCopy ( int is ) [inline]

Definition at line 251 of file MaterialBudgetData.h.

References theVolumeCopy.

Referenced by MaterialBudgetTree::fillEndTrack().

                                  {
    return theVolumeCopy[is];
  }

int MaterialBudgetData::getStepVolumeID ( int is ) [inline]

Definition at line 245 of file MaterialBudgetData.h.

References theVolumeID.

Referenced by MaterialBudgetTree::fillEndTrack().

                                {
    return theVolumeID[is];
  }

std::string MaterialBudgetData::getStepVolumeName ( int is ) [inline]

Definition at line 248 of file MaterialBudgetData.h.

References theVolumeName.

Referenced by MaterialBudgetTree::fillEndTrack().

                                        {
    return theVolumeName[is];
  }

float MaterialBudgetData::getStepVolumeX ( int is ) [inline]

Definition at line 254 of file MaterialBudgetData.h.

References theVolumeX.

Referenced by MaterialBudgetTree::fillEndTrack().

                                 {
    return theVolumeX[is];
  }

CLHEP::HepLorentzVector MaterialBudgetData::getStepVolumeXaxis ( int is ) [inline]

Definition at line 263 of file MaterialBudgetData.h.

References theVolumeXaxis1, theVolumeXaxis2, and theVolumeXaxis3.

Referenced by MaterialBudgetTree::fillEndTrack().

                                                     {
    return CLHEP::HepLorentzVector(theVolumeXaxis1[is],theVolumeXaxis2[is],theVolumeXaxis3[is]);
  }

float MaterialBudgetData::getStepVolumeY ( int is ) [inline]

Definition at line 257 of file MaterialBudgetData.h.

References theVolumeY.

Referenced by MaterialBudgetTree::fillEndTrack().

                                 {
    return theVolumeY[is];
  }

CLHEP::HepLorentzVector MaterialBudgetData::getStepVolumeYaxis ( int is ) [inline]

Definition at line 266 of file MaterialBudgetData.h.

References theVolumeYaxis1, theVolumeYaxis2, and theVolumeYaxis3.

Referenced by MaterialBudgetTree::fillEndTrack().

                                                     {
    return CLHEP::HepLorentzVector(theVolumeYaxis1[is],theVolumeYaxis2[is],theVolumeYaxis3[is]);
  }

float MaterialBudgetData::getStepVolumeZ ( int is ) [inline]

Definition at line 260 of file MaterialBudgetData.h.

References theVolumeZ.

Referenced by MaterialBudgetTree::fillEndTrack().

                                 {
    return theVolumeZ[is];
  }

CLHEP::HepLorentzVector MaterialBudgetData::getStepVolumeZaxis ( int is ) [inline]

Definition at line 269 of file MaterialBudgetData.h.

References theVolumeZaxis1, theVolumeZaxis2, and theVolumeZaxis3.

Referenced by MaterialBudgetTree::fillEndTrack().

                                                     {
    return CLHEP::HepLorentzVector(theVolumeZaxis1[is],theVolumeZaxis2[is],theVolumeZaxis3[is]);
  }

float MaterialBudgetData::getSupportDil ( int is ) const [inline]

Definition at line 142 of file MaterialBudgetData.h.

References theSupportDil.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                      {
    return theSupportDil[is]; }

float MaterialBudgetData::getSupportDmb ( int is ) const [inline]

Definition at line 124 of file MaterialBudgetData.h.

References theSupportDmb.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                                      {
    return theSupportDmb[is]; }

float MaterialBudgetData::getSupportFractionIL ( ) const [inline]

Definition at line 61 of file MaterialBudgetData.h.

References theSupportFractionIL.

                                     {
    return theSupportFractionIL; }

float MaterialBudgetData::getSupportFractionMB ( ) const [inline]

Definition at line 33 of file MaterialBudgetData.h.

References theSupportFractionMB.

                                     {
    return theSupportFractionMB; }

float MaterialBudgetData::getSupportIL ( ) const [inline]

Definition at line 77 of file MaterialBudgetData.h.

References theSupportIL.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                             {
    return theSupportIL; }

float MaterialBudgetData::getSupportMB ( ) const [inline]

Definition at line 47 of file MaterialBudgetData.h.

References theSupportMB.

Referenced by MaterialBudgetTrackerHistos::fillEndTrack().

                             {
    return theSupportMB; }

float MaterialBudgetData::getTotalIL ( ) const [inline]

Definition at line 75 of file MaterialBudgetData.h.

References theTotalIL.

Referenced by MaterialBudgetTree::fillEndTrack(), and MaterialBudgetTrackerHistos::fillEndTrack().

                           {
    return theTotalIL; }

float MaterialBudgetData::getTotalMB ( ) const [inline]

Definition at line 30 of file MaterialBudgetData.h.

References theTotalMB.

Referenced by MaterialBudgetTxt::fillEndTrack(), MaterialBudgetEcalHistos::fillEndTrack(), MaterialBudgetHistos::fillEndTrack(), MaterialBudgetTree::fillEndTrack(), MaterialBudgetTrackerHistos::fillEndTrack(), and MaterialBudgetTxt::fillPerStep().

                           {
    return theTotalMB; }

float MaterialBudgetData::getTrkLen ( ) const [inline]

Definition at line 110 of file MaterialBudgetData.h.

References theTrkLen.

Referenced by MaterialBudgetTxt::fillPerStep().

                          {
    return theTrkLen; }

void MaterialBudgetData::SetAllStepsToTree ( )

!!WARNING: this number is also hardcoded when booking the tree

Definition at line 23 of file MaterialBudgetData.cc.

References allStepsToTree, MAXNUMBERSTEPS, theAirDil, theAirDmb, theCablesDil, theCablesDmb, theCoolingDil, theCoolingDmb, theDil, theDmb, theElectronicsDil, theElectronicsDmb, theFinalX, theFinalY, theFinalZ, theInitialX, theInitialY, theInitialZ, theMaterialDensity, theMaterialID, theMaterialLambda0, theMaterialName, theMaterialX0, theOtherDil, theOtherDmb, theSensitiveDil, theSensitiveDmb, theStepFinalBeta, theStepFinalEnergy, theStepFinalEta, theStepFinalGamma, theStepFinalMass, theStepFinalPhi, theStepFinalPt, theStepFinalPx, theStepFinalPy, theStepFinalPz, theStepID, theStepInitialBeta, theStepInitialEnergy, theStepInitialEta, theStepInitialGamma, theStepInitialMass, theStepInitialPhi, theStepInitialPt, theStepInitialPx, theStepInitialPy, theStepInitialPz, theStepPostProcess, theStepPreProcess, theSupportDil, theSupportDmb, theVolumeCopy, theVolumeID, theVolumeName, theVolumeX, theVolumeXaxis1, theVolumeXaxis2, theVolumeXaxis3, theVolumeY, theVolumeYaxis1, theVolumeYaxis2, theVolumeYaxis3, theVolumeZ, theVolumeZaxis1, theVolumeZaxis2, and theVolumeZaxis3.

Referenced by MaterialBudgetAction::MaterialBudgetAction().

{
  allStepsToTree = true;
  MAXNUMBERSTEPS = 0;
  MAXNUMBERSTEPS = 10000; 
  theDmb = new float[MAXNUMBERSTEPS];
  theDil = new float[MAXNUMBERSTEPS];
  // rr
  theSupportDmb     = new float[MAXNUMBERSTEPS];
  theSensitiveDmb   = new float[MAXNUMBERSTEPS];
  theCablesDmb      = new float[MAXNUMBERSTEPS];
  theCoolingDmb     = new float[MAXNUMBERSTEPS];
  theElectronicsDmb = new float[MAXNUMBERSTEPS];
  theOtherDmb       = new float[MAXNUMBERSTEPS];
  theAirDmb         = new float[MAXNUMBERSTEPS];
  theSupportDil     = new float[MAXNUMBERSTEPS];
  theSensitiveDil   = new float[MAXNUMBERSTEPS];
  theCablesDil      = new float[MAXNUMBERSTEPS];
  theCoolingDil     = new float[MAXNUMBERSTEPS];
  theElectronicsDil = new float[MAXNUMBERSTEPS];
  theOtherDil       = new float[MAXNUMBERSTEPS];
  theAirDil         = new float[MAXNUMBERSTEPS];
  // rr
  theInitialX = new double[MAXNUMBERSTEPS];
  theInitialY = new double[MAXNUMBERSTEPS];
  theInitialZ = new double[MAXNUMBERSTEPS];
  theFinalX   = new double[MAXNUMBERSTEPS];
  theFinalY   = new double[MAXNUMBERSTEPS];
  theFinalZ   = new double[MAXNUMBERSTEPS];
  // rr
  theVolumeID     = new int[MAXNUMBERSTEPS];
  theVolumeName   = new std::string[MAXNUMBERSTEPS];
  theVolumeCopy   = new int[MAXNUMBERSTEPS];
  theVolumeX      = new float[MAXNUMBERSTEPS];
  theVolumeY      = new float[MAXNUMBERSTEPS];
  theVolumeZ      = new float[MAXNUMBERSTEPS];
  theVolumeXaxis1 = new float[MAXNUMBERSTEPS];
  theVolumeXaxis2 = new float[MAXNUMBERSTEPS];
  theVolumeXaxis3 = new float[MAXNUMBERSTEPS];
  theVolumeYaxis1 = new float[MAXNUMBERSTEPS];
  theVolumeYaxis2 = new float[MAXNUMBERSTEPS];
  theVolumeYaxis3 = new float[MAXNUMBERSTEPS];
  theVolumeZaxis1 = new float[MAXNUMBERSTEPS];
  theVolumeZaxis2 = new float[MAXNUMBERSTEPS];
  theVolumeZaxis3 = new float[MAXNUMBERSTEPS];
  theMaterialID      = new int[MAXNUMBERSTEPS];
  theMaterialName    = new std::string[MAXNUMBERSTEPS];
  theMaterialX0      = new float[MAXNUMBERSTEPS];
  theMaterialLambda0 = new float[MAXNUMBERSTEPS];
  theMaterialDensity = new float[MAXNUMBERSTEPS];
  theStepID             = new int[MAXNUMBERSTEPS];
  theStepInitialPt      = new float[MAXNUMBERSTEPS];
  theStepInitialEta     = new float[MAXNUMBERSTEPS];
  theStepInitialPhi     = new float[MAXNUMBERSTEPS];
  theStepInitialEnergy  = new float[MAXNUMBERSTEPS];
  theStepInitialPx      = new float[MAXNUMBERSTEPS];
  theStepInitialPy      = new float[MAXNUMBERSTEPS];
  theStepInitialPz      = new float[MAXNUMBERSTEPS];
  theStepInitialBeta    = new float[MAXNUMBERSTEPS];
  theStepInitialGamma   = new float[MAXNUMBERSTEPS];
  theStepInitialMass    = new float[MAXNUMBERSTEPS];
  theStepFinalPt        = new float[MAXNUMBERSTEPS];
  theStepFinalEta       = new float[MAXNUMBERSTEPS];
  theStepFinalPhi       = new float[MAXNUMBERSTEPS];
  theStepFinalEnergy    = new float[MAXNUMBERSTEPS];
  theStepFinalPx        = new float[MAXNUMBERSTEPS];
  theStepFinalPy        = new float[MAXNUMBERSTEPS];
  theStepFinalPz        = new float[MAXNUMBERSTEPS];
  theStepFinalBeta      = new float[MAXNUMBERSTEPS];
  theStepFinalGamma     = new float[MAXNUMBERSTEPS];
  theStepFinalMass      = new float[MAXNUMBERSTEPS];
  theStepPreProcess     = new int[MAXNUMBERSTEPS];
  theStepPostProcess    = new int[MAXNUMBERSTEPS];
  // rr
}