MuonBremsstrahlungSimulator.h

Go to the documentation of this file.

#ifndef MUONBREMSSTRAHLUNGSIMULATOR_H
#define MUONBREMSSTRAHLUNGSIMULATOR_H
 
#include "FastSimulation/MaterialEffects/interface/MaterialEffectsSimulator.h"
 
#include <string>
 
class TF1;
 
class ParticlePropagator;
class RandomEngineAndDistribution;
 
class MuonBremsstrahlungSimulator : public MaterialEffectsSimulator {
public:
  MuonBremsstrahlungSimulator(
      double A, double Z, double density, double radLen, double photonEnergyCut, double photonFractECut);
 
  ~MuonBremsstrahlungSimulator() override {}
 
  // Returns the actual Muon brem Energy
  inline const XYZTLorentzVector& deltaP_BremMuon() const { return deltaPMuon; }
 
  // Returns the actual photon brem Energy
  inline const XYZTLorentzVector& deltaP_BremPhoton() const { return brem_photon; }
 
private:
  TF1* f1;
 
  int npar;
 
  double photonEnergy;
  double bremProba;
  double photonFractE;
 
  double xmin, xmax, rand;
  double d;  //distance
 
  unsigned int poisson(double ymu);
 
  void compute(ParticlePropagator& Particle, RandomEngineAndDistribution const*) override;
 
  XYZTLorentzVector brem(ParticlePropagator& p, RandomEngineAndDistribution const*) const;
 
  double gbteth(const double ener, const double partm, const double efrac, RandomEngineAndDistribution const*) const;
 
  // The actual Muon Brem
  XYZTLorentzVector deltaPMuon;
  //The photon brem
  XYZTLorentzVector brem_photon;
};
#endif