nus17_fiducial_mass.C

Go to the documentation of this file.

#include "TFile.h"
#include "TF1.h"
#include "TH1.h"
#include "TEfficiency.h"
#include "TCanvas.h"
#include "TF1.h"

void nus17_fiducial_mass(std::string input)
{
  TH1::SetDefaultSumw2();
  
  TFile* fileIn = new TFile(input.c_str(), "READ");
 
  TH1D* hSigEffFull = new TH1D("hSigEff_FidMass","Nus17 NC Selection Efficiency;Diblocks;Efficiency",11,3.5,14.5);
  TH1D* hSigEffIdeal = new TH1D("hSigEff_SmallWindow","Nus17 NC Selection Efficiency;Diblocks;efficiency",11,3.5,14.5);

  TH1D* hNFull = (TH1D*)fileIn->Get("hN_FidMass");
  TH1D* hDFull = (TH1D*)fileIn->Get("hD_FidMass");
  TH1D* hNIdeal = (TH1D*)fileIn->Get("hN_SmallWindow");
  TH1D* hDIdeal = (TH1D*)fileIn->Get("hD_SmallWindow");

  hSigEffFull->Divide(hNFull, hDFull, 1.0, 1.0, "b");
  hSigEffIdeal->Divide(hNIdeal, hDIdeal, 1.0, 1.0, "b");


  TF1 *asymptote = new TF1("asymptote", "[0] + [1]/(0.993*x)", 0., 1.);
  asymptote->SetParameter(0, 1);
  asymptote->SetParameter(1, 1);

  TCanvas *c0 = new TCanvas("c0", "c0", 200, 10, 800, 600);
  c0->cd();
  hSigEffFull->Fit("asymptote");
  hSigEffFull->Draw();
  asymptote->Draw("same");
  c0->Update();
  c0->Print("Efficiency_vs_Diblocks.pdf");
  
  double a = asymptote->GetParameter(0);
  double del_a = asymptote->GetParError(0);
  double b = asymptote->GetParameter(1);
  double del_b = asymptote->GetParError(1);


  double ideal_eff = hSigEffIdeal->GetBinContent(11);

  std::cout << "Fiducial mass tables : " << std::endl;
  std::cout << "=======================" << std::endl;
  for(int i = 1; i <= hSigEffFull->GetNbinsX(); i++){
    double m_f = ((i+3)*0.993/ideal_eff)*(a+(b/((i+3)*0.993))); 

    std::cout << "N_diblock : " << i+3 << ", Fiducial Mass : " << m_f << ", Efficiency : " << hSigEffFull->GetBinContent(i) << std::endl;
  }

  std::cout << "=======================" << std::endl;
  std::cout << "Ideal Efficiency : " << ideal_eff << std::endl;

}