Geant4 Cross Reference

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


  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 // Hadrontherapy advanced example for Geant4
 27 // See more at: https://twiki.cern.ch/twiki/bi     27 // See more at: https://twiki.cern.ch/twiki/bin/view/Geant4/AdvancedExamplesHadrontherapy
 28                                                    28 
 29 #include "HadrontherapyAnalysisFileMessenger.h     29 #include "HadrontherapyAnalysisFileMessenger.hh"
                                                   >>  30 #include "HadrontherapyAnalysisManager.hh"
 30 #include "G4UIcmdWithAString.hh"                   31 #include "G4UIcmdWithAString.hh"
 31 #include "G4UIcmdWithABool.hh"                     32 #include "G4UIcmdWithABool.hh"
 32 #include "G4UIdirectory.hh"                        33 #include "G4UIdirectory.hh"
 33 #include "G4SystemOfUnits.hh"                      34 #include "G4SystemOfUnits.hh"
 34                                                    35 
 35 #include "HadrontherapyMatrix.hh"                  36 #include "HadrontherapyMatrix.hh"
 36 #include "HadrontherapyLet.hh"                     37 #include "HadrontherapyLet.hh"
 37                                                    38 
 38 //////////////////////////////////////////////     39 /////////////////////////////////////////////////////////////////////////////
 39 HadrontherapyAnalysisFileMessenger::Hadronther <<  40 HadrontherapyAnalysisFileMessenger::HadrontherapyAnalysisFileMessenger(HadrontherapyAnalysisManager* amgr)
 40 :AnalysisManager(amgr)                             41 :AnalysisManager(amgr)
 41 {                                              <<  42 {  
 42     secondaryCmd = new G4UIcmdWithABool("/anal <<  43   secondaryCmd = new G4UIcmdWithABool("/analysis/secondary",this);
 43     secondaryCmd -> SetParameterName("secondar <<  44   secondaryCmd -> SetParameterName("secondary", true);
 44     secondaryCmd -> SetDefaultValue("true");   <<  45   secondaryCmd -> SetDefaultValue("true");
 45     secondaryCmd -> SetGuidance("Set if dose/f <<  46   secondaryCmd -> SetGuidance("Set if dose/fluence for the secondary particles will be written" 
 46                                 "\n[usage]: /a <<  47         "\n[usage]: /analysis/secondary [true/false]"); 
 47     secondaryCmd -> AvailableForStates(G4State <<  48   secondaryCmd -> AvailableForStates(G4State_Idle, G4State_PreInit);
 48                                                <<  49 
 49     DoseMatrixCmd = new G4UIcmdWithAString("/a <<  50   DoseMatrixCmd = new G4UIcmdWithAString("/analysis/writeDoseFile",this);
 50     DoseMatrixCmd->SetGuidance("Write the dose <<  51   DoseMatrixCmd->SetGuidance("Write the dose/fluence to an ASCII file");
 51     DoseMatrixCmd->SetDefaultValue("Dose.out") <<  52   DoseMatrixCmd->SetDefaultValue("Dose.out");
 52     DoseMatrixCmd->SetParameterName("choice",t <<  53   DoseMatrixCmd->SetParameterName("choice",true); 
 53                                                <<  54 
 54     LetCmd = new G4UIcmdWithABool("/analysis/c <<  55   // With this messenger you can:
 55     LetCmd  -> SetParameterName("choice",true) <<  56   // give a name to the generated .root file
 56     LetCmd  -> SetDefaultValue(true);          <<  57   // One can use this messenger to define a different .root file name other then the default one 
 57     LetCmd  -> SetGuidance("Set if Let must be <<  58 #ifdef G4ANALYSIS_USE_ROOT
 58     LetCmd  -> AvailableForStates(G4State_Idle <<  59   FileNameCmd = new G4UIcmdWithAString("/analysis/setAnalysisFile",this);
 59                                                <<  60   FileNameCmd->SetGuidance("Set the .root filename for the root-output");
                                                   >>  61   FileNameCmd->SetDefaultValue("default.root");
                                                   >>  62   FileNameCmd->SetParameterName("choice",true); ///<doc did not say what second boolean really does
                                                   >>  63   FileNameCmd->AvailableForStates(G4State_Idle,G4State_PreInit);
                                                   >>  64 #endif
                                                   >>  65 
                                                   >>  66 
                                                   >>  67 LetCmd = new G4UIcmdWithABool("/analysis/computeLet",this);
                                                   >>  68   LetCmd  -> SetParameterName("choice",true); 
                                                   >>  69   LetCmd  -> SetDefaultValue(true);
                                                   >>  70   LetCmd  -> SetGuidance("Set if Let must be computed and write the ASCII filename for the Let");
                                                   >>  71   LetCmd  -> AvailableForStates(G4State_Idle, G4State_PreInit);
                                                   >>  72 
                                                   >>  73 
                                                   >>  74 
 60 }                                                  75 }
 61                                                    76 
 62 //////////////////////////////////////////////     77 /////////////////////////////////////////////////////////////////////////////
 63 HadrontherapyAnalysisFileMessenger::~Hadronthe     78 HadrontherapyAnalysisFileMessenger::~HadrontherapyAnalysisFileMessenger()
 64 {                                                  79 {
 65     delete secondaryCmd;                       <<  80   delete secondaryCmd; 
 66     delete DoseMatrixCmd;                      <<  81   delete DoseMatrixCmd; 
 67     delete LetCmd;                             <<  82   delete LetCmd;
                                                   >>  83 
                                                   >>  84 #ifdef G4ANALYSIS_USE_ROOT
                                                   >>  85   delete FileNameCmd;
                                                   >>  86 #endif
 68 }                                                  87 }
 69                                                    88 
 70 //////////////////////////////////////////////     89 /////////////////////////////////////////////////////////////////////////////
 71 void HadrontherapyAnalysisFileMessenger::SetNe     90 void HadrontherapyAnalysisFileMessenger::SetNewValue(G4UIcommand* command, G4String newValue)
 72 {                                              <<  91 { 
 73     if (command == secondaryCmd)                   92     if (command == secondaryCmd)
 74     {                                              93     {
 75         if (HadrontherapyMatrix::GetInstance() <<  94   if (HadrontherapyMatrix::GetInstance())
 76         {                                      <<  95   {
 77             HadrontherapyMatrix::GetInstance() <<  96       HadrontherapyMatrix::GetInstance() -> secondary = secondaryCmd -> GetNewBoolValue(newValue);
 78         }                                      <<  97   }
                                                   >>  98     }
                                                   >>  99 
                                                   >> 100     else if (command == DoseMatrixCmd) // Filename can be passed here TODO 
                                                   >> 101     { 
                                                   >> 102   if ( HadrontherapyMatrix * pMatrix = HadrontherapyMatrix::GetInstance() )
                                                   >> 103   {
                                                   >> 104       pMatrix -> TotalEnergyDeposit(); 
                                                   >> 105       pMatrix -> StoreDoseFluenceAscii(newValue);
                                                   >> 106 #ifdef G4ANALYSIS_USE_ROOT
                                                   >> 107       pMatrix -> StoreDoseFluenceRoot();
                                                   >> 108       HadrontherapyAnalysisManager::GetInstance() -> flush();     // Finalize & write the root file 
                                                   >> 109 #endif
                                                   >> 110   }
 79     }                                             111     }
 80                                                   112     
 81     else if (command == DoseMatrixCmd) // File << 113      else if (command == LetCmd)
 82     {                                             114     {
 83         if ( HadrontherapyMatrix * pMatrix = H << 115     if (HadrontherapyLet::GetInstance())
 84         {                                      << 116       HadrontherapyLet::GetInstance() -> doCalculation = LetCmd -> GetNewBoolValue(newValue);
 85             //pMatrix -> TotalEnergyDeposit(); << 
 86             pMatrix -> StoreDoseFluenceAscii(n << 
 87         }                                      << 
 88     }                                             117     }
 89                                                   118     
 90     else if (command == LetCmd)                << 119 #ifdef G4ANALYSIS_USE_ROOT
                                                   >> 120     else if (command == FileNameCmd)
 91     {                                             121     {
 92         if (HadrontherapyLet::GetInstance())   << 122   AnalysisManager->SetAnalysisFileName(newValue);
 93             HadrontherapyLet::GetInstance() -> << 123   HadrontherapyAnalysisManager::GetInstance() -> book(); // Book for a new ROOT TFile 
 94     }                                             124     }
                                                   >> 125 #endif
 95 }                                                 126 }
 96                                                   127 
 97                                                   128