EMShower.h

Go to the documentation of this file.

#ifndef EMShower_H
#define EMShower_H
 
#include "CommonTools/BaseParticlePropagator/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 RandomEngineAndDistribution;
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(RandomEngineAndDistribution const* engine,
           GammaFunctionGenerator* gamma,
           EMECALShowerParametrization* const myParam,
           std::vector<const RawParticle*>* const myPart,
           EcalHitMaker* const myGrid = nullptr,
           PreshowerHitMaker* const myPreshower = nullptr,
           bool bFixedLength = false);
 
  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 RandomEngineAndDistribution* random;
 
  // integer gamma function generator
  GammaFunctionGenerator* myGammaGenerator;
 
  bool bFixedLength_;
};
 
#endif