G4SPSEneDistribution Class Reference

#include <G4SPSEneDistribution.hh>

Public Member Functions
	G4SPSEneDistribution ()
	~G4SPSEneDistribution ()
void	SetEnergyDisType (G4String)
G4String	GetEnergyDisType ()
void	SetEmin (G4double)
G4double	GetEmin ()
G4double	GetArbEmin ()
void	SetEmax (G4double)
G4double	GetEmax ()
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 (G4ThreeVector)
void	ArbEnergyHisto (G4ThreeVector)
void	ArbEnergyHistoFile (G4String)
void	EpnEnergyHisto (G4ThreeVector)
void	InputEnergySpectra (G4bool)
void	InputDifferentialSpectra (G4bool)
void	ArbInterpolate (G4String)
G4String	GetIntType ()
void	Calculate ()
void	SetBiasRndm (G4SPSRandomGenerator *a)
void	ReSetHist (G4String)
void	SetVerbosity (G4int a)
G4double	GetWeight ()
G4double	GetMonoEnergy ()
G4double	GetSE ()
G4double	Getalpha ()
G4double	GetEzero ()
G4double	GetTemp ()
G4double	Getgrad ()
G4double	Getcept ()
G4PhysicsOrderedFreeVector	GetUserDefinedEnergyHisto ()
G4PhysicsOrderedFreeVector	GetArbEnergyHisto ()
G4double	GenerateOne (G4ParticleDefinition *)
G4double	GetProbability (G4double)

Detailed Description

Definition at line 146 of file G4SPSEneDistribution.hh.

Constructor & Destructor Documentation

G4SPSEneDistribution::G4SPSEneDistribution

(

)

Definition at line 53 of file G4SPSEneDistribution.cc.

References python.hepunit::MeV.

  : particle_definition(0), eneRndm(0), Splinetemp(0)
{
    //
    // 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.;
    Biased = false; // not biased
    EnergySpec = true; // true - energy spectra, false - momentum spectra
    DiffSpec = true; // true - differential spec, false integral spec
    IntType = "NULL"; // Interpolation type
    IPDFEnergyExist = false;
    IPDFArbExist = false;

    ArbEmin = 0.;
    ArbEmax = 1.e30;

    verbosityLevel = 0;

}

G4SPSEneDistribution::~G4SPSEneDistribution

(

)

Definition at line 87 of file G4SPSEneDistribution.cc.

                                            {
}

Member Function Documentation

void G4SPSEneDistribution::ArbEnergyHisto

(

G4ThreeVector

input

)

Definition at line 157 of file G4SPSEneDistribution.cc.

References G4cout, G4endl, G4PhysicsOrderedFreeVector::InsertValues(), CLHEP::Hep3Vector::x(), and CLHEP::Hep3Vector::y().

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

void G4SPSEneDistribution::ArbEnergyHistoFile

(

G4String

filename

)

Definition at line 168 of file G4SPSEneDistribution.cc.

References FatalException, G4Exception(), and G4PhysicsOrderedFreeVector::InsertValues().

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                                               {
    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);
    }
}

void G4SPSEneDistribution::ArbInterpolate

(

G4String

IType

)

Definition at line 295 of file G4SPSEneDistribution.cc.

References G4cout, G4endl, G4PhysicsOrderedFreeVector::GetMaxLowEdgeEnergy(), and G4PhysicsOrderedFreeVector::GetMinLowEdgeEnergy().

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                                        {
    if (EnergyDisType != "Arb")
        G4cout << "Error: this is for arbitrary distributions" << G4endl;
    IntType = IType;
    ArbEmax = ArbEnergyH.GetMaxLowEdgeEnergy();
    ArbEmin = ArbEnergyH.GetMinLowEdgeEnergy();

    // Now interpolate points
    if (IntType == "Lin")
        LinearInterpolation();
    if (IntType == "Log")
        LogInterpolation();
    if (IntType == "Exp")
        ExpInterpolation();
    if (IntType == "Spline")
        SplineInterpolation();
}

void G4SPSEneDistribution::Calculate

(

)

Definition at line 193 of file G4SPSEneDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                     {
    if (EnergyDisType == "Cdg")
        CalculateCdgSpectrum();
    else if (EnergyDisType == "Bbody")
        CalculateBbodySpectrum();
}

