Shielding.icc

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 // $Id$
00027 //
00028 //---------------------------------------------------------------------------
00029 //
00030 // ClassName:   
00031 //
00032 // Author: 2010 Tatsumi Koi, Gunter Folger
00033 //
00034 //   created from FTFP_BERT
00035 //
00036 // Modified:
00037 // 16.08.2010 H.Kurashige: Remove inclusion of G4ParticleWithCuts 
00038 // 26.04.2011 T.Koi: Add G4RadioactiveDecayPhysics
00039 // 16.10.2012 A.Ribon: Use new default stopping
00040 //
00041 //----------------------------------------------------------------------------
00042 //
00043 
00044 #include <iomanip>   
00045 
00046 #include "globals.hh"
00047 #include "G4ios.hh"
00048 #include "G4ProcessManager.hh"
00049 #include "G4ProcessVector.hh"
00050 #include "G4ParticleTypes.hh"
00051 #include "G4ParticleTable.hh"
00052 
00053 #include "G4Material.hh"
00054 #include "G4MaterialTable.hh"
00055 
00056 #include "G4DecayPhysics.hh"
00057 #include "G4RadioactiveDecayPhysics.hh"
00058 #include "G4EmStandardPhysics.hh"
00059 #include "G4EmExtraPhysics.hh"
00060 #include "G4IonQMDPhysics.hh"
00061 #include "G4StoppingPhysics.hh"
00062 #include "G4HadronElasticPhysicsHP.hh"
00063 #include "G4HadronElasticPhysicsLEND.hh"
00064 
00065 #include "G4DataQuestionaire.hh"
00066 #include "HadronPhysicsShielding.hh"
00067 
00068 //template<class T> TShielding<T>::TShielding(G4int ver):  T()
00069 template<class T> TShielding<T>::TShielding( G4int verbose, G4String LEN_model ):  T()
00070 {
00071   // default cut value  (1.0mm) 
00072   // defaultCutValue = 1.0*CLHEP::mm;
00073   G4DataQuestionaire it(photon, neutron, radioactive);
00074   G4cout << "<<< Geant4 Physics List simulation engine: Shielding 2.0"<<G4endl;
00075   G4cout <<G4endl;
00076   this->defaultCutValue = 0.7*CLHEP::mm;  
00077   this->SetVerboseLevel(verbose);
00078 
00079  // EM Physics
00080   this->RegisterPhysics( new G4EmStandardPhysics(verbose));
00081 
00082   // Synchroton Radiation & GN Physics
00083   this->RegisterPhysics( new G4EmExtraPhysics(verbose) );
00084 
00085   // Decays 
00086   this->RegisterPhysics( new G4DecayPhysics(verbose) );
00087   //if ( rad == true ) this->RegisterPhysics( new G4RadioactiveDecayPhysics(verbose) );
00088   this->RegisterPhysics( new G4RadioactiveDecayPhysics(verbose) );
00089 
00090   size_t find = LEN_model.find("LEND__");
00091 G4cout << LEN_model << G4endl;
00092   G4String evaluation;
00093   if ( find != G4String::npos )
00094   {
00095 G4cout << LEN_model << G4endl;
00096       evaluation=LEN_model;
00097       evaluation.erase(0,find+6);
00098       LEN_model="LEND";
00099   }
00100 
00101   // Hadron Elastic scattering
00102   if ( LEN_model == "HP" ) 
00103   {
00104      this->RegisterPhysics( new G4HadronElasticPhysicsHP(verbose) );
00105   }
00106   else if ( LEN_model == "LEND" ) 
00107   {
00108      this->RegisterPhysics( new G4HadronElasticPhysicsLEND(verbose,evaluation) );
00109      G4DataQuestionaire itt(lend);
00110   }
00111   else 
00112   {
00113      G4cout << "Shielding Physics List: Warning!" <<G4endl;
00114      G4cout << "\"" << LEN_model << "\" is not valid for the low energy neutorn model." <<G4endl;
00115      G4cout << "Neutron HP package will be used." <<G4endl;
00116      this->RegisterPhysics( new G4HadronElasticPhysicsHP(verbose) );
00117   } 
00118 
00119    // Hadron Physics
00120   HadronPhysicsShielding* hps = new HadronPhysicsShielding(verbose);
00121   if ( LEN_model == "HP" ) 
00122   {
00123      ;
00124   }
00125   else if ( LEN_model == "LEND" ) 
00126   {
00127      hps->UseLEND(evaluation); 
00128   }
00129   else 
00130   {
00131      //G4cout << "Shielding Physics List: Warning." <<G4endl;
00132      //G4cout << "Name of Low Energy Neutron model " << LEN_model << " is invalid." <<G4endl;
00133      //G4cout << "Will use neutron HP package." <<G4endl;
00134   }
00135   this->RegisterPhysics( hps );
00136   //this->RegisterPhysics( new HadronPhysicsShielding(verbose,lend));
00137 
00138 
00139   // Stopping Physics
00140   this->RegisterPhysics( new G4StoppingPhysics(verbose) );
00141 
00142   // Ion Physics
00143   this->RegisterPhysics( new G4IonQMDPhysics(verbose));
00144   
00145   // Neutron tracking cut --> not by default
00146   // this->RegisterPhysics( new G4NeutronTrackingCut(verbose));
00147 
00148 }
00149 
00150 template<class T> TShielding<T>::~TShielding()
00151 {
00152 }
00153 
00154 template<class T> void TShielding<T>::SetCuts()
00155 {
00156   if (this->verboseLevel >1){
00157     G4cout << "Shielding::SetCuts:";
00158   }  
00159   //  " G4VUserPhysicsList::SetCutsWithDefault" method sets 
00160   //   the default cut value for all particle types 
00161 
00162   this->SetCutsWithDefault();   
00163  
00164 //  if (this->verboseLevel > 0)
00165 //    G4VUserPhysicsList::DumpCutValuesTable();  
00166 }

Generated on Mon May 27 17:50:34 2013 for Geant4 by  doxygen 1.4.7