G4SPSRandomGenerator Class Reference

#include <G4SPSRandomGenerator.hh>

Collaboration diagram for G4SPSRandomGenerator:

Classes
struct	a_check
struct	bweights_t

Public Member Functions
	G4SPSRandomGenerator ()
	~G4SPSRandomGenerator ()
void	SetXBias (G4ThreeVector)
void	SetYBias (G4ThreeVector)
void	SetZBias (G4ThreeVector)
void	SetThetaBias (G4ThreeVector)
void	SetPhiBias (G4ThreeVector)
void	SetEnergyBias (G4ThreeVector)
void	SetPosThetaBias (G4ThreeVector)
void	SetPosPhiBias (G4ThreeVector)
G4double	GenRandX ()
G4double	GenRandY ()
G4double	GenRandZ ()
G4double	GenRandTheta ()
G4double	GenRandPhi ()
G4double	GenRandEnergy ()
G4double	GenRandPosTheta ()
G4double	GenRandPosPhi ()
void	SetIntensityWeight (G4double weight)
G4double	GetBiasWeight ()
void	ReSetHist (G4String)
void	SetVerbosity (G4int a)

Private Attributes
G4Cache< a_check >	local_IPDFXBias
G4bool	XBias
G4bool	IPDFXBias
G4PhysicsOrderedFreeVector	XBiasH
G4PhysicsOrderedFreeVector	IPDFXBiasH
G4Cache< a_check >	local_IPDFYBias
G4bool	YBias
G4bool	IPDFYBias
G4PhysicsOrderedFreeVector	YBiasH
G4PhysicsOrderedFreeVector	IPDFYBiasH
G4Cache< a_check >	local_IPDFZBias
G4bool	ZBias
G4bool	IPDFZBias
G4PhysicsOrderedFreeVector	ZBiasH
G4PhysicsOrderedFreeVector	IPDFZBiasH
G4Cache< a_check >	local_IPDFThetaBias
G4bool	ThetaBias
G4bool	IPDFThetaBias
G4PhysicsOrderedFreeVector	ThetaBiasH
G4PhysicsOrderedFreeVector	IPDFThetaBiasH
G4Cache< a_check >	local_IPDFPhiBias
G4bool	PhiBias
G4bool	IPDFPhiBias
G4PhysicsOrderedFreeVector	PhiBiasH
G4PhysicsOrderedFreeVector	IPDFPhiBiasH
G4Cache< a_check >	local_IPDFEnergyBias
G4bool	EnergyBias
G4bool	IPDFEnergyBias
G4PhysicsOrderedFreeVector	EnergyBiasH
G4PhysicsOrderedFreeVector	IPDFEnergyBiasH
G4Cache< a_check >	local_IPDFPosThetaBias
G4bool	PosThetaBias
G4bool	IPDFPosThetaBias
G4PhysicsOrderedFreeVector	PosThetaBiasH
G4PhysicsOrderedFreeVector	IPDFPosThetaBiasH
G4Cache< a_check >	local_IPDFPosPhiBias
G4bool	PosPhiBias
G4bool	IPDFPosPhiBias
G4PhysicsOrderedFreeVector	PosPhiBiasH
G4PhysicsOrderedFreeVector	IPDFPosPhiBiasH
G4Cache< bweights_t >	bweights
G4int	verbosityLevel
G4Mutex	mutex

Detailed Description

Andrea Dotti Feb 2015 Important: This is a shared class between threads. Only one thread should use the set-methods here. Note that this is exactly what is achieved using UI commands. If you use the set methods to set defaults in your application take care that only one thread is executing them. In addition take care of calling these methods before the run is started Do not use these setters during the event loop

Definition at line 156 of file G4SPSRandomGenerator.hh.

Constructor & Destructor Documentation

G4SPSRandomGenerator()

G4SPSRandomGenerator::G4SPSRandomGenerator

(

)

Definition at line 73 of file G4SPSRandomGenerator.cc.

{
    // Initialise all variables

    // Bias variables
    XBias = false;
    IPDFXBias = false;
    YBias = false;
    IPDFYBias = false;
    ZBias = false;
    IPDFZBias = false;
    ThetaBias = false;
    IPDFThetaBias = false;
    PhiBias = false;
    IPDFPhiBias = false;
    EnergyBias = false;
    IPDFEnergyBias = false;
    PosThetaBias = false;
    IPDFPosThetaBias = false;
    PosPhiBias = false;
    IPDFPosPhiBias = false;
    verbosityLevel = 0;
    G4MUTEXINIT(mutex);
}

~G4SPSRandomGenerator()

G4SPSRandomGenerator::~G4SPSRandomGenerator

(

)

Definition at line 98 of file G4SPSRandomGenerator.cc.

                                            {
    G4MUTEXDESTROY(mutex);
}

Member Function Documentation

GenRandEnergy()

G4double G4SPSRandomGenerator::GenRandEnergy

(

)

