Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/doiPET/doiPET.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 ]

  1 // ********************************************************************
  2 // * License and Disclaimer                                           *
  3 // *                                                                  *
  4 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  5 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  6 // * conditions of the Geant4 Software License,  included in the file *
  7 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  8 // * include a list of copyright holders.                             *
  9 // *                                                                  *
 10 // * Neither the authors of this software system, nor their employing *
 11 // * institutes,nor the agencies providing financial support for this *
 12 // * work  make  any representation or  warranty, express or implied, *
 13 // * regarding  this  software system or assume any liability for its *
 14 // * use.  Please see the license in the file  LICENSE  and URL above *
 15 // * for the full disclaimer and the limitation of liability.         *
 16 // *                                                                  *
 17 // * This  code  implementation is the result of  the  scientific and *
 18 // * technical work of the GEANT4 collaboration.                      *
 19 // * By using,  copying,  modifying or  distributing the software (or *
 20 // * any work based  on the software)  you  agree  to acknowledge its *
 21 // * use  in  resulting  scientific  publications,  and indicate your *
 22 // * acceptance of all terms of the Geant4 Software license.          *
 23 // ********************************************************************
 24 
 25 //GEANT4 - Depth-of-Interaction enabled Positron emission tomography (PET) advanced example 
 26 
 27 //Authors and contributors
 28 
 29 // Author list to be updated, with names of co-authors and contributors from National Institute of Radiological Sciences (NIRS)
 30 
 31 // Abdella M. Ahmed (1, 2), Andrew Chacon (1, 2), Harley Rutherford (1, 2),
 32 // Hideaki Tashima (3), Go Akamatsu (3), Akram Mohammadi (3), Eiji Yoshida (3), Taiga Yamaya (3)
 33 // Susanna Guatelli (2), and Mitra Safavi-Naeini (1, 2)
 34 
 35 // (1) Australian Nuclear Science and Technology Organisation, Australia
 36 // (2) University of Wollongong, Australia
 37 // (3) National Institute of Radiological Sciences, Japan
 38 
 39 #include "doiPETDetectorConstruction.hh"
 40 #include "doiPETPhysicsList.hh"
 41 #include "doiPETAnalysis.hh"
 42 #include "doiPETActionInitialization.hh"
 43 #include "Randomize.hh"
 44 #include "G4UImanager.hh"
 45 #include "G4VisExecutive.hh"
 46 #include "G4UIExecutive.hh"
 47 #include "G4SystemOfUnits.hh"
 48 #include "G4RunManagerFactory.hh"
 49 //
 50 //////////////////////////////////////////////////////////////////////////////
 51 
 52 int main(int argc,char** argv)
 53 {
 54   auto* runManager = G4RunManagerFactory::CreateRunManager();
 55   G4int nThreads = 4;
 56   runManager->SetNumberOfThreads(nThreads);
 57    
 58   G4UIExecutive* ui = 0;
 59   if ( argc == 1 ) {
 60     ui = new G4UIExecutive(argc, argv);
 61   }
 62   
 63   //Initialize analysis
 64   doiPETAnalysis* ptrAnalysis = doiPETAnalysis::GetInstance();
 65 
 66   runManager->SetUserInitialization(new doiPETDetectorConstruction());
 67 
 68   runManager->SetUserInitialization(new doiPETPhysicsList());
 69 
 70   // Set user action initialization
 71   runManager->SetUserInitialization(new doiPETActionInitialization(ptrAnalysis));
 72 
 73   G4double act  = 1000000 * becquerel;//Activity is set via run.mac file
 74   ptrAnalysis->SetActivity(act);
 75 
 76   G4double halfLife = 109.771 * 60 * s; //Halflife of F-18 as a default
 77   ptrAnalysis -> SetIsotopeHalfLife(halfLife);
 78 
 79   //Blurring specification of the scanner. see inputParameter.txt
 80   ptrAnalysis -> BlurringParameters();
 81 
 82   //Open file to write the output of the simulation
 83   ptrAnalysis->Open("result"); //file extention is affixed based on the type of the output (.root for root or .data for ascii)
 84 
 85 
 86   //
 87   ptrAnalysis -> PMTPosition();
 88   //Read reflector pattern from the inputParameter.txt file
 89   ptrAnalysis->ReadReflectorPattern();
 90 
 91   // Get the pointer to the User Interface manager
 92   G4UImanager* UImanager = G4UImanager::GetUIpointer();
 93 
 94   if ( ! ui ) { 
 95     // batch mode
 96     G4String command = "/control/execute ";
 97     G4String fileName = argv[1];
 98     UImanager->ApplyCommand(command+fileName);
 99   }
100   else { 
101     // interactive mode
102   G4VisManager* visManager = new G4VisExecutive;
103     visManager->Initialize();
104     UImanager->ApplyCommand("/control/execute init_vis.mac");
105     ui->SessionStart();
106     delete ui;
107   delete visManager;
108   }
109 
110   //close the file
111   ptrAnalysis->Close();
112   ptrAnalysis->Delete();
113   
114   delete runManager;
115   return 0;
116 
117 }
118 
119 //////////////////////////////////////////////////////////////////////////////
120 
121