EMShower.h

Go to the documentation of this file.

#ifndef EMShower_H
#define EMShower_H

#include "FastSimulation/Particle/interface/RawParticle.h"

//Famos Headers
#include "FastSimulation/ShowerDevelopment/interface/EMECALShowerParametrization.h"
#include "FastSimulation/ShowerDevelopment/interface/RadialInterval.h"
#include "CLHEP/GenericFunctions/IncompleteGamma.hh"
#include "DataFormats/Math/interface/Vector3D.h"

#include <vector>

class EcalHitMaker;
class PreshowerHitMaker;
class HcalHitMaker;
class GammaDistributionGenerator;
class RandomEngine;
class GammaFunctionGenerator;

class EMShower 
{

  typedef math::XYZVector XYZPoint;

  typedef std::pair<XYZPoint,double> Spot;
  typedef std::pair<unsigned int, double> Step;
  typedef std::vector<Step> Steps;
  typedef Steps::const_iterator step_iterator;

 public:

  EMShower(const RandomEngine* engine,
       GammaFunctionGenerator* gamma,
       EMECALShowerParametrization* const myParam,
       std::vector<const RawParticle*>* const myPart,
       EcalHitMaker  * const myGrid=NULL,
       PreshowerHitMaker * const myPreshower=NULL);

  virtual ~EMShower(){;}

  void prepareSteps();

  void compute();

  //  inline double getMeanDepth() const {return globalMeanDepth;};  

  inline double getMaximumOfShower() const {return globalMaximum;}

  void setGrid(EcalHitMaker * const myGrid) { theGrid=myGrid;}

  void setPreshower(PreshowerHitMaker * const myPresh ) ;

  void setHcal(HcalHitMaker * const myHcal);

 private:

  // The longitudinal development ersatzt.
  double gam(double x, double a) const;

  // Energy deposited in the layer t-dt-> t, in units of E0 (initial energy)
  double deposit(double t, double a, double b, double dt);

  // Energy deposited between 0 and t, in units of E0 (initial energy)
  double deposit(double a, double b, double t);

  // Set the intervals for the radial development
  void setIntervals(unsigned icomp,RadialInterval& rad);
  
  // The parametrization
  EMECALShowerParametrization* const theParam;

  // The Calorimeter properties
  const ECALProperties* theECAL;
  const HCALProperties* theHCAL;
  const PreshowerLayer1Properties* theLayer1;
  const PreshowerLayer2Properties* theLayer2;

  // The incident particle(s)
  std::vector<const RawParticle*>* const thePart;
  unsigned int nPart;

  // The basic quantities for the shower development.
  std::vector<double> theNumberOfSpots;
  std::vector<double> Etot;
  std::vector<double> E; 
  std::vector<double> photos;
  std::vector<double> T;
  std::vector<double> a; 
  std::vector<double> b;
  std::vector<double> Ti; 
  std::vector<double> TSpot;
  std::vector<double> aSpot; 
  std::vector<double> bSpot;

  // F.B : Use the maximum of the shower rather the center of gravity 
  //  std::vector<double> meanDepth;  
  //  double globalMeanDepth;
  std::vector<double> maximumOfShower;
  std::vector<std::vector<double> >depositedEnergy;
  std::vector<double> meanDepth;
  double innerDepth, outerDepth;
  
  double globalMaximum;

  double totalEnergy;

  // The steps for the longitudinal development
  Steps steps; 
  unsigned nSteps;
  bool stepsCalculated;

  // The crystal grid
  EcalHitMaker * theGrid;

  // The preshower 
  PreshowerHitMaker * thePreshower;

  // The HCAL hitmaker
  HcalHitMaker * theHcalHitMaker;

  // Is there a preshower ? 
  bool hasPreshower;
  // Histos
  //  Histos* myHistos;

  Genfun::IncompleteGamma myIncompleteGamma;

  // Random engine
  const RandomEngine* random;

  // integer gamma function generator
  GammaFunctionGenerator * myGammaGenerator;

  
  
};

#endif