Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/amsEcal/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/advanced/amsEcal/src/SteppingAction.cc (Version 11.3.0) and /examples/advanced/amsEcal/src/SteppingAction.cc (Version 10.2.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: SteppingAction.cc 83010 2014-07-24 14:53:07Z gcosmo $
 26 //                                                 27 //
 27 //....oooOO0OOooo........oooOO0OOooo........oo     28 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 28 //....oooOO0OOooo........oooOO0OOooo........oo     29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 29                                                    30 
 30 #include "SteppingAction.hh"                       31 #include "SteppingAction.hh"
 31                                                    32 
 32 #include "DetectorConstruction.hh"                 33 #include "DetectorConstruction.hh"
 33 #include "EventAction.hh"                          34 #include "EventAction.hh"
 34                                                    35 
 35 #include "G4Step.hh"                               36 #include "G4Step.hh"
 36                                                    37 
 37 //....oooOO0OOooo........oooOO0OOooo........oo     38 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 38                                                    39 
 39 SteppingAction::SteppingAction(DetectorConstru     40 SteppingAction::SteppingAction(DetectorConstruction* det, EventAction* evt)
 40 :G4UserSteppingAction(),detector(det),eventAct     41 :G4UserSteppingAction(),detector(det),eventAct(evt)
 41 {                                                  42 {
 42   first = true;                                    43   first = true;
 43   lvol_world = lvol_module = lvol_layer = lvol     44   lvol_world = lvol_module = lvol_layer = lvol_fiber = 0;
 44 }                                                  45 } 
 45                                                    46 
 46 //....oooOO0OOooo........oooOO0OOooo........oo     47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 47                                                    48 
 48 SteppingAction::~SteppingAction()                  49 SteppingAction::~SteppingAction()
 49 {}                                                 50 {}
 50                                                    51 
 51 //....oooOO0OOooo........oooOO0OOooo........oo     52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 52                                                    53 
 53 void SteppingAction::UserSteppingAction(const      54 void SteppingAction::UserSteppingAction(const G4Step* step )
 54 {                                                  55 { 
 55  //some initialisation                             56  //some initialisation
 56  //                                                57  // 
 57  if (first) {                                      58  if (first) {
 58    lvol_world  = detector->GetLvolWorld();         59    lvol_world  = detector->GetLvolWorld();
 59    lvol_module = detector->GetLvolModule();        60    lvol_module = detector->GetLvolModule();   
 60    lvol_layer  = detector->GetLvolLayer();         61    lvol_layer  = detector->GetLvolLayer();
 61    lvol_fiber  = detector->GetLvolFiber();         62    lvol_fiber  = detector->GetLvolFiber();      
 62    first = false;                                  63    first = false;   
 63  }                                                 64  }
 64                                                    65  
 65  //if no edep, return                              66  //if no edep, return
 66  //                                                67  //
 67  G4double edep = step->GetTotalEnergyDeposit()     68  G4double edep = step->GetTotalEnergyDeposit();
 68  ///if (edep == 0.) return;                        69  ///if (edep == 0.) return;
 69                                                    70  
 70  //locate point in geometry                        71  //locate point in geometry
 71  //                                                72  //
 72  G4int iModule = 0;                                73  G4int iModule = 0;
 73  G4int iLayer  = 0;                                74  G4int iLayer  = 0;
 74  G4int iFiber  = 0;                                75  G4int iFiber  = 0;
 75                                                    76    
 76  G4TouchableHandle touch1 = step->GetPreStepPo     77  G4TouchableHandle touch1 = step->GetPreStepPoint()->GetTouchableHandle(); 
 77  G4LogicalVolume* lvol = touch1->GetVolume()->     78  G4LogicalVolume* lvol = touch1->GetVolume()->GetLogicalVolume();
 78                                                    79   
 79       if (lvol == lvol_world) return;              80       if (lvol == lvol_world) return;
 80  else if (lvol == lvol_module) { iModule = tou     81  else if (lvol == lvol_module) { iModule = touch1->GetCopyNumber(0);}
 81  else if (lvol == lvol_layer)  { iLayer  = tou     82  else if (lvol == lvol_layer)  { iLayer  = touch1->GetCopyNumber(0);
 82                                  iModule = tou     83                                  iModule = touch1->GetCopyNumber(1);}
 83  else if (lvol == lvol_fiber)  { iFiber  = tou     84  else if (lvol == lvol_fiber)  { iFiber  = touch1->GetCopyNumber(0);
 84                                  iLayer  = tou     85                                  iLayer  = touch1->GetCopyNumber(1);
 85                                  iModule = tou     86                                  iModule = touch1->GetCopyNumber(2);}
 86                                                    87 
 87  // sum edep                                       88  // sum edep
 88  //                                                89  //
 89  eventAct->SumDeStep(iModule, iLayer, iFiber,      90  eventAct->SumDeStep(iModule, iLayer, iFiber, edep);         
 90 }                                                  91 }
 91                                                    92 
 92 //....oooOO0OOooo........oooOO0OOooo........oo     93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 93                                                    94 
 94 G4double SteppingAction::BirksAttenuation(cons     95 G4double SteppingAction::BirksAttenuation(const G4Step* aStep)
 95 {                                                  96 {
 96  //Example of Birk attenuation law in organic      97  //Example of Birk attenuation law in organic scintillators.
 97  //adapted from Geant3 PHYS337. See MIN 80 (19     98  //adapted from Geant3 PHYS337. See MIN 80 (1970) 239-244
 98  //                                                99  //
 99  const G4Material* material = aStep->GetTrack( << 100  G4Material* material = aStep->GetTrack()->GetMaterial();
100  G4double birk1       = material->GetIonisatio    101  G4double birk1       = material->GetIonisation()->GetBirksConstant();
101  G4double destep      = aStep->GetTotalEnergyD    102  G4double destep      = aStep->GetTotalEnergyDeposit();
102  G4double stepl       = aStep->GetStepLength()    103  G4double stepl       = aStep->GetStepLength();  
103  G4double charge      = aStep->GetTrack()->Get    104  G4double charge      = aStep->GetTrack()->GetDefinition()->GetPDGCharge();
104  //                                               105  //
105  G4double response = destep;                      106  G4double response = destep;
106  if (birk1*destep*stepl*charge != 0.)             107  if (birk1*destep*stepl*charge != 0.)
107    {                                              108    {
108      response = destep/(1. + birk1*destep/step    109      response = destep/(1. + birk1*destep/stepl);
109    }                                              110    }
110  return response;                                 111  return response;
111 }                                                 112 }
112                                                   113 
113 //....oooOO0OOooo........oooOO0OOooo........oo    114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
114                                                   115 
115                                                   116