Changeset 1313 for trunk/examples/advanced/hadrontherapy/include/HadrontherapyInteractionParameters.hh

Timestamp:

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

Author:

garnier

Message:

geant4.9.4 beta rc0

File:

• trunk/examples/advanced/hadrontherapy/include/HadrontherapyInteractionParameters.hh (modified) (2 diffs)

trunk/examples/advanced/hadrontherapy/include/HadrontherapyInteractionParameters.hh

@@ -34 +34 @@
 #include "G4NistElementBuilder.hh"
 
+#ifdef G4ANALYSIS_USE_ROOT 
+#include "TROOT.h"
+#include "TCanvas.h"
+#include "TFile.h"
+#include "TH1F.h" 
+#include "TH2F.h"
+#include "TGraph.h"
+#include "TLegend.h"
+#include "TLegendEntry.h"
+#include "TStyle.h"
+#endif
+
 class HadrontherapyDetectorConstruction;
 class HadrontherapyParameterMessenger; 
+class G4ParticleDefinition;
+class G4Material;
+
 class HadrontherapyInteractionParameters : public G4EmCalculator 
 {
 public:
 
-    HadrontherapyInteractionParameters();
-        ~HadrontherapyInteractionParameters();
+    HadrontherapyInteractionParameters(G4bool);
+    ~HadrontherapyInteractionParameters();
 
-// Get data for Mass SP (MeV*cm2/g)   
+// Get data for Mass SP    
 // G4NistMaterialBuilder class materials
 // User must provide: material kinetic energy lower limit, kinetic energy upper limit, number of points to retrieve,
 // [particle], [output filename].
 
-    bool GetStoppingTable (const G4String& vararg);
+    G4bool GetStoppingTable (const G4String& vararg);
+    G4double GetStopping (G4double energy,
+                          const G4ParticleDefinition*, 
+                          const G4Material*, 
+                          G4double density = 0.);
+#ifdef G4ANALYSIS_USE_ROOT 
+    void PlotStopping(const G4String&);
+#endif
     void ListOfNistMaterials (const G4String& vararg);
     void BeamOn();
@@ -55 +77 @@
 private:
     G4Material* GetNistMaterial(G4String material);
-
     G4NistElementBuilder* nistEle;
     G4NistMaterialBuilder* nistMat;
-    G4double kinEmin, kinEmax, npoints;
-    G4String particle, material, filename; 
     std::ofstream outfile;
     std::ostream data;
     G4Material* Pmaterial;
-    G4double density;
-    G4EmCalculator* emCal;
     HadrontherapyParameterMessenger* pMessenger; 
     bool beamFlag;
+
+#ifdef G4ANALYSIS_USE_ROOT 
+    TCanvas *theRootCanvas;
+    TGraph *theRootGraph;
+    TAxis *axisX, *axisY;
+#endif
+    G4double kinEmin, kinEmax, npoints;
+    G4String particle, material, filename; 
+    G4double dedxtot, density;
+    std::vector<G4double> energy;
+    std::vector<G4double> massDedx;
+
 };
 #endif