G4PhysicsListHelper.cc

Go to the documentation of this file.

//
// ********************************************************************
// * 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.          *
// ********************************************************************
//
//
// $Id: G4PhysicsListHelper.cc 97917 2016-06-23 07:12:01Z gcosmo $
//
// 
// ------------------------------------------------------------
//  GEANT 4 class header file 
//
// ------------------------------------------------------------
//  History
//       first version                   29 Apr 2011 by H.Kurashige
// ------------------------------------------------------------

#include "globals.hh"
#include "G4PhysicsListHelper.hh"
#include "G4ParticleDefinition.hh"
#include "G4ProcessManager.hh"
#include "G4ParticleTable.hh"

#include "G4ios.hh"
#include <iomanip>
#include <fstream>

G4ThreadLocal G4PhysicsListHelper* G4PhysicsListHelper::pPLHelper = 0;
 
G4PhysicsListHelper::G4PhysicsListHelper()
  :  useCoupledTransportation(false),
     theTransportationProcess(0),
     verboseLevel(1),
     theTable(0),
     sizeOfTable(0),
     ordParamFileName("")
{
  // pointer to the particle table
  theParticleTable = G4ParticleTable::GetParticleTable();
  aParticleIterator = theParticleTable->GetIterator();

  ReadOrdingParameterTable();

#ifdef G4VERBOSE
  if (verboseLevel >1){
    DumpOrdingParameterTable();
  }
#endif
}

G4PhysicsListHelper::~G4PhysicsListHelper()
{
  if (theTable !=0) {
    theTable->clear();
    delete theTable;
    theTable=0;
    sizeOfTable=0;
  }
  if (theTransportationProcess!=0) {
    delete theTransportationProcess;
    theTransportationProcess=0;
  }
}

G4PhysicsListHelper* G4PhysicsListHelper::GetPhysicsListHelper()  
{
  if (!pPLHelper)
  {
    pPLHelper = new G4PhysicsListHelper;
  }
  return pPLHelper;
}

void G4PhysicsListHelper::CheckParticleList() const
{
  bool isElectron = false;
  bool isPositron = false;
  bool isGamma    = false;
  bool isProton   = false;
  bool isGenericIon = false;
  bool isAnyIon   = false;
  bool isAnyChargedBaryon   = false;
  bool isEmProc   = false;

  // loop over all particles in G4ParticleTable
  aParticleIterator->reset();
  while( (*aParticleIterator)() ){
    G4ParticleDefinition* particle = aParticleIterator->value();
    G4String name = particle->GetParticleName();
    // check if any EM process exists
    if (!isEmProc) {
      G4ProcessVector* list = particle->GetProcessManager()->GetProcessList();
      for (int idx=0; idx<list->size(); idx++){
    isEmProc = ((*list)[idx])->GetProcessType() == fElectromagnetic;
    if (isEmProc) break;
      }
    }
    
    if      ( name == "e-") isElectron = true; 
    else if ( name == "e+") isPositron = true; 
    else if ( name == "gamma") isGamma = true; 
    else if ( name == "GenericIon") isGenericIon = true; 
    else if ( name == "proton") isProton = true; 
    else if ( particle->GetParticleType() == "nucleus") isAnyIon = true;
    else if ( particle->GetParticleType() == "baryon") {
       if ( particle->GetPDGCharge() != 0.0 ) isAnyChargedBaryon = true;
    }
  }

  if (!isEmProc) return;

  // RULE 1
  //  e+, e- and gamma should exist 
  //   if one of them exist
  bool isEmBasic =  isElectron || isPositron || isGamma;
  bool isMissingEmBasic =  !isElectron || !isPositron || !isGamma;
  if (isEmBasic && isMissingEmBasic) {
    G4String missingName="";
    if (!isElectron) missingName += "e- ";
    if (!isPositron) missingName += "e+ ";
    if (!isGamma) missingName += "gamma ";

#ifdef G4VERBOSE
    if (verboseLevel >0){
      G4cout << "G4PhysicsListHelper::CheckParticleList: "
         << missingName << " do not exist " << G4endl; 
      G4cout << " These particle are necessary for basic EM processes" 
         << G4endl;
    }
#endif
    G4Exception("G4PhysicsListHelper::CheckParticleList",
        "Run0101", FatalException,
        "Missing EM basic particle");
  }

  // RULE 2
  //  proton should exist 
  //   if any other charged baryon  exist
  if (!isProton && isAnyChargedBaryon) {
    G4String missingName="proton ";

#ifdef G4VERBOSE
    if (verboseLevel >0){
      G4cout << "G4PhysicsListHelper::CheckParticleList: "
         << missingName << " does not exist "<< G4endl; 
      G4cout << " Proton is necessary for EM baryon processes" << G4endl;
    }
#endif
    missingName += " should be created ";
    G4Exception("G4PhysicsListHelper::CheckParticleList",
        "Run0102", FatalException,
        "Missing Proton");
  }
   
  // RULE 3
  //  GenericIonn should exist 
  //   if any other ion  exist
  if (!isGenericIon && isAnyIon) {
    G4String missingName="GenericIon ";

#ifdef G4VERBOSE
    if (verboseLevel >0){
      G4cout << "G4PhysicsListHelper::CheckParticleList: "
         << missingName << " does not exist "<< G4endl; 
      G4cout << " GenericIon should be created if any ion is necessary" << G4endl;
    }
#endif
    G4Exception("G4PhysicsListHelper::CheckParticleList",
        "Run0103", FatalException,
        "Missing GenericIon");
  }
      
}


