Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/solidstate/phonon/src/G4PhononReflection.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/solidstate/phonon/src/G4PhononReflection.cc (Version 11.3.0) and /processes/solidstate/phonon/src/G4PhononReflection.cc (Version 10.4)


  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 /// \file processes/phonon/src/G4PhononReflect     26 /// \file processes/phonon/src/G4PhononReflection.cc
 27 /// \brief Implementation of the G4PhononRefle     27 /// \brief Implementation of the G4PhononReflection class
 28 //                                                 28 //
 29 // This process handles the interaction of pho     29 // This process handles the interaction of phonons with
 30 // boundaries. Implementation of this class is     30 // boundaries. Implementation of this class is highly 
 31 // geometry dependent.Currently, phonons are k     31 // geometry dependent.Currently, phonons are killed when
 32 // they reach a boundary. If the other side of     32 // they reach a boundary. If the other side of the 
 33 // boundary was Al, a hit is registered.           33 // boundary was Al, a hit is registered.
 34 //                                                 34 //  
                                                   >>  35 // $Id: G4PhononReflection.cc 76799 2013-11-15 20:30:53Z mkelsey $
 35 //                                                 36 //
 36 // 20131115  Throw exception if track's polari     37 // 20131115  Throw exception if track's polarization state is invalid.
 37                                                    38 
 38 #include "G4PhononReflection.hh"                   39 #include "G4PhononReflection.hh"
 39 #include "G4ExceptionSeverity.hh"                  40 #include "G4ExceptionSeverity.hh"
 40 #include "G4GeometryTolerance.hh"                  41 #include "G4GeometryTolerance.hh"
 41 #include "G4LatticePhysical.hh"                    42 #include "G4LatticePhysical.hh"
 42 #include "G4PhononLong.hh"                         43 #include "G4PhononLong.hh"
 43 #include "G4PhononTransFast.hh"                    44 #include "G4PhononTransFast.hh"
 44 #include "G4PhononTransSlow.hh"                    45 #include "G4PhononTransSlow.hh"
 45 #include "G4PhysicalConstants.hh"                  46 #include "G4PhysicalConstants.hh"
 46 #include "G4Step.hh"                               47 #include "G4Step.hh"
 47 #include "G4StepPoint.hh"                          48 #include "G4StepPoint.hh"
 48 #include "G4SystemOfUnits.hh"                      49 #include "G4SystemOfUnits.hh"
 49 #include "G4VParticleChange.hh"                    50 #include "G4VParticleChange.hh"
 50                                                    51 
 51                                                    52 
 52 G4PhononReflection::G4PhononReflection(const G     53 G4PhononReflection::G4PhononReflection(const G4String& aName)
 53   : G4VPhononProcess(aName),                       54   : G4VPhononProcess(aName),
 54     kCarTolerance(G4GeometryTolerance::GetInst     55     kCarTolerance(G4GeometryTolerance::GetInstance()->GetSurfaceTolerance()) {;}
 55                                                    56 
 56 G4PhononReflection::~G4PhononReflection() {;}      57 G4PhononReflection::~G4PhononReflection() {;}
 57                                                    58 
 58 //....oooOO0OOooo........oooOO0OOooo........oo     59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 59                                                    60 
 60 // Always return DBL_MAX and Forced. This ensu     61 // Always return DBL_MAX and Forced. This ensures that the process is
 61 // called at the end of every step. In PostSte     62 // called at the end of every step. In PostStepDoIt the process
 62 // decides whether the step encountered a volu     63 // decides whether the step encountered a volume boundary and a
 63 // reflection should be applied                    64 // reflection should be applied
 64                                                    65 
 65 G4double G4PhononReflection::GetMeanFreePath(c     66 G4double G4PhononReflection::GetMeanFreePath(const G4Track&, G4double,
 66                G4ForceCondition* condition) {      67                G4ForceCondition* condition) {
 67   *condition = Forced;                             68   *condition = Forced;
 68   return DBL_MAX;                                  69   return DBL_MAX;
 69 }                                                  70 }
 70                                                    71 
 71 //....oooOO0OOooo........oooOO0OOooo........oo     72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 72                                                    73 
 73 // This process handles the interaction of pho     74 // This process handles the interaction of phonons with
 74 // boundaries. Implementation of this class is     75 // boundaries. Implementation of this class is highly geometry
 75 // dependent.Currently, phonons are killed whe     76 // dependent.Currently, phonons are killed when they reach a
 76 // boundary. If the other side of the boundary     77 // boundary. If the other side of the boundary was Al, a hit is
 77 // registered.                                     78 // registered.
 78                                                    79   
 79 G4VParticleChange* G4PhononReflection::PostSte     80 G4VParticleChange* G4PhononReflection::PostStepDoIt(const G4Track& aTrack,
 80                 const G4Step& aStep) {             81                 const G4Step& aStep) { 
 81   aParticleChange.Initialize(aTrack);              82   aParticleChange.Initialize(aTrack);
 82                                                    83    
 83   //Check if current step is limited by a volu     84   //Check if current step is limited by a volume boundary
 84   G4StepPoint* postStepPoint = aStep.GetPostSt     85   G4StepPoint* postStepPoint = aStep.GetPostStepPoint();
 85   if (postStepPoint->GetStepStatus()!=fGeomBou     86   if (postStepPoint->GetStepStatus()!=fGeomBoundary) {
 86     //make sure that correct phonon velocity i     87     //make sure that correct phonon velocity is used after the step
 87     int pol = GetPolarization(aTrack);             88     int pol = GetPolarization(aTrack);
 88     if (pol < 0 || pol > 2) {                      89     if (pol < 0 || pol > 2) {
 89       G4Exception("G4PhononReflection::PostSte     90       G4Exception("G4PhononReflection::PostStepDoIt","Phonon001",
 90       EventMustBeAborted, "Track is not a phon     91       EventMustBeAborted, "Track is not a phonon");
 91       return &aParticleChange;    // NOTE: Wil     92       return &aParticleChange;    // NOTE: Will never get here
 92     }                                              93     }
 93                                                    94 
 94     // FIXME:  This should be using wave-vecto     95     // FIXME:  This should be using wave-vector, shouldn't it?
 95     G4double vg = theLattice->MapKtoV(pol, aTr     96     G4double vg = theLattice->MapKtoV(pol, aTrack.GetMomentumDirection());
 96                                                    97     
 97     //Since step was not a volume boundary, ju     98     //Since step was not a volume boundary, just set correct phonon velocity and return
 98     aParticleChange.ProposeVelocity(vg);           99     aParticleChange.ProposeVelocity(vg);
 99     return &aParticleChange;                      100     return &aParticleChange;
100   }                                               101   }
101                                                   102   
102   // do nothing but return is the step is too     103   // do nothing but return is the step is too short
103   // This is to allow actual reflection where     104   // This is to allow actual reflection where after
104   // the first boundary crossing a second, inf    105   // the first boundary crossing a second, infinitesimal
105   // step occurs crossing back into the origin    106   // step occurs crossing back into the original volume
106   if (aTrack.GetStepLength()<=kCarTolerance/2)    107   if (aTrack.GetStepLength()<=kCarTolerance/2) { 
107     return &aParticleChange;                      108     return &aParticleChange;
108   }                                               109   }
109                                                   110   
110   G4double eKin = aTrack.GetKineticEnergy();      111   G4double eKin = aTrack.GetKineticEnergy();     
111   aParticleChange.ProposeNonIonizingEnergyDepo    112   aParticleChange.ProposeNonIonizingEnergyDeposit(eKin);
112   aParticleChange.ProposeTrackStatus(fStopAndK    113   aParticleChange.ProposeTrackStatus(fStopAndKill);
113                                                   114   
114   return &aParticleChange;                        115   return &aParticleChange; 
115 }                                                 116 }
116                                                   117 
117 //....oooOO0OOooo........oooOO0OOooo........oo    118 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
118                                                   119 
119                                                   120 
120                                                   121