Ignore:
Timestamp:
Jan 8, 2010, 3:02:48 PM (16 years ago)
Author:
garnier
Message:

update to geant4.9.3

Location:
trunk/examples/extended/electromagnetic/TestEm11
Files:
41 edited

Legend:

Unmodified
Added
Removed

  • trunk/examples/extended/electromagnetic/TestEm11/GNUmakefile

    r807 r1230  
    1 # $Id: GNUmakefile,v 1.5 2006/09/20 09:29:37 maire Exp $
     1# $Id: GNUmakefile,v 1.10 2008/06/11 22:37:54 maire Exp $
    22# --------------------------------------------------------------
    33# GNUmakefile for examples module.  Gabriele Cosmo, 06/04/98.
    … …  
    1717#### G4ANALYSIS_USE := true
    1818
    19 ifdef G4ANALYSIS_USE
    20   CPPFLAGS += -DG4ANALYSIS_USE
    21 endif
    22 
    2319include $(G4INSTALL)/config/architecture.gmk
    24 
    25 ifdef G4ANALYSIS_USE
    26   # for the aida-config command see the README file
    27   CPPFLAGS += `aida-config --include`
    28   LOADLIBS += `aida-config --lib`
    29 endif
    3020
    3121include $(G4INSTALL)/config/binmake.gmk

  • trunk/examples/extended/electromagnetic/TestEm11/TestEm11.cc

    r807 r1230  
    2525//
    2626// $Id: TestEm11.cc,v 1.3 2006/06/29 16:39:18 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/DetectorConstruction.hh

    r807 r1230  
    2525//
    2626// $Id: DetectorConstruction.hh,v 1.2 2006/06/29 16:39:20 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/DetectorMessenger.hh

    r807 r1230  
    2525//
    2626// $Id: DetectorMessenger.hh,v 1.2 2006/06/29 16:39:22 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/EventAction.hh

    r807 r1230  
    2525//
    2626// $Id: EventAction.hh,v 1.2 2006/06/29 16:39:25 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/EventActionMessenger.hh

    r807 r1230  
    2525//
    2626// $Id: EventActionMessenger.hh,v 1.2 2006/06/29 16:39:27 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/HistoManager.hh

    r807 r1230  
    2525//
    2626// $Id: HistoManager.hh,v 1.5 2007/11/07 17:22:15 maire Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/HistoMessenger.hh

    r807 r1230  
    2424// ********************************************************************
    2525//
    26 // $Id: HistoMessenger.hh,v 1.4 2007/11/07 17:22:16 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26// $Id: HistoMessenger.hh,v 1.5 2008/01/11 13:01:57 jjacquem Exp $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    6363   G4UIcmdWithAString*        factoryCmd;
    6464   G4UIcmdWithAString*        typeCmd;
     65   G4UIcmdWithAString*        optionCmd;
    6566   G4UIcommand*               histoCmd;
    6667   G4UIcmdWithAnInteger*      prhistoCmd;
    6768   G4UIcmdWithAnInteger*      rmhistoCmd;   
    6869   G4UIcmdWithADoubleAndUnit* csdaCmd;
     70 
    6971};
    7072

  • trunk/examples/extended/electromagnetic/TestEm11/include/PhysListEmStandard.hh

    r807 r1230  
    2626//
    2727// $Id: PhysListEmStandard.hh,v 1.2 2006/06/29 16:39:41 gunter Exp $
    28 // GEANT4 tag $Name: $
     28// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/PhysicsList.hh

    r807 r1230  
    2626//
    2727// $Id: PhysicsList.hh,v 1.3 2006/06/29 16:39:43 gunter Exp $
    28 // GEANT4 tag $Name: $
     28// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/PhysicsListMessenger.hh

    r807 r1230  
    2525//
    2626// $Id: PhysicsListMessenger.hh,v 1.2 2006/06/29 16:39:46 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/PrimaryGeneratorAction.hh

    r807 r1230  
    2525//
    2626// $Id: PrimaryGeneratorAction.hh,v 1.2 2006/06/29 16:39:49 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/PrimaryGeneratorMessenger.hh

    r807 r1230  
    2525//
    2626// $Id: PrimaryGeneratorMessenger.hh,v 1.2 2006/06/29 16:39:52 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/RunAction.hh

    r807 r1230  
    2525//
    2626// $Id: RunAction.hh,v 1.3 2007/08/19 20:52:53 maire Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/StepMax.hh

    r807 r1230  
    2525//
    2626// $Id: StepMax.hh,v 1.3 2006/06/29 16:39:57 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/StepMaxMessenger.hh

    r807 r1230  
    2525//
    2626// $Id: StepMaxMessenger.hh,v 1.3 2006/06/29 16:40:00 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/SteppingAction.hh

    r807 r1230  
    2525//
    2626// $Id: SteppingAction.hh,v 1.2 2006/06/29 16:40:03 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/SteppingVerbose.hh

    r807 r1230  
    2525//
    2626// $Id: SteppingVerbose.hh,v 1.2 2006/06/29 16:40:06 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/include/TrackingAction.hh

    r807 r1230  
    2525//
    2626// $Id: TrackingAction.hh,v 1.2 2006/06/29 16:40:08 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/range.mac

    r807 r1230  
    1 # $Id: range.mac,v 1.4 2007/08/19 20:52:53 maire Exp $
     1# $Id: range.mac,v 1.5 2009/09/22 14:20:31 maire Exp $
    22#
    33# Macro file for "TestEm11.cc"
    … …  
    1313/testem/det/setSizeYZ 5 mm
    1414#
    15 /testem/phys/addPhysics  standard
    16 ###/testem/phys/addPhysics  penelope
    17 ###/testem/phys/addPhysics  livermore
     15/testem/phys/addPhysics  local
     16###/testem/phys/addPhysics  empenelope
     17###/testem/phys/addPhysics  emlivermore
    1818#
    1919/run/initialize

  • trunk/examples/extended/electromagnetic/TestEm11/run01.mac

    r807 r1230  
    1 # $Id: run01.mac,v 1.5 2006/07/06 15:56:38 maire Exp $
     1# $Id: run01.mac,v 1.6 2009/09/22 14:20:31 maire Exp $
    22#
    33# limit the step size from histo 1
    … …  
    1010/testem/det/setSizeYZ 1 mm
    1111#
    12 /testem/phys/addPhysics  standard       # em physics
    13 ###/testem/phys/addPhysics  penelope    # em physics
     12/testem/phys/addPhysics  local  # em physics
     13###/testem/phys/addPhysics  empenelope          # em physics
    1414#
    1515/run/initialize

  • trunk/examples/extended/electromagnetic/TestEm11/run02.mac

    r807 r1230  
    1 # $Id: run02.mac,v 1.4 2006/10/26 10:18:59 maire Exp $
     1# $Id: run02.mac,v 1.5 2009/09/22 14:20:31 maire Exp $
    22#
    33# limit the step size by geometry
    … …  
    1111/testem/det/setNbOfLayers 100
    1212#
    13 /testem/phys/addPhysics  standard       # em physics
    14 ###/testem/phys/addPhysics  penelope    # em physics
     13/testem/phys/addPhysics  local  # em physics
     14###/testem/phys/addPhysics  empenelope          # em physics
    1515#
    1616/run/initialize

  • trunk/examples/extended/electromagnetic/TestEm11/sandia.mac

    r807 r1230  
    1 # $Id: sandia.mac,v 1.1 2007/08/19 20:52:53 maire Exp $
     1# $Id: sandia.mac,v 1.9 2009/09/22 14:20:31 maire Exp $
    22#
    33# limit the step size from histo 8
    44#
    5 # G.J.Lockwood et al. Sandia report SAND79-0414.UC-34a, February 1987
     5# G.J.Lockwood et al.
     6#     Sandia report SAND79-0414.UC-34a, February 1987
    67# O.Kadri et al. NIM B 258 (2007) 381
    78#
    … …  
    1011#
    1112/testem/det/setMat G4_Al
     13###/testem/det/setMat G4_Mo
    1214###/testem/det/setMat G4_Ta
    1315/testem/det/setSizeX  1 cm
    1416/testem/det/setSizeYZ 1 cm
    1517#
    16 /testem/phys/addPhysics  standard       # em physics
     18/testem/phys/addPhysics  local  # em physics
    1719#
    1820/run/initialize
    1921#
    2022/gun/particle e-
    21 /gun/energy 500 keV
     23/gun/energy 314 keV
    2224#
    23 /testem/histo/setFileType hbook
    24 /testem/histo/setFileName Al.500keV
     25###/testem/histo/setFileType hbook
     26/testem/histo/setFileType root
     27/testem/histo/setFileName sandia/Al.314keV
    2528/testem/histo/setHisto 8 100 0.  1.0 none       #normalized Edep profile
    26 ###/testem/histo/setHisto 8 100 0.  0.5 none
     29###/testem/histo/setHisto 8 100 0.  0.5 none    #for Ta
    2730#
    2831/testem/event/printModulo 10000

  • trunk/examples/extended/electromagnetic/TestEm11/src/DetectorConstruction.cc

    r807 r1230  
    2424// ********************************************************************
    2525//
    26 // $Id: DetectorConstruction.cc,v 1.8 2007/11/07 17:22:16 maire Exp $
    27 // GEANT4 tag $Name: $
     26// $Id: DetectorConstruction.cc,v 1.9 2009/11/16 15:09:10 maire Exp $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828
    2929//
    … …  
    9191  man->FindOrBuildMaterial("G4_Si", isotopes);
    9292  man->FindOrBuildMaterial("G4_Fe", isotopes);
     93  man->FindOrBuildMaterial("G4_Cu", isotopes); 
    9394  man->FindOrBuildMaterial("G4_Ge", isotopes);
    9495  man->FindOrBuildMaterial("G4_Mo", isotopes);

  • trunk/examples/extended/electromagnetic/TestEm11/src/DetectorMessenger.cc

    r807 r1230  
    2525//
    2626// $Id: DetectorMessenger.cc,v 1.2 2006/06/29 16:40:12 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/EventAction.cc

    r807 r1230  
    2525//
    2626// $Id: EventAction.cc,v 1.2 2006/06/29 16:40:15 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/EventActionMessenger.cc

    r807 r1230  
    2525//
    2626// $Id: EventActionMessenger.cc,v 1.2 2006/06/29 16:40:18 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/HistoManager.cc

    r807 r1230  
    2424// ********************************************************************
    2525//
    26 // $Id: HistoManager.cc,v 1.6 2007/11/07 17:22:16 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26// $Id: HistoManager.cc,v 1.10 2008/09/22 15:11:34 maire Exp $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    5454 
    5555  fileName[0] = "dummy";
    56   fileType    = "hbook";
    57   fileOption  = "--noErrors uncompress"; 
     56  fileType    = "root";
     57  fileOption  = "--noErrors export=root uncompress"; 
    5858  // histograms
    5959  for (G4int k=0; k<MaxHisto; k++) {
    … …  
    223223{
    224224 if (ih > MaxHisto) {
    225     G4cout << "---> warning from HistoManager::RemoveHisto() : histo " << ih
    226            << "does not exist" << G4endl;
     225    G4cout << "---> warning from HistoManager::Scale() : histo " << ih
     226           << "does not exist  (fac = " << fac << ")" << G4endl;
    227227    return;
    228228  }
    … …  
    236236void HistoManager::PrintHisto(G4int ih)
    237237{
    238  if (ih < MaxHisto) ascii[ih] = true;
     238 if (ih < MaxHisto) { ascii[ih] = true; ascii[0] = true; }
    239239 else
    240240    G4cout << "---> warning from HistoManager::PrintHisto() : histo " << ih
    … …  
    249249{
    250250#ifdef G4ANALYSIS_USE
     251
     252 if (!ascii[0]) return;
    251253 
    252254 G4String name = fileName[0] + ".ascii";
    … …  
    262264      for (G4int iBin=0; iBin<Nbins[ih]; ++iBin) {
    263265         File << "  " << iBin << "\t"
    264               << histo[ih]->binMean(iBin) << "\t"
     266              << 0.5*(histo[ih]->axis().binLowerEdge(iBin) +
     267                      histo[ih]->axis().binUpperEdge(iBin)) << "\t"           
    265268              << histo[ih]->binHeight(iBin)
    266269              << G4endl;

  • trunk/examples/extended/electromagnetic/TestEm11/src/HistoMessenger.cc

    r807 r1230  
    2424// ********************************************************************
    2525//
    26 // $Id: HistoMessenger.cc,v 1.5 2007/11/07 17:22:16 maire Exp $
    27 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26// $Id: HistoMessenger.cc,v 1.6 2008/01/11 13:01:57 jjacquem Exp $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    5656  typeCmd->SetGuidance("set histograms file type: hbook, root, XML");
    5757  typeCmd->SetCandidates("hbook root XML");
     58
     59  optionCmd = new G4UIcmdWithAString("/testem/histo/setFileOption",this);
     60  optionCmd->SetGuidance("set option for the histograms file");
    5861
    5962  histoCmd = new G4UIcommand("/testem/histo/setHisto",this);
    … …  
    112115  delete factoryCmd;
    113116  delete histoDir;
     117  delete optionCmd;
    114118}
    115119
    … …  
    123127  if (command == typeCmd)
    124128    histoManager->SetFileType(newValues);
     129
     130  if (command == optionCmd)
     131    histoManager->SetFileOption(newValues);
    125132
    126133  if (command == histoCmd)

  • trunk/examples/extended/electromagnetic/TestEm11/src/PhysListEmStandard.cc

    r807 r1230  
    2424// ********************************************************************
    2525//
    26 //
    27 // $Id: PhysListEmStandard.cc,v 1.3 2006/06/29 16:40:34 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-01-patch-02 $
     26// $Id: PhysListEmStandard.cc,v 1.11 2009/11/16 15:09:10 maire Exp $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2928//
    3029//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    3938#include "G4PhotoElectricEffect.hh"
    4039
    41 #include "G4MultipleScattering.hh"
    42 
     40#include "G4eMultipleScattering.hh"
     41#include "G4UrbanMscModel93.hh"
    4342#include "G4eIonisation.hh"
    4443#include "G4eBremsstrahlung.hh"
    4544#include "G4eplusAnnihilation.hh"
    4645
     46#include "G4MuMultipleScattering.hh"
    4747#include "G4MuIonisation.hh"
    4848#include "G4MuBremsstrahlung.hh"
    4949#include "G4MuPairProduction.hh"
    5050
     51#include "G4hMultipleScattering.hh"
    5152#include "G4hIonisation.hh"
     53#include "G4hBremsstrahlung.hh"
     54#include "G4hPairProduction.hh"
     55
    5256#include "G4ionIonisation.hh"
     57#include "G4IonParametrisedLossModel.hh"
     58#include "G4NuclearStopping.hh"
     59
     60#include "G4EmProcessOptions.hh"
     61#include "G4MscStepLimitType.hh"
    5362
    5463//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    6877{
    6978  // Add standard EM Processes
    70  
     79  //
     80
    7181  theParticleIterator->reset();
    7282  while( (*theParticleIterator)() ){
    … …  
    7686     
    7787    if (particleName == "gamma") {
    78       // gamma         
     88      // gamma
    7989      pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
    8090      pmanager->AddDiscreteProcess(new G4ComptonScattering);
    … …  
    8393    } else if (particleName == "e-") {
    8494      //electron
    85       pmanager->AddProcess(new G4MultipleScattering, -1, 1,1);
    86       pmanager->AddProcess(new G4eIonisation,        -1, 2,2);
    87       pmanager->AddProcess(new G4eBremsstrahlung,    -1, 3,3);
    88            
     95      G4eMultipleScattering* msc = new G4eMultipleScattering();
     96      msc->AddEmModel(0, new G4UrbanMscModel93());
     97      pmanager->AddProcess(msc,                       -1, 1, 1);     
     98      pmanager->AddProcess(new G4eIonisation,         -1, 2, 2);
     99      pmanager->AddProcess(new G4eBremsstrahlung,     -1, 3, 3);
     100           
    89101    } else if (particleName == "e+") {
    90102      //positron
    91       pmanager->AddProcess(new G4MultipleScattering, -1, 1,1);
    92       pmanager->AddProcess(new G4eIonisation,        -1, 2,2);
    93       pmanager->AddProcess(new G4eBremsstrahlung,    -1, 3,3);
    94       pmanager->AddProcess(new G4eplusAnnihilation,   0,-1,4);
     103      G4eMultipleScattering* msc = new G4eMultipleScattering();
     104      msc->AddEmModel(0, new G4UrbanMscModel93());
     105      pmanager->AddProcess(msc,                       -1, 1, 1);           
     106      pmanager->AddProcess(new G4eIonisation,         -1, 2, 2);
     107      pmanager->AddProcess(new G4eBremsstrahlung,     -1, 3, 3);
     108      pmanager->AddProcess(new G4eplusAnnihilation,    0,-1, 4);
    95109           
    96     } else if( particleName == "mu+" ||
     110    } else if (particleName == "mu+" ||
    97111               particleName == "mu-"    ) {
    98112      //muon 
    99       pmanager->AddProcess(new G4MultipleScattering,-1, 1,1);
    100       pmanager->AddProcess(new G4MuIonisation,      -1, 2,2);
    101       pmanager->AddProcess(new G4MuBremsstrahlung,  -1, 3,3);
    102       pmanager->AddProcess(new G4MuPairProduction,  -1, 4,4);       
     113      pmanager->AddProcess(new G4MuMultipleScattering, -1, 1, 1);
     114      pmanager->AddProcess(new G4MuIonisation,         -1, 2, 2);
     115      pmanager->AddProcess(new G4MuBremsstrahlung,     -1, 3, 3);
     116      pmanager->AddProcess(new G4MuPairProduction,     -1, 4, 4);
     117             
     118    } else if( particleName == "proton" ||
     119               particleName == "pi-" ||
     120               particleName == "pi+"    ) {
     121      //proton 
     122      pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
     123      pmanager->AddProcess(new G4hIonisation,         -1, 2, 2);
     124      pmanager->AddProcess(new G4hBremsstrahlung,     -1, 3, 3);
     125      pmanager->AddProcess(new G4hPairProduction,     -1, 4, 4);       
    103126     
    104     } else if( particleName == "GenericIon" ) {
    105       pmanager->AddProcess(new G4MultipleScattering,-1, 1,1);
    106       pmanager->AddProcess(new G4ionIonisation,     -1, 2,2);
    107 
     127    } else if( particleName == "alpha" ||
     128               particleName == "He3"    ) {
     129      //alpha
     130      pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
     131      pmanager->AddProcess(new G4ionIonisation,       -1, 2, 2);
     132      pmanager->AddProcess(new G4NuclearStopping,     -1, 3,-1);
     133           
     134    } else if( particleName == "GenericIon" ) {
     135      //Ions
     136      pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
     137      G4ionIonisation* ionIoni = new G4ionIonisation();
     138      ionIoni->SetEmModel(new G4IonParametrisedLossModel());
     139      pmanager->AddProcess(ionIoni,                   -1, 2, 2);     
     140      pmanager->AddProcess(new G4NuclearStopping,     -1, 3,-1);     
     141     
    108142    } else if ((!particle->IsShortLived()) &&
    109143               (particle->GetPDGCharge() != 0.0) &&
    110144               (particle->GetParticleName() != "chargedgeantino")) {
    111145      //all others charged particles except geantino
    112       pmanager->AddProcess(new G4MultipleScattering,-1,1,1);
    113       pmanager->AddProcess(new G4hIonisation,       -1,2,2);
     146      pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
     147      pmanager->AddProcess(new G4hIonisation,         -1, 2, 2);
    114148    }
    115149  }
     150
     151  // Em options
     152  //
     153  // Main options and setting parameters are shown here.
     154  // Several of them have default values.
     155  //
     156  G4EmProcessOptions emOptions;
     157 
     158  //physics tables
     159  //
     160  emOptions.SetMinEnergy(100*eV);       //default   
     161  emOptions.SetMaxEnergy(100*TeV);      //default 
     162  emOptions.SetDEDXBinning(12*20);      //default=12*7
     163  emOptions.SetLambdaBinning(12*20);    //default=12*7
     164  emOptions.SetSplineFlag(true);        //default
     165     
     166  //multiple coulomb scattering
     167  //
     168  emOptions.SetMscStepLimitation(fUseDistanceToBoundary);  //default=fUseSafety
     169  emOptions.SetMscRangeFactor(0.04);    //default
     170  emOptions.SetMscGeomFactor (2.5);     //default       
     171  emOptions.SetSkin(3.);                //default
     172     
     173  //energy loss
     174  //
     175  emOptions.SetStepFunction(0.2, 100*um);       //default=(0.2, 1*mm)   
     176  emOptions.SetLinearLossLimit(1.e-2);          //default
     177   
     178  //ionization
     179  //
     180  emOptions.SetSubCutoff(false);        //default
    116181}
    117182

  • trunk/examples/extended/electromagnetic/TestEm11/src/PhysicsList.cc

    r807 r1230  
    2424// ********************************************************************
    2525//
    26 // $Id: PhysicsList.cc,v 1.5 2006/10/26 10:18:59 maire Exp $
    27 // GEANT4 tag $Name: $
     26// $Id: PhysicsList.cc,v 1.7 2009/09/22 14:20:31 maire Exp $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    … …  
    3434
    3535#include "PhysListEmStandard.hh"
    36 #include "PhysListEmLivermore.hh"
    37 #include "PhysListEmPenelope.hh"
     36
     37#include "G4EmStandardPhysics.hh"
     38#include "G4EmStandardPhysics_option1.hh"
     39#include "G4EmStandardPhysics_option2.hh"
     40#include "G4EmStandardPhysics_option3.hh"
     41#include "G4EmLivermorePhysics.hh"
     42#include "G4EmPenelopePhysics.hh"
    3843
    3944#include "G4LossTableManager.hh"
    … …  
    5762
    5863  // EM physics
    59   emPhysicsList = new PhysListEmStandard(emName = "standard");
     64  emPhysicsList = new PhysListEmStandard(emName = "local");
    6065
    6166}
    … …  
    186191  G4EmProcessOptions emOptions;
    187192  emOptions.SetBuildCSDARange(true);
    188  
     193  emOptions.SetDEDXBinningForCSDARange(8*20);
     194   
    189195  // decay process
    190196  //
    … …  
    206212  if (name == emName) return;
    207213
    208   if (name == "standard") {
     214  if (name == "local") {
    209215
    210216    emName = name;
    211217    delete emPhysicsList;
    212218    emPhysicsList = new PhysListEmStandard(name);
     219
     220  } else if (name == "emstandard_opt0") {
     221
     222    emName = name;
     223    delete emPhysicsList;
     224    emPhysicsList = new G4EmStandardPhysics();
    213225   
    214   } else if (name == "livermore") {
    215 
    216     emName = name;
    217     delete emPhysicsList;
    218     emPhysicsList = new PhysListEmLivermore(name);
     226  } else if (name == "emstandard_opt1") {
     227
     228    emName = name;
     229    delete emPhysicsList;
     230    emPhysicsList = new G4EmStandardPhysics_option1();
     231
     232  } else if (name == "emstandard_opt2") {
     233
     234    emName = name;
     235    delete emPhysicsList;
     236    emPhysicsList = new G4EmStandardPhysics_option2();
     237       
     238  } else if (name == "emstandard_opt3") {
     239
     240    emName = name;
     241    delete emPhysicsList;
     242    emPhysicsList = new G4EmStandardPhysics_option3();
     243           
     244  } else if (name == "emlivermore") {
     245    emName = name;
     246    delete emPhysicsList;
     247    emPhysicsList = new G4EmLivermorePhysics();
    219248   
    220   } else if (name == "penelope") {
    221 
    222     emName = name;
    223     delete emPhysicsList;
    224     emPhysicsList = new PhysListEmPenelope(name);
     249  } else if (name == "empenelope") {
     250    emName = name;
     251    delete emPhysicsList;
     252    emPhysicsList = new G4EmPenelopePhysics();
    225253       
    226254  } else {

  • trunk/examples/extended/electromagnetic/TestEm11/src/PhysicsListMessenger.cc

    r807 r1230  
    2525//
    2626// $Id: PhysicsListMessenger.cc,v 1.2 2006/06/29 16:40:38 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/PrimaryGeneratorAction.cc

    r807 r1230  
    2525//
    2626// $Id: PrimaryGeneratorAction.cc,v 1.2 2006/06/29 16:40:40 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/PrimaryGeneratorMessenger.cc

    r807 r1230  
    2525//
    2626// $Id: PrimaryGeneratorMessenger.cc,v 1.2 2006/06/29 16:40:43 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/RunAction.cc

    r807 r1230  
    2525//
    2626// $Id: RunAction.cc,v 1.7 2007/08/19 20:52:53 maire Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/StepMax.cc

    r807 r1230  
    2525//
    2626// $Id: StepMax.cc,v 1.3 2006/06/29 16:40:47 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/StepMaxMessenger.cc

    r807 r1230  
    2525//
    2626// $Id: StepMaxMessenger.cc,v 1.3 2006/06/29 16:40:53 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/SteppingAction.cc

    r807 r1230  
    2525//
    2626// $Id: SteppingAction.cc,v 1.5 2007/08/19 20:52:53 maire Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/SteppingVerbose.cc

    r807 r1230  
    2525//
    2626// $Id: SteppingVerbose.cc,v 1.2 2006/06/29 16:41:54 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/src/TrackingAction.cc

    r807 r1230  
    2525//
    2626// $Id: TrackingAction.cc,v 1.2 2006/06/29 16:41:56 gunter Exp $
    27 // GEANT4 tag $Name: $
     27// GEANT4 tag $Name: geant4-09-03-cand-01 $
    2828//
    2929//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

  • trunk/examples/extended/electromagnetic/TestEm11/water.mac

    r807 r1230  
    1 # $Id: water.mac,v 1.4 2006/09/18 17:26:20 maire Exp $
     1# $Id: water.mac,v 1.5 2009/09/22 14:20:31 maire Exp $
    22#
    33# no constraint on tracking step
    … …  
    1010/testem/det/setSizeYZ 3 cm
    1111#
    12 /testem/phys/addPhysics  standard       # em physics
    13 ###/testem/phys/addPhysics  g4v52       # em physics
     12/testem/phys/addPhysics  local  # em physics
    1413#
    1514/run/initialize
Note: See TracChangeset for help on using the changeset viewer.