GflashShowino Class Reference

#include <GflashShowino.h>

Public Member Functions
void	addEnergyDeposited (double energy)
double	getDepth ()
double	getEnergy ()
double	getEnergyDeposited ()
double	getGlobalTime ()
GflashTrajectory *	getHelix ()
double	getPathLength ()
double	getPathLengthAtShower ()
double	getPathLengthOnEcal ()
Gflash3Vector &	getPosition ()
Gflash3Vector &	getPositionAtShower ()
int	getShowerType ()
double	getStepLengthToHcal ()
double	getStepLengthToOut ()
	GflashShowino ()
void	initialize (int showerType, double energy, double globalTime, double charge, Gflash3Vector &position, Gflash3Vector &momentum, double magneticField)
void	setGlobalTime (double globalTime)
void	setPathLength (double pathLength)
void	setPosition (Gflash3Vector position)
void	updateShowino (double deltaStep)
	~GflashShowino ()

Private Member Functions
int	convertShowerType (int fastSimShowerType, const Gflash3Vector &pos)
void	evaluateLengths ()
Gflash3Vector &	simulateFirstInteractionPoint (int showType, Gflash3Vector &pos)

Private Attributes
double	theEnergy
double	theEnergyDeposited
double	theGlobalTime
GflashTrajectory *	theHelix
double	thePathLength
double	thePathLengthAtShower
double	thePathLengthOnEcal
Gflash3Vector	thePosition
Gflash3Vector	thePositionAtShower
int	theShowerType
double	theStepLengthToHcal
double	theStepLengthToOut

Detailed Description

Definition at line 8 of file GflashShowino.h.

Constructor & Destructor Documentation

GflashShowino::GflashShowino

(

)

Definition at line 4 of file GflashShowino.cc.

References theHelix.

  : theShowerType(-1), theEnergy(0), thePositionAtShower(0,0,0), thePathLengthAtShower(0), 
    thePathLengthOnEcal(0), theStepLengthToHcal(0), theStepLengthToOut(0),  
    thePathLength(0), theGlobalTime(0), thePosition(0,0,0), theEnergyDeposited(0)
{
  theHelix =  new GflashTrajectory;;
}

GflashShowino::~GflashShowino

(

)

Definition at line 12 of file GflashShowino.cc.

References theHelix.

{
  delete theHelix;
}

Member Function Documentation

void GflashShowino::addEnergyDeposited ( double  energy )

inline

Definition at line 41 of file GflashShowino.h.

References relval_parameters_module::energy, and theEnergyDeposited.

Referenced by GflashHadronShowerProfile::hadronicParameterization().

{ theEnergyDeposited += energy; }

int GflashShowino::convertShowerType ( int  fastSimShowerType, const Gflash3Vector &  pos  )

private

void GflashShowino::evaluateLengths ( )

private

Definition at line 61 of file GflashShowino.cc.

References eta(), Gflash::EtaMax, GflashTrajectory::getPathLengthAtRhoEquals(), GflashTrajectory::getPathLengthAtZ(), Gflash::kENCA, Gflash::kESPM, Gflash::kHB, Gflash::kHE, min, Gflash::RFrontCrystalEB, Gflash::Rmax, Gflash::Rmin, theHelix, thePathLength, thePathLengthAtShower, thePathLengthOnEcal, thePosition, theStepLengthToHcal, theStepLengthToOut, Gflash::ZFrontCrystalEE, Gflash::Zmax, and Gflash::Zmin.

