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