test_fid_boundary.C

void test_fid_boundary(){
  TFile *file = new TFile("processed_files/checkout_run1_intrinsic_nue_POT1.2E23.root");
  // TFile *file = new TFile("processed_files/checkout_run1_intrinsic_nue_LowE_POT6.1E23.root");
  TTree *T_eval = (TTree*)file->Get("wcpselection/T_eval");

  Float_t truth_vtxX, truth_vtxY, truth_vtxZ;
  T_eval->SetBranchStatus("*",0);
  T_eval->SetBranchStatus("truth_vtxX",1);
  T_eval->SetBranchStatus("truth_vtxY",1);
  T_eval->SetBranchStatus("truth_vtxZ",1);
  T_eval->SetBranchAddress("truth_vtxX",&truth_vtxX);
  T_eval->SetBranchAddress("truth_vtxY",&truth_vtxY);
  T_eval->SetBranchAddress("truth_vtxZ",&truth_vtxZ);

  float bx[2]={100,100};
  float by[2]={0,0};
  float bz[2]={500,500};
  for (Int_t i=0;i!=T_eval->GetEntries();i++){
    T_eval->GetEntry(i);
    if (truth_vtxX < bx[0]) bx[0] = truth_vtxX;
    if (truth_vtxX > bx[1]) bx[1] = truth_vtxX;
    if (truth_vtxY < by[0]) by[0] = truth_vtxY;
    if (truth_vtxY > by[1]) by[1] = truth_vtxY;
    if (truth_vtxZ < bz[0]) bz[0] = truth_vtxZ;
    if (truth_vtxZ > bz[1]) bz[1] = truth_vtxZ;
  }

  std::cout << bx[0] << " " << bx[1] << std::endl;
  std::cout << by[0] << " " << by[1] << std::endl;
  std::cout << bz[0] << " " << bz[1] << std::endl;


  TTree *T_BDTvars =  (TTree*)file->Get("wcpselection/T_BDTvars");
  T_eval->AddFriend(T_BDTvars);

  TCanvas *c1 = new TCanvas("c1","c1",1200,800);
  c1->Divide(2,2);
  c1->cd(1);
  TH1F *hx1 = new TH1F("hx1","hx1",100,0,256);
  TH1F *hx2 = new TH1F("hx2","hx2",100,0,256);

  T_eval->Project("hx1","truth_vtxX");
  T_eval->Project("hx2","truth_vtxX","T_BDTvars.nue_score>7.0");
  hx2->Divide(hx1);
  hx2->Draw();

  c1->cd(2);
  TH1F *hy1 = new TH1F("hy1","hy1",100,-116,118);
  TH1F *hy2 = new TH1F("hy2","hy2",100,-116,118);

  T_eval->Project("hy1","truth_vtxY");
  T_eval->Project("hy2","truth_vtxY","T_BDTvars.nue_score>7.0");
  hy2->Divide(hy1);
  hy2->Draw();


  c1->cd(3);
  TH1F *hz1 = new TH1F("hz1","hz1",100,0,1056);
  TH1F *hz2 = new TH1F("hz2","hz2",100,0,1056);

  T_eval->Project("hz1","truth_vtxZ");
  T_eval->Project("hz2","truth_vtxZ","T_BDTvars.nue_score>7.0");
  hz2->Divide(hz1);
  hz2->Draw();
  
  //  T_eval->Draw("truth_vtxX","T_BDTvars.nue_score>7.0");
}