ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4ContinuousGainOfEnergy.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4ContinuousGainOfEnergy.hh
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 //
28 // Class: G4ContinuousGainOfEnergy
29 // Author: L. Desorgher
30 // Organisation: SpaceIT GmbH
31 // Contract: ESA contract 21435/08/NL/AT
32 // Customer: ESA/ESTEC
34 //
35 // CHANGE HISTORY
36 // --------------
37 // ChangeHistory:
38 // -10 May 2007 creation by L. Desorgher
39 // -February-March 2009 Update for protons by L.Desorgher
40 // -July August 2009 Update for ion by L.Desorgher
41 //
42 //-------------------------------------------------------------
43 // Documentation:
44 // Continuous process acting on adjoint particles to compute the continuous gain of energy of charged particles when they are tracked back!
45 //
46 //
47 #ifndef G4ContinuousGainOfEnergy_h
48 #define G4ContinuousGainOfEnergy_h 1
49 
50 #include "G4VContinuousProcess.hh"
51 #include "globals.hh"
52 #include "G4Material.hh"
53 #include "G4MaterialCutsCouple.hh"
54 #include "G4Track.hh"
55 #include "G4UnitsTable.hh"
56 #include "G4ParticleChange.hh"
57 #include "G4VEnergyLossProcess.hh"
58 #include "G4ProductionCutsTable.hh"
59 
60 
61 class G4Step;
63 class G4VEmModel;
65 
66 
67 
69 {
70 public:
71 
72  G4ContinuousGainOfEnergy(const G4String& name = "EnergyGain",
74 
75  virtual ~G4ContinuousGainOfEnergy();
76 
77 
78 protected:
79 
80 
81  //------------------------------------------------------------------------
82  // Methods with standard implementation; may be overwritten if needed
83  //------------------------------------------------------------------------
84 protected:
85 
86 
88  G4double previousStepSize,
89  G4double currentMinimumStep,
90  G4double& currentSafety);
91 
92 
93  //------------------------------------------------------------------------
94  // Generic methods common to all processes
95  //------------------------------------------------------------------------
96 public:
97 
98 
99 
101 
103 
104 
106 
107 
108  void SetLossFluctuations(G4bool val);
109  inline void SetIsIntegral(G4bool val){is_integral= val;}
110 
112 
114 
115 protected:
116 
117 
118 
119 
120 private:
121 
122  void DefineMaterial(const G4MaterialCutsCouple* couple);
123  void SetDynamicMassCharge(const G4Track& track, G4double energy);
124 
125 
126  // hide assignment operator
127 
130 
131 
132 private:
133 
141 
142 
143 
147 
150 
151 
153 
154  //adding for Ions
155  //----------------
163 
164 
165 
166 
167 
168 };
169 
171 //
173  const G4MaterialCutsCouple* couple)
174 {
175  if(couple != currentCouple) {
176  currentCouple = couple;
177  currentMaterial = couple->GetMaterial();
178  currentCoupleIndex = couple->GetIndex();
180 
181  size_t idx=1;
182  const std::vector<G4double>* aVec = G4ProductionCutsTable::GetProductionCutsTable()->GetEnergyCutsVector(idx);
185  //G4cout<<"Define Material"<<G4endl;
186  //if(!meanFreePath) ResetNumberOfInteractionLengthLeft();
187  }
188 }
189 
190 #endif