// code for plotting Optical Beacon pulse widths graphs

// Marieke Navin - September 2004

{
  gROOT->Reset();                       // reset global variables
  
  gStyle->SetCanvasColor(10);
  gStyle->SetFrameFillColor(10);
  gStyle->SetOptStat(0);
  
  ifstream data;                        // declare a file to input from
  data.open("width_6face.dat");   // open different files containing values
 

//const Int_t n = 25;                        // number of points to be plotted
const Int_t n = 6; 
  //careful with no of points!

 Float_t x[n], y1[n], y2[n];

  
  for(Int_t i=0;i<n;i++){
    data>>x[i]>>y1[i]>>y2[i];
  }

  for (Int_t i=0; i<n; i++)
    {
      cout << x[i] << " " << y1[i] <<" " << y2[i] << endl;
    }
  
  TGraph* inputSpec = new TGraph(n, x, y1); //change for title as necessary
  TGraph* inputSpec2 = new TGraph(n, x, y2);

 
inputSpec->SetTitle("Amplitude measured on internal and external PMTs (6 Faces 15/465)");
  inputSpec->GetXaxis()->SetTitle("Angle (degrees)");
  inputSpec->GetYaxis()->SetTitle("Amplitude (mV)");

  inputSpec->SetName("6face");
  inputSpec->Draw("AC");
  
  inputSpec2->SetLineStyle(kDashDotted);
  inputSpec2->Draw("CP");
}