Changeset 1342 for trunk/examples/advanced/microbeam/src/MicrobeamRunAction.cc

Timestamp:

Nov 5, 2010, 4:08:39 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/examples/advanced/microbeam/src/MicrobeamRunAction.cc (modified) (5 diffs)

trunk/examples/advanced/microbeam/src/MicrobeamRunAction.cc

@@ -25 +25 @@
 //
 // -------------------------------------------------------------------
-// $Id: MicrobeamRunAction.cc,v 1.6 2006/06/29 16:05:37 gunter Exp $
+// $Id: MicrobeamRunAction.cc,v 1.7 2010/10/07 14:03:11 sincerti Exp $
 // -------------------------------------------------------------------
 
@@ -36 +36 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 
-MicrobeamRunAction::MicrobeamRunAction(MicrobeamDetectorConstruction* det)
-:Detector(det)
+MicrobeamRunAction::MicrobeamRunAction(MicrobeamDetectorConstruction* det,
+MicrobeamHistoManager * his)
+:Detector(det),Histo(his)
 {   
   saveRndm = 0;  
@@ -55 +56 @@
 {  
  
+  // Histograms
+  Histo->book();
+
   // save Rndm status
   if (saveRndm > 0)
@@ -99 +103 @@
   }   
   
-  FILE *myFile;
-  myFile=fopen("3DDose.txt","a");
   for (G4int i=0; i<nbOfPixels; i++) 
   {  
@@ -106 +108 @@
     v = mapVoxels[i];
     if ( (GetNumEvent()+1) !=0) 
-      fprintf (myFile,"%f %f %f %f \n", v.x(), v.y(), v.z(), dose3DDose[i]/(GetNumEvent()+1));
+      {
+          Histo->FillNtuple(4,0,v.x());
+          Histo->FillNtuple(4,1,v.y());
+          Histo->FillNtuple(4,2,v.z());
+          Histo->FillNtuple(4,3,dose3DDose[i]/(GetNumEvent()+1));
+          Histo->AddRowNtuple(4);                            
+      }
   }
-  fclose (myFile);                              
-    
+   
   G4cout << "-> Total number of particles detected by the gas detector : " << GetNbOfHitsGas() << G4endl;  
   G4cout << G4endl;    
+  
+  //save histograms      
+  Histo->save();
+
 }