G4SPSEneDistribution Class Reference

#include <G4SPSEneDistribution.hh>

Public Member Functions
	G4SPSEneDistribution ()
	~G4SPSEneDistribution ()
void	SetEnergyDisType (const G4String &)
const G4String &	GetEnergyDisType ()
void	SetEmin (G4double)
G4double	GetEmin () const
G4double	GetArbEmin ()
void	SetEmax (G4double)
G4double	GetEmax () const
G4double	GetArbEmax ()
void	SetMonoEnergy (G4double)
void	SetAlpha (G4double)
void	SetBiasAlpha (G4double)
void	SetTemp (G4double)
void	SetBeamSigmaInE (G4double)
void	SetEzero (G4double)
void	SetGradient (G4double)
void	SetInterCept (G4double)
void	UserEnergyHisto (const G4ThreeVector &)
void	ArbEnergyHisto (const G4ThreeVector &)
void	ArbEnergyHistoFile (const G4String &)
void	EpnEnergyHisto (const G4ThreeVector &)
void	InputEnergySpectra (G4bool)
void	InputDifferentialSpectra (G4bool)
void	ArbInterpolate (const G4String &)
const G4String &	GetIntType ()
void	Calculate ()
void	SetBiasRndm (G4SPSRandomGenerator *a)
void	ReSetHist (const G4String &)
void	SetVerbosity (G4int a)
G4double	GetWeight () const
G4double	GetMonoEnergy ()
G4double	GetSE ()
G4double	Getalpha () const
G4double	GetEzero () const
G4double	GetTemp ()
G4double	Getgrad () const
G4double	Getcept () const
G4PhysicsFreeVector	GetUserDefinedEnergyHisto ()
G4PhysicsFreeVector	GetArbEnergyHisto ()
G4double	GenerateOne (G4ParticleDefinition *)
G4double	GetProbability (G4double)
G4double	GetArbEneWeight (G4double)
void	ApplyEnergyWeight (G4bool val)
G4bool	IfApplyEnergyWeight () const

Detailed Description

Definition at line 70 of file G4SPSEneDistribution.hh.

Constructor & Destructor Documentation

G4SPSEneDistribution()

G4SPSEneDistribution::G4SPSEneDistribution

(

)

Definition at line 41 of file G4SPSEneDistribution.cc.

{
  G4MUTEXINIT(mutex);
 
  // Initialise all variables
 
  particle_energy = 1.0 * MeV;
  EnergyDisType = "Mono";
  weight=1.;
  MonoEnergy = 1 * MeV;
  Emin = 0.;
  Emax = 1.e30;
  alpha = 0.;
  biasalpha = 0.;
  prob_norm = 1.0;
  Ezero = 0.;
  SE = 0.;
  Temp = 0.;
  grad = 0.;
  cept = 0.;
  IntType = "NULL"; // Interpolation type
 
  ArbEmin = 0.;
  ArbEmax = 1.e30;
 
  verbosityLevel = 0;
    
  threadLocal_t& data = threadLocalData.Get();
  data.Emax = Emax;
  data.Emin = Emin;
  data.alpha =alpha;
  data.cept = cept;
  data.Ezero = Ezero;
  data.grad = grad;
  data.particle_energy = 0.;
  data.particle_definition = nullptr;
  data.weight = weight;
}

~G4SPSEneDistribution()

G4SPSEneDistribution::~G4SPSEneDistribution

(

)

Definition at line 80 of file G4SPSEneDistribution.cc.

{
  G4MUTEXDESTROY(mutex);
  if(Arb_grad_cept_flag)
  {
    delete [] Arb_grad;
    delete [] Arb_cept;
  }
    
  if(Arb_alpha_Const_flag)
  {
    delete [] Arb_alpha;
    delete [] Arb_Const;
  }
    
  if(Arb_ezero_flag)
  {
    delete [] Arb_ezero;
  }
  delete Bbody_x;
  delete BBHist;
  delete CP_x;
  delete CPHist;
  for (auto & it : SplineInt)
  {
    delete it;
    it = nullptr;
  }
  SplineInt.clear();
}

Member Function Documentation

ApplyEnergyWeight()

void G4SPSEneDistribution::ApplyEnergyWeight ( G4bool  val )

inline

Definition at line 196 of file G4SPSEneDistribution.hh.

{ applyEvergyWeight = val; }

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

ArbEnergyHisto()

void G4SPSEneDistribution::ArbEnergyHisto

(

const G4ThreeVector &

input

)

