HistoManager Class Reference

#include <HistoManager.hh>

Collaboration diagram for HistoManager:

Public Member Functions
	~HistoManager ()
void	BeginOfRun ()
void	EndOfRun ()
void	SetVerbose (G4int val)
void	SetParticleName (const G4String &)
void	SetElementName (const G4String &)
void	SetNumberOfBinsE (G4int val)
void	SetNumberOfBinsP (G4int val)
void	SetMinKinEnergy (G4double val)
void	SetMaxKinEnergy (G4double val)
void	SetMinMomentum (G4double val)
void	SetMaxMomentum (G4double val)
G4int	GetVerbose () const

Static Public Member Functions
static HistoManager *	GetPointer ()

Private Member Functions
	HistoManager ()

Private Attributes
Histo *	fHisto
const G4ParticleDefinition *	fNeutron
G4String	fParticleName
G4String	fElementName
G4double	fMinKinEnergy
G4double	fMaxKinEnergy
G4double	fMinMomentum
G4double	fMaxMomentum
G4int	fVerbose
G4int	fBinsE
G4int	fBinsP
G4bool	fIsInitialised

Static Private Attributes
static HistoManager *	fManager = 0

Detailed Description

Definition at line 60 of file HistoManager.hh.

Constructor & Destructor Documentation

HistoManager::HistoManager ( )

private

Definition at line 81 of file HistoManager.cc.

References fBinsE, fBinsP, fElementName, fHisto, fIsInitialised, fMaxKinEnergy, fMaxMomentum, fMinKinEnergy, fMinMomentum, fNeutron, fParticleName, fVerbose, and Histo::SetVerbose().

{
  fHisto     = new Histo();
  fNeutron   = G4Neutron::Neutron();
  fVerbose   = 1;

  fParticleName  = "proton";
  fElementName   = "Al";

  fMinKinEnergy  = 0.1*MeV;
  fMaxKinEnergy  = 10*TeV;
  fMinMomentum   = 1*MeV;
  fMaxMomentum   = 10*TeV;

  fBinsE    = 800;
  fBinsP    = 700;

  fHisto->SetVerbose(fVerbose);

  fIsInitialised = false;
}

Here is the call graph for this function:

HistoManager::~HistoManager

(

)

Definition at line 105 of file HistoManager.cc.

References fHisto.

{
  delete fHisto;
}

Member Function Documentation

void HistoManager::BeginOfRun

(

)

Definition at line 112 of file HistoManager.cc.

References Histo::Add1D(), Histo::Book(), fBinsE, fBinsP, fHisto, fIsInitialised, fMaxKinEnergy, fMaxMomentum, fMinKinEnergy, fMinMomentum, and Histo::SetHisto1D().

{
  G4double p1 = std::log10(fMinMomentum/GeV);
  G4double p2 = std::log10(fMaxMomentum/GeV);
  G4double e1 = std::log10(fMinKinEnergy/MeV);
  G4double e2 = std::log10(fMaxKinEnergy/MeV);

  //G4cout<<"e1= "<<e1<<" e2= "<<e2<<" p1= "<<p1<<" p2= "<<p2<<G4endl;

  if(fIsInitialised) {
    fHisto->SetHisto1D(0,fBinsP,p1,p2,1.0);
    fHisto->SetHisto1D(1,fBinsE,e1,e2,1.0);
    fHisto->SetHisto1D(2,fBinsP,p1,p2,1.0);
    fHisto->SetHisto1D(3,fBinsE,e1,e2,1.0);
    fHisto->SetHisto1D(4,fBinsE,e1,e2,1.0);
    fHisto->SetHisto1D(5,fBinsE,e1,e2,1.0);
    fHisto->SetHisto1D(6,fBinsE,e1,e2,1.0);
    fHisto->SetHisto1D(7,fBinsE,e1,e2,1.0);

  } else {
    fHisto->Add1D("h1","Elastic cross section (barn,1.0) as a functions of log10(p/GeV)",
                 fBinsP,p1,p2,1.0);
    fHisto->Add1D("h2","Elastic cross section (barn) as a functions of log10(E/MeV)",
                 fBinsE,e1,e2,1.0);
    fHisto->Add1D("h3","Inelastic cross section (barn) as a functions of log10(p/GeV)",
                 fBinsP,p1,p2,1.0);
    fHisto->Add1D("h4","Inelastic cross section (barn) as a functions of log10(E/MeV)",
                 fBinsE,e1,e2,1.0);
    fHisto->Add1D("h5","Capture cross section (barn) as a functions of log10(E/MeV)",
                 fBinsE,e1,e2,1.0);
    fHisto->Add1D("h6","Fission cross section (barn) as a functions of log10(E/MeV)",
                 fBinsE,e1,e2,1.0);
    fHisto->Add1D("h7","Charge exchange cross section (barn) as a functions of log10(E/MeV)",
                 fBinsE,e1,e2,1.0);
    fHisto->Add1D("h8","Total cross section (barn) as a functions of log10(E/MeV)",
                 fBinsE,e1,e2,1.0);
  }

  fIsInitialised = true;
  fHisto->Book();  
}

