G4EmConfigurator Class Reference

#include <G4EmConfigurator.hh>

Collaboration diagram for G4EmConfigurator:

Public Member Functions
	G4EmConfigurator (G4int verboseLevel=1)
	~G4EmConfigurator ()
void	SetExtraEmModel (const G4String &particleName, const G4String &processName, G4VEmModel *, const G4String &regionName="", G4double emin=0.0, G4double emax=DBL_MAX, G4VEmFluctuationModel *fm=0)
void	AddModels ()
void	PrepareModels (const G4ParticleDefinition *aParticle, G4VEnergyLossProcess *p)
void	PrepareModels (const G4ParticleDefinition *aParticle, G4VEmProcess *p)
void	PrepareModels (const G4ParticleDefinition *aParticle, G4VMultipleScattering *p)
void	Clear ()
void	SetVerbose (G4int value)

Private Member Functions
G4Region *	FindRegion (const G4String &)
void	SetModelForRegion (G4VEmModel *model, G4VEmFluctuationModel *fm, G4Region *reg, const G4String &particleName, const G4String &processName, G4double emin, G4double emax)
G4bool	UpdateModelEnergyRange (G4VEmModel *mod, G4double emin, G4double emax)
G4EmConfigurator &	operator= (const G4EmConfigurator &right)
	G4EmConfigurator (const G4EmConfigurator &)

Private Attributes
std::vector< G4VEmModel * >	models
std::vector< G4VEmFluctuationModel * >	flucModels
std::vector< G4String >	particles
std::vector< G4String >	processes
std::vector< G4String >	regions
std::vector< G4double >	lowEnergy
std::vector< G4double >	highEnergy
G4int	index
G4int	verbose

Detailed Description

Definition at line 63 of file G4EmConfigurator.hh.

Constructor & Destructor Documentation

G4EmConfigurator() [1/2]

G4EmConfigurator::G4EmConfigurator

(

G4int

verboseLevel = 1

)

Definition at line 65 of file G4EmConfigurator.cc.

                                           :verbose(val)
{
  index = -10;
}

~G4EmConfigurator()

G4EmConfigurator::~G4EmConfigurator

(

)

Definition at line 72 of file G4EmConfigurator.cc.

{}

G4EmConfigurator() [2/2]

G4EmConfigurator::G4EmConfigurator ( const G4EmConfigurator &  )

private

Member Function Documentation

AddModels()

void G4EmConfigurator::AddModels

(

)

Definition at line 113 of file G4EmConfigurator.cc.

{
  size_t n = models.size();
  if(0 < verbose) {
    G4cout << "### G4EmConfigurator::AddModels n= " << n << G4endl;
  }
  if(n > 0) {
    for(size_t i=0; i<n; ++i) {
      if(models[i]) {
        G4Region* reg = FindRegion(regions[i]);
        if(reg) {
          --index;
          SetModelForRegion(models[i],flucModels[i],reg,
                            particles[i],processes[i],
                            lowEnergy[i],highEnergy[i]);
        }
      }
    }
  }
  Clear();
}

Here is the call graph for this function:

Here is the caller graph for this function:

Clear()

void G4EmConfigurator::Clear

(

)

Definition at line 355 of file G4EmConfigurator.cc.

{
  particles.clear();
  processes.clear();
  models.clear();
  flucModels.clear();
  regions.clear();
  lowEnergy.clear();
  highEnergy.clear();
}

Here is the caller graph for this function:

FindRegion()

G4Region * G4EmConfigurator::FindRegion ( const G4String &  regionName )

private

Definition at line 368 of file G4EmConfigurator.cc.

