G4SPSRandomGenerator.hh

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4SPSRandomGenerator.hh

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//
// MODULE:        G4SPSRandomGenerator.hh
//
// Version:      1.0
// Date:         5/02/04
// Author:       Fan Lei 
// Organisation: QinetiQ ltd.
// Customer:     ESA/ESTEC
//
//
// CHANGE HISTORY
// --------------
// 06/06/2014 A Dotti
//    Note on thread safety: added a mutex to protect access to shared
//    resources (data members).
//    Getters and Setters are mutex'd but not the GetRand* methods,
//    because it is assumed these are called only during the event loop
//    during which the status of this class is invariant.
//
// 26/10/2004 F Lei
//    Created separated the theta, phi generators for position distributions.
//
// Version 1.0, 05/02/2004, Fan Lei, Created.
//    Based on the G4GeneralParticleSource class in Geant4 v6.0
//
//
// Class Description:
//
// Special random number generator used by G4GeneralParticleSource to allow 
// biasing applied at the lowest level for all distributions.
//
//
// MEMBER FUNCTIONS
// ----------------
//
// G4SPSRandomGenerator ()
//    Constructor: Initializes variables
//
// ~G4SPSRandomGenerator ()
//    Destructor: 
//
// void SetXBias(G4ThreeVector)
//    Allows the user to re-distribute the random
//    numbers used to generate x co-ordinates.
//
// void SetYBias(G4ThreeVector)
//    Allows the user to re-distribute the random
//    numbers used to generate y co-ordinates.
//
// void SetZBias(G4ThreeVector)
//    Allows the user to re-distribute the random
//    numbers used to generate z co-ordinates.
//
// void SetThetaBias(G4ThreeVector)
//    Allows the user to re-distribute the random
//    numbers used to generate values of theta.
//
// void SetPhiBias(G4ThreeVector)
//    Allows the user to re-distribute the random
//    numbers used to generate values of phi.
//
// void SetPosThetaBias(G4ThreeVector)
//    Allows the user to re-distribute the random
//    numbers used to generate values of theta for position distribution.
//
// void SetPosPhiBias(G4ThreeVector)
//    Allows the user to re-distribute the random
//    numbers used to generate values of phi for position distribution.
//
// void SetEnergyBias(G4ThreeVector)
//    Allows the user to re-distribute the random
//    numbers used to generate the energies.
//
// G4double GenRandX()
//    Generates the random number for x, with or without biasing.
//
// G4double GenRandY()
//    Generates the random number for y, with or without biasing.
//
// G4double GenRandZ()
//    Generates the random number for z, with or without biasing.
//
// G4double GenRandTheta()
//    Generates the random number for theta, with or without biasing.
//
// G4double GenRandPhi()
//    Generates the random number for phi, with or without biasing.
//
// G4double GenRandEnergy()
//    Generates the random number for energy, with or without biasing.
//
// G4double GenRandPosTheta()
//    Generates the random number for theta, with or without biasing for position distribution.
//
// G4double GenRandPosPhi()
//    Generates the random number for phi, with or without biasing for position distribution.
//
//  inline G4double GetBiasWeight()
//    Returns the weight change after biasing
// 
//  void ReSetHist(G4String);
//    Re-sets the histogram for user defined distribution
//
// void SetVerbosity(G4int)
//    Sets the verbosity level.
//
//
#ifndef G4SPSRandomGenerator_h
#define G4SPSRandomGenerator_h 1

#include "G4PhysicsOrderedFreeVector.hh"
#include "G4DataInterpolation.hh"
#include "G4ThreeVector.hh"
#include "G4Threading.hh"
#include "G4Cache.hh"

class G4SPSRandomGenerator {
public:
  G4SPSRandomGenerator();
  ~G4SPSRandomGenerator();

  //  static G4SPSRandomGenerator* getInstance ();

  // Biasing Methods
  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();

  // method to re-set the histograms
  void ReSetHist(G4String);

  // Set the verbosity level.
  void SetVerbosity(G4int a);

private:
  //Encapsulate in a struct
  //to gurantee that correct
  //initial state is set via constructor
        struct a_check {
          G4bool val;
          a_check() { val = false; }
        };
        //See .cc for an explanation of this
        //in method GenRandX()
  G4Cache<a_check> local_IPDFXBias;
  G4bool XBias, IPDFXBias;
  G4PhysicsOrderedFreeVector XBiasH;
  G4PhysicsOrderedFreeVector IPDFXBiasH;
        G4Cache<a_check> local_IPDFYBias;
  G4bool YBias, IPDFYBias;
  G4PhysicsOrderedFreeVector YBiasH;
  G4PhysicsOrderedFreeVector IPDFYBiasH;
        G4Cache<a_check> local_IPDFZBias;
  G4bool ZBias, IPDFZBias;
  G4PhysicsOrderedFreeVector ZBiasH;
  G4PhysicsOrderedFreeVector IPDFZBiasH;
        G4Cache<a_check> local_IPDFThetaBias;
  G4bool ThetaBias, IPDFThetaBias;
  G4PhysicsOrderedFreeVector ThetaBiasH;
  G4PhysicsOrderedFreeVector IPDFThetaBiasH;
        G4Cache<a_check> local_IPDFPhiBias;
  G4bool PhiBias, IPDFPhiBias;
  G4PhysicsOrderedFreeVector PhiBiasH;
  G4PhysicsOrderedFreeVector IPDFPhiBiasH;
        G4Cache<a_check> local_IPDFEnergyBias;
  G4bool EnergyBias, IPDFEnergyBias;
  G4PhysicsOrderedFreeVector EnergyBiasH;
  G4PhysicsOrderedFreeVector IPDFEnergyBiasH;
        G4Cache<a_check> local_IPDFPosThetaBias;
  G4bool PosThetaBias, IPDFPosThetaBias;
  G4PhysicsOrderedFreeVector PosThetaBiasH;
  G4PhysicsOrderedFreeVector IPDFPosThetaBiasH;
        G4Cache<a_check> local_IPDFPosPhiBias;
  G4bool PosPhiBias, IPDFPosPhiBias;
  G4PhysicsOrderedFreeVector PosPhiBiasH;
  G4PhysicsOrderedFreeVector IPDFPosPhiBiasH;

  //G4double alpha;   // for biasing energy
  struct bweights_t {
    G4double w[9];
    bweights_t();
    G4double& operator[] (const int i);
  };
  G4Cache<bweights_t> bweights;
  //G4double bweights[9]; //record x,y,z,theta,phi,energy,posThet,posPhi,intensity weights

  // Verbosity
  G4int verbosityLevel;

    G4Mutex mutex; //protect shared resources
};

#endif