{
TCanvas *c1 = new TCanvas("c1","residual studies 2", 500, 750); 
c1->Clear();
c1->Divide(2,4);
gStyle->SetOptFit(1011);
gStyle->SetOptStat(0);

c1_1->cd();
  r091->Fit("gaus");
  Double_t p091 = gaus->GetParameter(2);
c1_2->cd();
  r092->Fit("gaus");
  Double_t p092 = gaus->GetParameter(2);

c1_3->cd();
  r101->Fit("gaus");
  Double_t p101 = gaus->GetParameter(2);

c1_4->cd();
  r102->Fit("gaus");
  Double_t p102 = gaus->GetParameter(2);

c1_5->cd();
  r111->Fit("gaus");
  Double_t p111 = gaus->GetParameter(2);

c1_6->cd();
  r112->Fit("gaus");
  Double_t p112 = gaus->GetParameter(2);

c1_7->cd();
  r121->Fit("gaus");
  Double_t p121 = gaus->GetParameter(2);

c1_8->cd();
  r122->Fit("gaus");
  Double_t p122 = gaus->GetParameter(2);


ifstream infile ("cfmgmc.input");
Double_t cfm[10];
for (int ii=0; ii<10; ii++) {
  infile>>cfm[ii];
}
infile.close();

ofstream outfile ("res_out.txt", ios::app);
outfile<<" cfmgmc.input" ;
for (int ii=0; ii<10; ii++) {
  outfile<<cfm[ii]<<"  ";
}
outfile <<"\n\n layer 09: "<<p091<<" "<<p092<<endl;
outfile <<" layer 10: "<<p101<<" "<<p102<<endl;
outfile <<" layer 11: "<<p111<<" "<<p112<<endl;
outfile <<" layer 12: "<<p121<<" "<<p122<<endl;
outfile.close();

}