Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/muons/include/G4MuPairProduction.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/G4MuPairProduction.hh (Version 11.3.0) and /processes/electromagnetic/muons/include/G4MuPairProduction.hh (Version 5.2.p2)


  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: G4MuPairProduction.hh,v 1.14 2003/06/16 17:01:42 gunter Exp $
                                                   >>  25 // GEANT4 tag $Name: geant4-05-02-patch-01 $
 28 //                                                 26 //
 29 // GEANT4 Class header file                    <<  27 //--------------- G4MuPairProduction physics process ------------------
 30 //                                             <<  28 //                by Laszlo Urban, May 1998      
 31 //                                             <<  29 //------------------------------------------------------------------------------
 32 // File name:     G4MuPairProduction           << 
 33 //                                             << 
 34 // Author:        Laszlo Urban                 << 
 35 //                                             << 
 36 // Creation date: 30.05.1998                   << 
 37 //                                             << 
 38 // Modifications:                              << 
 39                                                << 
 40 // 10-02-00 modifications , new e.m. structure     30 // 10-02-00 modifications , new e.m. structure, L.Urban
 41 // 10-08-01 new methods Store/Retrieve Physics     31 // 10-08-01 new methods Store/Retrieve PhysicsTable (mma)
 42 // 29-10-01 all static functions no more inlin <<  32 // 29-10-01 all static functions no more inlined (mma) 
 43 // 10-05-02 V.Ivanchenko update to new design  <<  33 // 16-01-03 Migrade to cut per region (V.Ivanchenko)
 44 // 26-12-02 secondary production moved to deri <<  34 //------------------------------------------------------------------------------
 45 // 27-01-03 Make models region aware (V.Ivanch << 
 46 // 05-02-03 Fix compilation warnings (V.Ivanch << 
 47 // 13-02-03 SubCutoff regime is assigned to a  << 
 48 // 08-08-03 STD substitute standard  (V.Ivanch << 
 49 // 12-11-03 G4EnergyLossSTD -> G4EnergyLossPro << 
 50 // 21-01-04 Migrade to G4ParticleChangeForLoss << 
 51 // 28-04-04 Fix minor bug in energy balance (V << 
 52 // 17-08-04 Rename the process "Mu" -> "mu" (V << 
 53 // 08-11-04 Migration to new interface of Stor << 
 54 // 08-04-05 Major optimisation of internal int << 
 55 //                                             << 
 56 // Class Description:                          << 
 57 //                                             << 
 58 // This class manages the PairProduction proce << 
 59 // it inherites from G4VContinuousDiscreteProc << 
 60 //                                             << 
 61                                                    35 
 62 // ------------------------------------------- <<  36 #ifndef G4MuPairproduction_h
 63 //                                             <<  37 #define G4MuPairproduction_h 1
 64                                                << 
 65 #ifndef G4MuPairProduction_h                   << 
 66 #define G4MuPairProduction_h 1                 << 
 67                                                    38 
 68 //....oooOO0OOooo........oooOO0OOooo........oo     39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 69                                                    40 
                                                   >>  41 #include "G4ios.hh"
 70 #include "globals.hh"                              42 #include "globals.hh"
 71 #include "G4VEnergyLossProcess.hh"             <<  43 #include "Randomize.hh"
                                                   >>  44 #include "G4VMuEnergyLoss.hh"
                                                   >>  45 #include "G4Track.hh"
                                                   >>  46 #include "G4Step.hh"
                                                   >>  47 #include "G4OrderedTable.hh"
                                                   >>  48 #include "G4PhysicsTable.hh"
                                                   >>  49 #include "G4PhysicsLogVector.hh"
                                                   >>  50 #include "G4MaterialCutsCouple.hh"
                                                   >>  51 
                                                   >>  52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >>  53 
                                                   >>  54 class G4MuPairProduction : public G4VMuEnergyLoss
 72                                                    55 
 73 class G4MuPairProduction : public G4VEnergyLos << 
 74 {                                                  56 {
 75 public:                                        <<  57   public:
 76                                                    58 
 77   explicit G4MuPairProduction(const G4String&  <<  59      G4MuPairProduction(const G4String& processName = "MuPairProd");
 78                                                    60 
 79   ~G4MuPairProduction() override = default;    <<  61     ~G4MuPairProduction();
 80                                                    62 
 81   G4bool IsApplicable(const G4ParticleDefiniti <<  63      G4bool IsApplicable(const G4ParticleDefinition&);
 82                                                    64 
 83   G4double MinPrimaryEnergy(const G4ParticleDe <<  65      void BuildPhysicsTable(const G4ParticleDefinition& ParticleType);
 84           const G4Material*, G4double cut) ove << 
 85                                                    66 
 86   inline void SetLowestKineticEnergy(G4double  <<  67      void BuildLossTable(const G4ParticleDefinition& ParticleType);
 87                                                    68 
 88   // print description in html                 <<  69      void BuildLambdaTable(const G4ParticleDefinition& ParticleType);
 89   void ProcessDescription(std::ostream&) const << 
 90                                                    70 
 91   G4MuPairProduction & operator=(const G4MuPai <<  71      void PrintInfoDefinition();
 92   G4MuPairProduction(const G4MuPairProduction& << 
 93                                                    72 
 94 protected:                                     <<  73      G4double GetMeanFreePath( const G4Track& track,
                                                   >>  74                                G4double previousStepSize,
                                                   >>  75                                G4ForceCondition* condition );
 95                                                    76 
 96   // Print out of the class parameters         <<  77      G4VParticleChange *PostStepDoIt(const G4Track& track,
 97   void StreamProcessInfo(std::ostream& outFile <<  78                                      const G4Step& Step  );
 98                                                    79 
 99   void InitialiseEnergyLossProcess(const G4Par <<  80      G4double GetDMicroscopicCrossSection(
100                  const G4ParticleDefinition*)  <<  81                                       const G4ParticleDefinition* ParticleType,
                                                   >>  82                                       G4double KineticEnergy,
                                                   >>  83                                       G4double AtomicNumber,
                                                   >>  84                                       G4double PairEnergy);
101                                                    85 
102   const G4ParticleDefinition* theParticle = nu <<  86      G4bool StorePhysicsTable(G4ParticleDefinition* ,
103   G4double                    lowestKinEnergy; <<  87                     const G4String& directory, G4bool);
104   G4bool                      isInitialised =  <<  88       // store eLoss and MeanFreePath tables into an external file
                                                   >>  89       // specified by 'directory' (must exist before invokation)
105                                                    90 
106 };                                             <<  91      G4bool RetrievePhysicsTable(G4ParticleDefinition* ,
                                                   >>  92                const G4String& directory, G4bool);
                                                   >>  93       // retrieve eLoss and MeanFreePath tables from an external file
                                                   >>  94       // specified by 'directory'
107                                                    95 
108 //....oooOO0OOooo........oooOO0OOooo........oo <<  96   protected:
109                                                    97 
110 inline void G4MuPairProduction::SetLowestKinet <<  98      G4double ComputeMeanFreePath(const G4ParticleDefinition* ParticleType,
111 {                                              <<  99                                         G4double KineticEnergy,
112   lowestKinEnergy = e;                         << 100                                   const G4MaterialCutsCouple* couple);
113 }                                              << 101 
                                                   >> 102      void ComputePartialSumSigma(const G4ParticleDefinition* ParticleType,
                                                   >> 103                                        G4double KineticEnergy,
                                                   >> 104                                  const G4MaterialCutsCouple* couple);
                                                   >> 105 
                                                   >> 106      virtual G4double ComputeMicroscopicCrossSection(
                                                   >> 107                                       const G4ParticleDefinition* ParticleType,
                                                   >> 108                                       G4double KineticEnergy,
                                                   >> 109                                       G4double AtomicNumber,
                                                   >> 110                                       G4double ElectronEnergyCut,
                                                   >> 111                                       G4double PositronEnergyCut);
                                                   >> 112 
                                                   >> 113      G4double ComputeDDMicroscopicCrossSection(
                                                   >> 114                                       const G4ParticleDefinition* ParticleType,
                                                   >> 115                                       G4double KineticEnergy,
                                                   >> 116                                       G4double AtomicNumber,
                                                   >> 117                                       G4double PairEnergy,
                                                   >> 118                                       G4double asymmetry);
                                                   >> 119 
                                                   >> 120      G4double ComputeDMicroscopicCrossSection(
                                                   >> 121                                       const G4ParticleDefinition* ParticleType,
                                                   >> 122                                       G4double KineticEnergy,
                                                   >> 123                                       G4double AtomicNumber,
                                                   >> 124                                       G4double PairEnergy);
                                                   >> 125 
                                                   >> 126   private:
                                                   >> 127 
                                                   >> 128      G4MuPairProduction & operator=(const G4MuPairProduction &right);
                                                   >> 129      G4MuPairProduction(const G4MuPairProduction&);
                                                   >> 130 
                                                   >> 131      G4double ComputePairLoss( const G4ParticleDefinition* ParticleType,
                                                   >> 132                                G4double Z,G4double T,G4double ElectronCut,
                                                   >> 133                                      G4double PositronCut);
                                                   >> 134 
                                                   >> 135      G4Element* SelectRandomAtom(const G4MaterialCutsCouple* couple) const;
                                                   >> 136 
                                                   >> 137      void MakeSamplingTables( const G4ParticleDefinition* ParticleType );
                                                   >> 138 
                                                   >> 139   protected:
                                                   >> 140 
                                                   >> 141      virtual G4double SecondaryEnergyThreshold(size_t index);
                                                   >> 142   
                                                   >> 143   private:
                                                   >> 144 
                                                   >> 145      G4PhysicsTable* theMeanFreePathTable;
                                                   >> 146 
                                                   >> 147      G4OrderedTable PartialSumSigma;
                                                   >> 148 
                                                   >> 149      static G4double LowerBoundLambda; // bining for lambda table
                                                   >> 150      static G4double UpperBoundLambda;
                                                   >> 151      static G4int    NbinLambda;
                                                   >> 152      G4double LowestKineticEnergy,HighestKineticEnergy;
                                                   >> 153      G4int    TotBin;
                                                   >> 154 
                                                   >> 155      const std::vector<G4double>* electronEnergyCuts;
                                                   >> 156      const std::vector<G4double>* positronEnergyCuts;
                                                   >> 157 
                                                   >> 158      // tables for sampling
                                                   >> 159      static G4int nzdat,ntdat,NBIN;
                                                   >> 160      static G4double zdat[5],tdat[8];
                                                   >> 161      static G4double ya[1001],proba[5][8][1001];
                                                   >> 162      static G4double MinPairEnergy;
                                                   >> 163 
                                                   >> 164   public:
                                                   >> 165 
                                                   >> 166     static void SetLowerBoundLambda(G4double val);
                                                   >> 167     static void SetUpperBoundLambda(G4double val);
                                                   >> 168     static void SetNbinLambda(G4int n);
                                                   >> 169     static G4double GetLowerBoundLambda();
                                                   >> 170     static G4double GetUpperBoundLambda();
                                                   >> 171     static G4int GetNbinLambda();
                                                   >> 172 
                                                   >> 173 };
114                                                   174 
115 //....oooOO0OOooo........oooOO0OOooo........oo    175 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 176 
                                                   >> 177 #include "G4MuPairProduction.icc"
116                                                   178 
117 #endif                                            179 #endif
118                                                   180