Definition at line 315 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  G4double ehi = input.x(),
           val = input.y();
  if (verbosityLevel > 1)
  {
    G4cout << "In ArbEnergyHisto" << G4endl;
    G4cout << " " << ehi << " " << val << G4endl;
  }
  ArbEnergyH.InsertValues(ehi, val);
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

ArbEnergyHistoFile()

void G4SPSEneDistribution::ArbEnergyHistoFile

(

const G4String &

filename

)

Definition at line 328 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  std::ifstream infile(filename, std::ios::in);
  if (!infile)
  {
    G4Exception("G4SPSEneDistribution::ArbEnergyHistoFile", "Event0301",
                FatalException, "Unable to open the histo ASCII file");
  }
  G4double ehi, val;
  while (infile >> ehi >> val)
  {
    ArbEnergyH.InsertValues(ehi, val);
  }
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

ArbInterpolate()

void G4SPSEneDistribution::ArbInterpolate

(

const G4String &

IType

)

Definition at line 561 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
 
  IntType = IType;
  ArbEmax = ArbEnergyH.GetMaxEnergy();
  ArbEmin = ArbEnergyH.Energy(0);
 
  // Now interpolate points
 
  if (IntType == "Lin") LinearInterpolation();
  if (IntType == "Log") LogInterpolation();
  if (IntType == "Exp") ExpInterpolation();
  if (IntType == "Spline") SplineInterpolation();
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

Calculate()

void G4SPSEneDistribution::Calculate

(

)

Definition at line 360 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  if (EnergyDisType == "Cdg")
  {
    CalculateCdgSpectrum();
  }
  else if (EnergyDisType == "Bbody")
  {
    if(!BBhistInit)
    {
      BBInitHists();
    }
    CalculateBbodySpectrum();
  }
  else if (EnergyDisType == "CPow")
  {
    if(!CPhistInit)
    {
      CPInitHists();
    }
    CalculateCPowSpectrum();
  }
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

EpnEnergyHisto()

void G4SPSEneDistribution::EpnEnergyHisto

(

const G4ThreeVector &

input

)

Definition at line 344 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  G4double ehi = input.x(),
           val = input.y();
  if (verbosityLevel > 1)
  {
    G4cout << "In EpnEnergyHisto" << G4endl;
    G4cout << " " << ehi << " " << val << G4endl;
  }
  EpnEnergyH.InsertValues(ehi, val);
  Emax = ehi;
  threadLocalData.Get().Emax = Emax;
  Epnflag = true;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

GenerateOne()

G4double G4SPSEneDistribution::GenerateOne

(

G4ParticleDefinition *

a

)

Definition at line 1941 of file G4SPSEneDistribution.cc.

{
  // Copy global shared status to thread-local one
  //
  threadLocal_t& params = threadLocalData.Get();
  params.particle_definition=a;
  params.particle_energy=-1;
  if(applyEvergyWeight)
  {
    params.Emax = ArbEmax;
    params.Emin = ArbEmin;
  }
  else
  {
    params.Emax = Emax;
    params.Emin = Emin;
  }
  params.alpha = alpha;
  params.Ezero = Ezero;
  params.grad = grad;
  params.cept = cept;
  params.weight = weight;
  // particle_energy = -1.;
 
  if((EnergyDisType == "Mono") && ((MonoEnergy>Emax)||(MonoEnergy<Emin)))
  {
    G4ExceptionDescription ed;
    ed << "MonoEnergy " << G4BestUnit(MonoEnergy,"Energy")
       << " is outside of [Emin,Emax] = ["
       << G4BestUnit(Emin,"Energy") << ", "
       << G4BestUnit(Emax,"Energy") << ". MonoEnergy is used anyway.";
    G4Exception("G4SPSEneDistribution::GenerateOne()",
                "GPS0001", JustWarning, ed);
    params.particle_energy=MonoEnergy;
    return params.particle_energy;
  }
  while ( (EnergyDisType == "Arb")
        ?   (params.particle_energy < ArbEmin
          || params.particle_energy > ArbEmax)
        :   (params.particle_energy < params.Emin
          || params.particle_energy > params.Emax) )
  {
    if (Biased)
    {
      GenerateBiasPowEnergies();
    }
    else
    {
      if (EnergyDisType == "Mono")
      {
        GenerateMonoEnergetic();
      }
      else if (EnergyDisType == "Lin")
      {
        GenerateLinearEnergies(false);
      }
      else if (EnergyDisType == "Pow")
      {
        GeneratePowEnergies(false);
      }
      else if (EnergyDisType == "CPow")
      {
        GenerateCPowEnergies();
      }
      else if (EnergyDisType == "Exp")
      {
        GenerateExpEnergies(false);
      }
      else if (EnergyDisType == "Gauss")
      {
        GenerateGaussEnergies();
      }
      else if (EnergyDisType == "Brem")
      {
        GenerateBremEnergies();
      }
      else if (EnergyDisType == "Bbody")
      {
        GenerateBbodyEnergies();
      }
      else if (EnergyDisType == "Cdg")
      {
        GenerateCdgEnergies();
      }
      else if (EnergyDisType == "User")
      {
        GenUserHistEnergies();
      }
      else if (EnergyDisType == "Arb")
      {
        GenArbPointEnergies();
      }
      else if (EnergyDisType == "Epn")
      {
        GenEpnHistEnergies();
      }
      else
      {
        G4cout << "Error: EnergyDisType has unusual value" << G4endl;
      }
    }
  }
   return params.particle_energy;
}

