Geant4  10.03.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
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"
47 
48 #include "G4ParticleDefinition.hh"
49 #include "G4ParticleTable.hh"
50 #include "G4Gamma.hh"
51 #include "G4Electron.hh"
52 #include "G4Positron.hh"
53 #include "G4MuonPlus.hh"
54 #include "G4MuonMinus.hh"
55 
58 #include "G4MuonNuclearProcess.hh"
59 #include "G4MuonVDNuclearModel.hh"
60 
62 #include "G4AnnihiToMuPair.hh"
63 #include "G4eeToHadrons.hh"
64 
65 #include "G4PhysicsListHelper.hh"
66 #include "G4BuilderType.hh"
67 #include "G4AutoDelete.hh"
68 
69 // factory
71 //
73 
74 G4bool G4EmExtraPhysics::gnActivated = true;
75 G4bool G4EmExtraPhysics::munActivated = true;
76 G4bool G4EmExtraPhysics::synActivated = false;
77 G4bool G4EmExtraPhysics::synActivatedForAll = false;
78 G4bool G4EmExtraPhysics::gmumuActivated = false;
79 G4bool G4EmExtraPhysics::pmumuActivated = false;
80 G4bool G4EmExtraPhysics::phadActivated = false;
81 
82 G4ThreadLocal G4BertiniElectroNuclearBuilder* G4EmExtraPhysics::theGNPhysics = nullptr;
83 G4ThreadLocal G4SynchrotronRadiation* G4EmExtraPhysics::theSynchRad = nullptr;
84 G4ThreadLocal G4GammaConversionToMuons* G4EmExtraPhysics::theGammaToMuMu = nullptr;
85 G4ThreadLocal G4AnnihiToMuPair* G4EmExtraPhysics::thePosiToMuMu = nullptr;
86 G4ThreadLocal G4eeToHadrons* G4EmExtraPhysics::thePosiToHadrons = nullptr;
87 
89  G4VPhysicsConstructor("G4GammaLeptoNuclearPhys"),
90  verbose(ver)
91 {
92  theMessenger = new G4EmMessenger(this);
94  if(verbose > 1) G4cout << "### G4EmExtraPhysics" << G4endl;
95 }
96 
98  : G4EmExtraPhysics(1)
99 {}
100 
102 {
103  delete theMessenger;
104  theMessenger = nullptr;
105 }
106 
108 {
109  synActivated = val;
110 }
111 
113 {
114  synActivatedForAll = val;
115  if(synActivatedForAll) { synActivated = true; }
116 }
117 
119 {
120  gnActivated = val;
121 }
122 
124 {
125  munActivated = val;
126 }
127 
129 {
130  gmumuActivated = val;
131 }
132 
134 {
135  pmumuActivated = val;
136 }
137 
139 {
140  phadActivated = val;
141 }
142 
144 {
145  G4Gamma::Gamma();
150 }
151 
153 {
159 
161  if(gnActivated) {
162  theGNPhysics = new G4BertiniElectroNuclearBuilder();
163  theGNPhysics->Build();
164  //G4AutoDelete::Register(theGNPhysics);
165  }
166  if(munActivated) {
167  G4MuonNuclearProcess* muNucProcess = new G4MuonNuclearProcess();
168  G4MuonVDNuclearModel* muNucModel = new G4MuonVDNuclearModel();
169  muNucProcess->RegisterMe(muNucModel);
170  ph->RegisterProcess( muNucProcess, muonplus);
171  ph->RegisterProcess( muNucProcess, muonminus);
172  }
173  if(gmumuActivated) {
174  theGammaToMuMu = new G4GammaConversionToMuons();
175  ph->RegisterProcess(theGammaToMuMu, gamma);
176  }
177  if(pmumuActivated) {
178  thePosiToMuMu = new G4AnnihiToMuPair();
179  ph->RegisterProcess(thePosiToMuMu, positron);
180  }
181  if(phadActivated) {
182  thePosiToHadrons = new G4eeToHadrons();
183  ph->RegisterProcess(thePosiToHadrons, positron);
184  }
185  if(synActivated) {
186  theSynchRad = new G4SynchrotronRadiation();
187  ph->RegisterProcess( theSynchRad, electron);
188  ph->RegisterProcess( theSynchRad, positron);
189  //G4AutoDelete::Register(theSynchRad);
190  if(synActivatedForAll) {
191  auto myParticleIterator=GetParticleIterator();
192  myParticleIterator->reset();
193  G4ParticleDefinition* particle = nullptr;
194 
195  while( (*myParticleIterator)() ) {
196  particle = myParticleIterator->value();
197  if( particle->GetPDGStable() && particle->GetPDGCharge() != 0.0) {
198  if(verbose > 1) {
199  G4cout << "### G4SynchrotronRadiation for "
200  << particle->GetParticleName() << G4endl;
201  }
202  ph->RegisterProcess( theSynchRad, particle);
203  }
204  }
205  }
206  }
207 }
208 
virtual ~G4EmExtraPhysics()
static G4MuonPlus * MuonPlus()
Definition: G4MuonPlus.cc:99
void MuonNuclear(G4bool val)
G4bool GetPDGStable() const
#define G4ThreadLocal
Definition: tls.hh:89
int G4int
Definition: G4Types.hh:78
const G4String & GetParticleName() const
void RegisterMe(G4HadronicInteraction *a)
G4ParticleTable::G4PTblDicIterator * GetParticleIterator() const
G4GLOB_DLL std::ostream G4cout
bool G4bool
Definition: G4Types.hh:79
G4bool RegisterProcess(G4VProcess *process, G4ParticleDefinition *particle)
void GammaToMuMu(G4bool val)
G4EmExtraPhysics(G4int ver=1)
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
void PositronToMuMu(G4bool val)
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 PositronToHadrons(G4bool val)
void GammaNuclear(G4bool val)
void Synch(G4bool val)
G4double GetPDGCharge() const
#define G4_DECLARE_PHYSCONSTR_FACTORY(physics_constructor)