Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/optical/src/G4OpAbsorption.cc

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/optical/src/G4OpAbsorption.cc (Version 11.3.0) and /processes/optical/src/G4OpAbsorption.cc (Version 9.3.p2)


  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 // $Id: G4OpAbsorption.cc,v 1.9 2008/10/24 19:52:28 gum Exp $
                                                   >>  28 // GEANT4 tag $Name: geant4-09-03-patch-02 $
 27 //                                                 29 //
 28 //////////////////////////////////////////////     30 ////////////////////////////////////////////////////////////////////////
 29 // Optical Photon Absorption Class Implementat     31 // Optical Photon Absorption Class Implementation
 30 //////////////////////////////////////////////     32 ////////////////////////////////////////////////////////////////////////
 31 //                                                 33 //
 32 // File:        G4OpAbsorption.cc                  34 // File:        G4OpAbsorption.cc
 33 // Description: Discrete Process -- Absorption <<  35 // Description: Discrete Process -- Absorption of Optical Photons  
 34 // Version:     1.0                                36 // Version:     1.0
 35 // Created:     1996-05-21                         37 // Created:     1996-05-21
 36 // Author:      Juliet Armstrong                   38 // Author:      Juliet Armstrong
 37 // Updated:     2005-07-28 - add G4ProcessType     39 // Updated:     2005-07-28 - add G4ProcessType to constructor
 38 //              2000-09-18 by Peter Gumplinger     40 //              2000-09-18 by Peter Gumplinger
 39 //              > comment out warning - "No Ab <<  41 //              > comment out warning - "No Absorption length specified" 
 40 //              1997-04-09 by Peter Gumplinger     42 //              1997-04-09 by Peter Gumplinger
 41 //              > new physics/tracking scheme      43 //              > new physics/tracking scheme
 42 //              1998-08-25 by Stefano Magni        44 //              1998-08-25 by Stefano Magni
 43 //              > Change process to use G4Mate     45 //              > Change process to use G4MaterialPropertiesTables
 44 //              1998-09-03 by Peter Gumplinger     46 //              1998-09-03 by Peter Gumplinger
 45 //              > Protect G4MaterialPropertyVe     47 //              > Protect G4MaterialPropertyVector* AttenuationLengthVector
                                                   >>  48 // mail:        gum@triumf.ca
                                                   >>  49 //              magni@mi.infn.it
 46 //                                                 50 //
 47 //////////////////////////////////////////////     51 ////////////////////////////////////////////////////////////////////////
 48                                                    52 
 49 #include "G4ios.hh"                                53 #include "G4ios.hh"
 50 #include "G4OpProcessSubType.hh"                   54 #include "G4OpProcessSubType.hh"
 51 #include "G4OpticalParameters.hh"              << 
 52                                                    55 
 53 #include "G4OpAbsorption.hh"                       56 #include "G4OpAbsorption.hh"
 54                                                    57 
 55 //....oooOO0OOooo........oooOO0OOooo........oo <<  58 /////////////////////////
                                                   >>  59 // Class Implementation
                                                   >>  60 /////////////////////////
                                                   >>  61 
                                                   >>  62         //////////////
                                                   >>  63         // Operators
                                                   >>  64         //////////////
                                                   >>  65 
                                                   >>  66 // G4OpAbsorption::operator=(const G4OpAbsorption &right)
                                                   >>  67 // {
                                                   >>  68 // }
                                                   >>  69 
                                                   >>  70         /////////////////
                                                   >>  71         // Constructors
                                                   >>  72         /////////////////
                                                   >>  73 
 56 G4OpAbsorption::G4OpAbsorption(const G4String&     74 G4OpAbsorption::G4OpAbsorption(const G4String& processName, G4ProcessType type)
 57   : G4VDiscreteProcess(processName, type)      <<  75               : G4VDiscreteProcess(processName, type)
 58 {                                                  76 {
 59   Initialise();                                <<  77         if (verboseLevel>0) {
 60   if(verboseLevel > 0)                         <<  78            G4cout << GetProcessName() << " is created " << G4endl;
 61   {                                            <<  79         }
 62     G4cout << GetProcessName() << " is created << 
 63   }                                            << 
 64   SetProcessSubType(fOpAbsorption);            << 
 65 }                                              << 
 66                                                << 
 67 //....oooOO0OOooo........oooOO0OOooo........oo << 
 68 G4OpAbsorption::~G4OpAbsorption() = default;   << 
 69                                                    80 
 70 //....oooOO0OOooo........oooOO0OOooo........oo <<  81         SetProcessSubType(fOpAbsorption);
 71 void G4OpAbsorption::PreparePhysicsTable(const << 
 72 {                                              << 
 73   Initialise();                                << 
 74 }                                                  82 }
 75                                                    83 
 76 //....oooOO0OOooo........oooOO0OOooo........oo <<  84 // G4OpAbsorption::G4OpAbsorption(const G4OpAbsorpton &right)
 77 void G4OpAbsorption::Initialise()              <<  85 // {
 78 {                                              <<  86 // }
 79   SetVerboseLevel(G4OpticalParameters::Instanc <<  87 
 80 }                                              <<  88         ////////////////
                                                   >>  89         // Destructors
                                                   >>  90         ////////////////
                                                   >>  91 
                                                   >>  92 G4OpAbsorption::~G4OpAbsorption(){}
                                                   >>  93 
                                                   >>  94         ////////////
                                                   >>  95         // Methods
                                                   >>  96         ////////////
 81                                                    97 
 82 //....oooOO0OOooo........oooOO0OOooo........oo <<  98 // PostStepDoIt
 83 G4VParticleChange* G4OpAbsorption::PostStepDoI <<  99 // -------------
 84                                                << 100 //
                                                   >> 101 G4VParticleChange*
                                                   >> 102 G4OpAbsorption::PostStepDoIt(const G4Track& aTrack, const G4Step& aStep)
 85 {                                                 103 {
 86   aParticleChange.Initialize(aTrack);          << 104         aParticleChange.Initialize(aTrack);
 87                                                << 
 88   const G4DynamicParticle* aParticle = aTrack. << 
 89   G4double thePhotonMomentum         = aPartic << 
 90                                                   105 
 91   aParticleChange.ProposeLocalEnergyDeposit(th << 106         aParticleChange.ProposeTrackStatus(fStopAndKill);
 92   aParticleChange.ProposeTrackStatus(fStopAndK << 
 93                                                   107 
 94   if(verboseLevel > 1)                         << 108         if (verboseLevel>0) {
 95   {                                            << 109      G4cout << "\n** Photon absorbed! **" << G4endl;
 96     G4cout << "\n** OpAbsorption: Photon absor << 110         }
 97   }                                            << 111         return G4VDiscreteProcess::PostStepDoIt(aTrack, aStep);
 98   return G4VDiscreteProcess::PostStepDoIt(aTra << 
 99 }                                                 112 }
100                                                   113 
101 //....oooOO0OOooo........oooOO0OOooo........oo << 
102 G4double G4OpAbsorption::GetMeanFreePath(const << 
103                                          G4For << 
104 {                                              << 
105   const G4DynamicParticle* aParticle = aTrack. << 
106   G4MaterialPropertiesTable* MPT =             << 
107     aTrack.GetMaterial()->GetMaterialPropertie << 
108   G4double attLength = DBL_MAX;                << 
109                                                << 
110   if(MPT)                                      << 
111   {                                            << 
112     G4MaterialPropertyVector* attVector = MPT- << 
113     if(attVector)                              << 
114     {                                          << 
115       attLength =                              << 
116         attVector->Value(aParticle->GetTotalMo << 
117     }                                          << 
118   }                                            << 
119                                                   114 
120   return attLength;                            << 115 // GetMeanFreePath
121 }                                              << 116 // ---------------
                                                   >> 117 //
                                                   >> 118 G4double G4OpAbsorption::GetMeanFreePath(const G4Track& aTrack,
                                                   >> 119                  G4double ,
                                                   >> 120                  G4ForceCondition* )
                                                   >> 121 {
                                                   >> 122   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
                                                   >> 123         const G4Material* aMaterial = aTrack.GetMaterial();
                                                   >> 124 
                                                   >> 125   G4double thePhotonMomentum = aParticle->GetTotalMomentum();
                                                   >> 126 
                                                   >> 127   G4MaterialPropertiesTable* aMaterialPropertyTable;
                                                   >> 128   G4MaterialPropertyVector* AttenuationLengthVector;
                                                   >> 129   
                                                   >> 130         G4double AttenuationLength = DBL_MAX;
                                                   >> 131 
                                                   >> 132   aMaterialPropertyTable = aMaterial->GetMaterialPropertiesTable();
                                                   >> 133 
                                                   >> 134   if ( aMaterialPropertyTable ) {
                                                   >> 135      AttenuationLengthVector = aMaterialPropertyTable->
                                                   >> 136                                                 GetProperty("ABSLENGTH");
                                                   >> 137            if ( AttenuationLengthVector ){
                                                   >> 138              AttenuationLength = AttenuationLengthVector->
                                                   >> 139                                          GetProperty (thePhotonMomentum);
                                                   >> 140            }
                                                   >> 141            else {
                                                   >> 142 //             G4cout << "No Absorption length specified" << G4endl;
                                                   >> 143            }
                                                   >> 144         } 
                                                   >> 145         else {
                                                   >> 146 //           G4cout << "No Absorption length specified" << G4endl;
                                                   >> 147         }
122                                                   148 
123 //....oooOO0OOooo........oooOO0OOooo........oo << 149         return AttenuationLength;
124 void G4OpAbsorption::SetVerboseLevel(G4int ver << 
125 {                                              << 
126   verboseLevel = verbose;                      << 
127   G4OpticalParameters::Instance()->SetAbsorpti << 
128 }                                                 150 }
129                                                   151