Definition at line 568 of file G4SPSRandomGenerator.cc.

                                             {
    if (verbosityLevel >= 1)
        G4cout << "In GenRandEnergy" << G4endl;
    if (EnergyBias == false) {
        // Energy is not biased
        G4double rndm = G4UniformRand();
        return (rndm);
    } else {
        if ( local_IPDFEnergyBias.Get().val == false ) {
            local_IPDFEnergyBias.Get().val = true;
            // ENERGY is biased
            G4AutoLock l(&mutex);
        if (IPDFEnergyBias == false) {
            // IPDF has not been created, so create it
            G4double bins[1024], vals[1024], sum;
            G4int ii;
            G4int maxbin = G4int(EnergyBiasH.GetVectorLength());
            bins[0] = EnergyBiasH.GetLowEdgeEnergy(size_t(0));
            vals[0] = EnergyBiasH(size_t(0));
            sum = vals[0];
            for (ii = 1; ii < maxbin; ii++) {
                bins[ii] = EnergyBiasH.GetLowEdgeEnergy(size_t(ii));
                vals[ii] = EnergyBiasH(size_t(ii)) + vals[ii - 1];
                sum = sum + EnergyBiasH(size_t(ii));
            }
            IPDFEnergyBiasH = ZeroPhysVector;
            for (ii = 0; ii < maxbin; ii++) {
                vals[ii] = vals[ii] / sum;
                IPDFEnergyBiasH.InsertValues(bins[ii], vals[ii]);
            }
            // Make IPDFEnergyBias = true
            IPDFEnergyBias = true;
        }
        }
        // IPDF has been create so carry on
        G4double rndm = G4UniformRand();
        //  size_t weight_bin_no = IPDFEnergyBiasH.FindValueBinLocation(rndm);
        size_t numberOfBin = IPDFEnergyBiasH.GetVectorLength();
        G4int biasn1 = 0;
        G4int biasn2 = numberOfBin / 2;
        G4int biasn3 = numberOfBin - 1;
        while (biasn1 != biasn3 - 1) {
            if (rndm > IPDFEnergyBiasH(biasn2))
              biasn1 = biasn2;
            else
              biasn3 = biasn2;
            biasn2 = biasn1 + (biasn3 - biasn1 + 1) / 2;
        }
        bweights_t& w = bweights.Get();
        w[5] = IPDFEnergyBiasH(biasn2) - IPDFEnergyBiasH(biasn2 - 1);
        G4double xaxisl = IPDFEnergyBiasH.GetLowEdgeEnergy(size_t(biasn2 - 1));
        G4double xaxisu = IPDFEnergyBiasH.GetLowEdgeEnergy(size_t(biasn2));
        G4double NatProb = xaxisu - xaxisl;
        w[5] = NatProb / w[5];
        if (verbosityLevel >= 1)
          G4cout << "Energy bin weight " << w[5] << " " << rndm
          << G4endl;
        return (IPDFEnergyBiasH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenRandPhi()

G4double G4SPSRandomGenerator::GenRandPhi

(

)

Definition at line 508 of file G4SPSRandomGenerator.cc.

                                          {
    if (verbosityLevel >= 1)
        G4cout << "In GenRandPhi" << G4endl;
    if (PhiBias == false) {
        // Phi is not biased
        G4double rndm = G4UniformRand();
        return (rndm);
    } else {
        // Phi is biased
        if ( local_IPDFPhiBias.Get().val == false ) {
            local_IPDFPhiBias.Get().val = true;
            G4AutoLock l(&mutex);
            if (IPDFPhiBias == false) {
                // IPDF has not been created, so create it
                G4double bins[1024], vals[1024], sum;
                G4int ii;
                G4int maxbin = G4int(PhiBiasH.GetVectorLength());
                bins[0] = PhiBiasH.GetLowEdgeEnergy(size_t(0));
                vals[0] = PhiBiasH(size_t(0));
                sum = vals[0];
                for (ii = 1; ii < maxbin; ii++) {
                    bins[ii] = PhiBiasH.GetLowEdgeEnergy(size_t(ii));
                    vals[ii] = PhiBiasH(size_t(ii)) + vals[ii - 1];
                    sum = sum + PhiBiasH(size_t(ii));
                }

                for (ii = 0; ii < maxbin; ii++) {
                    vals[ii] = vals[ii] / sum;
                    IPDFPhiBiasH.InsertValues(bins[ii], vals[ii]);
                }
                // Make IPDFPhiBias = true
                IPDFPhiBias = true;
        }
        }
        // IPDF has been create so carry on
        G4double rndm = G4UniformRand();
        //      size_t weight_bin_no = IPDFPhiBiasH.FindValueBinLocation(rndm);
        size_t numberOfBin = IPDFPhiBiasH.GetVectorLength();
        G4int biasn1 = 0;
        G4int biasn2 = numberOfBin / 2;
        G4int biasn3 = numberOfBin - 1;
        while (biasn1 != biasn3 - 1) {
            if (rndm > IPDFPhiBiasH(biasn2))
              biasn1 = biasn2;
            else
              biasn3 = biasn2;
            biasn2 = biasn1 + (biasn3 - biasn1 + 1) / 2;
        }
        bweights_t& w = bweights.Get();
        w[4] = IPDFPhiBiasH(biasn2) - IPDFPhiBiasH(biasn2 - 1);
        G4double xaxisl = IPDFPhiBiasH.GetLowEdgeEnergy(size_t(biasn2 - 1));
        G4double xaxisu = IPDFPhiBiasH.GetLowEdgeEnergy(size_t(biasn2));
        G4double NatProb = xaxisu - xaxisl;
        w[4] = NatProb / w[4];
        if (verbosityLevel >= 1)
          G4cout << "Phi bin weight " << w[4] << " " << rndm << G4endl;
        return (IPDFPhiBiasH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenRandPosPhi()

G4double G4SPSRandomGenerator::GenRandPosPhi

(

)

Definition at line 693 of file G4SPSRandomGenerator.cc.

                                             {
    if (verbosityLevel >= 1)
        G4cout << "In GenRandPosPhi" << G4endl;
    if (PosPhiBias == false) {
        // PosPhi is not biased
        G4double rndm = G4UniformRand();
        return (rndm);
    } else {
        // PosPhi is biased
        if (local_IPDFPosPhiBias.Get().val == false ) {
            local_IPDFPosPhiBias.Get().val = true;
            G4AutoLock l(&mutex);
        if (IPDFPosPhiBias == false) {
            // IPDF has not been created, so create it
            G4double bins[1024], vals[1024], sum;
            G4int ii;
            G4int maxbin = G4int(PosPhiBiasH.GetVectorLength());
            bins[0] = PosPhiBiasH.GetLowEdgeEnergy(size_t(0));
            vals[0] = PosPhiBiasH(size_t(0));
            sum = vals[0];
            for (ii = 1; ii < maxbin; ii++) {
                bins[ii] = PosPhiBiasH.GetLowEdgeEnergy(size_t(ii));
                vals[ii] = PosPhiBiasH(size_t(ii)) + vals[ii - 1];
                sum = sum + PosPhiBiasH(size_t(ii));
            }

            for (ii = 0; ii < maxbin; ii++) {
                vals[ii] = vals[ii] / sum;
                IPDFPosPhiBiasH.InsertValues(bins[ii], vals[ii]);
            }
            // Make IPDFPosPhiBias = true
            IPDFPosPhiBias = true;
        }
        }
        // IPDF has been create so carry on
        G4double rndm = G4UniformRand();
        //      size_t weight_bin_no = IPDFPosPhiBiasH.FindValueBinLocation(rndm);
        size_t numberOfBin = IPDFPosPhiBiasH.GetVectorLength();
        G4int biasn1 = 0;
        G4int biasn2 = numberOfBin / 2;
        G4int biasn3 = numberOfBin - 1;
        while (biasn1 != biasn3 - 1) {
            if (rndm > IPDFPosPhiBiasH(biasn2))
              biasn1 = biasn2;
            else
              biasn3 = biasn2;
            biasn2 = biasn1 + (biasn3 - biasn1 + 1) / 2;
        }
        bweights_t& w = bweights.Get();
        w[7] = IPDFPosPhiBiasH(biasn2) - IPDFPosPhiBiasH(biasn2 - 1);
        G4double xaxisl = IPDFPosPhiBiasH.GetLowEdgeEnergy(size_t(biasn2 - 1));
        G4double xaxisu = IPDFPosPhiBiasH.GetLowEdgeEnergy(size_t(biasn2));
        G4double NatProb = xaxisu - xaxisl;
        w[7] = NatProb / w[7];
        if (verbosityLevel >= 1)
          G4cout << "PosPhi bin weight " << w[7] << " " << rndm
          << G4endl;
        return (IPDFPosPhiBiasH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenRandPosTheta()

G4double G4SPSRandomGenerator::GenRandPosTheta

(

)

Definition at line 629 of file G4SPSRandomGenerator.cc.

                                               {
    if (verbosityLevel >= 1) {
        G4cout << "In GenRandPosTheta" << G4endl;
        G4cout << "Verbosity " << verbosityLevel << G4endl;
    }
    if (PosThetaBias == false) {
        // Theta is not biased
        G4double rndm = G4UniformRand();
        return (rndm);
    } else {
        // Theta is biased
        if ( local_IPDFPosThetaBias.Get().val == false ) {
            local_IPDFPosThetaBias.Get().val = true;
            G4AutoLock l(&mutex);
            if (IPDFPosThetaBias == false) {
                // IPDF has not been created, so create it
                G4double bins[1024], vals[1024], sum;
                G4int ii;
                G4int maxbin = G4int(PosThetaBiasH.GetVectorLength());
                bins[0] = PosThetaBiasH.GetLowEdgeEnergy(size_t(0));
                vals[0] = PosThetaBiasH(size_t(0));
                sum = vals[0];
                for (ii = 1; ii < maxbin; ii++) {
                    bins[ii] = PosThetaBiasH.GetLowEdgeEnergy(size_t(ii));
                    vals[ii] = PosThetaBiasH(size_t(ii)) + vals[ii - 1];
                    sum = sum + PosThetaBiasH(size_t(ii));
                }

                for (ii = 0; ii < maxbin; ii++) {
                    vals[ii] = vals[ii] / sum;
                    IPDFPosThetaBiasH.InsertValues(bins[ii], vals[ii]);
                }
                // Make IPDFThetaBias = true
                IPDFPosThetaBias = true;
        }
        }
        // IPDF has been create so carry on
        G4double rndm = G4UniformRand();
        //      size_t weight_bin_no = IPDFThetaBiasH.FindValueBinLocation(rndm);
        size_t numberOfBin = IPDFPosThetaBiasH.GetVectorLength();
        G4int biasn1 = 0;
        G4int biasn2 = numberOfBin / 2;
        G4int biasn3 = numberOfBin - 1;
        while (biasn1 != biasn3 - 1) {
            if (rndm > IPDFPosThetaBiasH(biasn2))
              biasn1 = biasn2;
            else
              biasn3 = biasn2;
            biasn2 = biasn1 + (biasn3 - biasn1 + 1) / 2;
        }
        bweights_t& w = bweights.Get();
        w[6] = IPDFPosThetaBiasH(biasn2) - IPDFPosThetaBiasH(biasn2 - 1);
          G4double xaxisl =
              IPDFPosThetaBiasH.GetLowEdgeEnergy(size_t(biasn2 - 1));
          G4double xaxisu = IPDFPosThetaBiasH.GetLowEdgeEnergy(size_t(biasn2));
          G4double NatProb = xaxisu - xaxisl;
          w[6] = NatProb / w[6];
          if (verbosityLevel >= 1)
            G4cout << "PosTheta bin weight " << w[6] << " " << rndm
            << G4endl;
          return (IPDFPosThetaBiasH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenRandTheta()

G4double G4SPSRandomGenerator::GenRandTheta

(

)

Definition at line 445 of file G4SPSRandomGenerator.cc.

                                            {
    if (verbosityLevel >= 1) {
        G4cout << "In GenRandTheta" << G4endl;
        G4cout << "Verbosity " << verbosityLevel << G4endl;
    }
    if (ThetaBias == false) {
        // Theta is not biased
        G4double rndm = G4UniformRand();
        return (rndm);
    } else {
        // Theta is biased
        if ( local_IPDFThetaBias.Get().val == false ) {
            local_IPDFThetaBias.Get().val = true;
            G4AutoLock l(&mutex);
        if (IPDFThetaBias == false) {
            // IPDF has not been created, so create it
            G4double bins[1024], vals[1024], sum;
            G4int ii;
            G4int maxbin = G4int(ThetaBiasH.GetVectorLength());
            bins[0] = ThetaBiasH.GetLowEdgeEnergy(size_t(0));
            vals[0] = ThetaBiasH(size_t(0));
            sum = vals[0];
            for (ii = 1; ii < maxbin; ii++) {
                bins[ii] = ThetaBiasH.GetLowEdgeEnergy(size_t(ii));
                vals[ii] = ThetaBiasH(size_t(ii)) + vals[ii - 1];
                sum = sum + ThetaBiasH(size_t(ii));
            }

            for (ii = 0; ii < maxbin; ii++) {
                vals[ii] = vals[ii] / sum;
                IPDFThetaBiasH.InsertValues(bins[ii], vals[ii]);
            }
            // Make IPDFThetaBias = true
            IPDFThetaBias = true;
        }
        }
        // IPDF has been create so carry on
        G4double rndm = G4UniformRand();
        //      size_t weight_bin_no = IPDFThetaBiasH.FindValueBinLocation(rndm);
        size_t numberOfBin = IPDFThetaBiasH.GetVectorLength();
        G4int biasn1 = 0;
        G4int biasn2 = numberOfBin / 2;
        G4int biasn3 = numberOfBin - 1;
        while (biasn1 != biasn3 - 1) {
            if (rndm > IPDFThetaBiasH(biasn2))
              biasn1 = biasn2;
            else
              biasn3 = biasn2;
            biasn2 = biasn1 + (biasn3 - biasn1 + 1) / 2;
        }
        bweights_t& w = bweights.Get();
        w[3] = IPDFThetaBiasH(biasn2) - IPDFThetaBiasH(biasn2 - 1);
        G4double xaxisl = IPDFThetaBiasH.GetLowEdgeEnergy(size_t(biasn2 - 1));
        G4double xaxisu = IPDFThetaBiasH.GetLowEdgeEnergy(size_t(biasn2));
        G4double NatProb = xaxisu - xaxisl;
        w[3] = NatProb / w[3];
        if (verbosityLevel >= 1)
          G4cout << "Theta bin weight " << w[3] << " " << rndm
          << G4endl;
        return (IPDFThetaBiasH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenRandX()

G4double G4SPSRandomGenerator::GenRandX

(

)

Definition at line 250 of file G4SPSRandomGenerator.cc.

                                        {
    if (verbosityLevel >= 1)
        G4cout << "In GenRandX" << G4endl;
    if (XBias == false) {
        // X is not biased
        G4double rndm = G4UniformRand();
        return (rndm);
    } else {
        // X is biased
        //This is shared among threads, and we need to initialize
        //only once. Multiple instances of this class can exists
        //so we rely on a class-private, thread-private variable
        //to check if we need an initialiation. We do not use a lock here
        //because the boolean on which we check is thread private
        if ( local_IPDFXBias.Get().val == false ) {
            //For time that this thread arrived, here
            //Now two cases are possible: it is the first time
             //ANY thread has ever initialized this.
            //Now we need a lock. In any case, the thread local
            //variable can now be set to true
            local_IPDFXBias.Get().val = true;
            G4AutoLock l(&mutex);
            if (IPDFXBias == false) {
                // IPDF has not been created, so create it
                G4double bins[1024], vals[1024], sum;
                G4int ii;
                G4int maxbin = G4int(XBiasH.GetVectorLength());
                bins[0] = XBiasH.GetLowEdgeEnergy(size_t(0));
                vals[0] = XBiasH(size_t(0));
                sum = vals[0];
                for (ii = 1; ii < maxbin; ii++) {
                    bins[ii] = XBiasH.GetLowEdgeEnergy(size_t(ii));
                    vals[ii] = XBiasH(size_t(ii)) + vals[ii - 1];
                    sum = sum + XBiasH(size_t(ii));
                }

                for (ii = 0; ii < maxbin; ii++) {
                    vals[ii] = vals[ii] / sum;
                    IPDFXBiasH.InsertValues(bins[ii], vals[ii]);
                }
                // Make IPDFXBias = true
                IPDFXBias = true;
            }
        }
          // IPDF has been create so carry on
          G4double rndm = G4UniformRand();

          // Calculate the weighting: Find the bin that the determined
          // rndm is in and the weigthing will be the difference in the
          // natural probability (from the x-axis) divided by the
          // difference in the biased probability (the area).
          size_t numberOfBin = IPDFXBiasH.GetVectorLength();
          G4int biasn1 = 0;
          G4int biasn2 = numberOfBin / 2;
          G4int biasn3 = numberOfBin - 1;
          while (biasn1 != biasn3 - 1) {
              if (rndm > IPDFXBiasH(biasn2))
                biasn1 = biasn2;
              else
                biasn3 = biasn2;
              biasn2 = biasn1 + (biasn3 - biasn1 + 1) / 2;
          }
          // retrieve the areas and then the x-axis values
          bweights_t& w = bweights.Get();
          w[0] = IPDFXBiasH(biasn2) - IPDFXBiasH(biasn2 - 1);
          G4double xaxisl = IPDFXBiasH.GetLowEdgeEnergy(size_t(biasn2 - 1));
          G4double xaxisu = IPDFXBiasH.GetLowEdgeEnergy(size_t(biasn2));
          G4double NatProb = xaxisu - xaxisl;
          //G4cout << "X Bin weight " << bweights[0] << " " << rndm << G4endl;
          //G4cout << "lower and upper xaxis vals "<<xaxisl<<" "<<xaxisu<<G4endl;
          w[0] = NatProb / w[0];
          if (verbosityLevel >= 1)
            G4cout << "X bin weight " << w[0] << " " << rndm << G4endl;
          return (IPDFXBiasH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenRandY()

G4double G4SPSRandomGenerator::GenRandY

(

)

Definition at line 327 of file G4SPSRandomGenerator.cc.

                                        {
    if (verbosityLevel >= 1)
        G4cout << "In GenRandY" << G4endl;
    if (YBias == false) {
        // Y is not biased
        G4double rndm = G4UniformRand();
        return (rndm);
    } else {
        // Y is biased
        if ( local_IPDFYBias.Get().val == false ) {
            local_IPDFYBias.Get().val = true;
            G4AutoLock l(&mutex);
        if (IPDFYBias == false) {
            // IPDF has not been created, so create it
            G4double bins[1024], vals[1024], sum;
            G4int ii;
            G4int maxbin = G4int(YBiasH.GetVectorLength());
            bins[0] = YBiasH.GetLowEdgeEnergy(size_t(0));
            vals[0] = YBiasH(size_t(0));
            sum = vals[0];
            for (ii = 1; ii < maxbin; ii++) {
                bins[ii] = YBiasH.GetLowEdgeEnergy(size_t(ii));
                vals[ii] = YBiasH(size_t(ii)) + vals[ii - 1];
                sum = sum + YBiasH(size_t(ii));
            }

            for (ii = 0; ii < maxbin; ii++) {
                vals[ii] = vals[ii] / sum;
                IPDFYBiasH.InsertValues(bins[ii], vals[ii]);
            }
            // Make IPDFYBias = true
            IPDFYBias = true;
        }
        }       // IPDF has been create so carry on
        G4double rndm = G4UniformRand();
        size_t numberOfBin = IPDFYBiasH.GetVectorLength();
        G4int biasn1 = 0;
        G4int biasn2 = numberOfBin / 2;
        G4int biasn3 = numberOfBin - 1;
        while (biasn1 != biasn3 - 1) {
            if (rndm > IPDFYBiasH(biasn2))
              biasn1 = biasn2;
            else
              biasn3 = biasn2;
            biasn2 = biasn1 + (biasn3 - biasn1 + 1) / 2;
        }
        bweights_t& w = bweights.Get();
        w[1] = IPDFYBiasH(biasn2) - IPDFYBiasH(biasn2 - 1);
        G4double xaxisl = IPDFYBiasH.GetLowEdgeEnergy(size_t(biasn2 - 1));
        G4double xaxisu = IPDFYBiasH.GetLowEdgeEnergy(size_t(biasn2));
        G4double NatProb = xaxisu - xaxisl;
        w[1] = NatProb / w[1];
        if (verbosityLevel >= 1)
          G4cout << "Y bin weight " << w[1] << " " << rndm << G4endl;
        return (IPDFYBiasH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenRandZ()

G4double G4SPSRandomGenerator::GenRandZ

(

)

Definition at line 385 of file G4SPSRandomGenerator.cc.

                                        {
    if (verbosityLevel >= 1)
        G4cout << "In GenRandZ" << G4endl;
    if (ZBias == false) {
        // Z is not biased
        G4double rndm = G4UniformRand();
        return (rndm);
    } else {
        // Z is biased
        if (local_IPDFZBias.Get().val == false ) {
            local_IPDFZBias.Get().val = true;
            G4AutoLock l(&mutex);
        if (IPDFZBias == false) {
            // IPDF has not been created, so create it
            G4double bins[1024], vals[1024], sum;
            G4int ii;
            G4int maxbin = G4int(ZBiasH.GetVectorLength());
            bins[0] = ZBiasH.GetLowEdgeEnergy(size_t(0));
            vals[0] = ZBiasH(size_t(0));
            sum = vals[0];
            for (ii = 1; ii < maxbin; ii++) {
                bins[ii] = ZBiasH.GetLowEdgeEnergy(size_t(ii));
                vals[ii] = ZBiasH(size_t(ii)) + vals[ii - 1];
                sum = sum + ZBiasH(size_t(ii));
            }

            for (ii = 0; ii < maxbin; ii++) {
                vals[ii] = vals[ii] / sum;
                IPDFZBiasH.InsertValues(bins[ii], vals[ii]);
            }
            // Make IPDFZBias = true
            IPDFZBias = true;
        }
        }
        // IPDF has been create so carry on
        G4double rndm = G4UniformRand();
        //      size_t weight_bin_no = IPDFZBiasH.FindValueBinLocation(rndm);
        size_t numberOfBin = IPDFZBiasH.GetVectorLength();
        G4int biasn1 = 0;
        G4int biasn2 = numberOfBin / 2;
        G4int biasn3 = numberOfBin - 1;
        while (biasn1 != biasn3 - 1) {
            if (rndm > IPDFZBiasH(biasn2))
              biasn1 = biasn2;
            else
              biasn3 = biasn2;
            biasn2 = biasn1 + (biasn3 - biasn1 + 1) / 2;
        }
        bweights_t& w = bweights.Get();
        w[2] = IPDFZBiasH(biasn2) - IPDFZBiasH(biasn2 - 1);
        G4double xaxisl = IPDFZBiasH.GetLowEdgeEnergy(size_t(biasn2 - 1));
        G4double xaxisu = IPDFZBiasH.GetLowEdgeEnergy(size_t(biasn2));
        G4double NatProb = xaxisu - xaxisl;
        w[2] = NatProb / w[2];
        if (verbosityLevel >= 1)
          G4cout << "Z bin weight " << w[2] << " " << rndm << G4endl;
        return (IPDFZBiasH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetBiasWeight()

G4double G4SPSRandomGenerator::GetBiasWeight

(

)

Definition at line 186 of file G4SPSRandomGenerator.cc.

                                             {
  bweights_t& w = bweights.Get();
  return w[0] * w[1] * w[2] * w[3]
    * w[4] * w[5] * w[6] * w[7]
    * w[8];
}

Here is the caller graph for this function:

ReSetHist()

void G4SPSRandomGenerator::ReSetHist

(

G4String

atype

)

Definition at line 203 of file G4SPSRandomGenerator.cc.

                                                   {
    G4AutoLock l(&mutex);
    if (atype == "biasx") {
        XBias = false;
        IPDFXBias = false;
        local_IPDFXBias.Get().val = false;
        XBiasH = IPDFXBiasH = ZeroPhysVector;
    } else if (atype == "biasy") {
        YBias = false;
        IPDFYBias = false;
                local_IPDFYBias.Get().val = false;
        YBiasH = IPDFYBiasH = ZeroPhysVector;
    } else if (atype == "biasz") {
        ZBias = false;
        IPDFZBias = false;
                local_IPDFZBias.Get().val = false;
        ZBiasH = IPDFZBiasH = ZeroPhysVector;
    } else if (atype == "biast") {
        ThetaBias = false;
        IPDFThetaBias = false;
                local_IPDFThetaBias.Get().val = false;
        ThetaBiasH = IPDFThetaBiasH = ZeroPhysVector;
    } else if (atype == "biasp") {
        PhiBias = false;
        IPDFPhiBias = false;
                local_IPDFPhiBias.Get().val = false;
        PhiBiasH = IPDFPhiBiasH = ZeroPhysVector;
    } else if (atype == "biase") {
        EnergyBias = false;
        IPDFEnergyBias = false;
                local_IPDFEnergyBias.Get().val = false;
        EnergyBiasH = IPDFEnergyBiasH = ZeroPhysVector;
    } else if (atype == "biaspt") {
        PosThetaBias = false;
        IPDFPosThetaBias = false;
        local_IPDFPosThetaBias.Get().val = false;
        PosThetaBiasH = IPDFPosThetaBiasH = ZeroPhysVector;
    } else if (atype == "biaspp") {
        PosPhiBias = false;
        IPDFPosPhiBias = false;
        local_IPDFPosPhiBias.Get().val = false;
        PosPhiBiasH = IPDFPosPhiBiasH = ZeroPhysVector;
    } else {
        G4cout << "Error, histtype not accepted " << G4endl;
    }
}

Here is the caller graph for this function:

SetEnergyBias()

void G4SPSRandomGenerator::SetEnergyBias

(

G4ThreeVector

input

)

Definition at line 155 of file G4SPSRandomGenerator.cc.

                                                            {
    G4AutoLock l(&mutex);
    G4double ehi, val;
    ehi = input.x();
    val = input.y();
    EnergyBiasH.InsertValues(ehi, val);
    EnergyBias = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetIntensityWeight()

void G4SPSRandomGenerator::SetIntensityWeight

(

G4double

weight

)

Definition at line 182 of file G4SPSRandomGenerator.cc.

                                                             {
  bweights.Get()[8] = weight;
}

Here is the caller graph for this function:

SetPhiBias()

void G4SPSRandomGenerator::SetPhiBias

(

G4ThreeVector

input

)

Definition at line 146 of file G4SPSRandomGenerator.cc.

                                                         {
    G4AutoLock l(&mutex);
    G4double ehi, val;
    ehi = input.x();
    val = input.y();
    PhiBiasH.InsertValues(ehi, val);
    PhiBias = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetPosPhiBias()

void G4SPSRandomGenerator::SetPosPhiBias

(

G4ThreeVector

input

)

Definition at line 173 of file G4SPSRandomGenerator.cc.

                                                            {
    G4AutoLock l(&mutex);
    G4double ehi, val;
    ehi = input.x();
    val = input.y();
    PosPhiBiasH.InsertValues(ehi, val);
    PosPhiBias = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetPosThetaBias()

void G4SPSRandomGenerator::SetPosThetaBias

(

G4ThreeVector

input

)

Definition at line 164 of file G4SPSRandomGenerator.cc.

                                                              {
    G4AutoLock l(&mutex);
    G4double ehi, val;
    ehi = input.x();
    val = input.y();
    PosThetaBiasH.InsertValues(ehi, val);
    PosThetaBias = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetThetaBias()

void G4SPSRandomGenerator::SetThetaBias

(

G4ThreeVector

input

)

Definition at line 137 of file G4SPSRandomGenerator.cc.

                                                           {
    G4AutoLock l(&mutex);
    G4double ehi, val;
    ehi = input.x();
    val = input.y();
    ThetaBiasH.InsertValues(ehi, val);
    ThetaBias = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetVerbosity()

void G4SPSRandomGenerator::SetVerbosity

(

G4int

a

)

Definition at line 193 of file G4SPSRandomGenerator.cc.

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

SetXBias()

void G4SPSRandomGenerator::SetXBias

(

G4ThreeVector

input

)

Definition at line 110 of file G4SPSRandomGenerator.cc.

                                                       {
    G4AutoLock l(&mutex);
    G4double ehi, val;
    ehi = input.x();
    val = input.y();
    XBiasH.InsertValues(ehi, val);
    XBias = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetYBias()

void G4SPSRandomGenerator::SetYBias

(

G4ThreeVector

input

)

Definition at line 119 of file G4SPSRandomGenerator.cc.

                                                       {
    G4AutoLock l(&mutex);
    G4double ehi, val;
    ehi = input.x();
    val = input.y();
    YBiasH.InsertValues(ehi, val);
    YBias = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetZBias()

void G4SPSRandomGenerator::SetZBias

(

G4ThreeVector

input

)

Definition at line 128 of file G4SPSRandomGenerator.cc.

                                                       {
    G4AutoLock l(&mutex);
    G4double ehi, val;
    ehi = input.x();
    val = input.y();
    ZBiasH.InsertValues(ehi, val);
    ZBias = true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

bweights

G4Cache<bweights_t> G4SPSRandomGenerator::bweights

private

Definition at line 240 of file G4SPSRandomGenerator.hh.

EnergyBias

G4bool G4SPSRandomGenerator::EnergyBias

private

Definition at line 222 of file G4SPSRandomGenerator.hh.

EnergyBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::EnergyBiasH

private

Definition at line 223 of file G4SPSRandomGenerator.hh.

IPDFEnergyBias

G4bool G4SPSRandomGenerator::IPDFEnergyBias

private

Definition at line 222 of file G4SPSRandomGenerator.hh.

IPDFEnergyBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::IPDFEnergyBiasH

private

Definition at line 224 of file G4SPSRandomGenerator.hh.

IPDFPhiBias

G4bool G4SPSRandomGenerator::IPDFPhiBias

private

Definition at line 218 of file G4SPSRandomGenerator.hh.

IPDFPhiBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::IPDFPhiBiasH

private

Definition at line 220 of file G4SPSRandomGenerator.hh.

IPDFPosPhiBias

G4bool G4SPSRandomGenerator::IPDFPosPhiBias

private

Definition at line 230 of file G4SPSRandomGenerator.hh.

IPDFPosPhiBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::IPDFPosPhiBiasH

private

Definition at line 232 of file G4SPSRandomGenerator.hh.

IPDFPosThetaBias

G4bool G4SPSRandomGenerator::IPDFPosThetaBias

private

Definition at line 226 of file G4SPSRandomGenerator.hh.

IPDFPosThetaBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::IPDFPosThetaBiasH

private

Definition at line 228 of file G4SPSRandomGenerator.hh.

IPDFThetaBias

G4bool G4SPSRandomGenerator::IPDFThetaBias

private

Definition at line 214 of file G4SPSRandomGenerator.hh.

IPDFThetaBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::IPDFThetaBiasH

private

Definition at line 216 of file G4SPSRandomGenerator.hh.

IPDFXBias

G4bool G4SPSRandomGenerator::IPDFXBias

private

Definition at line 202 of file G4SPSRandomGenerator.hh.

IPDFXBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::IPDFXBiasH

private

Definition at line 204 of file G4SPSRandomGenerator.hh.

IPDFYBias

G4bool G4SPSRandomGenerator::IPDFYBias

private

Definition at line 206 of file G4SPSRandomGenerator.hh.

IPDFYBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::IPDFYBiasH

private

Definition at line 208 of file G4SPSRandomGenerator.hh.

IPDFZBias

G4bool G4SPSRandomGenerator::IPDFZBias

private

Definition at line 210 of file G4SPSRandomGenerator.hh.

IPDFZBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::IPDFZBiasH

private

Definition at line 212 of file G4SPSRandomGenerator.hh.

local_IPDFEnergyBias

G4Cache<a_check> G4SPSRandomGenerator::local_IPDFEnergyBias

private

Definition at line 221 of file G4SPSRandomGenerator.hh.

local_IPDFPhiBias

G4Cache<a_check> G4SPSRandomGenerator::local_IPDFPhiBias

private

Definition at line 217 of file G4SPSRandomGenerator.hh.

local_IPDFPosPhiBias

G4Cache<a_check> G4SPSRandomGenerator::local_IPDFPosPhiBias

private

Definition at line 229 of file G4SPSRandomGenerator.hh.

local_IPDFPosThetaBias

G4Cache<a_check> G4SPSRandomGenerator::local_IPDFPosThetaBias

private

Definition at line 225 of file G4SPSRandomGenerator.hh.

local_IPDFThetaBias

G4Cache<a_check> G4SPSRandomGenerator::local_IPDFThetaBias

private

Definition at line 213 of file G4SPSRandomGenerator.hh.

local_IPDFXBias

G4Cache<a_check> G4SPSRandomGenerator::local_IPDFXBias

private

Definition at line 201 of file G4SPSRandomGenerator.hh.

local_IPDFYBias

G4Cache<a_check> G4SPSRandomGenerator::local_IPDFYBias

private

Definition at line 205 of file G4SPSRandomGenerator.hh.

local_IPDFZBias

G4Cache<a_check> G4SPSRandomGenerator::local_IPDFZBias

private

Definition at line 209 of file G4SPSRandomGenerator.hh.

mutex

G4Mutex G4SPSRandomGenerator::mutex

private

Definition at line 246 of file G4SPSRandomGenerator.hh.

PhiBias

G4bool G4SPSRandomGenerator::PhiBias

private

Definition at line 218 of file G4SPSRandomGenerator.hh.

PhiBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::PhiBiasH

private

Definition at line 219 of file G4SPSRandomGenerator.hh.

PosPhiBias

G4bool G4SPSRandomGenerator::PosPhiBias

private

Definition at line 230 of file G4SPSRandomGenerator.hh.

PosPhiBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::PosPhiBiasH

private

Definition at line 231 of file G4SPSRandomGenerator.hh.

PosThetaBias

G4bool G4SPSRandomGenerator::PosThetaBias

private

Definition at line 226 of file G4SPSRandomGenerator.hh.

PosThetaBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::PosThetaBiasH

private

Definition at line 227 of file G4SPSRandomGenerator.hh.

ThetaBias

G4bool G4SPSRandomGenerator::ThetaBias

private

Definition at line 214 of file G4SPSRandomGenerator.hh.

ThetaBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::ThetaBiasH

private

Definition at line 215 of file G4SPSRandomGenerator.hh.

verbosityLevel

G4int G4SPSRandomGenerator::verbosityLevel

private

Definition at line 244 of file G4SPSRandomGenerator.hh.

XBias

G4bool G4SPSRandomGenerator::XBias

private

Definition at line 202 of file G4SPSRandomGenerator.hh.

XBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::XBiasH

private

Definition at line 203 of file G4SPSRandomGenerator.hh.

YBias

G4bool G4SPSRandomGenerator::YBias

private

Definition at line 206 of file G4SPSRandomGenerator.hh.

YBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::YBiasH

private

Definition at line 207 of file G4SPSRandomGenerator.hh.

ZBias

G4bool G4SPSRandomGenerator::ZBias

private

Definition at line 210 of file G4SPSRandomGenerator.hh.

ZBiasH

G4PhysicsOrderedFreeVector G4SPSRandomGenerator::ZBiasH

private

Definition at line 211 of file G4SPSRandomGenerator.hh.

---

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

• G4SPSRandomGenerator.hh
• G4SPSRandomGenerator.cc