#include "G4Transportation.hh"
#include "G4CoupledTransportation.hh"
#include "G4RunManagerKernel.hh"
#include "G4ScoringManager.hh"

void G4PhysicsListHelper::AddTransportation()
{
  G4int verboseLevelTransport = 0;

#ifdef G4VERBOSE
  if (verboseLevel >2){
    G4cout << "G4PhysicsListHelper::AddTransportation()  "<< G4endl;
  }
#endif

  G4int nParaWorld = 
    G4RunManagerKernel::GetRunManagerKernel()->GetNumberOfParallelWorld();
  
  if ( nParaWorld>0 || 
       useCoupledTransportation || 
       G4ScoringManager::GetScoringManagerIfExist()) {
#ifdef G4VERBOSE
    if (verboseLevel >0) {
      G4cout << " G4PhysicsListHelper::AddTransportation()"
         << "--- G4CoupledTransportation is used " 
         << G4endl;
    }
#endif
    theTransportationProcess = new G4CoupledTransportation(verboseLevelTransport);    
  } else {
    theTransportationProcess = new G4Transportation(verboseLevelTransport);
  }
 
  // loop over all particles in G4ParticleTable
  aParticleIterator->reset();
  while( (*aParticleIterator)() ){
    G4ParticleDefinition* particle = aParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    // Add transportation process for all particles 
    if ( pmanager == 0) {
      // Error !! no process manager
#ifdef G4VERBOSE
      if (verboseLevel>0){
    G4cout << "G4PhysicsListHelper::AddTransportation  "
           <<" : No Process Manager for "
           << particle->GetParticleName() << G4endl;
      }
#endif
      G4Exception("G4PhysicsListHelper::AddTransportation",
          "Run0104", FatalException,
          "No process manager");
      continue;
    } 
    // Molecule use different type transportation
    if(particle->GetParticleType() == "Molecule") continue; 

    // add transportation with ordering = ( -1, "first", "first" )
    pmanager ->AddProcess(theTransportationProcess);
    pmanager ->SetProcessOrderingToFirst(theTransportationProcess, idxAlongStep);
    pmanager ->SetProcessOrderingToFirst(theTransportationProcess, idxPostStep);
  }
}

#include "G4ProcessManager.hh"
  
