G4TablesForExtrapolator Class Reference

#include <G4TablesForExtrapolator.hh>

Public Member Functions
	G4TablesForExtrapolator (G4int verb, G4int bins, G4double e1, G4double e2)
	~G4TablesForExtrapolator ()
const G4PhysicsTable *	GetPhysicsTable (ExtTableType type) const
void	Initialisation ()
G4TablesForExtrapolator &	operator= (const G4TablesForExtrapolator &right)=delete
	G4TablesForExtrapolator (const G4TablesForExtrapolator &)=delete

Detailed Description

Definition at line 74 of file G4TablesForExtrapolator.hh.

Constructor & Destructor Documentation

G4TablesForExtrapolator() [1/2]

G4TablesForExtrapolator::G4TablesForExtrapolator ( G4int verb, G4int bins, G4double e1, G4double e2 )

explicit

Definition at line 74 of file G4TablesForExtrapolator.cc.

  : emin(e1), emax(e2), verbose(verb), nbins(bins)
{
  electron = G4Electron::Electron();
  positron = G4Positron::Positron();
  proton   = G4Proton::Proton();
  muonPlus = G4MuonPlus::MuonPlus();
  muonMinus= G4MuonMinus::MuonMinus();
}

Referenced by G4TablesForExtrapolator(), and operator=().

~G4TablesForExtrapolator()

G4TablesForExtrapolator::~G4TablesForExtrapolator

(

)

Definition at line 87 of file G4TablesForExtrapolator.cc.

{
  if(nullptr != dedxElectron) {
    dedxElectron->clearAndDestroy();
    delete dedxElectron;
  }
  if(nullptr != dedxPositron) {
    dedxPositron->clearAndDestroy();
    delete dedxPositron;
  }
  if(nullptr != dedxProton) {
    dedxProton->clearAndDestroy();
    delete dedxProton;
  }
  if(nullptr != dedxMuon) {
    dedxMuon->clearAndDestroy();
    delete dedxMuon;
  }
  if(nullptr != rangeElectron) {
    rangeElectron->clearAndDestroy();
    delete rangeElectron;
  }
  if(nullptr != rangePositron) {
    rangePositron->clearAndDestroy();
    delete rangePositron;
  }
  if(nullptr != rangeProton) {
    rangeProton->clearAndDestroy();
    delete rangeProton;
  }
  if(nullptr != rangeMuon) {
    rangeMuon->clearAndDestroy();
    delete rangeMuon;
  }
  if(nullptr != invRangeElectron) {
    invRangeElectron->clearAndDestroy();
    delete invRangeElectron;
  }
  if(nullptr != invRangePositron) {
    invRangePositron->clearAndDestroy();
    delete invRangePositron;
  }
  if(nullptr != invRangeProton) {
    invRangeProton->clearAndDestroy();
    delete invRangeProton;
  }
  if(nullptr != invRangeMuon) {
    invRangeMuon->clearAndDestroy();
    delete invRangeMuon;
  }
  if(nullptr != mscElectron) {
    mscElectron->clearAndDestroy();
    delete mscElectron;
  }
  delete pcuts;
  delete builder;
}

G4TablesForExtrapolator() [2/2]

G4TablesForExtrapolator::G4TablesForExtrapolator ( const G4TablesForExtrapolator & )

delete

Member Function Documentation

GetPhysicsTable()

const G4PhysicsTable * G4TablesForExtrapolator::GetPhysicsTable

(

ExtTableType

type

)

const

Definition at line 148 of file G4TablesForExtrapolator.cc.

{
  const G4PhysicsTable* table = nullptr;
  switch (type) 
    {
    case fDedxElectron:
      table = dedxElectron;
      break;
    case fDedxPositron:
      table = dedxPositron;
      break;
    case fDedxProton:
      table = dedxProton;
      break;
    case fDedxMuon:
      table = dedxMuon;
      break;
    case fRangeElectron:
      table = rangeElectron;
      break;
    case fRangePositron:
      table = rangePositron;
      break;
    case fRangeProton:
      table = rangeProton;
      break;
    case fRangeMuon:
      table = rangeMuon;
      break;
    case fInvRangeElectron:
      table = invRangeElectron;
      break;
    case fInvRangePositron:
      table = invRangePositron;
      break;
    case fInvRangeProton:
      table = invRangeProton;
      break;
    case fInvRangeMuon:
      table = invRangeMuon;
      break;
    case fMscElectron:
      table = mscElectron;
    }
  return table;
}

