Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/xray_fluorescence/src/XrayFluoSimulation.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/xray_fluorescence/src/XrayFluoSimulation.cc (Version 11.3.0) and /examples/advanced/xray_fluorescence/src/XrayFluoSimulation.cc (Version 10.0.p4)


  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 //                                                 26 //
                                                   >>  27 // $Id: XrayFluoSimulation.cc
 27 //                                                 28 //
 28 // Author: Elena Guardincerri                  <<  29 // Author: Elena Guardincerri 
 29 //                                                 30 //
 30 // History:                                        31 // History:
 31 // -----------                                     32 // -----------
 32 // 28 Nov 2001 Elena Guardincerri     Created      33 // 28 Nov 2001 Elena Guardincerri     Created
 33 // 24 Ago 2002 Splitted in a separet class  Al     34 // 24 Ago 2002 Splitted in a separet class  Alfonso Mantero
 34 // 18 Jan 2011 adapted to new deexcitation des     35 // 18 Jan 2011 adapted to new deexcitation design
 35 //                                                 36 //
 36 // -------------------------------------------     37 // -------------------------------------------------------------------
 37                                                <<  38 #include "G4RunManager.hh"
 38 #include "G4Types.hh"                          << 
 39 #include "G4RunManagerFactory.hh"              << 
 40 #include "G4UImanager.hh"                          39 #include "G4UImanager.hh"
 41 #include "Randomize.hh"                            40 #include "Randomize.hh"
                                                   >>  41 #ifdef G4VIS_USE
 42 #include "G4VisExecutive.hh"                       42 #include "G4VisExecutive.hh"
                                                   >>  43 #endif
                                                   >>  44 #ifdef G4UI_USE
 43 #include "G4UIExecutive.hh"                        45 #include "G4UIExecutive.hh"
                                                   >>  46 #endif
 44 #include "XrayFluoDetectorConstruction.hh"         47 #include "XrayFluoDetectorConstruction.hh"
 45 #include "XrayFluoPlaneDetectorConstruction.hh     48 #include "XrayFluoPlaneDetectorConstruction.hh"
 46 #include "XrayFluoMercuryDetectorConstruction.     49 #include "XrayFluoMercuryDetectorConstruction.hh"
 47 #include "XrayFluoPhysicsList.hh"                  50 #include "XrayFluoPhysicsList.hh"
                                                   >>  51 #include "XrayFluoPrimaryGeneratorAction.hh"
                                                   >>  52 #include "XrayFluoPlanePrimaryGeneratorAction.hh"
                                                   >>  53 #include "XrayFluoMercuryPrimaryGeneratorAction.hh"
                                                   >>  54 #include "XrayFluoRunAction.hh"
                                                   >>  55 #include "XrayFluoEventAction.hh"
                                                   >>  56 #include "XrayFluoSteppingAction.hh"
                                                   >>  57 #include "XrayFluoSteppingVerbose.hh"
 48 #include "XrayFluoSimulation.hh"                   58 #include "XrayFluoSimulation.hh"
 49                                                <<  59 #ifdef G4ANALYSIS_USE
 50 #include "XrayFluoActionInitializer.hh"        << 
 51 #include "XrayFluoAnalysisManager.hh"              60 #include "XrayFluoAnalysisManager.hh"
                                                   >>  61 #endif
                                                   >>  62 using namespace CLHEP;
                                                   >>  63 
                                                   >>  64 XrayFluoSimulation::XrayFluoSimulation(G4int seed):dir(seed)
                                                   >>  65 
                                                   >>  66 { }
                                                   >>  67 
                                                   >>  68 XrayFluoSimulation::~XrayFluoSimulation()
                                                   >>  69 
                                                   >>  70 { }
 52                                                    71 
 53 void XrayFluoSimulation::RunSimulation(int arg     72 void XrayFluoSimulation::RunSimulation(int argc,char* argv[])
 54 {                                                  73 {
 55   auto* runManager = G4RunManagerFactory::Crea <<  74 
 56   G4int nThreads = 4;                          <<  75   // choose the Random engine
 57   runManager->SetNumberOfThreads(nThreads);    <<  76   CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine);
                                                   >>  77   CLHEP::HepRandom::setTheSeed(dir);
                                                   >>  78 
                                                   >>  79   //XrayFluo Verbose output class
                                                   >>  80   G4VSteppingVerbose::SetInstance(new XrayFluoSteppingVerbose);
                                                   >>  81      
                                                   >>  82   // Construct the default run manager
                                                   >>  83   G4RunManager * runManager = new G4RunManager;
                                                   >>  84 
                                                   >>  85   // set mandatory initialization 
                                                   >>  86 
                                                   >>  87   XrayFluoPhysicsList* xrayList = new   XrayFluoPhysicsList();
 58                                                    88 
 59   // chosing Geometry setup                        89   // chosing Geometry setup
 60   G4int geometryNumber = 0;                    <<  90 
                                                   >>  91   G4int geometryNumber = 0; 
 61                                                    92 
 62   if (argc == 3){                                  93   if (argc == 3){
 63     geometryNumber =  atoi(argv[2]);               94     geometryNumber =  atoi(argv[2]);
 64                                                    95 
 65   }                                                96   }
 66   while ( (geometryNumber < 1) || (geometryNum <<  97   while ( (geometryNumber != 1) && (geometryNumber !=2) && (geometryNumber !=3) && (geometryNumber !=4)) {
 67     G4cout << "Please Select Simulation Geomet     98     G4cout << "Please Select Simulation Geometrical Set-Up: "<< G4endl;
 68     G4cout << "1 - Test Beam" << G4endl;           99     G4cout << "1 - Test Beam" << G4endl;
 69     G4cout << "2 - Infinite Plane" << G4endl;     100     G4cout << "2 - Infinite Plane" << G4endl;
 70     G4cout << "3 - Planet and Sun"<< G4endl;      101     G4cout << "3 - Planet and Sun"<< G4endl;
 71     G4cout << "4 - Phase-Space Production"<< G    102     G4cout << "4 - Phase-Space Production"<< G4endl;
 72                                                   103 
                                                   >> 104 
 73     G4cin >> geometryNumber;                      105     G4cin >> geometryNumber;
 74   }                                               106   }
 75                                                   107 
 76   // set analysis to have the messenger runnin << 108   XrayFluoDetectorConstruction* testBeamDetector = 0;
 77   XrayFluoAnalysisManager* analysis =          << 109   XrayFluoPlaneDetectorConstruction* planeDetector = 0;
 78     XrayFluoAnalysisManager::getInstance();    << 110   XrayFluoMercuryDetectorConstruction* mercuryDetector = 0;
                                                   >> 111 
 79                                                   112 
 80   // set mandatory initialization              << 
 81                                                   113 
 82   //Initialize geometry                        << 
 83   if (geometryNumber == 1 || geometryNumber ==    114   if (geometryNumber == 1 || geometryNumber == 4) {
 84     XrayFluoDetectorConstruction* testBeamDete << 115     testBeamDetector = XrayFluoDetectorConstruction::GetInstance();
 85       = XrayFluoDetectorConstruction::GetInsta << 
 86     if (geometryNumber == 4) {                    116     if (geometryNumber == 4) {
 87       testBeamDetector->PhaseSpaceOn();           117       testBeamDetector->PhaseSpaceOn();
 88       analysis->PhaseSpaceOn();                << 
 89     }                                             118     }
 90     runManager->SetUserInitialization(testBeam    119     runManager->SetUserInitialization(testBeamDetector);
 91   }                                               120   }
 92   else if (geometryNumber == 2) {                 121   else if (geometryNumber == 2) {
 93     XrayFluoPlaneDetectorConstruction*         << 122     planeDetector = XrayFluoPlaneDetectorConstruction::GetInstance();
 94       planeDetector = XrayFluoPlaneDetectorCon << 
 95     runManager->SetUserInitialization(planeDet    123     runManager->SetUserInitialization(planeDetector);
 96   }                                               124   }
 97   else if (geometryNumber == 3) {                 125   else if (geometryNumber == 3) {
 98     XrayFluoMercuryDetectorConstruction* mercu << 126     mercuryDetector = XrayFluoMercuryDetectorConstruction::GetInstance();
 99       XrayFluoMercuryDetectorConstruction::Get << 
100     runManager->SetUserInitialization(mercuryD    127     runManager->SetUserInitialization(mercuryDetector);
101   }                                               128   }
102                                                   129 
103                                                   130 
104   //Initialize physics                         << 131   runManager->SetUserInitialization(xrayList);
105   runManager->SetUserInitialization(new   Xray << 132   
106                                                << 133 #ifdef G4VIS_USE
107   //Initialize actions                         << 
108   runManager->SetUserInitialization            << 
109     (new XrayFluoActionInitializer(geometryNum << 
110                                                << 
111   //visualization manager                         134   //visualization manager
112   G4VisManager* visManager = new G4VisExecutiv    135   G4VisManager* visManager = new G4VisExecutive;
113   visManager->Initialize();                       136   visManager->Initialize();
                                                   >> 137 #endif
                                                   >> 138 
                                                   >> 139 #ifdef G4ANALYSIS_USE
                                                   >> 140   // set analysis to have the messenger running...
                                                   >> 141   XrayFluoAnalysisManager* analysis = XrayFluoAnalysisManager::getInstance();
                                                   >> 142 #endif
                                                   >> 143   XrayFluoEventAction* eventAction = 0;
                                                   >> 144   XrayFluoRunAction* runAction = new XrayFluoRunAction();
                                                   >> 145   XrayFluoSteppingAction* stepAction = new XrayFluoSteppingAction();
                                                   >> 146 
114                                                   147 
115   // get the pointer to the User Interface man << 148   //Selecting the PrimaryGenerator depending upon Geometry setup selected
                                                   >> 149 
                                                   >> 150   if (geometryNumber == 1 || geometryNumber == 4) {
                                                   >> 151     if (geometryNumber == 4) {
                                                   >> 152 #ifdef G4ANALYSIS_USE
                                                   >> 153      analysis->PhaseSpaceOn();
                                                   >> 154      analysis->CreatePersistency();
                                                   >> 155 #endif
                                                   >> 156     }
                                                   >> 157     eventAction = new XrayFluoEventAction(testBeamDetector);
                                                   >> 158     runManager->SetUserAction(new XrayFluoPrimaryGeneratorAction(testBeamDetector));
                                                   >> 159   }
                                                   >> 160 
                                                   >> 161   else if (geometryNumber == 2) {
                                                   >> 162     eventAction = new XrayFluoEventAction(planeDetector);
                                                   >> 163     runManager->SetUserAction(new XrayFluoPlanePrimaryGeneratorAction(planeDetector));
                                                   >> 164   }
                                                   >> 165 
                                                   >> 166   else if (geometryNumber == 3) {
                                                   >> 167     stepAction->SetMercuryFlag(true);
                                                   >> 168     eventAction = new XrayFluoEventAction(mercuryDetector);
                                                   >> 169     runManager->SetUserAction(new XrayFluoMercuryPrimaryGeneratorAction(mercuryDetector));
                                                   >> 170   }
                                                   >> 171  
                                                   >> 172   runManager->SetUserAction(eventAction); 
                                                   >> 173   runManager->SetUserAction(runAction);
                                                   >> 174   runManager->SetUserAction(stepAction);
                                                   >> 175   
                                                   >> 176   //Initialize G4 kernel
                                                   >> 177   //runManager->Initialize();
                                                   >> 178   
                                                   >> 179   // get the pointer to the User Interface manager 
116   G4UImanager* UImanager = G4UImanager::GetUIp    180   G4UImanager* UImanager = G4UImanager::GetUIpointer();
117                                                   181 
                                                   >> 182   if (getenv("G4VIS_USE")) {
                                                   >> 183     UImanager->ApplyCommand("/control/execute vis.mac");
                                                   >> 184   }     
                                                   >> 185 
118   if (argc == 1)   // Define UI session for in    186   if (argc == 1)   // Define UI session for interactive mode.
119     {                                             187     {
120       UImanager->ApplyCommand("/control/execut << 188       UImanager->ApplyCommand("/control/execute initInter.mac");     
121       UImanager->ApplyCommand("/control/execut << 189 #ifdef G4UI_USE
122       G4UIExecutive* ui = new G4UIExecutive(ar    190       G4UIExecutive* ui = new G4UIExecutive(argc, argv);
                                                   >> 191 
                                                   >> 192       if (ui->IsGUI())
                                                   >> 193         UImanager->ApplyCommand("/control/execute gui.mac");     
123       ui->SessionStart();                         194       ui->SessionStart();
124       delete ui;                                  195       delete ui;
                                                   >> 196 #endif
125     }                                             197     }
126   else           // Batch mode                    198   else           // Batch mode
127     {                                          << 199     { 
128       G4String command = "/control/execute ";     200       G4String command = "/control/execute ";
129       G4String fileName = argv[1];                201       G4String fileName = argv[1];
130       UImanager->ApplyCommand(command+fileName    202       UImanager->ApplyCommand(command+fileName);
131     }                                             203     }
132                                                << 204   
133   // job termination                              205   // job termination
                                                   >> 206   
                                                   >> 207 #ifdef G4VIS_USE
134   delete visManager;                              208   delete visManager;
135   G4cout << "visManager deleted"<< G4endl;        209   G4cout << "visManager deleted"<< G4endl;
                                                   >> 210 #endif
136                                                   211 
                                                   >> 212 //   if (testBeamDetector) delete testBeamDetector;
                                                   >> 213 //   if (planeDetector) delete planeDetector;
                                                   >> 214 //   if (mercuryDetector) delete mercuryDetector;
                                                   >> 215 
                                                   >> 216   
137   delete runManager;                              217   delete runManager;
                                                   >> 218   
                                                   >> 219    
                                                   >> 220 
138 }                                                 221 }
139                                                   222