{
gROOT->Reset();
gStyle->SetFillColor(0);
gStyle->SetCanvasColor(0);
gStyle->SetPadColor(0);
gStyle->SetTitleColor(0);
gStyle->SetStatColor(0);
gStyle->SetTextColor(1);
gStyle->SetPaperSize(21,28);
gStyle->SetOptTitle(0);
gStyle->SetMarkerColor(2);
gStyle->SetMarkerSize(0.7);
gStyle->SetMarkerStyle(22);

const Int_t MAX = 100;
const Float_t Total = 1000.0;

char readFile[20];

Int_t ntr[MAX], fntr[MAX], nhad[MAX];

Float_t x[MAX],ex[MAX],y[MAX],ey[MAX],fy[MAX],fey[MAX];
for (int i=0;i<MAX;i++) {
  ntr[i]=0;
  fntr[i]=0;
  nhad[i] = 0;
  x[i] = i+1;
  ex[i] = 0.0;
  y[i] = 0.0;
  ey[i] = 0.0;
  fy[i] = 0.0;
  fey[i] = 0.0;
}

for (int i=0; i<MAX; i++) {
  sprintf(readFile,"tauhist%d.root",i+1);
  TFile *f = new TFile(readFile,"READ");
  f->ls();
  gDirectory->cd("TauAnalysis");
  
  tau->Draw("ntr");
  TH1F *htemp = (TH1F*)gPad->GetPrimitive("htemp");
  if(htemp) ntr[i] = htemp->GetEntries();
  
  tau->Draw("fntr");
  htemp = (TH1F*)gPad->GetPrimitive("htemp");
  if(htemp) fntr[i] = htemp->GetEntries();
  
  nhad[i] = tau->GetMaximum("f_nhadronic");
  
  //    cout << "# of entries in TAU  : " << ntr[i] << endl;
  //    cout << "# of entries in FTAU : " << fntr[i] << endl;
  //    cout << "# of entries in FLEPTONIC : " << nhad[i] << endl;
  //    cout << "check point 0"<< endl;
  //    f->Close();
  //    cout << "check point 1"<< endl;
}

Float_t on,n,p,q;
for(int i=0; i<MAX; i++) {
  //  y[i] = Float_t(ntr[i])/Total;
  //  ey[i] = sqrt(Float_t(ntr[i]))/Total;
  n = Total;
  p = Float_t(ntr[i])/n;
  q = 1.0 - p;
  y[i] = p;
  ey[i] = sqrt( p * q / n );
  if (nhad[i]!=0){
    //    fy[i] = Float_t(fntr[i])/Float_t(nhad[i]);
    //    fey[i] = sqrt(Float_t(fntr[i]))/Float_t(nhad[i]);
    n = Float_t(nhad[i]);
    p = Float_t(fntr[i])/n;
    q = 1.0 - p;
    fy[i] = p;
    fey[i] = sqrt( p * q /n );
  }
  else {
    fy[i] = 0.0;
    fey[i] = 0.0;
  }
}
delete c1;

TCanvas *tc1;
TPaveText *pt;
TText *text, *t1;
TPad *c1;
TPostScript *ps;
TGraphErrors *graph1,*graph2;

ps = new TPostScript("fakeTauEff.eps",113);

ps->Range(20,24);
tc1 = new TCanvas("c1", "TauFinderModule Eff. vs Pt (Chirality -1)",600,800);
pt = new TPaveText(0.20,0.92,0.80,0.98);
text = pt->AddText(tc1->GetTitle());
pt->SetBorderSize(1);
pt->SetTextColor(28);
pt->Draw();

c1 = new TPad("c1",tc1->GetTitle(),0.0,0.0,1.0,0.9);
c1->Draw();
c1->Divide(1,2);
c1->cd(1);

graph1 = new TGraphErrors(MAX,x,y,ex,ey);
//graph1->SetMarkerColor(2);
//graph1->SetMarkerSize(1);
//graph1->SetMarkerStyle(21);
graph1->Draw("AP");
t1 = new TText(40,0.2,"Tau inclusive decay");
t1->SetTextSize(0.08);
t1->SetTextColor(4);
t1->Draw();

c1->cd(2);
graph2 = new TGraphErrors(MAX,x,fy,ex,fey);
//graph2->SetMarkerColor(2);
//graph1->SetMarkerSize(1);
//graph1->SetMarkerStyle(21);
graph2->Draw("AP");
t1 = new TText(40,0.2,"Tau hadronic decay");
t1->SetTextSize(0.08);
t1->SetTextColor(4);
t1->Draw();

ps->Close();

}