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