Geant4  10.03
G4AdjointBremsstrahlungModel.hh
Go to the documentation of this file.
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: G4AdjointBremsstrahlungModel.hh 100666 2016-10-31 10:27:00Z gcosmo $
27 //
29 // Class: G4AdjointBremsstrahlungModel
30 // Author: L. Desorgher
31 // Organisation: SpaceIT GmbH
32 // Contract: ESA contract 21435/08/NL/AT
33 // Customer: ESA/ESTEC
35 //
36 // CHANGE HISTORY
37 // --------------
38 // ChangeHistory:
39 // 15 June 2007 creation by L. Desorgher. Adapted from G4eBremsstrahlungModel
40 // 20-10-2009 Remove all the screening effect that are not considered in the direct models blow 10 GeV. L.Desorgher
41 // 4-11-2009 Implement the use of a simple biased differential cross section (C(Z)/Egamma) allowing a rapid computation of adjoint CS
42 // and rapid sampling of adjoint secondaries. By this way cross section matrices are not used anymore, avoiding a rather
43 // time consuming computation of adjoint brem cross section matrices for each material at initialisation. This mode is switch on/off
44 // by selecting SetUseMatrix(false)/ SetUseMatrix(true) in the constructor. L.Desorgher
45 //
46 //
47 //-------------------------------------------------------------
48 // Documentation:
49 // Adjoint Model for e- Bremsstrahlung
50 //
51 
52 
53 
54 #ifndef G4AdjointBremsstrahlungModel_h
55 #define G4AdjointBremsstrahlungModel_h 1
56 #include "globals.hh"
57 #include "G4VEmAdjointModel.hh"
59 #include "G4PhysicsTable.hh"
60 #include "G4EmModelManager.hh"
61 class G4Timer;
63 
64 {
65 public:
66 
70  virtual void SampleSecondaries(const G4Track& aTrack,
71  G4bool IsScatProjToProjCase,
72  G4ParticleChange* fParticleChange);
73  void RapidSampleSecondaries(const G4Track& aTrack,
74  G4bool IsScatProjToProjCase,
75  G4ParticleChange* fParticleChange);
77  const G4Material* aMaterial,
78  G4double kinEnergyProj, // kinetic energy of the primary particle before the interaction
79  G4double kinEnergyProd // kinetic energy of the secondary particle
80  );
82  const G4Material* aMaterial,
83  G4double kinEnergyProj, // kinetic energy of the primary particle before the interaction
84  G4double kinEnergyProd // kinetic energy of the secondary particle
85  );
87  const G4Material* aMaterial,
88  G4double kinEnergyProj, // kinetic energy of the primary particle before the interaction
89  G4double kinEnergyProd // kinetic energy of the secondary particle
90  );
91  virtual G4double AdjointCrossSection(const G4MaterialCutsCouple* aCouple,
92  G4double primEnergy,
93  G4bool IsScatProjToProjCase);
95  G4double primEnergy,
96  G4bool IsScatProjToProjCase);
97 
98 
99  // private void InitialiseFwdModels();
100 
101 
102 private:
106 
109  std::vector<G4DataVector*> partialSumSigma;
110  std::vector<float> SigmaPerAtom;
111 
112 
113 };
114 
115 
116 #endif
std::vector< G4DataVector * > partialSumSigma
virtual void SampleSecondaries(const G4Track &aTrack, G4bool IsScatProjToProjCase, G4ParticleChange *fParticleChange)
virtual G4double AdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
virtual G4double GetAdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
virtual G4double DiffCrossSectionPerVolumePrimToSecond(const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyProd)
bool G4bool
Definition: G4Types.hh:79
void RapidSampleSecondaries(const G4Track &aTrack, G4bool IsScatProjToProjCase, G4ParticleChange *fParticleChange)
G4double DiffCrossSectionPerVolumePrimToSecondApproximated2(const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyProd)
G4double DiffCrossSectionPerVolumePrimToSecondApproximated1(const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyProd)
double G4double
Definition: G4Types.hh:76