DirectElectronPairProduction.cxx

Go to the documentation of this file.

#include <cassert>

#include "DirectElectronPairProduction.h"
#include "KokoulinPetrukhinModel.h"

ClassImpT(DirectElectronPairProduction, T)

DirectElectronPairProduction<KokoulinPetrukhinModel> kokoulinpetrukhininstance;

//_________________________________________________________________________________________
template <class T> DirectElectronPairProduction<T>::DirectElectronPairProduction() :
EnergyLossProcess()
{
  fModel = new T();
}
//_________________________________________________________________________________________
template <class T>
DirectElectronPairProduction<T>::DirectElectronPairProduction(const Material & material) :
EnergyLossProcess()
{
  fModel = new T(material);

  Process::Process_t proc = fModel->ModeledProcess();

  assert( proc == Process::ePairProduction || proc == Process::eAll );
}
//_________________________________________________________________________________________
template <class T> DirectElectronPairProduction<T>::~DirectElectronPairProduction()
{
  if (fModel) delete fModel;
}
//_________________________________________________________________________________________
template <class T> double DirectElectronPairProduction<T>::dE_dx(double E) const
{
  if(E > Threshold() ) {

     ValidityRange_t vldrange = fModel->ValidityRange();

     if(E >= vldrange.Emin && E <= vldrange.Emax)  return fModel->dE_dx(E);
  }

  return 0;
}
//_________________________________________________________________________________________

---