Here is the call graph for this function:

void HistoManager::EndOfRun

(

)

Definition at line 156 of file HistoManager.cc.

References fBinsE, fBinsP, fElementName, fHisto, Histo::Fill(), fMaxKinEnergy, fMaxMomentum, fMinKinEnergy, fMinMomentum, fNeutron, fParticleName, fVerbose, and Histo::Save().

{
  if(fVerbose > 0) {
    G4cout << "HistoManager: End of run actions are started" << G4endl;
  }

  const G4Element* elm = 
    G4NistManager::Instance()->FindOrBuildElement(fElementName);
  const G4Material* mat = 
    G4NistManager::Instance()->FindOrBuildMaterial("G4_"+fElementName);
  const G4ParticleDefinition* particle = 
    G4ParticleTable::GetParticleTable()->FindParticle(fParticleName);

  G4cout << "### Fill Cross Sections for " << fParticleName 
         << " off " << fElementName
         << G4endl;
  if(fVerbose > 0) {
    G4cout << "------------------------------------------------------------------------" 
           << G4endl;
    G4cout << "    N     E(MeV)   Elastic(b)   Inelastic(b)";
    if(particle == fNeutron) { G4cout << " Capture(b)   Fission(b)"; }
    G4cout << "   Total(b)" << G4endl;     
    G4cout << "------------------------------------------------------------------------" 
           << G4endl;
  }
  if(!particle || !elm) {
    G4cout << "HistoManager WARNING Particle or element undefined" << G4endl;
    return;
  }

  G4int prec = G4cout.precision();
  G4cout.precision(4);

  G4HadronicProcessStore* store = G4HadronicProcessStore::Instance();
  G4double mass = particle->GetPDGMass();

  // Build histograms

  G4double p1 = std::log10(fMinMomentum/GeV);
  G4double p2 = std::log10(fMaxMomentum/GeV);
  G4double e1 = std::log10(fMinKinEnergy/MeV);
  G4double e2 = std::log10(fMaxKinEnergy/MeV);
  G4double de = (e2 - e1)/G4double(fBinsE);
  G4double dp = (p2 - p1)/G4double(fBinsP);

  G4double x  = e1 - de*0.5; 
  G4double e, p, xs, xtot;
  G4int i;
  for(i=0; i<fBinsE; i++) {
    x += de;
    e  = std::pow(10.,x)*MeV;
    if(fVerbose>0) G4cout << std::setw(5) << i << std::setw(12) << e;  
    xs = store->GetElasticCrossSectionPerAtom(particle,e,elm,mat);
    xtot = xs;
    if(fVerbose>0) G4cout << std::setw(12) << xs/barn;  
    fHisto->Fill(1, x, xs/barn);    
    xs = store->GetInelasticCrossSectionPerAtom(particle,e,elm,mat);
    xtot += xs;
    if(fVerbose>0) G4cout << " " << std::setw(12) << xs/barn;  
    fHisto->Fill(3, x, xs/barn);    
    if(particle == fNeutron) {
      xs = store->GetCaptureCrossSectionPerAtom(particle,e,elm,mat);
      xtot += xs;
      if(fVerbose>0) G4cout << " " << std::setw(12) << xs/barn;  
      fHisto->Fill(4, x, xs/barn);    
      xs = store->GetFissionCrossSectionPerAtom(particle,e,elm,mat);
      xtot += xs;
      if(fVerbose>0) G4cout << " " << std::setw(12) << xs/barn;  
      fHisto->Fill(5, x, xs/barn);    
    }
    xs = store->GetChargeExchangeCrossSectionPerAtom(particle,e,elm,mat);
    if(fVerbose>0) G4cout << " " << std::setw(12) << xtot/barn << G4endl;   
    fHisto->Fill(6, x, xs/barn);    
    fHisto->Fill(7, x, xtot/barn);    
  }

  x = p1 - dp*0.5; 
  for(i=0; i<fBinsP; i++) {
    x += dp;
    p  = std::pow(10.,x)*GeV;
    e  = std::sqrt(p*p + mass*mass) - mass;
    xs = store->GetElasticCrossSectionPerAtom(particle,e,elm,mat);
    fHisto->Fill(0, x, xs/barn);    
    xs = store->GetInelasticCrossSectionPerAtom(particle,e,elm,mat);
    fHisto->Fill(2, x, xs/barn); 
  }
  if(fVerbose > 0) {
    G4cout << "-----------------------------------------------------------------" 
           << G4endl;
  }
  G4cout.precision(prec);
  fHisto->Save();
}

