G4PhysicsTableHelper Class Reference

#include <G4PhysicsTableHelper.hh>

Static Public Member Functions
static G4PhysicsTable *	PreparePhysicsTable (G4PhysicsTable *physTable)
static G4bool	RetrievePhysicsTable (G4PhysicsTable *physTable, const G4String &fileName, G4bool ascii)
static void	SetPhysicsVector (G4PhysicsTable *physTable, size_t idx, G4PhysicsVector *vec)
static void	SetVerboseLevel (G4int value)
static G4int	GetVerboseLevel ()

Protected Member Functions
	G4PhysicsTableHelper ()
	~G4PhysicsTableHelper ()
	G4PhysicsTableHelper (const G4PhysicsTableHelper &right)
G4PhysicsTableHelper &	operator= (const G4PhysicsTableHelper &)

Static Protected Attributes
static G4ThreadLocal G4int	verboseLevel = 1

Detailed Description

Definition at line 52 of file G4PhysicsTableHelper.hh.

Constructor & Destructor Documentation

G4PhysicsTableHelper::G4PhysicsTableHelper ( )

protected

Definition at line 46 of file G4PhysicsTableHelper.cc.

{
}

G4PhysicsTableHelper::~G4PhysicsTableHelper ( )

protected

Definition at line 50 of file G4PhysicsTableHelper.cc.

{
}

G4PhysicsTableHelper::G4PhysicsTableHelper ( const G4PhysicsTableHelper &  right )

protected

Definition at line 54 of file G4PhysicsTableHelper.cc.

{
}

Member Function Documentation

G4int G4PhysicsTableHelper::GetVerboseLevel ( )

static

Definition at line 235 of file G4PhysicsTableHelper.cc.

{
  return verboseLevel;
}

G4PhysicsTableHelper & G4PhysicsTableHelper::operator= ( const G4PhysicsTableHelper &  )

protected

Definition at line 58 of file G4PhysicsTableHelper.cc.

{
  return *this;
}

G4PhysicsTable * G4PhysicsTableHelper::PreparePhysicsTable ( G4PhysicsTable *  physTable )

static

Definition at line 64 of file G4PhysicsTableHelper.cc.

