{
  gROOT->Reset();
  
#include <fstream.h>
  
  ifstream in1;
  in1.open("energydeposit.dat");
   
  Float_t x1,x2,x3,x4,x5;
  Float_t nlines = 0;
  Float_t etotal = 0;
  char id[10];
   
  TH1F *h1 = new TH1F("h1","",100,0,100);
  h1->SetLineColor(4);
  h1->SetLineWidth(2);
  h1->SetXTitle("Energy (MeV)");
  h1->SetYTitle("Counts");
                                                                                
 TH2F *h2 = new TH2F("h2","",100,-100,100, 100, 0.,10);
  h2->SetLineColor(4);
  h2->SetLineWidth(9);
  h2->SetXTitle("X (mm)");
  h2->SetYTitle("Energy Deposit (MeV)");

  TH2F *h3 = new TH2F("h3","",100,-100,100, 100, 0,10);
  h3->SetLineColor(25);
  h3->SetLineWidth(9);
  h3->SetXTitle("Y (mm)");
  h3->SetYTitle("Energy Deposit (MeV)");
  
  TH2F *h4 = new TH2F("h4","",100,-100,100, 100, 0.,5);
  h4->SetLineColor(15);
  h4->SetLineWidth(9);
  h4->SetXTitle("Z (mm)");
  h4->SetYTitle("Energy Deposit (MeV)");
 
  TH3F *h5 = new TH3F("h5","",100,-100,100,100,-100,100,100,0.,5);
  h5->SetLineColor(5);
  h5->SetLineWidth(9);
  h5->SetXTitle("X (mm)");
  h5->SetYTitle("Y (mm)");
  h5->SetZTitle("Z (mm)");
  
  TH2F *h6 = new TH2F("h6","",100,-100,100, 100, -100,100);
  h6->SetLineColor(15);
  h6->SetLineWidth(9);
  h6->SetXTitle("X (mm)");
  h6->SetYTitle("Y (mm)");                                                                          
 while (1)
  {
    in1 >> x1 >> x2 >> x3 >> x4 >> x5 >> id;
    if (!in1.good()) break;
   //cout << " x5 = " << x5  << "  x4 = " << x4 << endl;
                                                     
    if (x4 != 9999) etotal += x4;
    //cout << "Total E value = " << etotal << endl;
    if ( (x1 == 9999) && (x2 == 9999) && (x3 == 9999) && (x4 == 9999) && (x5 ==
9999))
   {
   //cout << etotal << endl;
    h1->Fill(etotal);
    etotal = 0;
   }
    h2->Fill(x1, x4);
    h3->Fill(x2, x4);
    h4->Fill(x3, x4);
    h5->Fill(x1,x2,x4);
    h6->Fill(x1, x2);


  }
                                                                                
  in1.close();
                                                                                
   TCanvas* c1 = new TCanvas("c1","Energy",200,10,400,300);
   h1->Draw("lego1");
   c1->Modified();
   c1->Update();
                                                                                
   TCanvas* c2 = new TCanvas("c2","EnergyLoss",200,10,400,300);
   h2->Draw("scat");
   //h2->Draw("surf,same");   
   c2->Modified();
   c2->Update();
  
   TCanvas* c3 = new TCanvas("c3","EnergyLoss",200,10,400,300);
   h3->Draw("scat");
   //h3->Draw("same");
   c3->Modified();
   c3->Update();
   
   TCanvas* c4 = new TCanvas("c4","EnergyLoss",200,10,400,300);
   h4->Draw("surf2z");
   c4->Modified();
   c4->Update();

   TCanvas* c5 = new TCanvas("c5","EnergyLoss",200,10,400,300);
   h5->Draw("scat");
   c5->Modified();
   c5->Update();
   
   TCanvas* c6 = new TCanvas("c6","EnergyLoss",200,10,400,300);
   h6->Draw("surf2z");
   c6->Modified();
   c6->Update();

}