Referenced by G4SingleParticleSource::GeneratePrimaryVertex().

Getalpha()

G4double G4SPSEneDistribution::Getalpha

(

)

const

Definition at line 262 of file G4SPSEneDistribution.cc.

{
  return threadLocalData.Get().alpha;
}

GetArbEmax()

G4double G4SPSEneDistribution::GetArbEmax

(

)

Definition at line 157 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return ArbEmax;
}

GetArbEmin()

G4double G4SPSEneDistribution::GetArbEmin

(

)

Definition at line 151 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return ArbEmin;
}

GetArbEnergyHisto()

G4PhysicsFreeVector G4SPSEneDistribution::GetArbEnergyHisto

(

)

Definition at line 294 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return ArbEnergyH;
}

GetArbEneWeight()

G4double G4SPSEneDistribution::GetArbEneWeight

(

G4double

ene

)

Definition at line 1679 of file G4SPSEneDistribution.cc.

{
  auto nbelow = IPDFArbEnergyH.FindBin(ene,(IPDFArbEnergyH.GetVectorLength())/2);
  G4double wei = 0.;
  if(IntType=="Lin")
  {
    // note: grad[i] and cept[i] are calculated with x[i-1] and x[i]
    auto gr = Arb_grad[nbelow + 1];
    auto ce = Arb_cept[nbelow + 1];
    wei = ene*gr + ce;
  }
  else if(IntType=="Log")
  {
    auto alp = Arb_alpha[nbelow + 1];
    auto cns = Arb_Const[nbelow + 1];
    wei = cns * std::pow(ene,alp);
  }
  else if(IntType=="Exp")
  {
    auto e0 = Arb_ezero[nbelow + 1];
    auto cns = Arb_Const[nbelow + 1];
    wei = cns * std::exp(-ene/e0);
  }
  else if(IntType=="Spline")
  {
    wei = SplineInt[nbelow+1]->CubicSplineInterpolation(ene);
  }
  return wei;
}

Referenced by G4SingleParticleSource::GeneratePrimaryVertex().

Getcept()

G4double G4SPSEneDistribution::Getcept

(

)

const

Definition at line 283 of file G4SPSEneDistribution.cc.

{
  return threadLocalData.Get().cept;
}

GetEmax()

G4double G4SPSEneDistribution::GetEmax

(

)

const

Definition at line 170 of file G4SPSEneDistribution.cc.

{
  return threadLocalData.Get().Emax;
}

GetEmin()

G4double G4SPSEneDistribution::GetEmin

(

)

const

Definition at line 146 of file G4SPSEneDistribution.cc.

{
  return threadLocalData.Get().Emin;
}

GetEnergyDisType()

const G4String & G4SPSEneDistribution::GetEnergyDisType

(

)

Definition at line 133 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return EnergyDisType;
}

Referenced by G4GeneralParticleSource::ListSource(), and G4GeneralParticleSourceMessenger::SetNewValue().

GetEzero()

G4double G4SPSEneDistribution::GetEzero

(

)

const

Definition at line 267 of file G4SPSEneDistribution.cc.

{
  return threadLocalData.Get().Ezero;
}

Getgrad()

G4double G4SPSEneDistribution::Getgrad

(

)

const

Definition at line 278 of file G4SPSEneDistribution.cc.

{
  return threadLocalData.Get().grad;
}

GetIntType()

const G4String & G4SPSEneDistribution::GetIntType

(

)

Definition at line 227 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return IntType;
}

GetMonoEnergy()

G4double G4SPSEneDistribution::GetMonoEnergy

(

)

Definition at line 250 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return MonoEnergy;
}

GetProbability()