{
  G4ProductionCutsTable* cutTable = G4ProductionCutsTable::GetProductionCutsTable();  
  size_t numberOfMCC = cutTable->GetTableSize(); 

  if ( physTable !=0) {
    // compare size of physics table and number of material-cuts-couple
    if ( physTable->size() < numberOfMCC) {
      // enlarge physcis table
      physTable->resize(numberOfMCC, (G4PhysicsVector*)(0));
#ifdef G4VERBOSE  
      if (verboseLevel>2) {
    G4cerr << "G4PhysicsTableHelper::PreparePhysicsTable  ";
    G4cerr << "Physics Table "<< physTable ;
    G4cerr << " is resized to " << numberOfMCC << G4endl;
      }
#endif 
    } else if ( physTable->size() > numberOfMCC){
      // ERROR: this situation should not occur  
      //  size of physics table is shorter than  number of material-cuts-couple
      physTable->resize(numberOfMCC);
#ifdef G4VERBOSE  
      if (verboseLevel>0) {
    G4cerr << "G4PhysicsTableHelper::PreparePhysicsTable  ";
    G4cerr << "Physics Table "<< physTable ;
    G4cerr << " is longer than number of material-cuts-couple " << G4endl;
      }
#endif 
      G4Exception( "G4PhysicsTableHelper::PreparePhysicsTable()",
           "ProcCuts001", FatalException, 
           "Physics Table is inconsistent with  material-cuts-couple");
    } 
  } else {
    // create PhysicsTable is given poitner is null
    physTable = new G4PhysicsTable(numberOfMCC);
    if (physTable!=0) {
      physTable->resize(numberOfMCC, (G4PhysicsVector*)(0));
    } else {
      G4Exception( "G4PhysicsTableHelper::PreparePhysicsTable()",
           "ProcCuts002", FatalException, 
           "Can't create Physics Table");
    }
  }

#ifdef G4VERBOSE  
  if (verboseLevel>2) {
    if ( physTable !=0) { 
      G4cerr << "Physics Table size "<< physTable->size();
    } else {
      G4cerr << "Physics Table does not exist   ";
    }
    G4cerr << ": number of material-cuts-couple " << numberOfMCC << G4endl;
  }
#endif 

  // Reset recal-needed flag for all physics vectors
  physTable->ResetFlagArray();

  for (size_t idx = 0; idx <numberOfMCC; idx +=1){
    const G4MaterialCutsCouple* mcc = cutTable->GetMaterialCutsCouple(idx);
    //check if re-calculation of the physics vector is needed 
    // MCC is not used
    if ( !mcc->IsUsed() ) physTable->ClearFlag(idx);

    // RecalcNeeded flag of MCC is not asserted 
    if ( !mcc->IsRecalcNeeded() ) physTable->ClearFlag(idx);
  }
  
  return physTable;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4PhysicsTableHelper::RetrievePhysicsTable ( G4PhysicsTable *  physTable, const G4String &  fileName, G4bool  ascii  )

static

Definition at line 137 of file G4PhysicsTableHelper.cc.

{
  if (physTable == 0) return false;
  
  // retrieve physics table from the given file
  G4PhysicsTable* tempTable = new G4PhysicsTable();
  if (! tempTable->RetrievePhysicsTable(fileName,ascii) ){
#ifdef G4VERBOSE  
    if (verboseLevel>1) {
      G4cerr << "G4PhysicsTableHelper::RetrievePhysicsTable  ";
      G4cerr << "Fail to retreive from "<< fileName << G4endl;
    }
#endif 
    G4Exception( "G4ProductionCutsTable::RetrievePhysicsTable()",
         "ProcCuts105",
         JustWarning, "Can not retrieve physics tables from file");
    delete tempTable;
    return false;
  } 

  G4ProductionCutsTable* cutTable = G4ProductionCutsTable::GetProductionCutsTable();  
  const G4MCCIndexConversionTable* converter = cutTable->GetMCCIndexConversionTable();

  // check physics table size
  if ( tempTable->size() != converter->size()){
#ifdef G4VERBOSE  
    if (verboseLevel>0) {
      G4cerr << "G4PhysicsTableHelper::RetrievePhysicsTable  ";
      G4cerr << "Size of the physics table in "<< fileName;
      G4cerr << "( size =" << tempTable->size() << ")";
      G4cerr << " is inconsistent with material-cut info";
      G4cerr << "( size =" << converter->size() << ")";
      G4cerr << G4endl;
    }
#endif
    G4Exception( "G4ProductionCutsTable::RetrievePhysicsTable()",
         "ProcCuts106",
         JustWarning, "Retrived file is inconsistent with current physics tables ");
    delete tempTable;
    return false;
  }
  
  // fill the given physics table with retrived physics vectors 
  for (size_t idx=0; idx<converter->size(); idx++){
    if (converter->IsUsed(idx)){
      if (converter->GetIndex(idx)<0) continue;
      size_t i = converter->GetIndex(idx);
      G4PhysicsVector* vec = (*physTable)[i];
       if (vec !=0 ) delete vec;
      (*physTable)[i] =  (*tempTable)[idx];
      physTable->ClearFlag(i);
    }
  }
  tempTable->clear();
  delete tempTable;

  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4PhysicsTableHelper::SetPhysicsVector ( G4PhysicsTable *  physTable, size_t  idx, G4PhysicsVector *  vec  )

static

Definition at line 199 of file G4PhysicsTableHelper.cc.

{
  if ( physTable ==0) {  return;  }

  if ( physTable->size() <= idx) {
#ifdef G4VERBOSE  
    if (verboseLevel>0) {
      G4cerr << "G4PhysicsTableHelper::SetPhysicsVector   ";
      G4cerr << "Given index (" << idx << ")  exceeds ";
      G4cerr << "size of the physics table ";
      G4cerr << "( size =" << physTable->size()<< ")";
      G4cerr << G4endl;
    }
#endif
    G4Exception( "G4ProductionCutsTable::SetPhysicsVector()",
         "ProcCuts107",
         JustWarning, "Illegal index ");
    return;
  } 

  // set physics vector 
  (*physTable)[idx] = vec;
  // clear flag
  physTable->ClearFlag(idx);
 

}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4PhysicsTableHelper::SetVerboseLevel ( G4int  value )

static

Definition at line 230 of file G4PhysicsTableHelper.cc.

{
  verboseLevel = value;
}

Member Data Documentation

G4ThreadLocal G4int G4PhysicsTableHelper::verboseLevel = 1

staticprotected

Definition at line 87 of file G4PhysicsTableHelper.hh.

---

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

• geant4.10.03.p01/source/processes/cuts/include/G4PhysicsTableHelper.hh
• geant4.10.03.p01/source/processes/cuts/src/G4PhysicsTableHelper.cc