void G4PhysicsListHelper::ReadOrdingParameterTable()
{
  G4bool readInFile = false;
  std::ifstream fIn;  

  if( getenv("G4ORDPARAMTABLE") ){
    ordParamFileName = getenv("G4ORDPARAMTABLE");
#ifdef G4VERBOSE
    if (verboseLevel >1){
      G4cout << "G4PhysicsListHelper::ReadOrdingParameterTable  :"
         << ordParamFileName << " is assigned to Ordering Parameter Table "
         << G4endl; 
    }
#endif
    // open input file //
    fIn.open(ordParamFileName, std::ios::in);
    // check if the file has been opened successfully 
    if (!fIn) {
#ifdef G4VERBOSE
      if (verboseLevel >0) {
    G4cout << "G4PhysicsListHelper::ReadOrdingParameterTable  "
           << " Can not open file " << ordParamFileName << G4endl;
      }
#endif
      G4Exception("G4PhysicsListHelper::ReadOrdingParameterTable",
          "Run0105", JustWarning, 
          "Fail to open ordering paramter table ");
    } else {
      readInFile = true;
    }
  }
 

  // create OrdParamTable   
  if (theTable !=0) {
    theTable->clear();
    delete theTable;
    theTable=0;
    sizeOfTable=0;
  }
  theTable = new G4OrdParamTable();
  sizeOfTable=0;

  if (readInFile){
    // read in the file and fill the table 
    while(!fIn.eof()) {
      G4PhysicsListOrderingParameter tmp;
      G4int flag;
      fIn >> tmp.processTypeName >>  tmp.processType >> tmp.processSubType
      >> tmp.ordering[0] >> tmp.ordering[1] >> tmp.ordering[2] >> flag;
      tmp.isDuplicable = (flag!=0);
      theTable->push_back(tmp);
      sizeOfTable +=1;  
    }
    fIn.close();
  } else {
    ReadInDefaultOrderingParameter();
  }

  if (sizeOfTable==0){
#ifdef G4VERBOSE
    if (verboseLevel >0) {
      G4cout << "G4PhysicsListHelper::ReadOrdingParameterTable "
         << " Empty file " << ordParamFileName << G4endl;
    }
#endif
    G4Exception("G4PhysicsListHelper::ReadOrdingParameterTable",
        "Run0106", JustWarning, 
        "The ordering parameter table is empty ");
    delete theTable;
    theTable=0;
    sizeOfTable=0;
  }
  return;  
}

void G4PhysicsListHelper::DumpOrdingParameterTable(G4int subType) const
{
 if (theTable==0) {
#ifdef G4VERBOSE
    if (verboseLevel >0) {
      G4cout << "G4PhysicsListHelper::DumpOrdingParameterTable   "
         << " No ordering parameter table  : " << ordParamFileName 
         << G4endl;
    }
#endif
    return;
  }
  G4cout << "G4PhysicsListHelper::DumpOrdingParameterTable  : "
     << ordParamFileName << G4endl; 
  G4cout << "          TypeName  " 
     << "    ProcessType" <<  "        SubType"
     << "         AtRest" <<  "      AlongStep" <<  "        PostStep"
     << "     Duplicable" << G4endl;
  for (int i=0; i<sizeOfTable ; i++){
    G4PhysicsListOrderingParameter* tmp=&(theTable->at(i));
    if ((subType>=0) && (subType!=tmp->processSubType)) continue;
    G4cout << std::setw(18)     << tmp->processTypeName 
       << std::setw(15)     << tmp->processType 
       << std::setw(15) << tmp->processSubType
       << std::setw(15) << tmp->ordering[0] 
       << std::setw(15) << tmp->ordering[1] 
       << std::setw(15) << tmp->ordering[2];
    if (tmp->isDuplicable) {
      G4cout << "  true";
    } else {
      G4cout << "  false";
    }
    G4cout <<G4endl;
  }  
}

G4PhysicsListOrderingParameter G4PhysicsListHelper::GetOrdingParameter(G4int subType) const
{
  G4PhysicsListOrderingParameter value;

 if (theTable==0) {
#ifdef G4VERBOSE
    if (verboseLevel >0) {
      G4cout << "G4PhysicsListHelper::GetOrderingParameter : "
         << " No ordering parameter table  : " << ordParamFileName 
         << G4endl;
    }
#endif
    return value;
  }

  for (int i=0; i<sizeOfTable ; i++){
    G4PhysicsListOrderingParameter* tmp=&(theTable->at(i));
    if (subType == tmp->processSubType){
      value.processTypeName = tmp->processTypeName; 
      value.processType     = tmp->processType; 
      value.processSubType  = tmp->processSubType; 
      value.ordering[0]     = tmp->ordering[0]; 
      value.ordering[1]     = tmp->ordering[1]; 
      value.ordering[2]     = tmp->ordering[2]; 
      value.isDuplicable    = tmp->isDuplicable;
    }
  }  
  return value;
}

