source: trunk/source/processes/electromagnetic/highenergy/src/G4eeCrossSections.cc @ 836

Last change on this file since 836 was 819, checked in by garnier, 16 years ago

import all except CVS
File size: 9.0 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: G4eeCrossSections.cc,v 1.6 2006/06/29 19:32:42 gunter Exp $
27// GEANT4 tag $Name:  $
28//
29// -------------------------------------------------------------------
30//
31// GEANT4 Class header file
32//
33//
34// File name:     G4eeCrossSections
35//
36// Author:        Vladimir Ivanchenko
37//
38// Creation date: 25.10.2003
39//
40// Modifications:
41//
42//
43// -------------------------------------------------------------------
44//
45
46
47//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
48//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
49
50#include "G4eeCrossSections.hh"
51#include "G4PionPlus.hh"
52#include "G4PionMinus.hh"
53#include "G4PionZero.hh"
54#include "G4KaonPlus.hh"
55#include "G4KaonMinus.hh"
56#include "G4KaonZeroLong.hh"
57#include "G4PhysicsLinearVector.hh"
58
59#include <iostream>
60#include <fstream>
61
62//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
63
64using namespace std;
65
66G4eeCrossSections::G4eeCrossSections()
67{
68  Initialise();
69}
70
71//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
72
73G4eeCrossSections::~G4eeCrossSections()
74{}
75
76//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
77
78void G4eeCrossSections::Initialise()
79{
80  MsPi = G4PionPlus::PionPlus()->GetPDGMass();
81  MsPi0= G4PionZero::PionZero()->GetPDGMass();
82  MsEta= 547.30*MeV;
83  MsEtap=957.78*MeV;
84  MsKs = G4KaonZeroLong::KaonZeroLong()->GetPDGMass();
85  MsKc=G4KaonPlus::KaonPlus()->GetPDGMass();
86  MsRho= 770.0*MeV;
87  MsOm = 781.94*MeV;
88  MsF0 = 980.0*MeV;
89  MsA0 = 983.4*MeV;
90  MsPhi= 1019.413*MeV;
91  MsK892 = 891.66*MeV;
92  MsK0892 = 896.10*MeV;
93  GRho = 150.7*MeV;
94  GOm = 8.41*MeV;
95  GPhi = 4.43*MeV;
96  GK892 = 50.8*MeV;
97  GK0892 = 50.5*MeV;
98  PhRho = 0.0;
99  PhOm = 0.0;
100  PhPhi = 155.0*degree;
101  PhRhoPi = 186.0*degree;
102     
103  BrRhoPiG = 4.5e-4;
104  BrRhoPi0G= 6.8e-4;
105  BrRhoEtaG= 2.4e-4;
106  BrRhoEe = 4.49e-5;
107  BrOm3Pi = 0.888;
108  BrOmPi0G= 0.085;
109  BrOmEtaG= 6.5e-4;
110  BrOm2Pi = 0.0221;
111  PhOm2Pi = 90.0;
112  BrOmEe = 7.07e-5;
113  BrPhi2Kc = 0.491;
114  BrPhiKsKl= 0.341;
115  BrPhi3Pi = 0.155;
116  BrPhiPi0G= 1.31e-3;
117  BrPhiEtaG= 1.26e-2;
118  BrPhi2Pi = 8.e-5;
119  PhPhi2Pi = -20.0*degree;
120  BrPhiEe = 2.99e-4;
121
122  MsRho3 = MsRho*MsRho*MsRho;
123  MsOm3  = MsOm*MsOm*MsOm;
124  MsPhi3 = MsPhi*MsPhi*MsPhi;
125
126  MeVnb = 3.8938e+11*nanobarn;
127  Alpha = 1.0/137.036;
128
129  AOmRho = 3.0;
130  ARhoPRho = 0.72;
131  cterm=0.;
132  mssig = 600.*MeV;
133  gsig = 500.*MeV;
134  brsigpipi = 1.;
135
136  msrho1450 = 1465.*MeV;
137  msrho1700 = 1700.*MeV;
138  grho1450 = 310.*MeV;
139  grho1700 = 240.*MeV;
140  arhoompi0 = 1.;
141  arho1450ompi0 = 1.;
142  arho1700ompi0 = 1.;
143  phrhoompi0 = 0.;
144  phrho1450ompi0 = pi;
145  phrho1700ompi0 = 0.;
146  aomrhopi0 = 1.;
147  phomrhopi0 = 0.;
148  arhopi0pi0g = 0.;
149  aompi0pi0g = 0.;
150  phrhopi0pi0g = 0.;
151  phompi0pi0g = 0.;
152  brrho1450ompi0 = 0.02;
153  brrho1450pipi = 0.50;
154  brrho1700ompi0 = 1.0;
155  brrho1700pipi = 0.02;
156  aphirhopi0 = 1.;
157  phphirhopi0 = pi;
158  arhosigg = 0.;
159  phrhosigg = 0.;
160  aomsigg = 0.;
161  phomsigg = 0.;
162
163  G4String w0, w1, w2;
164  ph3p = 0;
165
166  /*
167  G4double emin, emax;
168  G4int nbins;
169  const G4String fname = "wrhopi.wid";
170  ifstream  fi(fname.c_str());
171  fi >> w0 >> nbins >> w1 >> emin >> w2 >> emax;
172  emin *= MeV;
173  emax *= MeV;
174  ph3p = new G4PhysicsLinearVector(emin,emax,nbins);
175  G4int nlines = nbins/5;
176  G4double s0, s1, s2, s3, s4;
177  for(G4int i=0; i<nlines; i++) {
178    fi >> s0 >> s1 >> s2 >> s3 >> s4;
179    ph3p->PutValue(5*i, s0);
180    ph3p->PutValue(5*i + 1, s1);
181    ph3p->PutValue(5*i + 2, s2);
182    ph3p->PutValue(5*i + 3, s3);
183    ph3p->PutValue(5*i + 4, s4);
184  }
185  fi.close();
186  */
187}
188
189//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
190
191G4double G4eeCrossSections::CrossSection2pi(G4double e) 
192{
193 
194  complex<G4double> xr(cos(PhRho),sin(PhRho));
195  complex<G4double> xo(cos(PhOm2Pi),sin(PhOm2Pi));
196  complex<G4double> xf(cos(PhPhi2Pi),sin(PhPhi2Pi));
197
198  G4double s = e*e;
199  complex<G4double> drho = DpRho(e);
200  complex<G4double> dom  = DpOm(e);
201  complex<G4double> dphi = DpPhi(e);
202
203  complex<G4double> amp = 
204      sqrt(Width2p(s,MsRho,GRho,1.0,MsPi)*MsRho3*BrRhoEe*GRho)*xr/drho
205    + sqrt(Width2p(s,MsOm,GOm,BrOm2Pi,MsPi)*MsOm3*BrOmEe*GOm)*xo/dom
206    + sqrt(Width2p(s,MsPhi,GPhi,BrPhi2Pi,MsPi)*MsPhi3*BrPhiEe*GPhi)*xf/dphi;
207
208  G4double cross = 12.0*pi*MeVnb*norm(amp)/(e*s);
209
210  return cross;
211}
212
213//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
214
215G4double G4eeCrossSections::Width2p(G4double s, G4double mres, 
216                                    G4double g, G4double br, G4double m) 
217{
218  G4double m2    = 4.0*m*m;
219  G4double s0    = mres*mres;
220  G4double f     = (s - m2)/(s0 - m2);
221  if(f < 0.0)  f = 0.0;
222  return g*br*sqrt(f)*f*s0/s;
223}
224
225//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
226
227G4double G4eeCrossSections::Width3p(G4double s, G4double mres, 
228                                    G4double g, G4double br) 
229{
230  G4double w = PhaseSpace3p(sqrt(s));
231  G4double w0= PhaseSpace3p(mres);
232  G4double x = g*br*w/w0;
233  return x;
234}
235
236//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
237
238G4double G4eeCrossSections::PhaseSpace3p(G4double e) 
239{
240 
241  //  G4bool b;
242  //  G4double x = ph3p->GetValue(e, b);
243  G4double x = 1.0; 
244  G4double emev = e/MeV;
245  G4double y = 414.12/emev;
246  x *= pow(e/MsOm, 5.0) * pow(emev*0.1, 3.0)*(1.0 - y*y);
247  return x;
248}
249
250//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
251
252G4double G4eeCrossSections::WidthPg(G4double s, G4double mres, 
253                                    G4double g, G4double br, G4double m)
254{
255  G4double m2    = m*m;
256  G4double s0    = mres*mres;
257  G4double f     = (s - m2)*mres/((s0 - m2)*sqrt(s));
258  if(f < 0.0)  f = 0.0;
259  return g*br*f*f*f;
260}
261
262//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
263
264G4double G4eeCrossSections::WidthRho(G4double e) 
265{
266  G4double w = Width2p(e*e, MsRho, GRho, 1.0, MsPi);
267  return w;
268}
269
270//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
271
272G4double G4eeCrossSections::WidthOm(G4double e) 
273{
274  G4double s = e*e;
275  G4double w = (Width3p(s, MsOm, GOm, BrOm3Pi) +
276                WidthPg(s, MsOm, GOm, BrOmPi0G, MsPi0) +
277                WidthPg(s, MsOm, GOm, BrOmEtaG, MsEta) +
278                Width2p(s, MsOm, GOm, BrOm2Pi, MsPi)) /
279       (BrOm3Pi+BrOmPi0G+BrOmEtaG+BrOm2Pi);
280  return w;
281}
282
283//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
284
285G4double G4eeCrossSections::WidthPhi(G4double e) 
286{
287  G4double s = e*e;
288  G4double w = (Width3p(s, MsPhi, GPhi, BrPhi3Pi) +
289                WidthPg(s, MsPhi, GPhi, BrPhiPi0G, MsPi0) +
290                WidthPg(s, MsPhi, GPhi, BrPhiEtaG, MsEta) +
291                Width2p(s, MsPhi, GPhi, BrPhi2Kc, MsKc) +
292                Width2p(s, MsPhi, GPhi, BrPhiKsKl, MsKs)) /
293        (BrPhi3Pi+BrPhiPi0G+BrPhiEtaG+BrPhi2Kc+BrPhiKsKl);
294  return w;
295}
296
297//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
298
299complex<G4double> G4eeCrossSections::DpRho(G4double e) 
300{
301  complex<G4double> d(MsRho*MsRho - e*e, -e*WidthRho(e));
302  return d;
303}
304
305//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
306
307complex<G4double> G4eeCrossSections::DpOm(G4double e) 
308{
309  complex<G4double> d(MsOm*MsOm - e*e, -e*WidthOm(e));
310  return d;
311}
312
313//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
314
315complex<G4double> G4eeCrossSections::DpPhi(G4double e) 
316{
317  complex<G4double> d(MsPhi*MsPhi - e*e, -e*WidthPhi(e));
318  return d;
319}
320
321//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
Note: See TracBrowser for help on using the repository browser.