Initialisation()

void G4TablesForExtrapolator::Initialisation

(

)

Definition at line 197 of file G4TablesForExtrapolator.cc.

{
  if(verbose>1) {
    G4cout << "### G4TablesForExtrapolator::Initialisation" << G4endl;
  }
  G4int num = (G4int)G4Material::GetNumberOfMaterials();
  if(nmat == num) { return; }
  nmat = num;
  cuts.resize(nmat, DBL_MAX);
  couples.resize(nmat, nullptr);
 
  const G4MaterialTable* mtable = G4Material::GetMaterialTable();
  if(!pcuts) { pcuts = new G4ProductionCuts(); }
 
  for(G4int i=0; i<nmat; ++i) {
    couples[i] = new G4MaterialCutsCouple((*mtable)[i],pcuts);
  }
 
  dedxElectron     = PrepareTable(dedxElectron);
  dedxPositron     = PrepareTable(dedxPositron);
  dedxMuon         = PrepareTable(dedxMuon);
  dedxProton       = PrepareTable(dedxProton);
  rangeElectron    = PrepareTable(rangeElectron);
  rangePositron    = PrepareTable(rangePositron);
  rangeMuon        = PrepareTable(rangeMuon);
  rangeProton      = PrepareTable(rangeProton);
  invRangeElectron = PrepareTable(invRangeElectron);
  invRangePositron = PrepareTable(invRangePositron);
  invRangeMuon     = PrepareTable(invRangeMuon);
  invRangeProton   = PrepareTable(invRangeProton);
  mscElectron      = PrepareTable(mscElectron);
 
  builder = new G4LossTableBuilder(true);
  builder->SetBaseMaterialActive(false);
 
  if(verbose>1) {
    G4cout << "### G4TablesForExtrapolator Builds electron tables" 
       << G4endl;
  }
  ComputeElectronDEDX(electron, dedxElectron);
  builder->BuildRangeTable(dedxElectron,rangeElectron);  
  builder->BuildInverseRangeTable(rangeElectron, invRangeElectron);  
 
  if(verbose>1) {
    G4cout << "### G4TablesForExtrapolator Builds positron tables" 
       << G4endl;
  }
  ComputeElectronDEDX(positron, dedxPositron);
  builder->BuildRangeTable(dedxPositron, rangePositron);  
  builder->BuildInverseRangeTable(rangePositron, invRangePositron);  
 
  if(verbose>1) {
    G4cout << "### G4TablesForExtrapolator Builds muon tables" << G4endl;
  }
  ComputeMuonDEDX(muonPlus, dedxMuon);
  builder->BuildRangeTable(dedxMuon, rangeMuon);  
  builder->BuildInverseRangeTable(rangeMuon, invRangeMuon);  
  
  if(verbose>2) {
    G4cout << "DEDX MUON" << G4endl;
    G4cout << *dedxMuon << G4endl;
    G4cout << "RANGE MUON" << G4endl;
    G4cout << *rangeMuon << G4endl;
    G4cout << "INVRANGE MUON" << G4endl;
    G4cout << *invRangeMuon << G4endl;
  }
  if(verbose>1) {
    G4cout << "### G4TablesForExtrapolator Builds proton tables" 
       << G4endl;
  }
  ComputeProtonDEDX(proton, dedxProton);
  builder->BuildRangeTable(dedxProton, rangeProton);  
  builder->BuildInverseRangeTable(rangeProton, invRangeProton);  
 
  ComputeTrasportXS(electron, mscElectron);
}

operator=()

G4TablesForExtrapolator & G4TablesForExtrapolator::operator= ( const G4TablesForExtrapolator & right )

delete

---

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

• G4TablesForExtrapolator.hh
• G4TablesForExtrapolator.cc