Here is the call graph for this function:

HistoManager * HistoManager::GetPointer ( )

static

Definition at line 70 of file HistoManager.cc.

References fManager.

Referenced by RunAction::BeginOfRunAction(), RunAction::EndOfRunAction(), and DetectorMessenger::SetNewValue().

{
  if(!fManager) {
    static HistoManager manager;
    fManager = &manager;
  }
  return fManager;
}

Here is the caller graph for this function:

G4int HistoManager::GetVerbose ( ) const

inline

Definition at line 156 of file HistoManager.hh.

References fVerbose.

{
  return fVerbose;
}

void HistoManager::SetElementName ( const G4String &  name )

inline

Definition at line 121 of file HistoManager.hh.

References fElementName.

Referenced by DetectorMessenger::SetNewValue().

{
  fElementName = name;
}

Here is the caller graph for this function:

void HistoManager::SetMaxKinEnergy ( G4double  val )

inline

Definition at line 141 of file HistoManager.hh.

References fMaxKinEnergy, and fMinKinEnergy.

Referenced by DetectorMessenger::SetNewValue().

{
  if(val>fMinKinEnergy) { fMaxKinEnergy = val; }
}

Here is the caller graph for this function:

void HistoManager::SetMaxMomentum ( G4double  val )

inline

Definition at line 151 of file HistoManager.hh.

References fMaxMomentum, and fMinMomentum.

Referenced by DetectorMessenger::SetNewValue().

{
  if(val>fMinMomentum) { fMaxMomentum = val; }
}

Here is the caller graph for this function:

void HistoManager::SetMinKinEnergy ( G4double  val )

inline

Definition at line 136 of file HistoManager.hh.

References fMaxKinEnergy, and fMinKinEnergy.

Referenced by DetectorMessenger::SetNewValue().

{
  if(val>0 && val<fMaxKinEnergy) { fMinKinEnergy = val; }
}

Here is the caller graph for this function:

void HistoManager::SetMinMomentum ( G4double  val )

inline

Definition at line 146 of file HistoManager.hh.

References fMaxMomentum, and fMinMomentum.

