Geant4  10.02.p01
G4EmExtraPhysics.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4EmExtraPhysics.cc 66704 2013-01-10 18:20:17Z gunter $
27 //
28 //---------------------------------------------------------------------------
29 //
30 // ClassName: G4EmExtraPhysics
31 //
32 // Author: 2002 J.P. Wellisch
33 //
34 // Modified:
35 // 10.11.2005 V.Ivanchenko edit to provide a standard
36 // 19.06.2006 V.Ivanchenko add mu-nuclear process
37 // 16.10.2012 A.Ribon: renamed G4EmExtraBertiniPhysics as G4EmExtraPhysics
38 // 10.04.2014 A.Dotti: Add MT functionality for messenger
39 // 24.04.2014 A.Ribon: switched on muon-nuclear by default
40 //
41 //----------------------------------------------------------------------------
42 //
43 
44 #include "G4EmExtraPhysics.hh"
45 
46 #include "G4SystemOfUnits.hh"
48 
49 #include "G4ParticleDefinition.hh"
50 #include "G4ParticleTable.hh"
51 #include "G4Gamma.hh"
52 #include "G4Electron.hh"
53 #include "G4Positron.hh"
54 #include "G4MuonPlus.hh"
55 #include "G4MuonMinus.hh"
56 
59 #include "G4MuonNuclearProcess.hh"
60 #include "G4MuonVDNuclearModel.hh"
61 
62 #include "G4PhysicsListHelper.hh"
63 #include "G4BuilderType.hh"
64 #include "G4AutoDelete.hh"
65 
66 // factory
68 //
70 
75 
78 
80  G4VPhysicsConstructor("G4GammaLeptoNuclearPhys"),
81  verbose(ver)
82 {
83  theMessenger = new G4EmMessenger(this);
85  if(verbose > 1) G4cout << "### G4EmExtraPhysics" << G4endl;
86 }
87 
89  G4VPhysicsConstructor("G4GammaLeptoNuclearPhys"),
90  verbose(1)
91 {
92  theMessenger = new G4EmMessenger(this);
94  if(verbose > 1) G4cout << "### G4EmExtraPhysics" << G4endl;
95 }
96 
98 {
99  delete theMessenger;
100  theMessenger = 0;
101 }
102 
104 {
105  synActivated = val;
106 }
107 
109 {
110  synActivatedForAll = val;
111  if(synActivatedForAll) { synActivated = true; }
112 }
113 
115 {
116  gnActivated = val;
117 }
118 
120 {
121  munActivated = val;
122 }
123 
125 {
126  G4Gamma::Gamma();
131 }
132 
134 {
139 
141  if(gnActivated) {
143  theGNPhysics->Build();
144  //G4AutoDelete::Register(theGNPhysics);
145  }
146  if(munActivated) {
147  G4MuonNuclearProcess* muNucProcess = new G4MuonNuclearProcess();
148  G4MuonVDNuclearModel* muNucModel = new G4MuonVDNuclearModel();
149  muNucProcess->RegisterMe(muNucModel);
150  ph->RegisterProcess( muNucProcess, muonplus);
151  ph->RegisterProcess( muNucProcess, muonminus);
152  }
153  if(synActivated) {
155  ph->RegisterProcess( theSynchRad, electron);
156  ph->RegisterProcess( theSynchRad, positron);
157  //G4AutoDelete::Register(theSynchRad);
158  if(synActivatedForAll) {
159  aParticleIterator->reset();
160  G4ParticleDefinition* particle=0;
161 
162  while( (*aParticleIterator)() ) {
163  particle = aParticleIterator->value();
164  if( particle->GetPDGStable() && particle->GetPDGCharge() != 0.0) {
165  if(verbose > 1) {
166  G4cout << "### G4SynchrotronRadiation for "
167  << particle->GetParticleName() << G4endl;
168  }
169  ph->RegisterProcess( theSynchRad, particle);
170  }
171  }
172  }
173  }
174 }
175 
static G4bool munActivated
virtual ~G4EmExtraPhysics()
static G4MuonPlus * MuonPlus()
Definition: G4MuonPlus.cc:99
static G4ThreadLocal G4SynchrotronRadiation * theSynchRad
void MuonNuclear(G4bool val)
G4bool GetPDGStable() const
static G4bool synActivated
static G4ThreadLocal G4BertiniElectroNuclearBuilder * theGNPhysics
#define G4ThreadLocal
Definition: tls.hh:89
int G4int
Definition: G4Types.hh:78
const G4String & GetParticleName() const
void RegisterMe(G4HadronicInteraction *a)
G4GLOB_DLL std::ostream G4cout
bool G4bool
Definition: G4Types.hh:79
#define aParticleIterator
G4bool RegisterProcess(G4VProcess *process, G4ParticleDefinition *particle)
static G4bool gnActivated
G4EmExtraPhysics(G4int ver=1)
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
static G4Positron * Positron()
Definition: G4Positron.cc:94
static G4PhysicsListHelper * GetPhysicsListHelper()
void SynchAll(G4bool val)
static G4MuonMinus * MuonMinus()
Definition: G4MuonMinus.cc:100
static G4Electron * Electron()
Definition: G4Electron.cc:94
#define G4endl
Definition: G4ios.hh:61
void GammaNuclear(G4bool val)
void Synch(G4bool val)
G4double GetPDGCharge() const
G4EmMessenger * theMessenger
static G4bool synActivatedForAll
G4_DECLARE_PHYSCONSTR_FACTORY(G4EmExtraPhysics)