G4bool G4PhysicsListHelper::RegisterProcess(G4VProcess*            process,
                        G4ParticleDefinition*  particle)
{
  if (theTable==0) {
#ifdef G4VERBOSE
    if (verboseLevel >0) {
      G4cout << "G4PhysicsListHelper::RegisterProcess :"
         << " No ordering parameter table  : " << ordParamFileName 
         << G4endl;
    }
#endif
    G4Exception("G4PhysicsListHelper::RegisterPorcess",
        "Run0107", FatalException, 
        "No Ordering Parameter Table");
    return false;
  }

  const G4String pName = process->GetProcessName(); 
  const G4int pType    = process->GetProcessType();
  const G4int pSubType = process->GetProcessSubType();
  
#ifdef G4VERBOSE
  if (verboseLevel >2) {
    G4cout << "G4PhysicsListHelper::RegisterProcess :"
       << pName << " Process Type = " << pType 
       << " SubType = "<< pSubType  
       << " to " << particle->GetParticleName()
       << G4endl;
  }
#endif

  // Check Process Type/SubType
  if ((pType <1)||(pSubType<1)) {
#ifdef G4VERBOSE
    if (verboseLevel >0) {
      G4cout << "G4PhysicsListHelper::RegisterProcess :"
         << pName << " for " << particle->GetParticleName()
         << " has illegal Process Type = " << pType 
         << " SubType = "<< pSubType << G4endl;
    }
#endif
    G4Exception("G4PhysicsListHelper::RegisterPorcess",
        "Run0108", FatalException, 
        "No Matching process Type/SubType");
    return false;
  }
  
  G4bool isFound = false;
  G4int  ord[3];
  G4bool duplicable = false;
  for (int i=0; i<sizeOfTable ; i++){
    G4PhysicsListOrderingParameter* tmp=&(theTable->at(i));
    if ((tmp->processType==pType)&&(tmp->processSubType==pSubType)){
      ord[0] = tmp->ordering[0]; 
      ord[1] = tmp->ordering[1]; 
      ord[2] = tmp->ordering[2];
      duplicable = tmp->isDuplicable;
      isFound = true;
      break;
    }
  } 
  if (!isFound) {
#ifdef G4VERBOSE
    if (verboseLevel >0) {
      G4cout << "G4PhysicsListHelper::RegisterProcess :"
         << pName << " for " << particle->GetParticleName()
         << " with  type/subtype =" 
         << pType << "/" << pSubType 
         << "  is not reigstered in OrdingParameterTable  "
         << G4endl;
    }
#endif
    G4Exception("G4PhysicsListHelper::RegisterPorcess",
        "Run0109", FatalException, 
        "No Matching process Type/SubType");
    return false;
  }

  // Check Process Manager
  G4ProcessManager* pManager = particle->GetProcessManager();
  if ( pManager == 0) {
      // Error !! no process manager
#ifdef G4VERBOSE
    if (verboseLevel>0){
      G4cout << "G4PhysicsListHelper::RegisterProcess "
         <<" : No Process Manager for "
         << particle->GetParticleName() << G4endl;
    }
#endif
    G4Exception("G4PhysicsListHelper::RegisterProcess   ",
        "Riun0110", FatalException,
        "No process manager");
    return false;
  }

  // Check Duplication
  if (!duplicable){
    G4bool duplicated = false;
    G4ProcessVector* pList = pManager->GetProcessList();
    for (G4int idx=0; idx<pList->size(); idx++) {
      const G4VProcess* p = (*pList)[idx];
      if ((p->GetProcessType()== pType)  && 
      (p->GetProcessSubType()== pSubType)){
    duplicated = true;
#ifdef G4VERBOSE
    if (verboseLevel >0) {
      G4cout << "G4PhysicsListHelper::RegisterProcess :"
         << pName << " for " << particle->GetParticleName()
         << " with  type/subtype =" 
         << pType << "/" << pSubType 
         << "  is has same subType as "
         << p->GetProcessName()
         << " for " << particle->GetParticleName()
         << G4endl;
      G4cout << "It will not be added !!" << G4endl;
    }
#endif
    G4Exception("G4PhysicsListHelper::RegisterPorcess",
            "Run0111", JustWarning, 
            "Duplication of processes");
      }
    }
    if (duplicated) return false;
  }

  // Add Process
  G4int code = pManager ->AddProcess(process);
  if (code <0) return false;

  // Set Ordering Parameter
  for(G4int idx=0; idx<3; idx++){
    G4ProcessVectorDoItIndex idxOrd = static_cast<G4ProcessVectorDoItIndex>(idx);
    if (ord[idx]<0) {
      // Do Nothing because NO DOIT
    } else if (ord[idx]==0) {
      pManager->SetProcessOrderingToFirst( process, idxOrd );
    } else if (ord[idx]<9999) {
      pManager->SetProcessOrdering( process, idxOrd , ord[idx]);
    } else {
      pManager->SetProcessOrderingToLast( process, idxOrd );
    } 
  } 
#ifdef G4VERBOSE
  if (verboseLevel >1) {
    G4cout << "G4PhysicsListHelper::RegisterProcess :"
       << pName << " for " << particle->GetParticleName()
       << " with  type/subtype =" 
       << pType << "/" << pSubType 
       << " is sucessfully registered with ordering parameters "
       << ord[0] << ":" << ord[1] << ":" << ord[2]
       << G4endl;
    }
#endif
  return true;
}

