source: trunk/source/physics_lists/builders/src/G4HadronElasticPhysics.cc @ 1340

Last change on this file since 1340 was 1340, checked in by garnier, 14 years ago

update ti head

File size: 6.9 KB
Line 
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: G4HadronElasticPhysics.cc,v 1.16 2010/09/24 12:54:13 vnivanch Exp $
27// GEANT4 tag $Name: phys-lists-V09-03-34 $
28//
29//---------------------------------------------------------------------------
30//
31// ClassName:   G4HadronElasticPhysics
32//
33// Author: 23 November 2006 V. Ivanchenko
34//
35// Modified:
36// 21.03.2007 V.Ivanchenko Use G4BGGNucleonElasticXS and G4BGGPionElasticXS;
37//                         Reduce thresholds for HE and Q-models to zero
38// 03.06.2010 V.Ivanchenko cleanup constructors and ConstructProcess method
39// 29.07.2010 V.Ivanchenko rename this class from G4HadronHElasticPhysics to
40//                         G4HadronElasticPhysics, old version of the class
41//                         is renamed to G4HadronElasticPhysics93
42//
43//----------------------------------------------------------------------------
44//
45// CHIPS for sampling scattering for p and n
46// Glauber model for samplimg of high energy pi+- (E > 1GeV)
47// LHEP sampling model for the other particle
48// BBG cross sections for p, n and pi+-
49// LHEP cross sections for other particles
50
51#include "G4HadronElasticPhysics.hh"
52
53#include "G4ParticleDefinition.hh"
54#include "G4ProcessManager.hh"
55
56#include "G4MesonConstructor.hh"
57#include "G4BaryonConstructor.hh"
58#include "G4IonConstructor.hh"
59#include "G4Neutron.hh"
60
61#include "G4WHadronElasticProcess.hh"
62#include "G4VHadronElastic.hh"
63#include "G4CHIPSElastic.hh"
64#include "G4ElasticHadrNucleusHE.hh"
65#include "G4BGGNucleonElasticXS.hh"
66#include "G4BGGPionElasticXS.hh"
67#include "G4NeutronElasticXS.hh"
68#include "G4CHIPSElasticXS.hh"
69
70G4HadronElasticPhysics::G4HadronElasticPhysics(G4int ver)
71  : G4VPhysicsConstructor("hElasticWEL_CHIPS"), verbose(ver), 
72    wasActivated(false)
73{
74  if(verbose > 1) { 
75    G4cout << "### G4HadronElasticPhysics: " << GetPhysicsName() 
76           << G4endl; 
77  }
78}
79
80G4HadronElasticPhysics::G4HadronElasticPhysics(const G4String&,
81    G4int ver, G4bool, const G4String&)
82  : G4VPhysicsConstructor("hElasticWEL_CHIPS"), verbose(ver), 
83    wasActivated(false)
84{
85  if(verbose > 1) { 
86    G4cout << "### G4HadronElasticPhysics: " << GetPhysicsName() 
87           << G4endl; 
88  }
89}
90
91G4HadronElasticPhysics::~G4HadronElasticPhysics()
92{}
93
94void G4HadronElasticPhysics::ConstructParticle()
95{
96  // G4cout << "G4HadronElasticPhysics::ConstructParticle" << G4endl;
97  G4MesonConstructor pMesonConstructor;
98  pMesonConstructor.ConstructParticle();
99
100  G4BaryonConstructor pBaryonConstructor;
101  pBaryonConstructor.ConstructParticle();
102
103  //  Construct light ions
104  G4IonConstructor pConstructor;
105  pConstructor.ConstructParticle(); 
106}
107
108void G4HadronElasticPhysics::ConstructProcess()
109{
110  if(wasActivated) { return; }
111  wasActivated = true;
112
113  G4double elimit = 1.0*GeV;
114  if(verbose > 1) {
115    G4cout << "### HadronElasticPhysics Construct Processes with HE limit " 
116           << elimit << " MeV" << G4endl;
117  }
118
119  G4VHadronElastic* plep0 = new G4VHadronElastic();
120  G4VHadronElastic* plep1 = new G4VHadronElastic();
121  plep1->SetMaxEnergy(elimit);
122
123  G4CHIPSElastic* chipsp = new G4CHIPSElastic();
124  G4CHIPSElastic* chipsn = new G4CHIPSElastic();
125
126  G4ElasticHadrNucleusHE* he = new G4ElasticHadrNucleusHE(); 
127  he->SetMinEnergy(elimit);
128
129  theParticleIterator->reset();
130  while( (*theParticleIterator)() )
131  {
132    G4ParticleDefinition* particle = theParticleIterator->value();
133    G4String pname = particle->GetParticleName();
134    if(pname == "anti_lambda"  ||
135       pname == "anti_neutron" ||
136       pname == "anti_omega-"  || 
137       pname == "anti_proton"  || 
138       pname == "anti_sigma-"  || 
139       pname == "anti_sigma+"  || 
140       pname == "anti_xi-"  || 
141       pname == "anti_xi0"  || 
142       pname == "kaon-"     || 
143       pname == "kaon+"     || 
144       pname == "kaon0S"    || 
145       pname == "kaon0L"    || 
146       pname == "lambda"    || 
147       pname == "omega-"    || 
148       pname == "pi-"       || 
149       pname == "pi+"       || 
150       pname == "proton"    || 
151       pname == "sigma-"    || 
152       pname == "sigma+"    || 
153       pname == "xi-"       || 
154       pname == "alpha"     ||
155       pname == "deuteron"  ||
156       pname == "triton") {
157     
158      G4ProcessManager* pmanager = particle->GetProcessManager();
159      G4WHadronElasticProcess* hel = new G4WHadronElasticProcess();
160      if(pname == "proton") { 
161        //hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
162        hel->AddDataSet(new G4CHIPSElasticXS());
163        hel->RegisterMe(chipsp);
164      } else if (pname == "pi+" || pname == "pi-") { 
165        hel->AddDataSet(new G4BGGPionElasticXS(particle));
166        hel->RegisterMe(plep1);
167        hel->RegisterMe(he);
168      } else {
169        hel->RegisterMe(plep0);
170      }
171      pmanager->AddDiscreteProcess(hel);
172      if(verbose > 1) {
173        G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
174               << " added for " << particle->GetParticleName() << G4endl;
175      }
176
177      // neutron case
178    } else if(pname == "neutron") {   
179
180      G4ProcessManager* pmanager = particle->GetProcessManager();
181      G4WHadronElasticProcess* hel = new G4WHadronElasticProcess();
182      //hel->AddDataSet(new G4BGGNucleonElasticXS(particle));
183      hel->AddDataSet(new G4CHIPSElasticXS());
184      hel->RegisterMe(chipsn);
185
186      pmanager->AddDiscreteProcess(hel);
187
188      if(verbose > 1) {
189        G4cout << "### HadronElasticPhysics: " << hel->GetProcessName()
190               << " added for " << particle->GetParticleName() << G4endl;
191      }
192    }
193  }
194}
195
196
Note: See TracBrowser for help on using the repository browser.