Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/muons/include/G4MuIonisation.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/electromagnetic/muons/include/G4MuIonisation.hh (Version 11.3.0) and /processes/electromagnetic/muons/include/G4MuIonisation.hh (Version 5.2.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
 26 //                                                 23 //
 27 // ------------------------------------------- <<  24 // $Id: G4MuIonisation.hh,v 1.15 2003/06/16 17:01:40 gunter Exp $
                                                   >>  25 // GEANT4 tag $Name: geant4-05-02-patch-01 $
 28 //                                                 26 //
 29 // GEANT4 Class header file                    <<  27 // --------------- G4MuIonisation physics process ------------------------------
                                                   >>  28 //                 by Laszlo Urban, September 1997 
                                                   >>  29 // -----------------------------------------------------------------------------
 30 //                                                 30 //
 31 //                                             <<  31 // 10-02-00 modifications , new e.m. structure, L.Urban
 32 // File name:     G4MuIonisation               << 
 33 //                                             << 
 34 // Author:        Laszlo Urban                 << 
 35 //                                             << 
 36 // Creation date: 30.05.1997                   << 
 37 //                                             << 
 38 // Modifications:                              << 
 39 //                                             << 
 40 // corrected by L.Urban on 24/09/97            << 
 41 // corrected by L.Urban on 13/01/98            << 
 42 // bugs fixed by L.Urban on 02/02/99           << 
 43 // 10/02/00 modifications , new e.m. structure << 
 44 // 10-08-01 new methods Store/Retrieve Physics     32 // 10-08-01 new methods Store/Retrieve PhysicsTable (mma)
 45 // 14-08-01 new function ComputeRestrictedMean <<  33 // 29-08-01 new function ComputeRestrictedMeandEdx() + 'cleanup' (mma)
 46 // 19-09-01 come back to previous process name <<  34 // 19-09-01 come back to the old process name 'MuIoni'
 47 // 29-10-01 all static functions no more inlin <<  35 // 29-10-01 all static functions no more inlined (mma)
 48 // 10-05-02 V.Ivanchenko update to new design  <<  36 // 16-01-03 Migrade to cut per region (V.Ivanchenko)
 49 // 09-12-02 V.Ivanchenko remove warning        << 
 50 // 26-12-02 Secondary production moved to deri << 
 51 // 24-01-03 Make models region aware (V.Ivanch << 
 52 // 05-02-03 Fix compilation warnings (V.Ivanch << 
 53 // 13-02-03 SubCutoff regime is assigned to a  << 
 54 // 23-05-03 Add fluctuation model as a member  << 
 55 // 03-06-03 Add SetIntegral method to choose f << 
 56 // 03-06-03 Fix initialisation problem for STD << 
 57 // 08-08-03 STD substitute standard  (V.Ivanch << 
 58 // 12-11-03 G4EnergyLossSTD -> G4EnergyLossPro << 
 59 // 21-01-04 Migrade to G4ParticleChangeForLoss << 
 60 // 17-08-04 Rename the process "Mu" -> "mu" (V << 
 61 // 08-11-04 Migration to new interface of Stor << 
 62 // 08-04-05 Major optimisation of internal int << 
 63 //                                                 37 //
 64 // Class Description:                          <<  38 // -----------------------------------------------------------------------------
                                                   >>  39 
                                                   >>  40 // Class description
 65 //                                                 41 //
 66 // This class manages the ionisation process f     42 // This class manages the ionisation process for muons.
 67 // it inherites from G4VContinuousDiscreteProc <<  43 // it inherites from G4VContinuousDiscreteProcess via G4VMuEnergyLoss.
 68 //                                                 44 //
                                                   >>  45 // Class description - end
 69                                                    46 
 70 // ------------------------------------------- <<  47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 71 //                                             <<  48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 72                                                    49 
 73 #ifndef G4MuIonisation_h                           50 #ifndef G4MuIonisation_h
 74 #define G4MuIonisation_h 1                         51 #define G4MuIonisation_h 1
 75                                                    52 
 76 #include "G4VEnergyLossProcess.hh"             <<  53 #include "G4VMuEnergyLoss.hh"
 77 #include "globals.hh"                          <<  54 #include "G4MaterialCutsCouple.hh"
 78 #include "G4VEmModel.hh"                       << 
 79                                                    55 
 80 class G4Material;                              <<  56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 81                                                    57 
 82 class G4MuIonisation : public G4VEnergyLossPro <<  58 class G4MuIonisation : public G4VMuEnergyLoss
 83 {                                              << 
 84                                                    59 
 85 public:                                        <<  60 {
                                                   >>  61   public:   // with description
 86                                                    62 
 87   explicit G4MuIonisation(const G4String& name <<  63     G4MuIonisation(const G4String& processName = "MuIoni");
 88                                                    64 
 89   ~G4MuIonisation() override = default;        <<  65    ~G4MuIonisation();
 90                                                    66 
 91   G4bool IsApplicable(const G4ParticleDefiniti <<  67     G4bool IsApplicable(const G4ParticleDefinition&);
                                                   >>  68       // return true for charged particles, false otherwise
 92                                                    69 
 93   G4double MinPrimaryEnergy(const G4ParticleDe <<  70     void BuildPhysicsTable(const G4ParticleDefinition& aParticleType);
 94                             const G4Material*, <<  71       // this function overloads a virtual function of the base class.
                                                   >>  72       // It is invoked by the G4ParticleWithCuts::SetCut() method.
                                                   >>  73       // It invokes BuildLambdaTable(), BuildLossTable(), BuildDEDXTable()
                                                   >>  74 
                                                   >>  75     void BuildLossTable(const G4ParticleDefinition& aParticleType);
                                                   >>  76       // build the dE/dx tables due to the ionisation, for every materials.
                                                   >>  77       // (restricted stopping power, Bethe-Bloch formula)
                                                   >>  78 
                                                   >>  79     void BuildLambdaTable(const G4ParticleDefinition& aParticleType);
                                                   >>  80       // build mean free path tables for the delta rays production.
                                                   >>  81       // the tables are built for every materials.
                                                   >>  82 
                                                   >>  83     G4bool StorePhysicsTable(G4ParticleDefinition* ,
                                                   >>  84                  const G4String& directory, G4bool);
                                                   >>  85       // store eLoss and MeanFreePath tables into an external file
                                                   >>  86       // specified by 'directory' (must exist before invokation)
                                                   >>  87 
                                                   >>  88     G4bool RetrievePhysicsTable(G4ParticleDefinition* ,
                                                   >>  89               const G4String& directory, G4bool);
                                                   >>  90       // retrieve eLoss and MeanFreePath tables from an external file
                                                   >>  91       // specified by 'directory'
                                                   >>  92 
                                                   >>  93     virtual void PrintInfoDefinition();
                                                   >>  94       // Print few lines of informations about the process: validity range,
                                                   >>  95       // origine ..etc..
                                                   >>  96       // Invoked by BuildPhysicsTable().
                                                   >>  97 
                                                   >>  98     G4double GetMeanFreePath(const G4Track& track,
                                                   >>  99                              G4double previousStepSize,
                                                   >> 100                              G4ForceCondition* condition );
                                                   >> 101       // It returns the MeanFreePath of the process for the current track :
                                                   >> 102       // (energy, material)
                                                   >> 103       // The previousStepSize and G4ForceCondition* are not used.
                                                   >> 104       // This function overloads a virtual function of the base class.
                                                   >> 105       // It is invoked by the ProcessManager of the Particle.
                                                   >> 106 
                                                   >> 107     G4VParticleChange* PostStepDoIt(const G4Track& track,
                                                   >> 108                                           const G4Step& Step  );
                                                   >> 109       // It computes the final state of the process (at end of step),
                                                   >> 110       // returned as a ParticleChange object.
                                                   >> 111       // This function overloads a virtual function of the base class.
                                                   >> 112       // It is invoked by the ProcessManager of the Particle.
                                                   >> 113 
                                                   >> 114 
                                                   >> 115   protected:   // with description
                                                   >> 116 
                                                   >> 117     virtual G4double ComputeRestrictedMeandEdx(
                                                   >> 118                             const G4ParticleDefinition& aParticleType,
                                                   >> 119                             G4double KineticEnergy,
                                                   >> 120                             const G4Material* material,
                                                   >> 121                             G4double DeltaThreshold);
                                                   >> 122       // computes restricted mean dE/dx in Geant4 internal units.
                                                   >> 123 
                                                   >> 124     virtual G4double ComputeCrossSectionPerAtom(
                                                   >> 125                             const G4ParticleDefinition& aParticleType,
                                                   >> 126                             G4double KineticEnergy,
                                                   >> 127                             G4double AtomicNumber,
                                                   >> 128           G4double DeltaThreshold);
                                                   >> 129       // computes total cross section per atom in Geant4 internal units.
                                                   >> 130 
                                                   >> 131     virtual G4double ComputeDifCrossSectionPerAtom(
                                                   >> 132                             const G4ParticleDefinition& aParticleType,
                                                   >> 133                             G4double KineticEnergy,
                                                   >> 134                             G4double AtomicNumber,
                                                   >> 135           G4double KnockonEnergy);
                                                   >> 136       // computes differential cross section per atom.
                                                   >> 137 
                                                   >> 138   protected:
                                                   >> 139 
                                                   >> 140     virtual G4double SecondaryEnergyThreshold(size_t index);
                                                   >> 141 
                                                   >> 142     G4PhysicsTable* theMeanFreePathTable;
                                                   >> 143 
                                                   >> 144   private:
                                                   >> 145 
                                                   >> 146     // hide assignment operator
                                                   >> 147     G4MuIonisation & operator=(const G4MuIonisation &right);
                                                   >> 148     G4MuIonisation(const G4MuIonisation&);
                                                   >> 149 
                                                   >> 150   private:
                                                   >> 151 
                                                   >> 152     static G4double LowerBoundLambda;      // bining for lambda table
                                                   >> 153     static G4double UpperBoundLambda;
                                                   >> 154     static G4int    NbinLambda;
                                                   >> 155 
                                                   >> 156     G4double LowestKineticEnergy;         // binning for dE/dx table
                                                   >> 157     G4double HighestKineticEnergy;
                                                   >> 158     G4int    TotBin;
                                                   >> 159 
                                                   >> 160     const std::vector<G4double>* secondaryEnergyCuts;
                                                   >> 161 
                                                   >> 162   public:  // with description
                                                   >> 163 
                                                   >> 164     static void SetLowerBoundLambda(G4double val);
                                                   >> 165     static void SetUpperBoundLambda(G4double val);
                                                   >> 166     static void SetNbinLambda(G4int n);
                                                   >> 167         // set the parameters of the mean free path table.
                                                   >> 168 
                                                   >> 169     static G4double GetLowerBoundLambda();
                                                   >> 170     static G4double GetUpperBoundLambda();
                                                   >> 171     static G4int GetNbinLambda();
                                                   >> 172       // get the parameters of the mean free path table.
                                                   >> 173 };
 95                                                   174 
                                                   >> 175 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 96                                                   176 
 97   // print description in html                 << 177 #include "G4MuIonisation.icc"
 98   void ProcessDescription(std::ostream&) const << 
 99                                                   178 
100   // hide assignment operator                  << 179 #endif
101   G4MuIonisation & operator=(const G4MuIonisat << 
102   G4MuIonisation(const G4MuIonisation&) = dele << 
103                                                   180 
104 protected:                                     << 
105                                                   181 
106   void InitialiseEnergyLossProcess(const G4Par << 
107                                    const G4Par << 
108                                                   182 
109 private:                                       << 
110                                                   183 
111   const G4ParticleDefinition* theParticle = nu << 
112   const G4ParticleDefinition* theBaseParticle  << 
113                                                   184 
114   G4double mass = 0.0;                         << 
115   G4double ratio = 0.0;                        << 
116                                                   185 
117   G4bool isInitialised = false;                << 
118 };                                             << 
119                                                   186 
120 //....oooOO0OOooo........oooOO0OOooo........oo << 
121                                                   187 
122 #endif                                         << 
123                                                   188