Changeset 1313 for trunk/examples/advanced/hadrontherapy/src/HadrontherapyRunAction.cc

Timestamp:

Jun 14, 2010, 3:54:58 PM (14 years ago)

Author:

garnier

Message:

geant4.9.4 beta rc0

File:

• trunk/examples/advanced/hadrontherapy/src/HadrontherapyRunAction.cc (modified) (4 diffs)

trunk/examples/advanced/hadrontherapy/src/HadrontherapyRunAction.cc

@@ -27 +27 @@
 // See more at: http://g4advancedexamples.lngs.infn.it/Examples/hadrontherapy
 
+
 #include "HadrontherapyRunAction.hh"
 #include "HadrontherapyEventAction.hh"
@@ -38 +39 @@
 #include "HadrontherapyRunAction.hh"
 
+
+#include "HadrontherapyAnalysisManager.hh"
+#include "HadrontherapyMatrix.hh"
+
 HadrontherapyRunAction::HadrontherapyRunAction()
 {
@@ -48 +53 @@
 void HadrontherapyRunAction::BeginOfRunAction(const G4Run* aRun)
 {       
-   G4RunManager::GetRunManager()->SetRandomNumberStore(true);
-   G4cout << "Run " << aRun -> GetRunID() << " starts ..." << G4endl;
+    G4RunManager::GetRunManager()-> SetRandomNumberStore(true);
+    G4cout << "Run " << aRun -> GetRunID() << " starts ..." << G4endl;
 
-   electromagnetic = 0;
-   hadronic = 0;
+    // Warning! any beamOn will reset all data (dose, fluence, histograms, etc)!
+    //
+
+    // Initialize matrix with energy of primaries clearing data inside
+    if (HadrontherapyMatrix::GetInstance()) 
+    {
+        HadrontherapyMatrix::GetInstance() -> Initialize();
+    }
+
+#ifdef G4ANALYSIS_USE_ROOT
+    HadrontherapyAnalysisManager::GetInstance() -> flush();     // Finalize the root file 
+    // Initialize root analysis ----> book
+    HadrontherapyAnalysisManager::GetInstance() -> book();
+#endif
+    electromagnetic = 0;
+    hadronic = 0;
 }
 
 void HadrontherapyRunAction::EndOfRunAction(const G4Run*)
 {
-  //   G4cout << " Summary of Run " << aRun -> GetRunID() <<" :"<< G4endl;
-  //G4cout << "Number of electromagnetic processes of primary particles in the phantom:"
-  //     << electromagnetic << G4endl;
-  //G4cout << "Number of hadronic processes of primary particles in the phantom:"
-  //     << hadronic << G4endl;
+    // Store dose & fluence data to ASCII & ROOT files 
+    if ( HadrontherapyMatrix::GetInstance() )
+    {
+        HadrontherapyMatrix::GetInstance() -> TotalEnergyDeposit(); 
+        HadrontherapyMatrix::GetInstance() -> StoreDoseFluenceAscii();
+
+#ifdef G4ANALYSIS_USE_ROOT
+        if (HadrontherapyAnalysisManager::GetInstance())
+        { 
+            HadrontherapyMatrix::GetInstance() -> StoreDoseFluenceRoot();
+        }
+#endif
+    }
+
+    //G4cout << " Summary of Run " << aRun -> GetRunID() <<" :"<< G4endl;
+    //G4cout << "Number of electromagnetic processes of primary particles in the phantom:"
+    //     << electromagnetic << G4endl;
+    //G4cout << "Number of hadronic processes of primary particles in the phantom:"
+    //     << hadronic << G4endl;
 }
 void HadrontherapyRunAction::AddEMProcess()
@@ -72 +105 @@
 }
 
-
-