Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/neuron/src/SteppingAction.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 /examples/extended/medical/dna/neuron/src/SteppingAction.cc (Version 11.3.0) and /examples/extended/medical/dna/neuron/src/SteppingAction.cc (Version 10.4.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 // This example is provided by the Geant4-DNA      26 // This example is provided by the Geant4-DNA collaboration
 27 // Any report or published results obtained us     27 // Any report or published results obtained using the Geant4-DNA software
 28 // shall cite the following Geant4-DNA collabo     28 // shall cite the following Geant4-DNA collaboration publication:
 29 // Med. Phys. 37 (2010) 4692-4708                  29 // Med. Phys. 37 (2010) 4692-4708
 30 // and papers                                      30 // and papers
 31 // M. Batmunkh et al. J Radiat Res Appl Sci 8      31 // M. Batmunkh et al. J Radiat Res Appl Sci 8 (2015) 498-507
 32 // O. Belov et al. Physica Medica 32 (2016) 15     32 // O. Belov et al. Physica Medica 32 (2016) 1510-1520
 33 // The Geant4-DNA web site is available at htt     33 // The Geant4-DNA web site is available at http://geant4-dna.org
 34 //                                             <<  34 // 
 35 // -------------------------------------------     35 // -------------------------------------------------------------------
 36 // November 2016                                   36 // November 2016
 37 // -------------------------------------------     37 // -------------------------------------------------------------------
 38 //                                                 38 //
 39 /// \file SteppingAction.cc                    <<  39 // $ID$
                                                   >>  40 /// \file SteppingAction.cc 
 40 /// \brief Implementation of the SteppingActio     41 /// \brief Implementation of the SteppingAction class
 41                                                    42 
                                                   >>  43 #include "Analysis.hh"
 42 #include "SteppingAction.hh"                       44 #include "SteppingAction.hh"
 43                                                << 
 44 #include "PrimaryGeneratorAction.hh"           << 
 45 #include "Run.hh"                              << 
 46 #include "RunAction.hh"                            45 #include "RunAction.hh"
 47                                                <<  46 #include "PrimaryGeneratorAction.hh"
 48 #include "G4AnalysisManager.hh"                <<  47 #include "G4SystemOfUnits.hh"
                                                   >>  48 #include "G4SteppingManager.hh"
                                                   >>  49 #include "G4VTouchable.hh"
                                                   >>  50 //
                                                   >>  51 //#include "NeuronLoadDataFile.hh"
                                                   >>  52 #include "G4ios.hh"  
 49 #include "G4LogicalVolume.hh"                      53 #include "G4LogicalVolume.hh"
 50 #include "G4LogicalVolumeStore.hh"                 54 #include "G4LogicalVolumeStore.hh"
                                                   >>  55 #include "G4VPhysicalVolume.hh"
 51 #include "G4PhysicalVolumeStore.hh"                56 #include "G4PhysicalVolumeStore.hh"
                                                   >>  57 #include "Run.hh"
 52 #include "G4RunManager.hh"                         58 #include "G4RunManager.hh"
 53 #include "G4SteppingManager.hh"                << 
 54 #include "G4SystemOfUnits.hh"                  << 
 55 #include "G4VPhysicalVolume.hh"                << 
 56 #include "G4VTouchable.hh"                     << 
 57 #include "G4ios.hh"                            << 
 58                                                    59 
 59 //....oooOO0OOooo........oooOO0OOooo........oo     60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 60                                                    61 
 61 SteppingAction::SteppingAction(RunAction* run) <<  62 SteppingAction::SteppingAction(RunAction* run)  
                                                   >>  63 : G4UserSteppingAction(), fRunAction(run)
                                                   >>  64 {
                                                   >>  65  //fNeuronLoadParamz = new NeuronLoadDataFile(); 
                                                   >>  66  //fpNeuronHits = new NeuronHitCompartments();
                                                   >>  67 }
 62                                                    68 
 63 //....oooOO0OOooo........oooOO0OOooo........oo     69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 64                                                    70 
 65 void SteppingAction::UserSteppingAction(const  <<  71 SteppingAction::~SteppingAction()
 66 {                                                  72 {
 67   if (step->GetPostStepPoint()->GetProcessDefi << 
 68     G4double edepStep = step->GetTotalEnergyDe << 
 69     G4VPhysicalVolume* volumeStep = step->GetP << 
 70     G4TouchableHandle touchStep = step->GetPre << 
 71     G4VPhysicalVolume* volumeMedium = G4Physic << 
 72     G4VPhysicalVolume* volumeSlice =           << 
 73       G4PhysicalVolumeStore::GetInstance()->Ge << 
 74                                                << 
 75     // count processes                         << 
 76     //                                         << 
 77     const G4StepPoint* endPoint = step->GetPos << 
 78     const G4VProcess* process = endPoint->GetP << 
 79     Run* run = static_cast<Run*>(G4RunManager: << 
 80     run->CountProcesses(process);              << 
 81                                                << 
 82     // Edep in all volume                      << 
 83     if (edepStep > 0.) {                       << 
 84       fRunAction->AddEdepALL(edepStep);        << 
 85                                                << 
 86       // Edep in a target volume (Bounding sli << 
 87       if (volumeStep == volumeSlice) {         << 
 88         fRunAction->AddEdepSlice(edepStep);    << 
 89       }                                        << 
 90                                                << 
 91       // Edep in Soma                          << 
 92       if (volumeStep->GetName() == "Soma") {   << 
 93         fRunAction->AddEdepSoma(edepStep);     << 
 94         run->AddSomaCompart(touchStep->GetCopy << 
 95       }                                        << 
 96                                                << 
 97       // Edep in Dendrites                     << 
 98       if (volumeStep->GetName() == "Dendrites" << 
 99         fRunAction->AddEdepDend(edepStep);     << 
100         run->AddDendCompart(touchStep->GetCopy << 
101       }                                        << 
102                                                << 
103       // Edep in Axon                          << 
104       if (volumeStep->GetName() == "Axon") {   << 
105         fRunAction->AddEdepAxon(edepStep);     << 
106         run->AddAxonCompart(touchStep->GetCopy << 
107       }                                        << 
108                                                << 
109       // Edep in whole Neuron                  << 
110       if (volumeStep != volumeMedium && volume << 
111         fRunAction->AddEdepNeuron(edepStep);   << 
112       }                                        << 
113                                                << 
114       // Edep outside bounding scice           << 
115       if (volumeStep == volumeMedium) {        << 
116         fRunAction->AddEdepMedium(edepStep);   << 
117       }                                        << 
118     }  // end edep>0                           << 
119   }                                            << 
120 }                                                  73 }
                                                   >>  74 
                                                   >>  75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >>  76 
                                                   >>  77 void SteppingAction::UserSteppingAction(const G4Step* step)
                                                   >>  78 { 
                                                   >>  79 
                                                   >>  80   if (step->GetPostStepPoint()->GetProcessDefinedStep()->
                                                   >>  81       GetProcessName()!="Transportation")
                                                   >>  82   {
                                                   >>  83    G4double edepStep = step->GetTotalEnergyDeposit() ;
                                                   >>  84    G4VPhysicalVolume* volumeStep = step->GetPreStepPoint()->
                                                   >>  85                                    GetPhysicalVolume();
                                                   >>  86    G4TouchableHandle touchStep = step->GetPreStepPoint()->GetTouchableHandle();
                                                   >>  87    G4VPhysicalVolume* volumeMedium = G4PhysicalVolumeStore::GetInstance()
                                                   >>  88                                    ->GetVolume("Medium");
                                                   >>  89    G4VPhysicalVolume* volumeSlice = G4PhysicalVolumeStore::GetInstance()
                                                   >>  90                                    ->GetVolume("BoundingSlice");
                                                   >>  91  
                                                   >>  92    // count processes
                                                   >>  93    // 
                                                   >>  94    const G4StepPoint* endPoint = step->GetPostStepPoint();
                                                   >>  95    const G4VProcess* process   = endPoint->GetProcessDefinedStep();
                                                   >>  96    Run* run = static_cast<Run*>(
                                                   >>  97         G4RunManager::GetRunManager()->GetNonConstCurrentRun());
                                                   >>  98    run->CountProcesses(process);
                                                   >>  99 
                                                   >> 100  // Edep in all volume
                                                   >> 101  if (edepStep > 0. ) 
                                                   >> 102  {   
                                                   >> 103    fRunAction->AddEdepALL(edepStep);
                                                   >> 104    //fpNeuronHits->RegisterEdep(volumeStep, edepStep);
                                                   >> 105 
                                                   >> 106    // Edep in a target volume (Bounding slice)
                                                   >> 107    if (volumeStep == volumeSlice) 
                                                   >> 108    {
                                                   >> 109      fRunAction->AddEdepSlice(edepStep);
                                                   >> 110 
                                                   >> 111    }
                                                   >> 112 
                                                   >> 113    // Edep in Soma
                                                   >> 114    if (volumeStep->GetName() == "Soma")
                                                   >> 115    {
                                                   >> 116      fRunAction->AddEdepSoma(edepStep);  
                                                   >> 117   // fRunAction->AddSomaCompart(touchStep->GetCopyNumber(), edepStep);
                                                   >> 118      run->AddSomaCompart(touchStep->GetCopyNumber(), edepStep);    
                                                   >> 119    }
                                                   >> 120 
                                                   >> 121    // Edep in Dendrites
                                                   >> 122    if (volumeStep->GetName() == "Dendrites")
                                                   >> 123    {
                                                   >> 124      fRunAction->AddEdepDend(edepStep);
                                                   >> 125      run->AddDendCompart(touchStep->GetCopyNumber(), edepStep);
                                                   >> 126    }
                                                   >> 127 
                                                   >> 128    // Edep in Axon
                                                   >> 129    if (volumeStep->GetName() == "Axon")
                                                   >> 130    {   
                                                   >> 131      fRunAction->AddEdepAxon(edepStep); 
                                                   >> 132      run->AddAxonCompart(touchStep->GetCopyNumber(), edepStep);
                                                   >> 133    } 
                                                   >> 134   
                                                   >> 135    // Edep in whole Neuron
                                                   >> 136    if (volumeStep != volumeMedium && volumeStep != volumeSlice)
                                                   >> 137    {
                                                   >> 138      fRunAction->AddEdepNeuron(edepStep); 
                                                   >> 139    }
                                                   >> 140 
                                                   >> 141    // Edep outside bounding scice
                                                   >> 142    if (volumeStep == volumeMedium)
                                                   >> 143    {
                                                   >> 144      fRunAction->AddEdepMedium(edepStep);
                                                   >> 145    } 
                                                   >> 146  
                                                   >> 147  } // end edep>0  
                                                   >> 148   }
                                                   >> 149 
                                                   >> 150 }    
121                                                   151