Referenced by initialize().

                                    {
  //thePathLengthAtShower: path Length from the origin to the shower starting point in cm
  //theStepLengthToOut: the total path length from the starting point of
  //                    shower to the maximum distance inside paramerized envelopes
  double eta = thePosition.getEta();
  
  if (std::fabs(eta) < Gflash::EtaMax[Gflash::kESPM]) {
    thePathLengthOnEcal   = theHelix->getPathLengthAtRhoEquals(Gflash::RFrontCrystalEB);
    thePathLengthAtShower = theHelix->getPathLengthAtRhoEquals(thePosition.getRho());
    theStepLengthToOut    = std::min(300.,theHelix->getPathLengthAtRhoEquals(Gflash::Rmax[Gflash::kHB]) - thePathLengthAtShower);
    theStepLengthToHcal   = theHelix->getPathLengthAtRhoEquals(Gflash::Rmin[Gflash::kHB]) - thePathLengthAtShower;
  }
  else if (std::fabs(eta) < Gflash::EtaMax[Gflash::kENCA] ) {
    double zsign = (eta > 0) ? 1.0 : -1.0 ; 
    thePathLengthOnEcal   = theHelix->getPathLengthAtZ(zsign*Gflash::ZFrontCrystalEE);
    thePathLengthAtShower = theHelix->getPathLengthAtZ(thePosition.getZ());
    theStepLengthToOut    = std::min(300.,theHelix->getPathLengthAtZ(zsign*Gflash::Zmax[Gflash::kHE]) - thePathLengthAtShower);
    theStepLengthToHcal   = theHelix->getPathLengthAtZ(zsign*Gflash::Zmin[Gflash::kHE]) - thePathLengthAtShower;
  }
  else { 
    //@@@extend for HF later
    theStepLengthToOut = 200.0;
  }

  thePathLength = thePathLengthAtShower;

}

double GflashShowino::getDepth ( void  )

inline

Definition at line 36 of file GflashShowino.h.

References thePathLength, and thePathLengthAtShower.

Referenced by GflashHadronShowerProfile::longitudinalProfile().

{ return (thePathLength-thePathLengthAtShower); }

double GflashShowino::getEnergy ( )

inline

Definition at line 24 of file GflashShowino.h.

References theEnergy.

Referenced by GflashHadronShowerProfile::getNumberOfSpots(), GflashHadronShowerProfile::hadronicParameterization(), GflashPiKShowerProfile::loadParameters(), GflashKaonPlusShowerProfile::loadParameters(), GflashKaonMinusShowerProfile::loadParameters(), GflashProtonShowerProfile::loadParameters(), GflashAntiProtonShowerProfile::loadParameters(), and GflashEMShowerProfile::parameterization().

{ return theEnergy; }

double GflashShowino::getEnergyDeposited ( )

inline

Definition at line 35 of file GflashShowino.h.

References theEnergyDeposited.

Referenced by GflashHadronShowerProfile::hadronicParameterization().

{ return theEnergyDeposited ; }

double GflashShowino::getGlobalTime ( )

inline

Definition at line 32 of file GflashShowino.h.

References theGlobalTime.

Referenced by GflashHadronShowerProfile::hadronicParameterization(), and GflashEMShowerProfile::parameterization().

{ return theGlobalTime; }

GflashTrajectory* GflashShowino::getHelix ( )

inline

Definition at line 30 of file GflashShowino.h.

References theHelix.

Referenced by GflashEMShowerProfile::getDistanceToOut(), GflashHadronShowerProfile::hadronicParameterization(), CalorimetryManager::HDShowerSimulation(), GflashHadronShowerProfile::hoProfile(), and GflashEMShowerProfile::parameterization().

{ return theHelix; }

double GflashShowino::getPathLength ( )

inline

Definition at line 33 of file GflashShowino.h.

References thePathLength.

Referenced by GflashHadronShowerProfile::hadronicParameterization().

{ return thePathLength; }

double GflashShowino::getPathLengthAtShower ( )

inline

Definition at line 26 of file GflashShowino.h.

References thePathLengthAtShower.

Referenced by GflashEMShowerProfile::getDistanceToOut(), GflashHadronShowerProfile::hadronicParameterization(), CalorimetryManager::HDShowerSimulation(), and GflashEMShowerProfile::parameterization().

{ return thePathLengthAtShower; }

double GflashShowino::getPathLengthOnEcal ( )

inline

Definition at line 25 of file GflashShowino.h.

References thePathLengthOnEcal.

Referenced by CalorimetryManager::HDShowerSimulation().

{ return thePathLengthOnEcal; }

Gflash3Vector& GflashShowino::getPosition ( )

inline

Definition at line 34 of file GflashShowino.h.

References thePosition.

Referenced by GflashEMShowerProfile::getDistanceToOut(), GflashHadronShowerProfile::hadronicParameterization(), and GflashHadronShowerProfile::longitudinalProfile().

{ return thePosition; }

Gflash3Vector& GflashShowino::getPositionAtShower ( )

inline

Definition at line 27 of file GflashShowino.h.

References thePositionAtShower.

