diff --git a/NPLib/Physics/NPInelasticBreakup.cxx b/NPLib/Physics/NPInelasticBreakup.cxx
index 30ca3bae55e9a760d289b0e04cbd0b218ece09b1..21e2636eaa5a4b4c71ad3a49741fa7c4f1b1f677 100644
--- a/NPLib/Physics/NPInelasticBreakup.cxx
+++ b/NPLib/Physics/NPInelasticBreakup.cxx
@@ -212,7 +212,7 @@ void InelasticBreakup::ReadConfigurationFile(NPL::InputParser parser) {
 
     if (blocks[i]->HasToken("CrossSectionPath")) {
       vector<string> file = blocks[i]->GetVectorString("CrossSectionPath");
-      TH1F* CStemp = Read1DProfile(file[0], file[1]);
+      TH1D* CStemp = Read1DProfile(file[0], file[1]);
 
       // multiply CStemp by sin(theta)
       TF1* fsin = new TF1("sin", Form("1/(sin(x*%f/180.))", M_PI), 0, 180);
@@ -225,7 +225,7 @@ void InelasticBreakup::ReadConfigurationFile(NPL::InputParser parser) {
       fLabCrossSection = true;
 
       vector<string> file = blocks[i]->GetVectorString("LabCrossSectionPath");
-      TH1F* CStemp = Read1DProfile(file[0], file[1]);
+      TH1D* CStemp = Read1DProfile(file[0], file[1]);
 
       // multiply CStemp by sin(theta)
       TF1* fsin = new TF1("sin", Form("1/(sin(x*%f/180.))", M_PI), 0, 180);
diff --git a/NPLib/Physics/NPInelasticBreakup.h b/NPLib/Physics/NPInelasticBreakup.h
index 82f89e480b6d38e41ebcbce2f58cbe5e2beae2ef..0a77961d864e0c8cf3dd95670c869eb400e4e380 100644
--- a/NPLib/Physics/NPInelasticBreakup.h
+++ b/NPLib/Physics/NPInelasticBreakup.h
@@ -35,7 +35,7 @@ using namespace NPL;
 #include "NPReaction.h"
 #include "TCanvas.h"
 #include "TGraph.h"
-#include "TH1F.h"
+#include "TH1.h"
 #include "TLorentzRotation.h"
 #include "TLorentzVector.h"
 #include "TRandom.h"
@@ -81,9 +81,9 @@ namespace NPL {
     double fQValue;                            // Q-value in MeV
     double fBeamEnergy;                        // Beam energy in MeV
     double fThetaCM;                           // Center-of-mass angle in radian
-    TH1F* fCrossSectionHist;                   // Differential cross section in CM frame
+    TH1D* fCrossSectionHist;                   // Differential cross section in CM frame
     TH2F* fDoubleDifferentialCrossSectionHist; // Diff. CS CM frame vs Beam E
-    TH1F* fExcitationEnergyHist;               // Distribution of Excitation energy
+    TH1D* fExcitationEnergyHist;               // Distribution of Excitation energy
    public:
     // Getters and Setters
     void SetBeamEnergy(const double& eBeam) { fBeamEnergy = eBeam; }
@@ -92,7 +92,7 @@ namespace NPL {
     void SetExcitationLight(const double& exci) { fExcitationLight = exci; }
     void SetExcitationHeavy(const double& exci) { fExcitationHeavy = exci; }
     void SetVerboseLevel(const int& verbose) { fVerboseLevel = verbose; }
-    void SetCrossSectionHist(TH1F* CrossSectionHist) {
+    void SetCrossSectionHist(TH1D* CrossSectionHist) {
       delete fCrossSectionHist;
       fCrossSectionHist = CrossSectionHist;
     }
@@ -113,7 +113,7 @@ namespace NPL {
      Particle* GetNucleus3() { return GetParticle3(); }
      Particle* GetNucleus4() { return GetParticle4(); }
  */
-    TH1F* GetCrossSectionHist() const { return fCrossSectionHist; }
+    TH1D* GetCrossSectionHist() const { return fCrossSectionHist; }
     int GetVerboseLevel() const { return fVerboseLevel; }
 
    public: