MaterialEffectsSimulator.h

Go to the documentation of this file.

#ifndef MATERIALEFFECTSSIMULATOR_H
#define MATERIALEFFECTSSIMULATOR_H

#include "DataFormats/GeometryVector/interface/GlobalVector.h"

#include "FastSimulation/ParticlePropagator/interface/ParticlePropagator.h"

#include <vector>

class RandomEngineAndDistribution;

class MaterialEffectsSimulator
{
 public:

  typedef std::vector<RawParticle>::const_iterator RHEP_const_iter;

  // Constructor : default values are for Silicon
  MaterialEffectsSimulator(double A = 28.0855,
               double Z = 14.0000,
               double density = 2.329,
               double radLen = 9.360);

  virtual ~MaterialEffectsSimulator();


  inline double theA() const { return A; }
  inline double theZ() const { return Z; }
  inline double rho() const { return density; }
  inline double radLenIncm() const { return radLen; }
  inline double excitE() const { return 12.5E-9*theZ(); }
  inline double eMass() const { return 0.000510998902; }

  void updateState(ParticlePropagator& myTrack, double radlen, RandomEngineAndDistribution const*);
  
  inline RHEP_const_iter beginDaughters() const {return _theUpdatedState.begin();}

  inline RHEP_const_iter endDaughters() const {return _theUpdatedState.end();}

  inline unsigned nDaughters() const {return _theUpdatedState.size();}

  inline void setNormalVector(const GlobalVector& normal) { theNormalVector = normal; }
  
  XYZVector orthogonal(const XYZVector&) const; 

  inline int closestDaughterId() { return theClosestChargedDaughterId; } 

  virtual void save() {};

 private:

  virtual void compute(ParticlePropagator& Particle, RandomEngineAndDistribution const*) = 0;

  inline double radiationLength() const {return radLengths;}

 protected:

  std::vector<RawParticle> _theUpdatedState;

  double radLengths;

  // Material properties
  double A;
  double Z;
  double density;
  double radLen;

  GlobalVector theNormalVector;

  int theClosestChargedDaughterId;

};

#endif