Referenced by GflashHadronShowerProfile::hadronicParameterization(), GflashProtonShowerProfile::loadParameters(), GflashPiKShowerProfile::loadParameters(), GflashKaonPlusShowerProfile::loadParameters(), GflashKaonMinusShowerProfile::loadParameters(), GflashAntiProtonShowerProfile::loadParameters(), and GflashEMShowerProfile::parameterization().

{ return thePositionAtShower; }

int GflashShowino::getShowerType ( )

inline

Definition at line 23 of file GflashShowino.h.

References theShowerType.

Referenced by GflashHadronShowerProfile::getNumberOfSpots(), GflashProtonShowerProfile::loadParameters(), GflashPiKShowerProfile::loadParameters(), GflashKaonPlusShowerProfile::loadParameters(), GflashKaonMinusShowerProfile::loadParameters(), GflashAntiProtonShowerProfile::loadParameters(), and GflashHadronShowerProfile::longitudinalProfile().

{ return theShowerType; }

double GflashShowino::getStepLengthToHcal ( )

inline

Definition at line 28 of file GflashShowino.h.

References theStepLengthToHcal.

Referenced by GflashHadronShowerProfile::longitudinalProfile().

{ return theStepLengthToHcal; }

double GflashShowino::getStepLengthToOut ( )

inline

Definition at line 29 of file GflashShowino.h.

References theStepLengthToOut.

Referenced by GflashHadronShowerProfile::hadronicParameterization().

{ return theStepLengthToOut; }

void GflashShowino::initialize	(	int	showerType,
		double	energy,
		double	globalTime,
		double	charge,
		Gflash3Vector &	position,
		Gflash3Vector &	momentum,
		double	magneticField
	)

Definition at line 17 of file GflashShowino.cc.

References relval_parameters_module::energy, evaluateLengths(), Gflash::findShowerType(), GflashTrajectory::initializeTrajectory(), position, simulateFirstInteractionPoint(), theEnergy, theEnergyDeposited, theGlobalTime, theHelix, thePosition, thePositionAtShower, and theShowerType.

Referenced by GflashEMShowerProfile::initialize(), and GflashHadronShowerProfile::initialize().

                                         {

  theEnergy = energy;
  theGlobalTime = globalTime;
  theEnergyDeposited = 0.0;

  // inside the magnetic field (tesla unit);
  theHelix->initializeTrajectory(momentum,position,charge,magneticField);

  if(showerType<100) {
    thePositionAtShower = position;  
    thePosition = thePositionAtShower;
    theShowerType = showerType;

  }
  else {
    //this input is from FastSimulation
    //1. simulate the shower starting position
    thePositionAtShower = simulateFirstInteractionPoint(showerType,position);
    thePosition = thePositionAtShower;

    //2. find shower type depending on where is the shower starting point
    theShowerType = Gflash::findShowerType(thePositionAtShower);

  }  

  evaluateLengths();

}

void GflashShowino::setGlobalTime ( double  globalTime )

inline

Definition at line 38 of file GflashShowino.h.

References theGlobalTime.

{ theGlobalTime = globalTime; }

void GflashShowino::setPathLength ( double  pathLength )

inline

Definition at line 39 of file GflashShowino.h.

References thePathLength.

Referenced by GflashHadronShowerProfile::hadronicParameterization().

{ thePathLength = pathLength; }

void GflashShowino::setPosition ( Gflash3Vector  position )

inline

Definition at line 40 of file GflashShowino.h.

References position, and thePosition.

{ thePosition = position; }

Gflash3Vector & GflashShowino::simulateFirstInteractionPoint ( int  showType, Gflash3Vector &  pos  )

private

fraction before the crystal, but inside Ecal

Definition at line 90 of file GflashShowino.cc.

References GflashTrajectory::getGflashTrajectoryPoint(), GflashTrajectory::getPathLengthAtRhoEquals(), GflashTrajectory::getPathLengthAtZ(), GflashTrajectoryPoint::getPosition(), Gflash::kHB, Gflash::kHE, Gflash::LengthCrystalEB, Gflash::LengthCrystalEE, funct::log(), Gflash::RFrontCrystalEB, Gflash::Rmax, Gflash::Rmin, funct::sin(), theEnergy, theHelix, thePathLengthOnEcal, Gflash::ZFrontCrystalEE, Gflash::Zmax, and Gflash::Zmin.

