LoadPDF.cxx

Go to the documentation of this file.

#include <iostream>
#include "TSystem.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TFile.h"
#include "LoadPDF.h"
#include <fstream>
#include <cassert>

using std::cout;
using std::endl;

// ******************* CONSTRUCTOR & DESTRUCTOR ***********************

LoadPDF::LoadPDF() : fHist(6),fHist2(2){

  for (unsigned int i=0;i<6;i++) fHist[i]=0;
  for (unsigned int i=0;i<2;i++) fHist2[i]=0;

}

LoadPDF::LoadPDF(const char* name) : fHist(6),fHist2(2){

  for (unsigned int i=0;i<6;i++) fHist[i]=0;
  for (unsigned int i=0;i<2;i++) fHist2[i]=0;

  //Get file
  TFile file_pdf(name,"READ");

  //Get PDFs 
  fHist2[0] = dynamic_cast<TH2*>(file_pdf.Get("hist_qpeqp_planes_nu"));
  fHist2[1] = dynamic_cast<TH2*>(file_pdf.Get("hist_qpeqp_planes_anu"));
  fHist[0] = dynamic_cast<TH1*>(file_pdf.Get("hist_planes_nu"));
  fHist[1] = dynamic_cast<TH1*>(file_pdf.Get("hist_planes_anu"));
  fHist[2] = dynamic_cast<TH1*>(file_pdf.Get("hist_y_nu"));
  fHist[3] = dynamic_cast<TH1*>(file_pdf.Get("hist_y_anu"));
  fHist[4] = dynamic_cast<TH1*>(file_pdf.Get("hist_cosz_nu"));
  fHist[5] = dynamic_cast<TH1*>(file_pdf.Get("hist_cosz_anu"));
    
  for (unsigned int i=0;i<6;i++){
    assert(fHist[i]);
    fHist[i]->SetDirectory(0);
    minPDFx_1D[i]=fHist[i]->GetXaxis()->GetXmin();
    maxPDFx_1D[i]=fHist[i]->GetXaxis()->GetXmax();
  }
  for (unsigned int i=0;i<2;i++){
    assert(fHist2[i]);
    fHist2[i]->SetDirectory(0);
    minPDFx_2D[i]=fHist2[i]->GetXaxis()->GetXmin();
    maxPDFx_2D[i]=fHist2[i]->GetXaxis()->GetXmax();
    minPDFy_2D[i]=fHist2[i]->GetYaxis()->GetXmin();
    maxPDFy_2D[i]=fHist2[i]->GetYaxis()->GetXmax();
    
  }
  

}//end of constructor

LoadPDF::~LoadPDF(){

}

// ******************** PUBLIC METHODS ***********************
//---------------------------------------------
Double_t LoadPDF::GetProb1_nu(Double_t qpeqp_in, Int_t planes_in){

  Double_t proba = fHist2[0]->GetBinContent(fHist2[0]->FindBin(qpeqp_in,planes_in));

  if(proba==0 || qpeqp_in < minPDFx_2D[0] || qpeqp_in > maxPDFx_2D[0] || planes_in < minPDFy_2D[0] || planes_in > maxPDFy_2D[0]){
    proba=0.000001;
  }
  return proba;
}

//---------------------------------------------
Double_t LoadPDF::GetProb1_anu(Double_t qpeqp_in, Int_t planes_in){
  Double_t proba = fHist2[1]->GetBinContent(fHist2[1]->FindBin(qpeqp_in,planes_in));

  if(proba==0 || qpeqp_in < minPDFx_2D[1] || qpeqp_in > maxPDFx_2D[1] || planes_in < minPDFy_2D[1] || planes_in > maxPDFy_2D[1]){
    proba=0.000001;
  }
  return proba;
}

//---------------------------------------------
Double_t LoadPDF::GetProb2_nu(Int_t planes_in){
  Double_t proba = fHist[0]->GetBinContent(fHist[0]->FindBin(planes_in)); 

  if(proba==0 || planes_in < minPDFx_1D[0] || planes_in > maxPDFx_1D[0]){
    proba=0.000001;
  }
  return proba;
}

//---------------------------------------------
Double_t LoadPDF::GetProb2_anu(Int_t planes_in){
  Double_t proba = fHist[1]->GetBinContent(fHist[1]->FindBin(planes_in));

  if(proba==0 || planes_in < minPDFx_1D[1] || planes_in > maxPDFx_1D[1]){
    proba=0.000001;
  }
  return proba;
}

//---------------------------------------------
Double_t LoadPDF::GetProb3_nu(Double_t y_in){
  Double_t proba = fHist[2]->GetBinContent(fHist[2]->FindBin(y_in)); 

  if(proba==0 || y_in < minPDFx_1D[2] || y_in > maxPDFx_1D[2]){
    proba=0.000001;
  }
  return proba;
}

//---------------------------------------------
Double_t LoadPDF::GetProb3_anu(Double_t y_in){  
  Double_t proba = fHist[3]->GetBinContent(fHist[3]->FindBin(y_in));

  if(proba==0 || y_in < minPDFx_1D[3] || y_in > maxPDFx_1D[3]){
    proba=0.000001;
  }
  return proba;
}

//---------------------------------------------
Double_t LoadPDF::GetProb4_nu(Double_t cosz_in){
  Double_t proba = fHist[4]->GetBinContent(fHist[4]->FindBin(cosz_in)); 

  if(proba==0 || cosz_in < minPDFx_1D[4] || cosz_in > maxPDFx_1D[4]){
    proba=0.000001;
  }
  return proba;
}

//---------------------------------------------
Double_t LoadPDF::GetProb4_anu(Double_t cosz_in){
  Double_t proba = fHist[5]->GetBinContent(fHist[5]->FindBin(cosz_in));

  if(proba==0 || cosz_in < minPDFx_1D[5] || cosz_in > maxPDFx_1D[5]){
    proba=0.000001;
  }
  return proba;
}

---