source: trunk/source/processes/electromagnetic/highenergy/include/G4eeToHadronsModel.hh @ 846

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

import all except CVS
File size: 5.4 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: G4eeToHadronsModel.hh,v 1.6 2007/05/22 17:37:30 vnivanch Exp $
27// GEANT4 tag $Name:  $
28//
29// -------------------------------------------------------------------
30//
31// GEANT4 Class header file
32//
33//
34// File name:     G4eeToHadronsModel
35//
36// Author:        Vladimir Ivanchenko
37//
38// Creation date: 25.10.2003
39//
40// Modifications:
41// 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko)
42// 18-05-05 Use optimized interfaces (V.Ivantchenko)
43//
44
45//
46// Class Description:
47//
48
49// -------------------------------------------------------------------
50//
51
52#ifndef G4eeToHadronsModel_h
53#define G4eeToHadronsModel_h 1
54
55#include "G4VEmModel.hh"
56
57class G4PhysicsVector;
58class G4Vee2hadrons;
59
60class G4eeToHadronsModel : public G4VEmModel
61{
62
63public:
64
65  G4eeToHadronsModel(const G4Vee2hadrons*, G4int ver=0,
66                     const G4String& nam = "eeToHadrons");
67
68  virtual ~G4eeToHadronsModel();
69
70  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
71
72  G4double PeakEnergy() const;
73
74  virtual G4double CrossSectionPerVolume(const G4Material*,
75                                         const G4ParticleDefinition*,
76                                         G4double kineticEnergy,
77                                         G4double cutEnergy,
78                                         G4double maxEnergy);
79
80  virtual G4double ComputeCrossSectionPerAtom(
81                                         const G4ParticleDefinition*,
82                                         G4double kineticEnergy,
83                                         G4double Z, G4double A,
84                                         G4double cutEnergy = 0.0,
85                                         G4double maxEnergy = DBL_MAX);
86
87  virtual G4double ComputeCrossSectionPerElectron(
88                                         const G4ParticleDefinition*,
89                                         G4double kineticEnergy,
90                                         G4double cutEnergy = 0.0,
91                                         G4double maxEnergy = DBL_MAX);
92
93  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
94                                 const G4MaterialCutsCouple*,
95                                 const G4DynamicParticle*,
96                                 G4double tmin = 0.0,
97                                 G4double maxEnergy = DBL_MAX);
98
99  G4DynamicParticle* GenerateCMPhoton(G4double);
100
101private:
102
103  void ComputeCMCrossSectionPerElectron();
104
105  // hide assignment operator
106  G4eeToHadronsModel & operator=(const  G4eeToHadronsModel &right);
107  G4eeToHadronsModel(const  G4eeToHadronsModel&);
108
109  const G4Vee2hadrons*  model;
110  G4ParticleDefinition* theGamma;
111  G4PhysicsVector*      crossPerElectron;
112  G4PhysicsVector*      crossBornPerElectron;
113  G4bool                isInitialised;
114  G4int                 nbins;
115  G4int                 verbose;
116
117  G4double              lowKinEnergy;
118  G4double              peakKinEnergy;
119  G4double              highKinEnergy;
120
121  G4double              emin;
122  G4double              epeak;
123  G4double              emax;
124};
125
126//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
127
128inline G4double G4eeToHadronsModel::PeakEnergy() const
129{
130  return peakKinEnergy;
131}
132
133//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
134
135inline G4double G4eeToHadronsModel::CrossSectionPerVolume(
136                                      const G4Material* mat,
137                                      const G4ParticleDefinition* p,
138                                      G4double kineticEnergy,
139                                      G4double, G4double)
140{
141  return mat->GetElectronDensity()*
142    ComputeCrossSectionPerElectron(p, kineticEnergy);
143}
144
145//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
146
147inline G4double G4eeToHadronsModel::ComputeCrossSectionPerAtom(
148                                      const G4ParticleDefinition* p,
149                                      G4double kineticEnergy,
150                                      G4double Z, G4double,
151                                      G4double, G4double)
152{
153  return Z*ComputeCrossSectionPerElectron(p, kineticEnergy);
154}
155
156//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
157
158#endif
Note: See TracBrowser for help on using the repository browser.