| Geant4 Cross Reference |
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 //---------------------------------------------------------------------------
27 //
28 // ClassName: G4StoppingPhysicsFritiofWithBinaryCascade
29 //
30 // Author: Alberto Ribon
31 //
32 // Date: July 2019
33 //
34 // Modified:
35 //
36 //----------------------------------------------------------------------------
37
38 #include "G4StoppingPhysicsFritiofWithBinaryCascade.hh"
39 #include "G4SystemOfUnits.hh"
40 #include "G4HadronicAbsorptionBertini.hh"
41 #include "G4HadronicAbsorptionFritiof.hh"
42 #include "G4HadronicAbsorptionFritiofWithBinaryCascade.hh"
43 #include "G4MuonMinusCapture.hh"
44 #include "G4ParticleDefinition.hh"
45 #include "G4ProcessManager.hh"
46 #include "G4LeptonConstructor.hh"
47 #include "G4MesonConstructor.hh"
48 #include "G4BaryonConstructor.hh"
49 #include "G4MuonMinus.hh"
50 #include "G4PionMinus.hh"
51 #include "G4BuilderType.hh"
52
53 // factory
54 #include "G4PhysicsConstructorFactory.hh"
55 //
56 G4_DECLARE_PHYSCONSTR_FACTORY( G4StoppingPhysicsFritiofWithBinaryCascade );
57
58 G4StoppingPhysicsFritiofWithBinaryCascade::
59 G4StoppingPhysicsFritiofWithBinaryCascade( G4int ver ) :
60 G4StoppingPhysicsFritiofWithBinaryCascade("stopping", ver)
61 {}
62
63 G4StoppingPhysicsFritiofWithBinaryCascade::
64 G4StoppingPhysicsFritiofWithBinaryCascade( const G4String& name, G4int ver,
65 G4bool UseMuonMinusCapture ) :
66 G4VPhysicsConstructor( name ), verbose( ver ),
67 useMuonMinusCapture( UseMuonMinusCapture )
68 {
69 SetPhysicsType(bStopping);
70 if(verbose > 1) G4cout << "### G4StoppingPhysicsFritiofWithBinaryCascade"
71 << G4endl;
72 }
73
74 G4StoppingPhysicsFritiofWithBinaryCascade::~G4StoppingPhysicsFritiofWithBinaryCascade() {}
75
76
77 void G4StoppingPhysicsFritiofWithBinaryCascade::ConstructParticle() {
78 // G4cout << "G4StoppingPhysicsFritiofWithBinaryCascade::ConstructParticle" << G4endl;
79 G4LeptonConstructor pLeptonConstructor;
80 pLeptonConstructor.ConstructParticle();
81
82 G4MesonConstructor pMesonConstructor;
83 pMesonConstructor.ConstructParticle();
84
85 G4BaryonConstructor pBaryonConstructor;
86 pBaryonConstructor.ConstructParticle();
87 }
88
89 void G4StoppingPhysicsFritiofWithBinaryCascade::ConstructProcess() {
90 if ( verbose > 1 ) G4cout << "### G4StoppingPhysicsFritiofWithBinaryCascade::ConstructProcess " << G4endl;
91
92 G4MuonMinusCapture* muProcess = nullptr;
93 if ( useMuonMinusCapture ) {
94 muProcess = new G4MuonMinusCapture();
95 }
96
97 G4HadronicAbsorptionBertini* hBertiniProcess = new G4HadronicAbsorptionBertini;
98 G4HadronicAbsorptionFritiof* hFritiofProcess = new G4HadronicAbsorptionFritiof;
99 G4HadronicAbsorptionFritiofWithBinaryCascade* hFritiofWithBinaryCascadeProcess =
100 new G4HadronicAbsorptionFritiofWithBinaryCascade;
101
102 G4double mThreshold = 130.0*MeV;
103
104 // Add Stopping Process
105 G4ParticleDefinition* particle = 0;
106 G4ProcessManager* pmanager = 0;
107
108 auto myParticleIterator=GetParticleIterator();
109 myParticleIterator->reset();
110
111 while ( (*myParticleIterator)() ) {
112
113 particle = myParticleIterator->value();
114 pmanager = particle->GetProcessManager();
115
116 if (useMuonMinusCapture && particle == G4MuonMinus::MuonMinus() ) {
117 pmanager->AddRestProcess( muProcess );
118 if ( verbose > 1 ) {
119 G4cout << "### G4MuonMinusCapture added for "
120 << particle->GetParticleName() << G4endl;
121 }
122 }
123
124 if ( particle->GetPDGCharge() <= 0.0 &&
125 particle->GetPDGMass() > mThreshold &&
126 ( ! particle->IsShortLived() ) ) {
127
128 // Use Fritiof/BinaryCascade for: anti-proton and anti-neutron
129 if ( particle == G4AntiProton::Definition() ||
130 particle == G4AntiNeutron::Definition() ) {
131 if ( hFritiofWithBinaryCascadeProcess->IsApplicable( *particle ) ) {
132 pmanager->AddRestProcess( hFritiofWithBinaryCascadeProcess );
133 if ( verbose > 1 ) {
134 G4cout << "### G4HadronicAbsorptionFritiofWithBinaryCascade added for "
135 << particle->GetParticleName() << G4endl;
136 }
137 }
138
139 // Use Fritiof/Precompound for:
140 // anti-lambda, anti-sigma0, anti-sigma+, anti-xi0 and anti-nuclei
141 } else if ( particle == G4AntiLambda::Definition() ||
142 particle == G4AntiSigmaZero::Definition() ||
143 particle == G4AntiSigmaPlus::Definition() ||
144 particle == G4AntiXiZero::Definition() ||
145 particle->GetBaryonNumber() < -1 ) { // Anti-nuclei
146 if ( hFritiofProcess->IsApplicable( *particle ) ) {
147 pmanager->AddRestProcess( hFritiofProcess );
148 if ( verbose > 1 ) {
149 G4cout << "### G4HadronicAbsorptionFritiof added for "
150 << particle->GetParticleName() << G4endl;
151 }
152 }
153
154 // Use Bertini for pi-, K-, Sigma-, Xi-, and Omega-
155 } else if ( particle == G4PionMinus::Definition() ||
156 particle == G4KaonMinus::Definition() ||
157 particle == G4SigmaMinus::Definition() ||
158 particle == G4XiMinus::Definition() ||
159 particle == G4OmegaMinus::Definition() ) {
160 if ( hBertiniProcess->IsApplicable( *particle ) ) {
161 pmanager->AddRestProcess( hBertiniProcess );
162 if ( verbose > 1 ) {
163 G4cout << "### G4HadronicAbsorptionBertini added for "
164 << particle->GetParticleName() << G4endl;
165 }
166 }
167
168 } else {
169 if ( verbose > 1 ) {
170 G4cout << "WARNING in G4StoppingPhysicsFritiofWithBinaryCascade::ConstructProcess: \
171 not able to deal with nuclear stopping of "
172 << particle->GetParticleName() << G4endl;
173 }
174 }
175 }
176 } // end of while loop
177 }
178