Referenced by initialize().

                                                                                                          {

  //determine the shower starting point (ssp):
  //the position at the entrance + the mean free path till the inelastic interaction inside calo

  double depthAtShower = 0.0;

  //set thePathLengthOnEcal, the pathLength at the reference (r=123.8 for barrel and z=304.5 for endcap)
    
  //effective interaction length fitter to ssp from Geant4
  double effectiveLambda = 0.0;
  if(theEnergy > 0.0 &&  theEnergy < 15) {
    effectiveLambda = 24.6+2.6*std::tanh(3.0*(std::log(theEnergy)-1.43));
  }
  else {
    effectiveLambda = 28.4+1.20*std::tanh(1.5*(std::log(theEnergy)-4.3));
  }
  //fraction before the crystal, but inside Ecal
  //  double frac_ssp1 = 1.5196e-01+1.3300e-01*tanh(-4.6971e-01*(std::log(theEnergy)+2.4162e+00));
  //fraction after the crystal, but before Hcal
  double frac_ssp2 = 2.8310e+00+2.6766e+00*tanh(-4.8068e-01*(std::log(theEnergy)+3.4857e+00));

  if(fastSimShowerType == 100 ) { //fastTrack.onEcal() == 1

    //    double rhoTemp = Gflash::ROffCrystalEB + Gflash::LengthCrystalEB*std::sin(position.getTheta());
    double rhoTemp = Gflash::LengthCrystalEB*std::sin(position.getTheta());
    thePathLengthOnEcal  = theHelix->getPathLengthAtRhoEquals(Gflash::RFrontCrystalEB);
    double pathLengthAt2 = theHelix->getPathLengthAtRhoEquals(Gflash::RFrontCrystalEB + rhoTemp );
    double pathLengthAt3 = theHelix->getPathLengthAtRhoEquals(Gflash::Rmin[Gflash::kHB]);

    /*
    if(CLHEP::HepUniformRand() < frac_ssp1 ) {
      depthAtShower = (pathLengthAt1-thePathLengthOnEcal)*CLHEP::HepUniformRand();
    }
    else {
    */
    //inside the crystal
    //      depthAtShower = (pathLengthAt1-thePathLengthOnEcal) - effectiveLambda*log(CLHEP::HepUniformRand());
      depthAtShower = - effectiveLambda*log(CLHEP::HepUniformRand());
      //after the crystal
      if(depthAtShower > (pathLengthAt2 - thePathLengthOnEcal) ) {
    //before Hcal
    if(CLHEP::HepUniformRand() < frac_ssp2 ) {
      depthAtShower = (pathLengthAt2 - thePathLengthOnEcal) + (pathLengthAt3-pathLengthAt2)*CLHEP::HepUniformRand();
    }
    //inside Hcal
    else {
      depthAtShower = (pathLengthAt3 - thePathLengthOnEcal) - effectiveLambda*log(CLHEP::HepUniformRand());
      //check whether the shower starts beyond HB
      double pathLengthAt4 = theHelix->getPathLengthAtRhoEquals(Gflash::Rmax[Gflash::kHB]);
      if(depthAtShower > (pathLengthAt4 - thePathLengthOnEcal)) {
        depthAtShower = (pathLengthAt4 - pathLengthAt3)*CLHEP::HepUniformRand();
      }
    }
      }
      //    }
  }
  else if(fastSimShowerType == 101 ) { //fastTrack.onEcal() == 2 

    double zTemp = Gflash::LengthCrystalEE;
    double zsign = (position.getEta() > 0) ? 1.0 : -1.0 ; 

    thePathLengthOnEcal  = theHelix->getPathLengthAtZ(zsign*Gflash::ZFrontCrystalEE);
    double pathLengthAt2 = theHelix->getPathLengthAtZ(zsign*(Gflash::ZFrontCrystalEE + zTemp));
    double pathLengthAt3 = theHelix->getPathLengthAtZ(zsign*Gflash::Zmin[Gflash::kHE]);
      
    /*
    if(CLHEP::HepUniformRand() <  frac_ssp1 ) {
      depthAtShower = (pathLengthAt1-thePathLengthOnEcal)*CLHEP::HepUniformRand();
    }
    else {
    */
    //      depthAtShower = (pathLengthAt1-thePathLengthOnEcal)-effectiveLambda*std::log(CLHEP::HepUniformRand());
      depthAtShower = -effectiveLambda*std::log(CLHEP::HepUniformRand());

      if(depthAtShower > (pathLengthAt2 - thePathLengthOnEcal) ) {
    if(CLHEP::HepUniformRand() < frac_ssp2 ) {
      depthAtShower = (pathLengthAt2 - thePathLengthOnEcal) +(pathLengthAt3 - pathLengthAt2)*CLHEP::HepUniformRand();
    }
    else {
      depthAtShower = (pathLengthAt3 - thePathLengthOnEcal)-effectiveLambda*std::log(CLHEP::HepUniformRand());
      //shower starts beyond HE
      double pathLengthAt4 = theHelix->getPathLengthAtZ(zsign*Gflash::Zmax[Gflash::kHE]);
      if(depthAtShower > (pathLengthAt4 - thePathLengthOnEcal)) {
        depthAtShower = (pathLengthAt4 - pathLengthAt3)*CLHEP::HepUniformRand();
      }
    }
      }
      //    }
  }
  else {
    depthAtShower = 0.0; 
  }

  double pathLength  = thePathLengthOnEcal + depthAtShower;
  GflashTrajectoryPoint trajectoryPoint;
  theHelix->getGflashTrajectoryPoint(trajectoryPoint,pathLength);

  //return the initial showino position at the shower starting position
  return trajectoryPoint.getPosition();

}