G4double G4SPSEneDistribution::GetProbability

(

G4double

ene

)

Definition at line 2046 of file G4SPSEneDistribution.cc.

{
  G4double prob = 1.;
 
  threadLocal_t& params = threadLocalData.Get();
  if (EnergyDisType == "Lin")
  {
    if (prob_norm == 1.)
    {
      prob_norm = 0.5*params.grad*params.Emax*params.Emax
                + params.cept*params.Emax
                - 0.5*params.grad*params.Emin*params.Emin
                - params.cept*params.Emin;
    }
    prob = params.cept + params.grad * ene;
    prob /= prob_norm;
  }
  else if (EnergyDisType == "Pow")
  {
    if (prob_norm == 1.)
    {
      if (alpha != -1.)
      {
        G4double emina = std::pow(params.Emin, params.alpha + 1);
        G4double emaxa = std::pow(params.Emax, params.alpha + 1);
        prob_norm = 1./(1.+alpha) * (emaxa - emina);
      }
      else
      {
        prob_norm = std::log(params.Emax) - std::log(params.Emin) ;
      }
    }
    prob = std::pow(ene, params.alpha)/prob_norm;
  }
  else if (EnergyDisType == "Exp")
  {
    if (prob_norm == 1.)
    {
      prob_norm = -params.Ezero*(std::exp(-params.Emax/params.Ezero)
                               - std::exp(params.Emin/params.Ezero));
    }  
    prob = std::exp(-ene / params.Ezero);
    prob /= prob_norm;
  }
  else if (EnergyDisType == "Arb")
  {
    prob = ArbEnergyH.Value(ene);
 
    if (prob <= 0.)
    {
      G4cout << " Warning:G4SPSEneDistribution::GetProbability: prob<= 0. "
             << prob << " " << ene << G4endl;
      prob = 1e-30;
    }
  }
  else
  {
    G4cout << "Error: EnergyDisType not supported" << G4endl;
  }
 
  return prob;
}

GetSE()

G4double G4SPSEneDistribution::GetSE

(

)

Definition at line 256 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return SE;
}

GetTemp()

G4double G4SPSEneDistribution::GetTemp

(

)

Definition at line 272 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return Temp;
}

GetUserDefinedEnergyHisto()

G4PhysicsFreeVector G4SPSEneDistribution::GetUserDefinedEnergyHisto

(

)

Definition at line 288 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  return UDefEnergyH;
}

GetWeight()

G4double G4SPSEneDistribution::GetWeight

(

)

const

Definition at line 245 of file G4SPSEneDistribution.cc.

{
  return threadLocalData.Get().weight;
}

Referenced by G4SingleParticleSource::GeneratePrimaryVertex().

IfApplyEnergyWeight()

G4bool G4SPSEneDistribution::IfApplyEnergyWeight ( ) const

inline

Definition at line 197 of file G4SPSEneDistribution.hh.

{ return applyEvergyWeight; }

Referenced by G4SingleParticleSource::GeneratePrimaryVertex().

InputDifferentialSpectra()

void G4SPSEneDistribution::InputDifferentialSpectra

(

G4bool

value

)

Definition at line 548 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
 
  // Allows user to specify integral or differential spectra
  //
  DiffSpec = value; // true = differential, false = integral
  if (verbosityLevel > 1)
  {
    G4cout << "Diffspec has value " << DiffSpec << G4endl;
  }
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

InputEnergySpectra()

void G4SPSEneDistribution::InputEnergySpectra

(

G4bool

value

)

Definition at line 535 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
 
  // Allows user to specify spectrum is momentum
  //
  EnergySpec = value; // false if momentum
  if (verbosityLevel > 1)
  {
    G4cout << "EnergySpec has value " << EnergySpec << G4endl;
  }
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

ReSetHist()

void G4SPSEneDistribution::ReSetHist

(

const G4String &

atype

)

Definition at line 1914 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  if (atype == "energy")
  {
    UDefEnergyH = IPDFEnergyH = ZeroPhysVector;
    IPDFEnergyExist = false;
    Emin = 0.;
    Emax = 1e30;
  }
  else if (atype == "arb")
  {
    ArbEnergyH = IPDFArbEnergyH = ZeroPhysVector;
    IPDFArbExist = false;
  }
  else if (atype == "epn")
  {
    UDefEnergyH = IPDFEnergyH = ZeroPhysVector;
    IPDFEnergyExist = false;
    EpnEnergyH = ZeroPhysVector;
  }
  else
  {
    G4cout << "Error, histtype not accepted " << G4endl;
  }
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetAlpha()

