Changeset 1313 for trunk/examples/advanced/hadrontherapy/src/HadrontherapyEventAction.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/HadrontherapyEventAction.cc (modified) (4 diffs)

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

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: HadrontherapyEventAction.cc; 
-// See more at: http://g4advancedexamples.lngs.infn.it/Examples/hadrontherapy
-
+// $Id: HadrontherapyEventAction.cc;
+//
+// See more at: http://workgroup.lngs.infn.it/geant4lns
+//
+// ----------------------------------------------------------------------------
+//                 GEANT 4 - Hadrontherapy example
+// ----------------------------------------------------------------------------
+// Code developed by:
+//
+// G.A.P. Cirrone(a)*, G. Candiano, F. Di Rosa(a), S. Guatelli(b), G. Russo(a)
+// 
+// (a) Laboratori Nazionali del Sud 
+//     of the National Institute for Nuclear Physics, Catania, Italy
+// (b) National Institute for Nuclear Physics Section of Genova, genova, Italy
+// 
+// * cirrone@lns.infn.it
+// --------------------------------------------------------------
 #include "G4Event.hh"
 #include "G4EventManager.hh"
 #include "G4HCofThisEvent.hh"
 #include "G4VHitsCollection.hh"
-#include "G4TrajectoryContainer.hh"
-#include "G4Trajectory.hh"
-#include "G4VVisManager.hh"
 #include "G4SDManager.hh"
 #include "G4VVisManager.hh"
+
 #include "HadrontherapyEventAction.hh"
 #include "HadrontherapyDetectorHit.hh"
@@ -64 +76 @@
   //printing survey
   if (evtNb%printModulo == 0)
-    G4cout << "\n---> Begin of Event: " << evtNb << G4endl;
+     G4cout << "\n---> Begin of Event: " << evtNb << G4endl;
    
   G4SDManager* pSDManager = G4SDManager::GetSDMpointer();
   if(hitsCollectionID == -1)
     hitsCollectionID = pSDManager -> GetCollectionID("HadrontherapyDetectorHitsCollection");
+  
 }
 
 /////////////////////////////////////////////////////////////////////////////
 void HadrontherapyEventAction::EndOfEventAction(const G4Event* evt)
-{  
+{ 
   if(hitsCollectionID < 0)
   return;
   G4HCofThisEvent* HCE = evt -> GetHCofThisEvent();
- 
+
+  // Clear voxels hit list 
+  HadrontherapyMatrix* matrix = HadrontherapyMatrix::GetInstance();
+  if (matrix) matrix -> ClearHitTrack(); 
+
   if(HCE)
   {
@@ -83 +100 @@
     if(CHC)
      {
-           matrix = HadrontherapyMatrix::getInstance();  
        if(matrix)
           { 
@@ -101 +117 @@
     }
   }
-  // Extract the trajectories and draw them in the visualisation
-
-  if (G4VVisManager::GetConcreteInstance())
-    {
-      G4TrajectoryContainer * trajectoryContainer = evt -> GetTrajectoryContainer();
-      G4int n_trajectories = 0;
-      if (trajectoryContainer) n_trajectories = trajectoryContainer -> entries();
-      
-      for (G4int i = 0; i < n_trajectories; i++) 
-        {
-          G4Trajectory* trj = (G4Trajectory*)
-            ((*(evt -> GetTrajectoryContainer()))[i]);
-          if(drawFlag == "all") trj -> DrawTrajectory(50);
-          else if((drawFlag == "charged")&&(trj -> GetCharge() != 0.))
-            trj -> DrawTrajectory(50);
-          else if ((drawFlag == "neutral")&&(trj -> GetCharge() == 0.))
-            trj -> DrawTrajectory(50);               
-        }
-    }
 }