source: trunk/source/physics_lists/builders/src/G4HadronElasticPhysicsXS.cc @ 1316

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

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

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