void fscaltest()
{
  gStyle->SetOptStat(kFALSE);   

   const Int_t NN =25;
 
  ifstream infile, infile1, infile2, infile3, infile4, infile5, infile6;
  infile.open("FqeStdError.dat",std::ifstream::in);
  infile1.open("syswo32.dat",std::ifstream::in);
  
  vector<double> M,eM, Fq2e, Fq3e, Fq4e,Fq5e, eFq2e, eFq3e, eFq4e, eFq5e;
  vector<double> Fq2e1, Fq3e1, Fq4e1,Fq5e1, Fq6e1, Fq7e1, Fq8e1;
  Double_t m, avbin, fq2e, efq2e, fq3e, efq3e, fq4e, efq4e, fq5e, efq5e, fq6e, efq6e, fq7e, efq7e, fq8e, efq8e;
    
  while(infile >> m >> avbin >> fq2e >> efq2e >> fq3e >> efq3e >> fq4e >> efq4e >> fq5e >> efq5e>> fq6e >> efq6e >> fq7e >> efq7e){
    M.push_back(log(m*m));
    Fq2e.push_back(log(fq2e));
    Fq3e.push_back(log(fq3e));
    Fq4e.push_back(log(fq4e));
    Fq5e.push_back(log(fq5e));
   
   
    eM.push_back(0.0);
    eFq2e.push_back( efq2e );
    eFq3e.push_back( efq3e );
    eFq4e.push_back( efq4e );
    eFq5e.push_back( efq5e );

  }

 while(infile1 >> m >> fq2e >> fq3e >> fq4e >> fq5e ){
    Fq2e1.push_back((fq2e));
    Fq3e1.push_back((fq3e));
    Fq4e1.push_back((fq4e));
    Fq5e1.push_back((fq5e));

  }
   Int_t n = M.size()-NN;   
   cout << "size of the array="<<n<<endl;
 
   TCanvas *c2   = new TCanvas("c2","Factorial Moments",20,10,500,500);
   
   TGraphMultiErrors* gme1 = new TGraphMultiErrors("gme1", "", n, &Fq2e[NN], &Fq3e[NN],&eFq2e[NN],&eFq2e[NN],&eFq3e[NN], &eFq3e[NN]);
   gme1->AddYError(n, &Fq2e1[NN], &Fq2e1[NN]);
   gme1->SetMarkerStyle(8);
   gme1->SetMarkerColor(4);
   gme1->SetLineColor(4);
   gme1->GetAttLine(0)->SetLineColor(4);
   gme1->GetAttLine(1)->SetLineColor(kGray);
   gme1->GetAttFill(1)->SetFillStyle(0);
 
   gme1->Draw("a p s ; ; 5 s=0.2");
     
     TF1 *f2 = new TF1("f2","pol1",0.045,.075);
     gme1->Fit("f2","R"); 
        
    }