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 8.1.p1)


  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.7 2006/06/29 21:08:50 gunter Exp $
                                                   >>  28 // GEANT4 tag $Name: geant4-08-01-patch-01 $
 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"               << 
 51 #include "G4OpticalParameters.hh"              << 
 52                                                << 
 53 #include "G4OpAbsorption.hh"                       54 #include "G4OpAbsorption.hh"
 54                                                    55 
 55 //....oooOO0OOooo........oooOO0OOooo........oo <<  56 /////////////////////////
 56 G4OpAbsorption::G4OpAbsorption(const G4String& <<  57 // Class Implementation
 57   : G4VDiscreteProcess(processName, type)      <<  58 /////////////////////////
 58 {                                              <<  59 
 59   Initialise();                                <<  60         //////////////
 60   if(verboseLevel > 0)                         <<  61         // Operators
 61   {                                            <<  62         //////////////
 62     G4cout << GetProcessName() << " is created <<  63 
 63   }                                            <<  64 // G4OpAbsorption::operator=(const G4OpAbsorption &right)
 64   SetProcessSubType(fOpAbsorption);            <<  65 // {
 65 }                                              <<  66 // }
                                                   >>  67 
                                                   >>  68         /////////////////
                                                   >>  69         // Constructors
                                                   >>  70         /////////////////
 66                                                    71 
 67 //....oooOO0OOooo........oooOO0OOooo........oo <<  72 G4OpAbsorption::G4OpAbsorption(const G4String& processName, G4ProcessType type)
 68 G4OpAbsorption::~G4OpAbsorption() = default;   <<  73               : G4VDiscreteProcess(processName, type)
 69                                                << 
 70 //....oooOO0OOooo........oooOO0OOooo........oo << 
 71 void G4OpAbsorption::PreparePhysicsTable(const << 
 72 {                                                  74 {
 73   Initialise();                                <<  75         if (verboseLevel>0) {
                                                   >>  76            G4cout << GetProcessName() << " is created " << G4endl;
                                                   >>  77         }
 74 }                                                  78 }
 75                                                    79 
 76 //....oooOO0OOooo........oooOO0OOooo........oo <<  80 // G4OpAbsorption::G4OpAbsorption(const G4OpAbsorpton &right)
 77 void G4OpAbsorption::Initialise()              <<  81 // {
 78 {                                              <<  82 // }
 79   SetVerboseLevel(G4OpticalParameters::Instanc <<  83 
 80 }                                              <<  84         ////////////////
                                                   >>  85         // Destructors
                                                   >>  86         ////////////////
                                                   >>  87 
                                                   >>  88 G4OpAbsorption::~G4OpAbsorption(){}
                                                   >>  89 
                                                   >>  90         ////////////
                                                   >>  91         // Methods
                                                   >>  92         ////////////
 81                                                    93 
 82 //....oooOO0OOooo........oooOO0OOooo........oo <<  94 // PostStepDoIt
 83 G4VParticleChange* G4OpAbsorption::PostStepDoI <<  95 // -------------
 84                                                <<  96 //
                                                   >>  97 G4VParticleChange*
                                                   >>  98 G4OpAbsorption::PostStepDoIt(const G4Track& aTrack, const G4Step& aStep)
 85 {                                                  99 {
 86   aParticleChange.Initialize(aTrack);          << 100         aParticleChange.Initialize(aTrack);
 87                                                   101 
 88   const G4DynamicParticle* aParticle = aTrack. << 102         aParticleChange.ProposeTrackStatus(fStopAndKill);
 89   G4double thePhotonMomentum         = aPartic << 
 90                                                   103 
 91   aParticleChange.ProposeLocalEnergyDeposit(th << 104         if (verboseLevel>0) {
 92   aParticleChange.ProposeTrackStatus(fStopAndK << 105      G4cout << "\n** Photon absorbed! **" << G4endl;
 93                                                << 106         }
 94   if(verboseLevel > 1)                         << 107         return G4VDiscreteProcess::PostStepDoIt(aTrack, aStep);
 95   {                                            << 
 96     G4cout << "\n** OpAbsorption: Photon absor << 
 97   }                                            << 
 98   return G4VDiscreteProcess::PostStepDoIt(aTra << 
 99 }                                                 108 }
100                                                   109 
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                                                   110 
120   return attLength;                            << 111 // GetMeanFreePath
121 }                                              << 112 // ---------------
                                                   >> 113 //
                                                   >> 114 G4double G4OpAbsorption::GetMeanFreePath(const G4Track& aTrack,
                                                   >> 115                  G4double ,
                                                   >> 116                  G4ForceCondition* )
                                                   >> 117 {
                                                   >> 118   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
                                                   >> 119         const G4Material* aMaterial = aTrack.GetMaterial();
                                                   >> 120 
                                                   >> 121   G4double thePhotonMomentum = aParticle->GetTotalMomentum();
                                                   >> 122 
                                                   >> 123   G4MaterialPropertiesTable* aMaterialPropertyTable;
                                                   >> 124   G4MaterialPropertyVector* AttenuationLengthVector;
                                                   >> 125   
                                                   >> 126         G4double AttenuationLength = DBL_MAX;
                                                   >> 127 
                                                   >> 128   aMaterialPropertyTable = aMaterial->GetMaterialPropertiesTable();
                                                   >> 129 
                                                   >> 130   if ( aMaterialPropertyTable ) {
                                                   >> 131      AttenuationLengthVector = aMaterialPropertyTable->
                                                   >> 132                                                 GetProperty("ABSLENGTH");
                                                   >> 133            if ( AttenuationLengthVector ){
                                                   >> 134              AttenuationLength = AttenuationLengthVector->
                                                   >> 135                                          GetProperty (thePhotonMomentum);
                                                   >> 136            }
                                                   >> 137            else {
                                                   >> 138 //             G4cout << "No Absorption length specified" << G4endl;
                                                   >> 139            }
                                                   >> 140         } 
                                                   >> 141         else {
                                                   >> 142 //           G4cout << "No Absorption length specified" << G4endl;
                                                   >> 143         }
122                                                   144 
123 //....oooOO0OOooo........oooOO0OOooo........oo << 145         return AttenuationLength;
124 void G4OpAbsorption::SetVerboseLevel(G4int ver << 
125 {                                              << 
126   verboseLevel = verbose;                      << 
127   G4OpticalParameters::Instance()->SetAbsorpti << 
128 }                                                 146 }
129                                                   147