{
Int_t i,j,k;
ifstream fin, fin2, fin3, fin4;
fin.open("FinalMasterListNOCSB.txt");
Int_t n=0;
while (!fin.eof())
{
  string blah;
  fin >> blah;
  n= n+1;
}
n= n/10;
fin.close();

Double_t CS[n], Ratio[n], eRatio[n], Pressure[n];

string Energy;
fin2.open("FinalMasterListNOCSB.txt");
for (i= 0; i<n; i++)
{ 
  string blah;
  fin2 >> blah;
  fin2 >> Energy;
  fin2 >> CS[i];
  fin2 >> blah;
  fin2 >> blah;
  fin2 >> blah;
  fin2 >> blah;
  fin2 >> Ratio[i];
  fin2 >> eRatio[i];
  fin2 >> Pressure[i];
}
fin2.close();


Int_t CS9n=0;
for(i=0;i<n;i++)
{
  if ( CS[i] == 9 ) {CS9n = CS9n+1;} 
}
Double_t CS9Pressure[CS9n], CS9ePressure[CS9n], CS9Ratio[CS9n], CS9eRatio[CS9n];
CS9n=0;
for(i=0;i<n;i++)
{
  if ( CS[i] == 9 ) 
  {CS9Pressure[CS9n]=Pressure[i];CS9ePressure[CS9n]=0;CS9Ratio[CS9n]=Ratio[i];CS9eRatio[CS9n]=eRatio[i];CS9n=CS9n+1;} 
}
Double_t *X9T = CS9Pressure;
Double_t *eX9T = CS9ePressure;
Double_t *Y9T = CS9Ratio;
Double_t *eY9T = CS9eRatio;
TGraphErrors *graph9 = new TGraphErrors(CS9n,X9T,Y9T,eX9T,eY9T);

Int_t CS10n=0;
for(i=0;i<n;i++)
{
  if ( CS[i] == 10 ) {CS10n = CS10n+1;} 
}
Double_t CS10Pressure[CS10n], CS10ePressure[CS10n], CS10Ratio[CS10n], CS10eRatio[CS10n];
CS10n=0;
for(i=0;i<n;i++)
{
  if ( CS[i] == 10 ) 
  {CS10Pressure[CS10n]=Pressure[i];CS10ePressure[CS10n]=0;CS10Ratio[CS10n]=Ratio[i];CS10eRatio[CS10n]=eRatio[i];
	  CS10n=CS10n+1;} 
}
Double_t *X10T = CS10Pressure;
Double_t *eX10T = CS10ePressure;
Double_t *Y10T = CS10Ratio;
Double_t *eY10T = CS10eRatio;
TGraphErrors *graph10 = new TGraphErrors(CS10n,X10T,Y10T,eX10T,eY10T);

Int_t CS11n=0;
for(i=0;i<n;i++)
{
  if ( CS[i] == 11 ) {CS11n = CS11n+1;} 
}
Double_t CS11Pressure[CS11n], CS11ePressure[CS11n], CS11Ratio[CS11n], CS11eRatio[CS11n];
CS11n=0;
for(i=0;i<n;i++)
{
  if ( CS[i] == 11 ) 
  {CS11Pressure[CS11n]=Pressure[i];CS11ePressure[CS11n]=0;CS11Ratio[CS11n]=Ratio[i];CS11eRatio[CS11n]=eRatio[i];
	  CS11n=CS11n+1;} 
}
Double_t *X11T = CS11Pressure;
Double_t *eX11T = CS11ePressure;
Double_t *Y11T = CS11Ratio;
Double_t *eY11T = CS11eRatio;
TGraphErrors *graph11 = new TGraphErrors(CS11n,X11T,Y11T,eX11T,eY11T);

string names[100];
names.append("Charge State Fractions of ^{48}Ti (CSB OUT) at ");
names.append(Energy);
names.append(" keV/u vs Pressure");

graph9->SetTitle(names.c_str());
graph9->GetXaxis()->SetTitle("Pressure in Torr");
graph9->GetXaxis()->CenterTitle();
graph9->GetYaxis()->SetTitle("FCCH/FC4");
graph9->GetYaxis()->SetTitleOffset(1.2);
graph9->GetYaxis()->CenterTitle();
graph9->SetMarkerStyle(21);
graph9->SetMarkerSize(1);
graph9->SetLineColor(2); // 2,4,3
graph9->SetMarkerColor(2);
graph9->Draw("ALP");

graph10->SetMarkerStyle(21);
graph10->SetMarkerSize(1);
graph10->SetLineColor(7); // 2,4,3
graph10->SetMarkerColor(7);
graph10->Draw("LP");

graph11->SetLineColor(1); // 2,4,3
graph11->SetMarkerStyle(21);
graph11->SetMarkerSize(1);
graph11->SetMarkerColor(1);
graph11->Draw("LP");



leg = new TLegend(0.65,0.65,0.89,0.89);
leg->SetHeader("Charge State");
leg->AddEntry(graph9,"9 +","l");
leg->AddEntry(graph10,"10 +","l");
leg->AddEntry(graph11,"11 +","l");
                                                                                
leg->Draw();
}