Referenced by DetectorMessenger::SetNewValue().

{
  if(val>0 && val<fMaxMomentum) { fMinMomentum = val; }
}

Here is the caller graph for this function:

void HistoManager::SetNumberOfBinsE ( G4int  val )

inline

Definition at line 126 of file HistoManager.hh.

References fBinsE.

Referenced by DetectorMessenger::SetNewValue().

{
  if(val>0) { fBinsE = val; } 
}

Here is the caller graph for this function:

void HistoManager::SetNumberOfBinsP ( G4int  val )

inline

Definition at line 131 of file HistoManager.hh.

References fBinsP.

Referenced by DetectorMessenger::SetNewValue().

{
  if(val>0) { fBinsP = val; } 
}

Here is the caller graph for this function:

void HistoManager::SetParticleName ( const G4String &  name )

inline

Definition at line 116 of file HistoManager.hh.

References fParticleName.

Referenced by DetectorMessenger::SetNewValue().

{
  fParticleName = name;
}

Here is the caller graph for this function:

void HistoManager::SetVerbose

(

G4int

val

)

Definition at line 252 of file HistoManager.cc.

References fHisto, fVerbose, and Histo::SetVerbose().

Referenced by DetectorMessenger::SetNewValue().

{
  fVerbose = val; 
  fHisto->SetVerbose(val);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Field Documentation

G4int HistoManager::fBinsE

private

Definition at line 110 of file HistoManager.hh.

Referenced by BeginOfRun(), EndOfRun(), HistoManager(), and SetNumberOfBinsE().

G4int HistoManager::fBinsP

private

Definition at line 111 of file HistoManager.hh.

Referenced by BeginOfRun(), EndOfRun(), HistoManager(), and SetNumberOfBinsP().

G4String HistoManager::fElementName

private

Definition at line 102 of file HistoManager.hh.

Referenced by EndOfRun(), HistoManager(), and SetElementName().

Histo* HistoManager::fHisto

private

Definition at line 97 of file HistoManager.hh.

Referenced by BeginOfRun(), EndOfRun(), HistoManager(), SetVerbose(), and ~HistoManager().

G4bool HistoManager::fIsInitialised

private

Definition at line 113 of file HistoManager.hh.

Referenced by BeginOfRun(), and HistoManager().

HistoManager * HistoManager::fManager = 0

staticprivate

Definition at line 95 of file HistoManager.hh.

Referenced by GetPointer().

G4double HistoManager::fMaxKinEnergy

private

Definition at line 105 of file HistoManager.hh.

Referenced by BeginOfRun(), EndOfRun(), HistoManager(), SetMaxKinEnergy(), and SetMinKinEnergy().

G4double HistoManager::fMaxMomentum

private

Definition at line 107 of file HistoManager.hh.

Referenced by BeginOfRun(), EndOfRun(), HistoManager(), SetMaxMomentum(), and SetMinMomentum().

G4double HistoManager::fMinKinEnergy

private

Definition at line 104 of file HistoManager.hh.

Referenced by BeginOfRun(), EndOfRun(), HistoManager(), SetMaxKinEnergy(), and SetMinKinEnergy().

G4double HistoManager::fMinMomentum

private

Definition at line 106 of file HistoManager.hh.

Referenced by BeginOfRun(), EndOfRun(), HistoManager(), SetMaxMomentum(), and SetMinMomentum().

const G4ParticleDefinition* HistoManager::fNeutron

private

Definition at line 99 of file HistoManager.hh.

Referenced by EndOfRun(), and HistoManager().

G4String HistoManager::fParticleName

private

Definition at line 101 of file HistoManager.hh.

Referenced by EndOfRun(), HistoManager(), and SetParticleName().

G4int HistoManager::fVerbose

private

Definition at line 109 of file HistoManager.hh.

Referenced by EndOfRun(), GetVerbose(), HistoManager(), and SetVerbose().

---

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

• HistoManager.hh
• HistoManager.cc