source: huonglan/ANA2011/ExpectedNumber.C

{
  TFile *f  = TFile::Open("resSignal11Jul_Nobtag.root");
  TH1F *hloose  = (TH1F*)f->Get("loose");
  TH1F *hlooseE = (TH1F*)f->Get("looseE");
  TH1F *hlooseM = (TH1F*)f->Get("looseM");
  TH1F *hnjets  = (TH1F*)f->Get("njets");
  TH1F *hnjetsE = (TH1F*)f->Get("njetsE");
  TH1F *hnjetsM = (TH1F*)f->Get("njetsM");

  double  nevtIn1 = 0., nevtIn2  = 0., nevtIn3  = 0.;
  double nevtOut1 = 0., nevtOut2 = 0., nevtOut3 = 0.;
  int nbin  = hloose->GetNbinsX();
  for (int ib = 1; ib < (nbin+1); ib++){
    double bv = hloose->GetBinContent(ib);
    double bc = hloose->GetBinCenter(ib) - 0.5;
    if (bc==0)
      nevtIn1 = bv;
    if (bc==15)
      nevtOut1 = bv;
  }

  double n[4], nE[4], nM[4];
  double Nexpt1[4] ;
  double Nexpt1E[4];
  double Nexpt1M[4];
  for (int i = 0; i < 4; i++){
    n[i]  = 0;
    nE[i] = 0;
    nM[i] = 0;
    Nexpt1[i]  = 0;
    Nexpt1E[i] = 0;
    Nexpt1M[i] = 0;
  }

  
  for (int i = 1; i < 29; i++){
    int ib     = i - 1;
    double iv  = hnjets->GetBinContent(ib);
    double ivE = hnjetsE->GetBinContent(ib);
    double ivM = hnjetsM->GetBinContent(ib);

    cout << "ib = " << ib << " hnjetsE->GetBinContent(ib) is " << hnjetsE->GetBinContent(ib) << endl;

    if (ib >= 4)
      cout << "ivE = " << ivE << endl;

    if (ib == 1){
      n[0]  = iv;
      nE[0] = ivE;
      nM[0] = ivM;
    }
    if (ib == 2){
      n[1]  = iv;
      nE[1] = ivE;
      nM[1] = ivM;
    }
    if (ib == 3){
      n[2]  = iv;
      nE[2] = ivE;
      nM[2] = ivM;
    }
    if (ib >= 4){
      n[3]  += iv;
      nE[3] += ivE;
      nM[3] += ivM;
      cout << "nE[3] = " << nE[3] << endl;
    }
  }

  for (int i = 0; i < 4; i++){
    Nexpt1[i]  = 1.43*690*0.685*n[i]/nevtIn1;
    Nexpt1E[i] = 1.43*690*0.685*nE[i]/nevtIn1;
    Nexpt1M[i] = 1.43*690*0.685*nM[i]/nevtIn1;
    if (i==3)
      cout << "Nexpt1E[3] = " << Nexpt1E[3] << endl;
  }

  cout << "sans btag : " << endl;
  cout << "1 jet : " << Nexpt1[0] << endl;
  cout << "2 jet : " << Nexpt1[1] << endl;
  cout << "3 jet : " << Nexpt1[2] << endl;
  cout << "4 jet : " << Nexpt1[3] << endl;

  cout << "echan : " << endl;
  cout << "1 jet : " << Nexpt1E[0] << endl;
  cout << "2 jet : " << Nexpt1E[1] << endl;
  cout << "3 jet : " << Nexpt1E[2] << endl;
  cout << "4 jet : " << Nexpt1E[3] << endl;

  cout << "mchan : " << endl;
  cout << "1 jet : " << Nexpt1M[0] << endl;
  cout << "2 jet : " << Nexpt1M[1] << endl;
  cout << "3 jet : " << Nexpt1M[2] << endl;
  cout << "4 jet : " << Nexpt1M[3] << endl;

}