| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 ////////////////////////////////////////////// 26 ////////////////////////////////////////////////////////////////////////////////
27 // Class: G4AdjointhIonisationModel 27 // Class: G4AdjointhIonisationModel
28 // Author: L. Desorgher 28 // Author: L. Desorgher
29 // Organisation: SpaceIT GmbH 29 // Organisation: SpaceIT GmbH
30 // 30 //
31 // Adjoint EM model for discrete reverse hadr 31 // Adjoint EM model for discrete reverse hadron ionisation.
32 // Tested only for protons. 32 // Tested only for protons.
33 ////////////////////////////////////////////// 33 ////////////////////////////////////////////////////////////////////////////////
34 34
35 #ifndef G4AdjointhIonisationModel_h 35 #ifndef G4AdjointhIonisationModel_h
36 #define G4AdjointhIonisationModel_h 1 36 #define G4AdjointhIonisationModel_h 1
37 37
38 #include "globals.hh" 38 #include "globals.hh"
39 #include "G4VEmAdjointModel.hh" 39 #include "G4VEmAdjointModel.hh"
40 40
41 class G4MaterialCutsCouple; 41 class G4MaterialCutsCouple;
42 class G4ParticleChange; 42 class G4ParticleChange;
43 class G4ParticleDefinition; 43 class G4ParticleDefinition;
44 class G4Track; 44 class G4Track;
45 class G4VEmModel; 45 class G4VEmModel;
46 46
47 class G4AdjointhIonisationModel : public G4VEm 47 class G4AdjointhIonisationModel : public G4VEmAdjointModel
48 { 48 {
49 public: 49 public:
50 explicit G4AdjointhIonisationModel(G4Particl 50 explicit G4AdjointhIonisationModel(G4ParticleDefinition* pDef);
51 51
52 ~G4AdjointhIonisationModel() override; 52 ~G4AdjointhIonisationModel() override;
53 53
54 void SampleSecondaries(const G4Track& aTrack 54 void SampleSecondaries(const G4Track& aTrack, G4bool isScatProjToProj,
55 G4ParticleChange* fPa 55 G4ParticleChange* fParticleChange) override;
56 56
57 void RapidSampleSecondaries(const G4Track& a 57 void RapidSampleSecondaries(const G4Track& aTrack, G4bool isScatProjToProj,
58 G4ParticleChange 58 G4ParticleChange* fParticleChange);
59 59
60 G4double DiffCrossSectionPerAtomPrimToSecond 60 G4double DiffCrossSectionPerAtomPrimToSecond(
61 G4double kinEnergyProj, // kin energy of 61 G4double kinEnergyProj, // kin energy of primary before interaction
62 G4double kinEnergyProd, // kinetic energy 62 G4double kinEnergyProd, // kinetic energy of the secondary particle
63 G4double Z, G4double A = 0.) override; 63 G4double Z, G4double A = 0.) override;
64 64
65 G4double AdjointCrossSection(const G4Materia 65 G4double AdjointCrossSection(const G4MaterialCutsCouple* aCouple,
66 G4double primEn 66 G4double primEnergy,
67 G4bool isScatPr 67 G4bool isScatProjToProj) override;
68 68
69 G4double GetSecondAdjEnergyMaxForScatProjToP 69 G4double GetSecondAdjEnergyMaxForScatProjToProj(
70 G4double primAdjEnergy) override; 70 G4double primAdjEnergy) override;
71 71
72 G4double GetSecondAdjEnergyMinForScatProjToP 72 G4double GetSecondAdjEnergyMinForScatProjToProj(G4double primAdjEnergy,
73 73 G4double tcut = 0.) override;
74 74
75 G4double GetSecondAdjEnergyMaxForProdToProj( 75 G4double GetSecondAdjEnergyMaxForProdToProj(G4double primAdjEnergy) override;
76 76
77 G4double GetSecondAdjEnergyMinForProdToProj( 77 G4double GetSecondAdjEnergyMinForProdToProj(G4double primAdjEnergy) override;
78 78
79 G4AdjointhIonisationModel(G4AdjointhIonisati 79 G4AdjointhIonisationModel(G4AdjointhIonisationModel&) = delete;
80 G4AdjointhIonisationModel& operator=(const G 80 G4AdjointhIonisationModel& operator=(const G4AdjointhIonisationModel& right) =
81 delete; 81 delete;
82 82
83 private: 83 private:
84 void DefineProjectileProperty(); 84 void DefineProjectileProperty();
85 85
86 G4VEmModel* fBraggDirectEMModel; 86 G4VEmModel* fBraggDirectEMModel;
87 87
88 // projectile properties 88 // projectile properties
89 G4double fMass = 0.; 89 G4double fMass = 0.;
90 G4double fSpin = 0.; 90 G4double fSpin = 0.;
91 G4double fMagMoment2 = 0.; 91 G4double fMagMoment2 = 0.;
92 G4double fMassRatio = 0.; 92 G4double fMassRatio = 0.;
93 G4double fFormFact = 0.; 93 G4double fFormFact = 0.;
94 G4double fOnePlusRatio2 = 0.; 94 G4double fOnePlusRatio2 = 0.;
95 G4double fOneMinusRatio2 = 0.; 95 G4double fOneMinusRatio2 = 0.;
96 }; 96 };
97 97
98 #endif 98 #endif
99 99