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