void G4PhysicsListHelper::ReadInDefaultOrderingParameter()
{
  
  G4PhysicsListOrderingParameter tmp;
  
  tmp.processTypeName = "Transportation";
  tmp.processType     = 1;
  tmp.processSubType  = 91;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     =  0;
  tmp.ordering[2]     =  0;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "CoupleTrans";
  tmp.processType     = 1;
  tmp.processSubType  = 92;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     =  0;
  tmp.ordering[2]     =  0;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "CoulombScat";
  tmp.processType     = 2;
  tmp.processSubType  =  1;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     = 1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Ionisation";
  tmp.processType     = 2;
  tmp.processSubType  =  2;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     =  2;
  tmp.ordering[2]     =  2;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Brems";
  tmp.processType     = 2;
  tmp.processSubType  =  3;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  3;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "PairProdCharged";
  tmp.processType     = 2;
  tmp.processSubType  =  4;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  4;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Annih";
  tmp.processType     = 2;
  tmp.processSubType  =  5;
  tmp.ordering[0]     =  5;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  5;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "AnnihToMuMu";
  tmp.processType     = 2;
  tmp.processSubType  =  6;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  6;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "AnnihToHad";
  tmp.processType     = 2;
  tmp.processSubType  =  7;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  7;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "NuclearStopp";
  tmp.processType     = 2;
  tmp.processSubType  =  8;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     =  8;
  tmp.ordering[2]     = -1;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Msc";
  tmp.processType     = 2;
  tmp.processSubType  = 10;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     =  1;
  tmp.ordering[2]     =  1;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Rayleigh";
  tmp.processType     = 2;
  tmp.processSubType  = 11;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "PhotoElectric";
  tmp.processType     = 2;
  tmp.processSubType  = 12;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Compton";
  tmp.processType     = 2;
  tmp.processSubType  = 13;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Conv";
  tmp.processType     = 2;
  tmp.processSubType  = 14;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "ConvToMuMu";
  tmp.processType     = 2;
  tmp.processSubType  = 15;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Cerenkov";
  tmp.processType     = 2;
  tmp.processSubType  = 21;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "Scintillation";
  tmp.processType     = 2;
  tmp.processSubType  = 22;
  tmp.ordering[0]     =  9999;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  9999;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "SynchRad";
  tmp.processType     = 2;
  tmp.processSubType  = 23;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "TransRad";
  tmp.processType     = 2;
  tmp.processSubType  = 24;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "OpAbsorb";
  tmp.processType     = 3;
  tmp.processSubType  = 31;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "OpBoundary";
  tmp.processType     = 3;
  tmp.processSubType  = 32;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "OpRayleigh";
  tmp.processType     = 3;
  tmp.processSubType  = 33;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "OpWLS";
  tmp.processType     = 3;
  tmp.processSubType  = 34;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "OpMieHG";
  tmp.processType     = 3;
  tmp.processSubType  = 35;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "DNAElastic";
  tmp.processType     = 2;
  tmp.processSubType  = 51;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "DNAExcit";
  tmp.processType     = 2;
  tmp.processSubType  = 52;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "DNAIonisation";
  tmp.processType     = 2;
  tmp.processSubType  = 53;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "DNAVibExcit";
  tmp.processType     = 2;
  tmp.processSubType  = 54;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "DNAAttachment";
  tmp.processType     = 2;
  tmp.processSubType  = 55;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "DNAChargeDec";
  tmp.processType     = 2;
  tmp.processSubType  = 56;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "DNAChargeInc";
  tmp.processType     = 2;
  tmp.processSubType  = 57;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "DNAElectronSolvatation";
  tmp.processType     = 2;
  tmp.processSubType  = 58;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;

  tmp.processTypeName = "DNAMolecularDecay";
  tmp.processType     = 6;
  tmp.processSubType  = 59;
  tmp.ordering[0]     =  1000;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     = -1;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;

  tmp.processTypeName = "ITTransportation";
  tmp.processType     = 1;
  tmp.processSubType  = 60;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     =  0;
  tmp.ordering[2]     =  0;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;

  tmp.processTypeName = "DNABrownianTransportation";
  tmp.processType     = 1;
  tmp.processSubType  = 61;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = 0; 
  tmp.ordering[2]     = 0;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;

  tmp.processTypeName = "DNADoubleIonisation";
  tmp.processType     = 2;
  tmp.processSubType  = 62;
  tmp.ordering[0]     =  -1;
  tmp.ordering[1]     =  -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;

  tmp.processTypeName = "DNADoubleCapture";
  tmp.processType     = 2;
  tmp.processSubType  = 63;
  tmp.ordering[0]     =  -1;
  tmp.ordering[1]     =  -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;

  tmp.processTypeName = "DNAIonisingTransfer";
  tmp.processType     = 2;
  tmp.processSubType  = 64;
  tmp.ordering[0]     =  -1;
  tmp.ordering[1]     =  -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;

  tmp.processTypeName = "HadElastic";
  tmp.processType     = 4;
  tmp.processSubType  = 111;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "HadInElastic";
  tmp.processType     = 4;
  tmp.processSubType  = 121;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "HadCapture";
  tmp.processType     = 4;
  tmp.processSubType  = 131;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "HadFission";
  tmp.processType     = 4;
  tmp.processSubType  = 141;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName = "HadAtRest";
  tmp.processType     = 4;
  tmp.processSubType  = 151;
  tmp.ordering[0]     =  1000;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     = -1;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "HadCEX";
  tmp.processType     = 4;
  tmp.processSubType  = 161;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "Decay";
  tmp.processType     = 6;
  tmp.processSubType  = 201;
  tmp.ordering[0]     =  1000;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "DecayWSpin";
  tmp.processType     = 6;
  tmp.processSubType  = 202;
  tmp.ordering[0]     =  1000;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "DecayPiSpin";
  tmp.processType     = 6;
  tmp.processSubType  = 203;
  tmp.ordering[0]     =  1000;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "DecayRadio";
  tmp.processType     = 6;
  tmp.processSubType  = 210;
  tmp.ordering[0]     =  1000;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "DecayUnKnown";
  tmp.processType     = 6;
  tmp.processSubType  = 211;
  tmp.ordering[0]     =  1000;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "DecayExt";
  tmp.processType     = 6;
  tmp.processSubType  = 231;
  tmp.ordering[0]     =  1000;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "StepLimiter";
  tmp.processType     = 7;
  tmp.processSubType  = 401;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "UsrSepcCuts";
  tmp.processType     = 7;
  tmp.processSubType  = 402;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "NeutronKiller";
  tmp.processType     = 7;
  tmp.processSubType  = 403;
  tmp.ordering[0]     = -1;
  tmp.ordering[1]     = -1;
  tmp.ordering[2]     =  1000;
  tmp.isDuplicable =  false;
  theTable->push_back(tmp);
  sizeOfTable +=1;  

  tmp.processTypeName =  "ParallelWorld";
  tmp.processType     = 10;
  tmp.processSubType  = 491;
  tmp.ordering[0]     = 9900;
  tmp.ordering[1]     = 1;
  tmp.ordering[2]     =  9900;
  tmp.isDuplicable =  true;
  theTable->push_back(tmp);
  sizeOfTable +=1;  
}