{
  // search for region
  G4Region* reg = 0;
  G4RegionStore* regStore = G4RegionStore::GetInstance();
  G4String r = regionName;
  if(r == "" || r == "world" || r == "World") { 
    r = "DefaultRegionForTheWorld";
  }
  reg = regStore->GetRegion(r, true); 
  if(!reg) {
    G4cout << "### G4EmConfigurator WARNING: fails to find a region <"
           << r << G4endl;
  } else if(verbose > 1) {
    G4cout << "### G4EmConfigurator finds out G4Region <" << r << ">" 
           << G4endl;
  } 
  return reg;  
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

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

private

PrepareModels() [1/3]

void G4EmConfigurator::PrepareModels	(	const G4ParticleDefinition *	aParticle,
		G4VEnergyLossProcess *	p
	)

Definition at line 232 of file G4EmConfigurator.cc.

{
  size_t n = particles.size();
  if(1 < verbose) {
    G4cout << " G4EmConfigurator::PrepareModels for EnergyLoss n= " 
           << n << G4endl;
  }
  if(n > 0) {
    G4String particleName = aParticle->GetParticleName(); 
    G4String processName  = p->GetProcessName(); 
    //G4cout <<  particleName << "  " <<  processName << G4endl;
    for(size_t i=0; i<n; ++i) {
      //G4cout <<  particles[i] << "  " <<  processes[i] << G4endl;
      if(processName == processes[i]) {
        if((particleName == particles[i]) ||
           (particles[i] == "all") ||
           (particles[i] == "charged" && aParticle->GetPDGCharge() != 0.0)) {
          G4Region* reg = FindRegion(regions[i]);
          //G4cout << "Region " << reg << G4endl;
          if(reg) {
            --index;
            G4VEmModel* mod = models[i];
            G4VEmFluctuationModel* fm = flucModels[i];
            if(mod) {
              if(UpdateModelEnergyRange(mod, lowEnergy[i], highEnergy[i])) {
                p->AddEmModel(index,mod,fm,reg);
                if(1 < verbose) {
                  G4cout << "### Added eloss model order= " << index << " for " 
                         << particleName << " and " << processName << G4endl;
                }
              }
            } else if(fm) {
              p->SetFluctModel(fm);
            }
          }
        }
      }
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

PrepareModels() [2/3]

void G4EmConfigurator::PrepareModels	(	const G4ParticleDefinition *	aParticle,
		G4VEmProcess *	p
	)

Definition at line 277 of file G4EmConfigurator.cc.

{
  size_t n = particles.size();
  if(1 < verbose) {
    G4cout << " G4EmConfigurator::PrepareModels for EM process n= " 
           << n << G4endl;
  }
  if(n > 0) {
    G4String particleName = aParticle->GetParticleName(); 
    G4String processName  = p->GetProcessName(); 
    //G4cout <<  particleName << "  " <<  particleName << G4endl;
    for(size_t i=0; i<n; ++i) {
      if(processName == processes[i]) {
        if((particleName == particles[i]) ||
           (particles[i] == "all") ||
           (particles[i] == "charged" && aParticle->GetPDGCharge() != 0.0)) {
          G4Region* reg = FindRegion(regions[i]);
          //G4cout << "Region " << reg << G4endl;
          if(reg) {
            --index;
            G4VEmModel* mod = models[i];
            if(mod) {
              if(UpdateModelEnergyRange(mod, lowEnergy[i], highEnergy[i])) { 
                p->AddEmModel(index,mod,reg);
                if(1 < verbose) {
                  G4cout << "### Added em model order= " << index << " for " 
                         << particleName << " and " << processName << G4endl;
                }
              }
            }
          }
        }
      }
    }
  }
}

Here is the call graph for this function:

PrepareModels() [3/3]

void G4EmConfigurator::PrepareModels	(	const G4ParticleDefinition *	aParticle,
		G4VMultipleScattering *	p
	)

Definition at line 318 of file G4EmConfigurator.cc.

{
  size_t n = particles.size();
  if(1 < verbose) {
    G4cout << " G4EmConfigurator::PrepareModels for MSC process n= " 
           << n << G4endl;
  }

  if(n > 0) {
    G4String particleName = aParticle->GetParticleName(); 
    G4String processName  = p->GetProcessName(); 
    for(size_t i=0; i<n; ++i) {
      if(processName == processes[i]) {
        if((particleName == particles[i]) ||
           (particles[i] == "all") ||
           (particles[i] == "charged" && aParticle->GetPDGCharge() != 0.0)) {
          G4Region* reg = FindRegion(regions[i]);
          if(reg) {
            --index;
            G4VEmModel* mod = models[i];
            if(mod) {
              if(UpdateModelEnergyRange(mod, lowEnergy[i], highEnergy[i])) { 
                p->AddEmModel(index,mod,reg);
                //G4cout << "### Added msc model order= " << index << " for " 
                //     << particleName << " and " << processName << G4endl;
              }
            }
          }
        }
      }
    }
  }
}

Here is the call graph for this function:

SetExtraEmModel()

void G4EmConfigurator::SetExtraEmModel	(	const G4String &	particleName,
		const G4String &	processName,
		G4VEmModel *	mod,
		const G4String &	regionName = "",
		G4double	emin = 0.0,
		G4double	emax = DBL_MAX,
		G4VEmFluctuationModel *	fm = 0
	)

Definition at line 77 of file G4EmConfigurator.cc.

{
  if(1 < verbose) {
    G4cout << " G4EmConfigurator::SetExtraEmModel " << mod->GetName()
           << " for " << particleName 
           << " and " << processName 
           << " in the region <" << regionName
           << "> Emin(MeV)= " << emin/MeV
           << " Emax(MeV)= " << emax/MeV
           << G4endl;
  }
  if(mod) {
    models.push_back(mod);
    flucModels.push_back(fm);
    emax = std::min(emax, mod->HighEnergyLimit());
    mod->SetActivationHighEnergyLimit(emax);
  } else {
    models.push_back(new G4DummyModel());
    flucModels.push_back(0);
  }

  particles.push_back(particleName);
  processes.push_back(processName);
  regions.push_back(regionName);
  lowEnergy.push_back(emin);
  highEnergy.push_back(emax);
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetModelForRegion()

void G4EmConfigurator::SetModelForRegion ( G4VEmModel *  model, G4VEmFluctuationModel *  fm, G4Region *  reg, const G4String &  particleName, const G4String &  processName, G4double  emin, G4double  emax  )

private

Definition at line 137 of file G4EmConfigurator.cc.

{
  if(1 < verbose) {
    G4cout << " G4EmConfigurator::SetModelForRegion: " << mod->GetName() 
           << G4endl;
    G4cout << " For " << particleName 
           << " and " << processName 
           << " in the region <" << reg->GetName()
           << " Emin(MeV)= " << emin/MeV
           << " Emax(MeV)= " << emax/MeV;
    if(fm) { G4cout << " FLmodel " << fm->GetName(); }
    G4cout << G4endl;
  }
  G4ParticleTable::G4PTblDicIterator* theParticleIterator = 
    G4ParticleTable::GetParticleTable()->GetIterator(); 

  // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  theParticleIterator->reset();
  while( (*theParticleIterator)() ) {
    const G4ParticleDefinition* part = theParticleIterator->value();

    //G4cout << particleName << " " << part->GetParticleName() << G4endl;

    if((part->GetParticleName() == particleName) ||
       (particleName == "all") ||
       (particleName == "charged" && part->GetPDGCharge() != 0.0)) {

      // search for process
      G4ProcessManager* pmanager = part->GetProcessManager();
      G4ProcessVector* plist = pmanager->GetProcessList();
      G4int np = pmanager->GetProcessListLength();
  
      //G4cout << processName << " in list of " << np << G4endl;

      G4VProcess* proc = 0;
      for(G4int i=0; i<np; ++i) {
        if(processName == (*plist)[i]->GetProcessName()) {
          proc = (*plist)[i];
          break;
        }
      }
      if(!proc) {
        G4cout << "### G4EmConfigurator WARNING: fails to find a process <"
               << processName << "> for " << particleName << G4endl;
        return;        
      } 

      if(mod) {
        if(!UpdateModelEnergyRange(mod, emin, emax)) { return; }
      }
      // classify process
      G4int ii = proc->GetProcessSubType();
      if(10 == ii && mod) {
        G4VMultipleScattering* p = static_cast<G4VMultipleScattering*>(proc);
        p->AddEmModel(index,mod,reg);
        if(1 < verbose) {
          G4cout << "### Added msc model order= " << index << " for " 
                 << particleName << " and " << processName << G4endl;
        }
        return;
      } else if(2 <= ii && 4 >= ii) {
        G4VEnergyLossProcess* p = static_cast<G4VEnergyLossProcess*>(proc);
        if(!mod && fm) {
          p->SetFluctModel(fm);
        } else { 
          p->AddEmModel(index,mod,fm,reg);
          if(1 < verbose) {
            G4cout << "### Added eloss model order= " << index << " for " 
                   << particleName << " and " << processName << G4endl;
          }
        }
        return;
      } else if(mod) {
        G4VEmProcess* p = static_cast<G4VEmProcess*>(proc);
        p->AddEmModel(index,mod,reg);
        if(1 < verbose) {
          G4cout << "### Added em model order= " << index << " for " 
                 << particleName << " and " << processName << G4endl;
        }
        return;
      } else {
        return;
      }
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetVerbose()

void G4EmConfigurator::SetVerbose ( G4int  value )

inline

Definition at line 136 of file G4EmConfigurator.hh.

{
  verbose = value;
}

Here is the caller graph for this function:

UpdateModelEnergyRange()

G4bool G4EmConfigurator::UpdateModelEnergyRange ( G4VEmModel *  mod, G4double  emin, G4double  emax  )

private

Definition at line 390 of file G4EmConfigurator.cc.

{
  // energy limits
  G4double e1 = std::max(emin,mod->LowEnergyLimit());
  G4double e2 = std::min(emax,mod->HighEnergyLimit());
  if(e2 <= e1) {
    G4cout << "### G4EmConfigurator WARNING: empty energy interval"
           << " for <" << mod->GetName() 
           << ">  Emin(MeV)= " << e1/CLHEP::MeV 
           << ">  Emax(MeV)= " << e2/CLHEP::MeV 
           << G4endl;
    return false;        
  }
  mod->SetLowEnergyLimit(e1);
  mod->SetHighEnergyLimit(e2);
  if(verbose > 1) {
    G4cout << "### G4EmConfigurator for " << mod->GetName() 
           << " Emin(MeV)= " << e1/MeV << " Emax(MeV)= " << e2/MeV 
           << G4endl;
  } 
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

flucModels

std::vector<G4VEmFluctuationModel*> G4EmConfigurator::flucModels

private

Definition at line 123 of file G4EmConfigurator.hh.

highEnergy

std::vector<G4double> G4EmConfigurator::highEnergy

private

Definition at line 128 of file G4EmConfigurator.hh.

index

G4int G4EmConfigurator::index

private

Definition at line 130 of file G4EmConfigurator.hh.

lowEnergy

std::vector<G4double> G4EmConfigurator::lowEnergy

private

Definition at line 127 of file G4EmConfigurator.hh.

models

std::vector<G4VEmModel*> G4EmConfigurator::models

private

Definition at line 122 of file G4EmConfigurator.hh.

particles

std::vector<G4String> G4EmConfigurator::particles

private

Definition at line 124 of file G4EmConfigurator.hh.

processes

std::vector<G4String> G4EmConfigurator::processes

private

Definition at line 125 of file G4EmConfigurator.hh.

regions

std::vector<G4String> G4EmConfigurator::regions

private

Definition at line 126 of file G4EmConfigurator.hh.

verbose

G4int G4EmConfigurator::verbose

private

Definition at line 131 of file G4EmConfigurator.hh.

---

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

• G4EmConfigurator.hh
• G4EmConfigurator.cc