void G4SPSEneDistribution::SetAlpha

(

G4double

alp

)

Definition at line 186 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  alpha = alp;
  threadLocalData.Get().alpha = alpha;
}

Referenced by G4AdjointPrimaryGenerator::G4AdjointPrimaryGenerator(), and G4GeneralParticleSourceMessenger::SetNewValue().

SetBeamSigmaInE()

void G4SPSEneDistribution::SetBeamSigmaInE

(

G4double

e

)

Definition at line 181 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  SE = e;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetBiasAlpha()

void G4SPSEneDistribution::SetBiasAlpha

(

G4double

alp

)

Definition at line 193 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  biasalpha = alp;
  Biased = true;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetBiasRndm()

void G4SPSEneDistribution::SetBiasRndm

(

G4SPSRandomGenerator *

a

)

Definition at line 233 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  eneRndm = a;
}

Referenced by G4SingleParticleSource::G4SingleParticleSource().

SetEmax()

void G4SPSEneDistribution::SetEmax

(

G4double

ema

)

Definition at line 163 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  Emax = ema;
  threadLocalData.Get().Emax = Emax;
}

Referenced by G4AdjointPrimaryGenerator::GenerateAdjointPrimaryVertex(), G4AdjointPrimaryGenerator::GenerateFwdPrimaryVertex(), and G4GeneralParticleSourceMessenger::SetNewValue().

SetEmin()

void G4SPSEneDistribution::SetEmin

(

G4double

emi

)

Definition at line 139 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  Emin = emi;
  threadLocalData.Get().Emin = Emin;
}

Referenced by G4AdjointPrimaryGenerator::GenerateAdjointPrimaryVertex(), G4AdjointPrimaryGenerator::GenerateFwdPrimaryVertex(), and G4GeneralParticleSourceMessenger::SetNewValue().

SetEnergyDisType()

void G4SPSEneDistribution::SetEnergyDisType

(

const G4String &

DisType

)

Definition at line 111 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  EnergyDisType = DisType;
  if (EnergyDisType == "User")
  {
    UDefEnergyH = IPDFEnergyH = ZeroPhysVector;
    IPDFEnergyExist = false;
  }
  else if (EnergyDisType == "Arb")
  {
    ArbEnergyH = IPDFArbEnergyH = ZeroPhysVector;
    IPDFArbExist = false;
  }
  else if (EnergyDisType == "Epn")
  {
    UDefEnergyH = IPDFEnergyH = ZeroPhysVector;
    IPDFEnergyExist = false;
    EpnEnergyH = ZeroPhysVector;
  }
}

Referenced by G4AdjointPrimaryGenerator::G4AdjointPrimaryGenerator(), and G4GeneralParticleSourceMessenger::SetNewValue().

SetEzero()

void G4SPSEneDistribution::SetEzero

(

G4double

eze

)

Definition at line 206 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  Ezero = eze;
  threadLocalData.Get().Ezero = Ezero;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetGradient()

void G4SPSEneDistribution::SetGradient

(

G4double

gr

)

Definition at line 213 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  grad = gr;
  threadLocalData.Get().grad = grad;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetInterCept()

void G4SPSEneDistribution::SetInterCept

(

G4double

c

)

Definition at line 220 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  cept = c;
  threadLocalData.Get().cept = cept;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetMonoEnergy()

void G4SPSEneDistribution::SetMonoEnergy

(

G4double

menergy

)

Definition at line 175 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  MonoEnergy = menergy;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetTemp()

void G4SPSEneDistribution::SetTemp

(

G4double

tem

)

Definition at line 200 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  Temp = tem;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetVerbosity()

void G4SPSEneDistribution::SetVerbosity

(

G4int

a

)

Definition at line 239 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  verbosityLevel = a;
}

Referenced by G4SingleParticleSource::SetVerbosity().

UserEnergyHisto()

void G4SPSEneDistribution::UserEnergyHisto

(

const G4ThreeVector &

input

)

Definition at line 300 of file G4SPSEneDistribution.cc.

{
  G4AutoLock l(&mutex);
  G4double ehi = input.x(),
           val = input.y();
  if (verbosityLevel > 1)
  {
    G4cout << "In UserEnergyHisto" << G4endl;
    G4cout << " " << ehi << " " << val << G4endl;
  }
  UDefEnergyH.InsertValues(ehi, val);
  Emax = ehi;
  threadLocalData.Get().Emax = Emax;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

---

The documentation for this class was generated from the following files:

• G4SPSEneDistribution.hh
• G4SPSEneDistribution.cc