Ignore:
Timestamp:
Nov 5, 2010, 4:08:39 PM (14 years ago)
Author:
garnier
Message:

update ti head

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/examples/novice/N06/src/ExN06PhysicsList.cc

    r1337 r1342  
    2525//
    2626//
    27 // $Id: ExN06PhysicsList.cc,v 1.17 2009/11/10 05:16:23 gum Exp $
    28 // GEANT4 tag $Name: geant4-09-04-beta-01 $
     27// $Id: ExN06PhysicsList.cc,v 1.19 2010/10/23 19:14:03 gum Exp $
     28// GEANT4 tag $Name: examples-V09-03-09 $
    2929//
    3030//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    3131//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
    32 
    33 #include "G4ios.hh"
    34 #include <iomanip>
    3532
    3633#include "globals.hh"
    … …  
    4239#include "G4ParticleTable.hh"
    4340
    44 #include "G4Material.hh"
    45 #include "G4MaterialTable.hh"
    46 
    4741#include "G4ProcessManager.hh"
    48 #include "G4ProcessVector.hh"
    4942
    5043#include "G4Cerenkov.hh"
    … …  
    5245#include "G4OpAbsorption.hh"
    5346#include "G4OpRayleigh.hh"
     47#include "G4OpMieHG.hh"
    5448#include "G4OpBoundaryProcess.hh"
    5549
    … …  
    6155ExN06PhysicsList::ExN06PhysicsList() :  G4VUserPhysicsList()
    6256{
    63   theCerenkovProcess           = 0;
    64   theScintillationProcess      = 0;
    65   theAbsorptionProcess         = 0;
    66   theRayleighScatteringProcess = 0;
    67   theBoundaryProcess           = 0;
     57  theCerenkovProcess           = NULL;
     58  theScintillationProcess      = NULL;
     59  theAbsorptionProcess         = NULL;
     60  theRayleighScatteringProcess = NULL;
     61  theMieHGScatteringProcess    = NULL;
     62  theBoundaryProcess           = NULL;
    6863 
    6964  pMessenger = new ExN06PhysicsListMessenger(this); 
    … …  
    110105{
    111106  // leptons
     107  //  e+/-
    112108  G4Electron::ElectronDefinition();
    113109  G4Positron::PositronDefinition();
     110  // mu+/-
     111  G4MuonPlus::MuonPlusDefinition();
     112  G4MuonMinus::MuonMinusDefinition();
     113  // nu_e
    114114  G4NeutrinoE::NeutrinoEDefinition();
    115115  G4AntiNeutrinoE::AntiNeutrinoEDefinition();
    116   G4MuonPlus::MuonPlusDefinition();
    117   G4MuonMinus::MuonMinusDefinition();
     116  // nu_mu
    118117  G4NeutrinoMu::NeutrinoMuDefinition();
    119118  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
    … …  
    124123void ExN06PhysicsList::ConstructMesons()
    125124{
    126  //  mesons
     125  //  mesons
    127126  G4PionPlus::PionPlusDefinition();
    128127  G4PionMinus::PionMinusDefinition();
    … …  
    134133void ExN06PhysicsList::ConstructBaryons()
    135134{
    136 //  barions
     135  //  barions
    137136  G4Proton::ProtonDefinition();
    138137  G4AntiProton::AntiProtonDefinition();
     138
    139139  G4Neutron::NeutronDefinition();
    140140  G4AntiNeutron::AntiNeutronDefinition();
    … …  
    207207    // gamma
    208208      // Construct processes for gamma
    209 //      pmanager->AddDiscreteProcess(new G4GammaConversion());
    210 //      pmanager->AddDiscreteProcess(new G4ComptonScattering());
    211 //      pmanager->AddDiscreteProcess(new G4PhotoElectricEffect());
     209      pmanager->AddDiscreteProcess(new G4GammaConversion());
     210      pmanager->AddDiscreteProcess(new G4ComptonScattering());
     211      pmanager->AddDiscreteProcess(new G4PhotoElectricEffect());
    212212
    213213    } else if (particleName == "e-") {
    … …  
    238238      if ((particle->GetPDGCharge() != 0.0) &&
    239239          (particle->GetParticleName() != "chargedgeantino")) {
    240      // all others charged particles except geantino
     240       // all others charged particles except geantino
    241241       pmanager->AddProcess(new G4hMultipleScattering(),-1,1,1);
    242242       pmanager->AddProcess(new G4hIonisation(),       -1,2,2);
    … …  
    251251{
    252252  theCerenkovProcess           = new G4Cerenkov("Cerenkov");
    253   theScintillationProcess = new G4Scintillation("Scintillation");
    254   theAbsorptionProcess     = new G4OpAbsorption();
     253  theScintillationProcess      = new G4Scintillation("Scintillation");
     254  theAbsorptionProcess         = new G4OpAbsorption();
    255255  theRayleighScatteringProcess = new G4OpRayleigh();
    256   theBoundaryProcess  = new G4OpBoundaryProcess();
     256  theMieHGScatteringProcess    = new G4OpMieHG();
     257  theBoundaryProcess           = new G4OpBoundaryProcess();
    257258
    258259//  theCerenkovProcess->DumpPhysicsTable();
    259260//  theScintillationProcess->DumpPhysicsTable();
    260 //  theAbsorptionProcess->DumpPhysicsTable();
    261261//  theRayleighScatteringProcess->DumpPhysicsTable();
    262262
    … …  
    296296      pmanager->AddDiscreteProcess(theAbsorptionProcess);
    297297      pmanager->AddDiscreteProcess(theRayleighScatteringProcess);
     298      pmanager->AddDiscreteProcess(theMieHGScatteringProcess);
    298299      pmanager->AddDiscreteProcess(theBoundaryProcess);
    299300    }
    … …  
    309310  theAbsorptionProcess->SetVerboseLevel(verbose);
    310311  theRayleighScatteringProcess->SetVerboseLevel(verbose);
     312  theMieHGScatteringProcess->SetVerboseLevel(verbose);
    311313  theBoundaryProcess->SetVerboseLevel(verbose); 
    312314}
Note: See TracChangeset for help on using the changeset viewer.