Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/transportation/include/G4NeutronKiller.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/transportation/include/G4NeutronKiller.hh (Version 11.3.0) and /processes/transportation/include/G4NeutronKiller.hh (Version 8.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 // $Id: G4NeutronKiller.hh,v 1.1 2006/11/17 15:53:04 vnivanch Exp $
                                                   >>  27 // GEANT4 tag $Name: geant4-08-03-patch-02 $
 26 //                                                 28 //
 27 //--------------------------------------------     29 //---------------------------------------------------------------------------
 28 //                                                 30 //
 29 // ClassName:   G4NeutronKiller                    31 // ClassName:   G4NeutronKiller
 30 //                                                 32 //
 31 // Description: The process to kill neutrons t     33 // Description: The process to kill neutrons to save CPU
 32 //                                                 34 //
 33 // Author:      V.Ivanchenko 26/09/00 for HARP     35 // Author:      V.Ivanchenko 26/09/00 for HARP software
 34 //                                                 36 //
 35 //--------------------------------------------     37 //----------------------------------------------------------------------------
 36 //                                                 38 //
 37 // Class description:                              39 // Class description:
 38 //                                                 40 //
 39 // G4NeutronKiller allows to remove unwanted n     41 // G4NeutronKiller allows to remove unwanted neutrons from simulation in 
 40 // order to improve CPU performance. There are     42 // order to improve CPU performance. There are two parameters:
 41 //                                                 43 //                 
 42 // 1) low energy threshold for neutron kinetic     44 // 1) low energy threshold for neutron kinetic energy (default 0)
 43 // 2) time limit for neutron track (default DB     45 // 2) time limit for neutron track (default DBL_MAX) 
 44 //                                                 46 // 
 45 // These parameters can be changed via Set met     47 // These parameters can be changed via Set methods or by UI commands:
 46 //   /physics_engine/neutron/energyCut             48 //   /physics_engine/neutron/energyCut
 47 //   /physics_engine/neutron/timeLimit             49 //   /physics_engine/neutron/timeLimit
 48 //                                                 50 //
 49 // If a neutron track is killed no energy depo     51 // If a neutron track is killed no energy deposition is added to the step 
 50 //                                                 52 //
 51                                                    53 
 52 //....oooOO0OOooo........oooOO0OOooo........oo     54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 53 //....oooOO0OOooo........oooOO0OOooo........oo     55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 54                                                    56 
 55 #ifndef NeutronKiller_h                            57 #ifndef NeutronKiller_h
 56 #define NeutronKiller_h 1                          58 #define NeutronKiller_h 1
 57                                                    59 
 58 #include "globals.hh"                              60 #include "globals.hh"
 59 #include "G4VDiscreteProcess.hh"                   61 #include "G4VDiscreteProcess.hh"
 60 #include "G4ParticleDefinition.hh"                 62 #include "G4ParticleDefinition.hh"
 61 #include "G4Step.hh"                               63 #include "G4Step.hh"
 62 #include "G4Track.hh"                              64 #include "G4Track.hh"
 63                                                    65 
 64 class G4NeutronKillerMessenger;                    66 class G4NeutronKillerMessenger;
 65                                                    67 
 66 //....oooOO0OOooo........oooOO0OOooo........oo     68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 67                                                    69 
 68 class G4NeutronKiller : public G4VDiscreteProc     70 class G4NeutronKiller : public G4VDiscreteProcess
 69 {                                                  71 {
 70 public:                                            72 public:
 71                                                    73 
 72   G4NeutronKiller(const G4String& processName      74   G4NeutronKiller(const G4String& processName = "nKiller",
 73       G4ProcessType   aType =  fGeneral );     <<  75       G4ProcessType   aType = fNotDefined);
 74                                                    76 
 75   virtual ~G4NeutronKiller();                      77   virtual ~G4NeutronKiller();
 76                                                    78 
 77   virtual G4bool IsApplicable(const G4Particle <<  79   G4bool IsApplicable(const G4ParticleDefinition&);
 78                                                    80 
 79   void  SetTimeLimit(G4double);                    81   void  SetTimeLimit(G4double);
 80                                                    82 
 81   void  SetKinEnergyLimit(G4double);               83   void  SetKinEnergyLimit(G4double);
 82                                                    84 
 83   virtual G4double PostStepGetPhysicalInteract <<  85   G4double PostStepGetPhysicalInteractionLength( const G4Track& track,
 84                G4double previousStepSize,      <<  86              G4double previousStepSize,
 85                G4ForceCondition* condition);   <<  87              G4ForceCondition* condition);
 86                                                    88 
 87   virtual G4VParticleChange* PostStepDoIt(cons <<  89   G4VParticleChange* PostStepDoIt(const G4Track&, const G4Step&);
 88                                                    90 
 89   virtual G4double GetMeanFreePath(const G4Tra <<  91   G4double GetMeanFreePath(const G4Track&, G4double,G4ForceCondition*);
 90                                                    92 
 91 private:                                           93 private:
 92                                                    94 
 93   // hide assignment operator as private           95   // hide assignment operator as private
 94   G4NeutronKiller(const G4NeutronKiller&);         96   G4NeutronKiller(const G4NeutronKiller&);
 95   G4NeutronKiller& operator = (const G4Neutron     97   G4NeutronKiller& operator = (const G4NeutronKiller &right);
 96                                                    98 
 97   G4double kinEnergyThreshold;                     99   G4double kinEnergyThreshold;
 98   G4double timeThreshold;                         100   G4double timeThreshold;
 99                                                   101      
100   G4NeutronKillerMessenger* pMess;                102   G4NeutronKillerMessenger* pMess;
101 };                                                103 };
                                                   >> 104 
                                                   >> 105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 106 
                                                   >> 107 inline G4double G4NeutronKiller::PostStepGetPhysicalInteractionLength( 
                                                   >> 108          const G4Track& aTrack,
                                                   >> 109          G4double, G4ForceCondition* condition)
                                                   >> 110 {
                                                   >> 111   // condition is set to "Not Forced"
                                                   >> 112   *condition = NotForced;
                                                   >> 113   
                                                   >> 114   G4double limit = DBL_MAX; 
                                                   >> 115   if(aTrack.GetGlobalTime() > timeThreshold || 
                                                   >> 116      aTrack.GetKineticEnergy() < kinEnergyThreshold) limit = 0.0;
                                                   >> 117   return limit;
                                                   >> 118 }
                                                   >> 119 
                                                   >> 120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 121 
                                                   >> 122 inline G4double G4NeutronKiller::GetMeanFreePath(const G4Track&,G4double,
                                                   >> 123              G4ForceCondition*)
                                                   >> 124 {
                                                   >> 125   return DBL_MAX;
                                                   >> 126 }    
                                                   >> 127 
                                                   >> 128 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 129 
                                                   >> 130 inline G4VParticleChange* G4NeutronKiller::PostStepDoIt(const G4Track& aTrack, 
                                                   >> 131               const G4Step&)
                                                   >> 132 {
                                                   >> 133   pParticleChange->Initialize(aTrack);
                                                   >> 134   pParticleChange->ProposeTrackStatus(fStopAndKill);
                                                   >> 135   return pParticleChange;
                                                   >> 136 }
102                                                   137 
103 //....oooOO0OOooo........oooOO0OOooo........oo    138 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
104                                                   139 
105 #endif                                            140 #endif
106                                                   141 
107                                                   142