void G4SPSEneDistribution::EpnEnergyHisto

(

G4ThreeVector

input

)

Definition at line 180 of file G4SPSEneDistribution.cc.

References G4cout, G4endl, G4PhysicsOrderedFreeVector::InsertValues(), CLHEP::Hep3Vector::x(), and CLHEP::Hep3Vector::y().

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

G4double G4SPSEneDistribution::GenerateOne

(

G4ParticleDefinition *

a

)

Definition at line 1217 of file G4SPSEneDistribution.cc.

References test::a, G4cout, and G4endl.

Referenced by G4SingleParticleSource::GeneratePrimaryVertex().

                                                                  {
    particle_definition = a;
    particle_energy = -1.;

    while ((EnergyDisType == "Arb") ? (particle_energy < ArbEmin
            || particle_energy > ArbEmax) : (particle_energy < Emin
            || particle_energy > Emax)) {
        if (Biased) {
            GenerateBiasPowEnergies();
        } else {
            if (EnergyDisType == "Mono")
                GenerateMonoEnergetic();
            else if (EnergyDisType == "Lin")
                GenerateLinearEnergies();
            else if (EnergyDisType == "Pow")
                GeneratePowEnergies();
            else if (EnergyDisType == "Exp")
                GenerateExpEnergies();
            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 particle_energy;
}

G4double G4SPSEneDistribution::Getalpha ( )

inline

Definition at line 220 of file G4SPSEneDistribution.hh.

                        {
        return alpha;
    }

G4double G4SPSEneDistribution::GetArbEmax ( )

inline

Definition at line 170 of file G4SPSEneDistribution.hh.

                                 {
        return ArbEmax;
    }

G4double G4SPSEneDistribution::GetArbEmin ( )

inline

Definition at line 161 of file G4SPSEneDistribution.hh.

                                 {
        return ArbEmin;
    }

G4PhysicsOrderedFreeVector G4SPSEneDistribution::GetArbEnergyHisto ( )

inline

Definition at line 245 of file G4SPSEneDistribution.hh.

                                                          {
        return ArbEnergyH;
    }

G4double G4SPSEneDistribution::Getcept ( )

inline

Definition at line 236 of file G4SPSEneDistribution.hh.

                       {
        return cept;
    }

G4double G4SPSEneDistribution::GetEmax ( )

inline

Definition at line 166 of file G4SPSEneDistribution.hh.

                              {
        return Emax;
    }

G4double G4SPSEneDistribution::GetEmin ( )

inline

Definition at line 157 of file G4SPSEneDistribution.hh.

                              {
        return Emin;
    }

G4String G4SPSEneDistribution::GetEnergyDisType ( )

inline

Definition at line 152 of file G4SPSEneDistribution.hh.

Referenced by UltraPrimaryGeneratorAction::GeneratePrimaries().

                                       {
        return EnergyDisType;
    }

G4double G4SPSEneDistribution::GetEzero ( )

inline

Definition at line 224 of file G4SPSEneDistribution.hh.

                        {
        return Ezero;
    }

G4double G4SPSEneDistribution::Getgrad ( )

inline

Definition at line 232 of file G4SPSEneDistribution.hh.

                       {
        return grad;
    }

G4String G4SPSEneDistribution::GetIntType ( )

inline

Definition at line 190 of file G4SPSEneDistribution.hh.

                                 {
        return IntType;
    }

G4double G4SPSEneDistribution::GetMonoEnergy ( )

inline

Definition at line 212 of file G4SPSEneDistribution.hh.

                             {
        return MonoEnergy;
    }

G4double G4SPSEneDistribution::GetProbability

(

G4double

ene

)

Definition at line 1256 of file G4SPSEneDistribution.cc.

References G4cout, G4endl, and G4PhysicsVector::Value().

                                                          {
    G4double prob = 1.;

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

G4double G4SPSEneDistribution::GetSE ( )

inline

Definition at line 216 of file G4SPSEneDistribution.hh.

                     {
        return SE;
    }

G4double G4SPSEneDistribution::GetTemp ( )

inline

Definition at line 228 of file G4SPSEneDistribution.hh.

                       {
        return Temp;
    }

G4PhysicsOrderedFreeVector G4SPSEneDistribution::GetUserDefinedEnergyHisto ( )

inline

Definition at line 241 of file G4SPSEneDistribution.hh.

                                                                  {
        return UDefEnergyH;
    }

G4double G4SPSEneDistribution::GetWeight ( )

inline

Definition at line 208 of file G4SPSEneDistribution.hh.

Referenced by G4SingleParticleSource::GeneratePrimaryVertex().

                         {
        return weight;
    }

void G4SPSEneDistribution::InputDifferentialSpectra

(

G4bool

value

)

Definition at line 288 of file G4SPSEneDistribution.cc.

References G4cout, and G4endl.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

void G4SPSEneDistribution::InputEnergySpectra

(

G4bool

value

)

Definition at line 281 of file G4SPSEneDistribution.cc.

References G4cout, and G4endl.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

void G4SPSEneDistribution::ReSetHist

(

G4String

atype

)

Definition at line 1199 of file G4SPSEneDistribution.cc.

References G4cout, and G4endl.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                                   {
    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;
    }
}

void G4SPSEneDistribution::SetAlpha

(

G4double

alp

)

Definition at line 120 of file G4SPSEneDistribution.cc.

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

                                                {
    alpha = alp;
}

void G4SPSEneDistribution::SetBeamSigmaInE

(

G4double

e

)

Definition at line 117 of file G4SPSEneDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                                     {
    SE = e;
}

void G4SPSEneDistribution::SetBiasAlpha

(

G4double

alp

)

Definition at line 124 of file G4SPSEneDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                                    {
    biasalpha = alp;
    Biased = true;
}

void G4SPSEneDistribution::SetBiasRndm ( G4SPSRandomGenerator *  a )

inline

Definition at line 196 of file G4SPSEneDistribution.hh.

References test::a.

Referenced by G4SingleParticleSource::G4SingleParticleSource().

                                              {
        eneRndm = a;
    }

void G4SPSEneDistribution::SetEmax

(

G4double

ema

)

Definition at line 109 of file G4SPSEneDistribution.cc.

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

                                               {
    Emax = ema;
}

void G4SPSEneDistribution::SetEmin

(

G4double

emi

)

Definition at line 105 of file G4SPSEneDistribution.cc.

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

                                               {
    Emin = emi;
}

void G4SPSEneDistribution::SetEnergyDisType

(

G4String

DisType

)

Definition at line 90 of file G4SPSEneDistribution.cc.

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

                                                            {
    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;
    }
}

void G4SPSEneDistribution::SetEzero

(

G4double

eze

)

Definition at line 133 of file G4SPSEneDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                                {
    Ezero = eze;
}

void G4SPSEneDistribution::SetGradient

(

G4double

gr

)

Definition at line 137 of file G4SPSEneDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                                  {
    grad = gr;
}

void G4SPSEneDistribution::SetInterCept

(

G4double

c

)

Definition at line 141 of file G4SPSEneDistribution.cc.

References test::c.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                                  {
    cept = c;
}

void G4SPSEneDistribution::SetMonoEnergy

(

G4double

menergy

)

Definition at line 113 of file G4SPSEneDistribution.cc.

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

                                                         {
    MonoEnergy = menergy;
}

void G4SPSEneDistribution::SetTemp

(

G4double

tem

)

Definition at line 129 of file G4SPSEneDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

                                               {
    Temp = tem;
}

void G4SPSEneDistribution::SetVerbosity ( G4int  a )

inline

Definition at line 203 of file G4SPSEneDistribution.hh.

References test::a.

Referenced by G4SingleParticleSource::SetVerbosity().

                               {
        verbosityLevel = a;
    }

void G4SPSEneDistribution::UserEnergyHisto

(

G4ThreeVector

input

)

Definition at line 145 of file G4SPSEneDistribution.cc.

References G4cout, G4endl, G4PhysicsOrderedFreeVector::InsertValues(), CLHEP::Hep3Vector::x(), and CLHEP::Hep3Vector::y().

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

---

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

• G4SPSEneDistribution.hh
• G4SPSEneDistribution.cc