Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/hadrontherapy/src/HadrontherapyEventAction.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/hadrontherapy/src/HadrontherapyEventAction.cc (Version 11.3.0) and /examples/advanced/hadrontherapy/src/HadrontherapyEventAction.cc (Version 9.4.p3)


  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 // Hadrontherapy advanced example for Geant4   <<  26 // This is the *BASIC* version of Hadrontherapy, a Geant4-based application
 27 // See more at: https://twiki.cern.ch/twiki/bi <<  27 // See more at: http://g4advancedexamples.lngs.infn.it/Examples/hadrontherapy
                                                   >>  28 //
                                                   >>  29 // Visit the Hadrontherapy web site (http://www.lns.infn.it/link/Hadrontherapy) to request 
                                                   >>  30 // the *COMPLETE* version of this program, together with its documentation;
                                                   >>  31 // Hadrontherapy (both basic and full version) are supported by the Italian INFN
                                                   >>  32 // Institute in the framework of the MC-INFN Group
                                                   >>  33 //
 28                                                    34 
 29 #include "G4SystemOfUnits.hh"                  << 
 30 #include "G4Event.hh"                              35 #include "G4Event.hh"
 31 #include "G4EventManager.hh"                       36 #include "G4EventManager.hh"
 32 #include "G4HCofThisEvent.hh"                      37 #include "G4HCofThisEvent.hh"
 33 #include "G4VHitsCollection.hh"                    38 #include "G4VHitsCollection.hh"
 34 #include "G4SDManager.hh"                          39 #include "G4SDManager.hh"
 35 #include "G4VVisManager.hh"                        40 #include "G4VVisManager.hh"
 36                                                    41 
 37 #include "HadrontherapyEventAction.hh"             42 #include "HadrontherapyEventAction.hh"
 38 #include "HadrontherapyDetectorHit.hh"             43 #include "HadrontherapyDetectorHit.hh"
 39 #include "HadrontherapyDetectorSD.hh"              44 #include "HadrontherapyDetectorSD.hh"
 40 #include "HadrontherapyDetectorConstruction.hh     45 #include "HadrontherapyDetectorConstruction.hh"
 41 #include "HadrontherapyMatrix.hh"                  46 #include "HadrontherapyMatrix.hh"
 42 #include "HadrontherapyEventActionMessenger.hh     47 #include "HadrontherapyEventActionMessenger.hh"
 43                                                    48 
 44 //////////////////////////////////////////////     49 /////////////////////////////////////////////////////////////////////////////
 45 HadrontherapyEventAction::HadrontherapyEventAc     50 HadrontherapyEventAction::HadrontherapyEventAction() :
 46 drawFlag("all" ),printModulo(10), pointerEvent <<  51   drawFlag("all" ),printModulo(1000), pointerEventMessenger(0)
 47 {                                              <<  52 { 
 48     hitsCollectionID = -1;                     <<  53   hitsCollectionID = -1;
 49     pointerEventMessenger = new HadrontherapyE <<  54   pointerEventMessenger = new HadrontherapyEventActionMessenger(this);
 50 }                                                  55 }
 51                                                    56 
 52 //////////////////////////////////////////////     57 /////////////////////////////////////////////////////////////////////////////
 53 HadrontherapyEventAction::~HadrontherapyEventA     58 HadrontherapyEventAction::~HadrontherapyEventAction()
 54 {                                                  59 {
 55     delete pointerEventMessenger;              <<  60  delete pointerEventMessenger;
 56 }                                                  61 }
 57                                                    62 
 58 //////////////////////////////////////////////     63 /////////////////////////////////////////////////////////////////////////////
 59 void HadrontherapyEventAction::BeginOfEventAct <<  64 void HadrontherapyEventAction::BeginOfEventAction(const G4Event* evt)
 60 {                                              <<  65 { 
 61     G4SDManager* pSDManager = G4SDManager::Get <<  66   G4int evtNb = evt->GetEventID();
 62     if(hitsCollectionID == -1)                 <<  67   
 63         hitsCollectionID = pSDManager -> GetCo <<  68   //printing survey
                                                   >>  69   if (evtNb%printModulo == 0)
                                                   >>  70      G4cout << "\n---> Begin of Event: " << evtNb << G4endl;
                                                   >>  71    
                                                   >>  72   G4SDManager* pSDManager = G4SDManager::GetSDMpointer();
                                                   >>  73   if(hitsCollectionID == -1)
                                                   >>  74     hitsCollectionID = pSDManager -> GetCollectionID("HadrontherapyDetectorHitsCollection");
                                                   >>  75   
 64 }                                                  76 }
 65                                                    77 
 66 //////////////////////////////////////////////     78 /////////////////////////////////////////////////////////////////////////////
 67 void HadrontherapyEventAction::EndOfEventActio     79 void HadrontherapyEventAction::EndOfEventAction(const G4Event* evt)
 68 {                                              <<  80 { 
 69     if(hitsCollectionID < 0)                   <<  81   if(hitsCollectionID < 0)
 70         return;                                <<  82   return;
 71     G4HCofThisEvent* HCE = evt -> GetHCofThisE <<  83   G4HCofThisEvent* HCE = evt -> GetHCofThisEvent();
 72                                                <<  84 
 73     // Clear voxels hit list                   <<  85   // Clear voxels hit list 
 74     HadrontherapyMatrix* matrix = Hadrontherap <<  86   HadrontherapyMatrix* matrix = HadrontherapyMatrix::GetInstance();
 75     if (matrix) matrix -> ClearHitTrack();     <<  87   if (matrix) matrix -> ClearHitTrack(); 
 76                                                <<  88 
 77     if(HCE)                                    <<  89   if(HCE)
 78     {                                          <<  90   {
 79         HadrontherapyDetectorHitsCollection* C <<  91     HadrontherapyDetectorHitsCollection* CHC = (HadrontherapyDetectorHitsCollection*)(HCE -> GetHC(hitsCollectionID));
 80         if(CHC)                                <<  92     if(CHC)
 81         {                                      <<  93      {
 82             if(matrix)                         <<  94        if(matrix)
 83             {                                  <<  95     { 
 84                 // Fill the matrix with the in <<  96         // Fill the matrix with the information: voxel and associated energy deposit 
 85                 // in the detector at the end  <<  97           // in the detector at the end of the event
 86                                                <<  98 
 87                 size_t HitCount = CHC -> entri <<  99     G4int HitCount = CHC -> entries();
 88                 for (size_t h=0; h<HitCount; h << 100     for (G4int h=0; h<HitCount; h++)
 89                 {                              << 101       {
 90                     G4int i = ((*CHC)[h]) -> G << 102         G4int i = ((*CHC)[h]) -> GetXID();
 91                     G4int j = ((*CHC)[h]) -> G << 103         G4int j = ((*CHC)[h]) -> GetYID();
 92                     G4int k = ((*CHC)[h]) -> G << 104         G4int k = ((*CHC)[h]) -> GetZID();
 93                     G4double energyDeposit = ( << 105               G4double energyDeposit = ((*CHC)[h]) -> GetEdep();
 94                     matrix -> Fill(i, j, k, en << 106               matrix -> Fill(i, j, k, energyDeposit/MeV);              
 95                 }                              << 107       }
 96             }                                  << 108     }
 97         }                                      << 
 98     }                                             109     }
                                                   >> 110   }
 99 }                                                 111 }
                                                   >> 112 
100                                                   113