{
Int_t i;
ifstream fin7, fin8;
fin7.open("OmegaGammas.txt");
Int_t p=0;
while (!fin7.eof()) { string blah; fin7 >> blah; p= p+1;}
p=(p-1)/2;
fin7.close();
Double_t OmegaGamma[p], Energy[p];
fin8.open("OmegaGammas.txt");
for(i=0;i<p;i++){ fin8 >> OmegaGamma[i]; fin8 >> Energy[i];}
fin8.close();
Int_t n = p*2;
Double_t Param[n];
for(i=0;i<n;i++)
{
if( i%2 == 0) Param[i] = OmegaGamma[i/2];
if( i%2 == 1) Param[i] = Energy[i/2];
}

ofstream fout; 
fout.open("dumbfix.txt");
for(i=0;i<n;i++){ fout << i << endl;}
fout.close();
ifstream fin;
fin.open("dumbfix.txt");
string names; 
names = "0";
for(i=0;i<p;i++)
{
 names.append(" + 1.54*pow(10,11)*pow(x/pow(10,9),-3/2)*[");
 string blah;
 fin >> blah;
 names.append(blah);
 names.append("]*exp(-11.605*pow(10,9)*(40*[");
 fin >> blah;
 names.append(blah);
 names.append("]*0.091+5.069)/x)");
}
fin.close();
cout << names << endl;

TF1 *rate = new TF1("rate", names.c_str(),1*pow(10,9),6*pow(10,9));
for(i=0;i<n;i++)
{
rate->SetParameter(i,i);
}
rate->Draw();

}