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