29 typedef std::pair<XYZPoint,double>
Spot;
30 typedef std::pair<unsigned int, double>
Step;
31 typedef std::vector<Step>
Steps;
49 double gam(
double x,
double a)
const {
return pow(x,a-1.)*
exp(-x); }
55 double fsq = factor * factor;
56 return ((fsq + 1.)/fsq) * r * r / (r*r + R*
R) ;
61 int indexFinder(
double x,
const std::vector<double> & Fhist);
std::vector< double > lamdepth
std::vector< double > eStep
std::vector< double > rlamStep
Steps::const_iterator step_iterator
std::vector< int > nspots
std::vector< double > lamstep
int indexFinder(double x, const std::vector< double > &Fhist)
std::pair< unsigned int, double > Step
std::vector< int > detector
const ECALProperties * theECALproperties
bool compute()
Compute the shower longitudinal and lateral development.
double gam(double x, double a) const
HcalHitMaker * theHcalHitMaker
std::vector< double > x0curr
double transProb(double factor, double R, double r)
std::vector< double > lamcurr
const HCALProperties * theHCALproperties
XYZVectorD XYZVector
spatial vector with cartesian internal representation
void makeSteps(int nsteps)
std::vector< double > lamtotal
HFShower(const RandomEngine *engine, HDShowerParametrization *myParam, EcalHitMaker *myGrid, HcalHitMaker *myHcalHitMaker, int onECAL, double epart)
std::pair< XYZPoint, double > Spot
std::vector< Step > Steps
const RandomEngine * random
HDShowerParametrization * theParam
Power< A, B >::type pow(const A &a, const B &b)
std::vector< double > x0depth