resolution_event_tree.C File Reference

#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncVars.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncCuts.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncExtra.h"
#include "CAFAna/Core/EventList.h"
#include "TSystem.h"
#include "CAFAna/Core/Var.h"
#include <fstream>
#include "NDAna/muonid/NDXSecMuonPID.h"

Go to the source code of this file.

Namespaces
	nuebarccinc

Functions
void	load_libs_muonid ()
std::vector< std::string >	get_files (std::string file_list)
void	resolution_event_tree (std::string file_list="", std::string eventlist="resolution_event_tree")

Function Documentation

std::vector<std::string> get_files

(

std::string

file_list

)

Definition at line 32 of file resolution_event_tree.C.

References file, compareCafs::files, in, and string.

Referenced by build_tools::merge_extraps(), and resolution_event_tree().

{
  std::vector<std::string> files;
  std::ifstream in(file_list);
  
  std::string file;
  while(std::getline(in, file)) {
    if(file.size() > 0) files.push_back(file);
  }
  return files;
}

void load_libs_muonid

(

)

Definition at line 24 of file resolution_event_tree.C.

References exit(), and runNovaSAM::ret.

Referenced by resolution_event_tree().

{
  const int ret = gSystem->Load("libNDAnamuonid");
  if(ret != 0) exit(1);
}

void resolution_event_tree	(	std::string	file_list = "",
		std::string	eventlist = "resolution_event_tree"
	)

Definition at line 44 of file resolution_event_tree.C.

References abs(), caf::Proxy< caf::SRVertexBranch >::elastic, make_root_from_grid_output::file_list, compareCafs::files, caf::Proxy< caf::SRElastic >::fuzzyk, get_files(), if(), iscc, caf::Proxy< caf::SRElastic >::IsValid, ana::kCycle, ana::kElectronPDG, nuebarccinc::kHasReconstructedElectronProng, ana::xsec::numubarcc::kIsCC, nuebarccinc::kMuonIDProd4Cut, nuebarccinc::kNuebarCC, nuebarccinc::kNueCC, nuebarccinc::kPreselectionLoose, kRun, kSlice, ana::kSubrun, nuebarccinc::kTrueElectronCosTheta, nuebarccinc::kTrueElectronE, nuebarccinc::kTrueElectronRecoCosTheta, nuebarccinc::kTrueElectronRecoEFHCTuned, nuebarccinc::kTrueElectronRecoEUntuned, load_libs_muonid(), ana::MakeEventTTreeFile(), caf::Proxy< caf::StandardRecord >::mc, caf::Proxy< caf::SRFuzzyK >::npng, caf::Proxy< caf::SRTruthBranch >::nu, make_root_from_grid_output::pdg, caf::Proxy< caf::SRFuzzyK >::png, nuebarccinc::PROD5_MC_RHC_NOMINAL, SIMPLENUTRUTHVAR, SIMPLEVAR, sr, nuebarccinc::true_electron_png(), ana::VarFromNuTruthVar(), and caf::Proxy< caf::StandardRecord >::vtx.

{
  
  std::vector<std::string> files;
  if(file_list != "") files= get_files(file_list);

  const Var kNuPDG = VarFromNuTruthVar(SIMPLENUTRUTHVAR(pdg));
  const Var kElectronPDG([](const caf::SRProxy * sr)->int
                         {
                           if ( !sr->vtx.elastic.IsValid ) return -1;
                           if (sr->vtx.elastic.fuzzyk.npng < 1) return -1;
                           return sr->vtx.elastic.fuzzyk.png[nuebarccinc::true_electron_png(sr)].truth.pdg;
                         });
  const Var kElectronProngIdx([](const caf::SRProxy * sr)->int
                              {
                                return nuebarccinc::true_electron_png(sr);
                              });

  const Var kIsCC = VarFromNuTruthVar(SIMPLENUTRUTHVAR(iscc));
  const Var kHasElectronFuzzyK([](const caf::SRProxy * sr)->int
                               {
                                 if ( !sr->vtx.elastic.IsValid ) return -1;
                                 if (sr->vtx.elastic.fuzzyk.npng < 1) return -1;
                                 int npngs = sr->vtx.elastic.fuzzyk.npng;
                                 bool haselectron = false;
                                 for(auto ipng = 0; ipng < npngs; ipng++) {
                                   if(abs(sr->vtx.elastic.fuzzyk.png[ipng].truth.pdg) == 11) {
                                     haselectron = true;
                                     break;
                                   }
                                 }
                                 return haselectron;
                               });

  const Var kHasElectronPrim([](const caf::SRProxy * sr)->int
                             {
                               if ( !sr->vtx.elastic.IsValid ) return -1;
                               if (sr->vtx.elastic.fuzzyk.npng < 1) return -1;
                               int npngs = sr->mc.nu[0].prim.size();
                               bool haselectron = false;
                               for(auto ipng = 0; ipng < npngs; ipng++) {
                                 if(abs(sr->mc.nu[0].prim[ipng].pdg) == 11) {
                                   haselectron = true;
                                   break;
                                 }
                               }
                               return haselectron;
                             });
  
  
  const Var kRun = SIMPLEVAR(hdr.run);
  const Var kSubrun = SIMPLEVAR(hdr.subrun);
  const Var kCycle = SIMPLEVAR(hdr.cycle);
  const Var kBatch = SIMPLEVAR(hdr.batch);
  const Var kEvent = SIMPLEVAR(hdr.evt);
  const Var kSlice = SIMPLEVAR(hdr.subevt);


  load_libs_muonid();
  const Cut kSelectedNueNuebarCC = 
    nuebarccinc::kPreselectionLoose && 
    nuebarccinc::kMuonIDProd4Cut    && 
    (nuebarccinc::kNueCC || nuebarccinc::kNuebarCC);

  std::vector<std::pair<std::string, Var> > floatvars 
    = { {"reco_cos_theta", nuebarccinc::kTrueElectronRecoCosTheta},
        {"true_cos_theta", nuebarccinc::kTrueElectronCosTheta},
        {"reco_e_untuned", nuebarccinc::kTrueElectronRecoEUntuned},
        {"reco_e_fhctuned", nuebarccinc::kTrueElectronRecoEFHCTuned},
        {"true_e", nuebarccinc::kTrueElectronE} };

  std::vector<std::pair<std::string, Var> > intvars
    = { {"nupdg", kNuPDG},
        {"electron_pdg", kElectronPDG},
        {"iscc", kIsCC},
        {"has_electron_fuzzyk", kHasElectronFuzzyK},
        {"has_electron_prim", kHasElectronPrim},
        {"electron_png_idx", kElectronProngIdx},
        {"run", kRun},
        {"subrun", kSubrun},
        {"cycle", kCycle},
        {"batch", kBatch},
        {"event", kEvent},
        {"slice", kSlice} };
         

  if(file_list != "")
    MakeEventTTreeFile(files,
                       eventlist + ".root",
                       {{"selected", nuebarccinc::kHasReconstructedElectronProng && kSelectedNueNuebarCC}},
                       floatvars,
                       intvars);                     
  else
    MakeEventTTreeFile(nuebarccinc::PROD5_MC_RHC_NOMINAL,
                       eventlist + ".root",
                       //{{"selected", nuebarccinc::kHasReconstructedElectronProng && kSelectedNueNuebarCC && kBackwardsElectrons}},
                       {{"selected", nuebarccinc::kHasReconstructedElectronProng && kSelectedNueNuebarCC}},
                       floatvars,
                       intvars);                     
                    
}