#include <TH2D.h>
#include <TFile.h>
#include "myHist.C"
#include <string>
void CompareDataMC()
{
  int ymin(46), ymax(60) ;
  double xmin(0.75) , xmax(1.25) ;
  TH2D* dp = (TH2D*) GetFromFile( "PtVsEtaGoodPlus", "OutputHistFileDataNoDofCut.root")  ;
  TH2D* dm = (TH2D*) GetFromFile( "PtVsEtaGoodMinus", "OutputHistFileDataNoDofCut.root")  ;
  TH1D* drat = dynamic_cast <TH1D*> ((dp->ProjectionY())->Clone("drat")) ;
  //  TH1D* drat = dynamic_cast <TH1D*> ((dp->ProjectionX("a",ymin,ymax))->Clone("drat")) ;

  zone(1,2) ;
  drat->Sumw2() ;
  drat->Divide( dm->ProjectionY() ) ;
  //  drat->Divide( dm->ProjectionX("b",ymin,ymax) ) ;
  drat->SetLineColor(36) ;
  drat->SetMaximum(xmax) ;
  drat->SetMinimum(xmin) ;
  drat->SetTitle("Eta, Plus Over Minus" );
  drat->SetStats(0) ;
  drat->DrawCopy() ;

  TH2D* mp = (TH2D*) GetFromFile( "PtVsEtaGoodPlus", "OutputHistFileMCNoDofCut.root")  ;
  TH2D* mm = (TH2D*) GetFromFile( "PtVsEtaGoodMinus", "OutputHistFileMCNoDofCut.root")  ;
  TH1D* mrat = dynamic_cast <TH1D*> ((mp->ProjectionY())->Clone("drat")) ;
  //  TH1D* mrat = dynamic_cast <TH1D*> ((mp->ProjectionX("c",ymin,ymax))->Clone("drat")) ;
  mrat->Sumw2() ;
  mrat->Divide( mm->ProjectionY() ) ;
  //  mrat->Divide( mm->ProjectionX("d",ymin,ymax) ) ;
  mrat->SetLineColor(46) ;
  mrat->SetTitle("Plus Over Minus" );
  mrat->DrawCopy("same") ;

  drat->SetTitle("Eta, Data Over MC" );
  drat->Divide( mrat ) ;
  NewC->cd(2) ;
  drat->SetMaximum(0.9*xmax) ;
  drat->SetMinimum(1.1*xmin) ;
  drat->SetStats(1) ;
  drat->Draw() ;
  drat->Fit("pol0") ;
}