/data/refman/pasoursint/CMSSW_6_1_2_SLHC2/src/FastSimulation/MaterialEffects/interface/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 RandomEngine;

class MuonBremsstrahlungSimulator : public MaterialEffectsSimulator
{
 public:

  MuonBremsstrahlungSimulator(const RandomEngine* engine,
                              double A, double Z, double density,
                              double radLen,double photonEnergyCut,double photonFractECut); 

  ~MuonBremsstrahlungSimulator() {}


// 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);

  XYZTLorentzVector brem(ParticlePropagator& p)const;

  double gbteth(const double ener,
                const double partm,
                const double efrac) const;

  // The actual Muon Brem
   XYZTLorentzVector deltaPMuon;
   //The photon brem
   XYZTLorentzVector brem_photon;



};
#endif