void GflashShowino::updateShowino

(

double

deltaStep

)

Definition at line 49 of file GflashShowino.cc.

References GflashTrajectory::getGflashTrajectoryPoint(), GflashTrajectoryPoint::getPosition(), theEnergy, theGlobalTime, theHelix, thePathLength, and thePosition.

Referenced by GflashHadronShowerProfile::hadronicParameterization().

{
  thePathLength += deltaStep;
  //trajectory point of showino along the shower depth at the pathLength
  GflashTrajectoryPoint trajectoryShowino;
  theHelix->getGflashTrajectoryPoint(trajectoryShowino,thePathLength);

  thePosition = trajectoryShowino.getPosition();

  theGlobalTime +=  (theEnergy/100.0)*deltaStep/30.0; //@@@calculate exact time change in nsec
}

Member Data Documentation

double GflashShowino::theEnergy

private

Definition at line 52 of file GflashShowino.h.

Referenced by getEnergy(), initialize(), simulateFirstInteractionPoint(), and updateShowino().

double GflashShowino::theEnergyDeposited

private

Definition at line 63 of file GflashShowino.h.

Referenced by addEnergyDeposited(), getEnergyDeposited(), and initialize().

double GflashShowino::theGlobalTime

private

Definition at line 61 of file GflashShowino.h.

Referenced by getGlobalTime(), initialize(), setGlobalTime(), and updateShowino().

GflashTrajectory* GflashShowino::theHelix

private

Definition at line 65 of file GflashShowino.h.

Referenced by evaluateLengths(), getHelix(), GflashShowino(), initialize(), simulateFirstInteractionPoint(), updateShowino(), and ~GflashShowino().

double GflashShowino::thePathLength

private

Definition at line 60 of file GflashShowino.h.

Referenced by evaluateLengths(), getDepth(), getPathLength(), setPathLength(), and updateShowino().

double GflashShowino::thePathLengthAtShower

private

Definition at line 54 of file GflashShowino.h.

Referenced by evaluateLengths(), getDepth(), and getPathLengthAtShower().

double GflashShowino::thePathLengthOnEcal

private

Definition at line 55 of file GflashShowino.h.

Referenced by evaluateLengths(), getPathLengthOnEcal(), and simulateFirstInteractionPoint().

Gflash3Vector GflashShowino::thePosition

private

Definition at line 62 of file GflashShowino.h.

Referenced by evaluateLengths(), getPosition(), initialize(), setPosition(), and updateShowino().

Gflash3Vector GflashShowino::thePositionAtShower

private

Definition at line 53 of file GflashShowino.h.

Referenced by getPositionAtShower(), and initialize().

int GflashShowino::theShowerType

private

Definition at line 51 of file GflashShowino.h.

Referenced by getShowerType(), and initialize().

double GflashShowino::theStepLengthToHcal

private

Definition at line 56 of file GflashShowino.h.

Referenced by evaluateLengths(), and getStepLengthToHcal().

double GflashShowino::theStepLengthToOut

private

Definition at line 57 of file GflashShowino.h.

Referenced by evaluateLengths(), and getStepLengthToOut().