dunemva::CAFMaker Class Reference

Inheritance diagram for dunemva::CAFMaker:

Public Member Functions
	CAFMaker (fhicl::ParameterSet const &pset)
virtual	~CAFMaker ()
void	beginJob () override
void	endJob () override
void	beginSubRun (const art::SubRun &sr) override
void	endSubRun (const art::SubRun &sr) override
void	reconfigure (fhicl::ParameterSet const &pset)
void	analyze (art::Event const &evt) override
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
std::string	fMVASelectLabel
std::string	fMVASelectNueLabel
std::string	fMVASelectNumuLabel
std::string	fCVNLabel
std::string	fRegCNNLabel
std::string	fEnergyRecoNueLabel
std::string	fEnergyRecoNumuLabel
std::string	fMVAMethod
float	fOscPro
double	fWeight
TTree *	fTree
TTree *	fMetaTree
int	fNwgt [knShifts]
double	fCvWgts [knShifts]
double	fWgts [knShifts][kmaxRwgts]
int	fIsFD
int	fIsFHC
int	fRun
int	fSubrun
int	fEvent
int	fIsCC
int	fNuPDG
int	fNuPDGunosc
int	fMode
int	fLepPDG
double	fEv
double	fQ2
double	fW
double	fX
double	fY
double	fNuMomX
double	fNuMomY
double	fNuMomZ
double	fLepMomX
double	fLepMomY
double	fLepMomZ
double	fLepE
double	fLepNuAngle
int	nP
int	nN
int	nPip
int	nPim
int	nPi0
int	nKp
int	nKm
int	nK0
int	nEM
int	nOtherHad
int	nNucleus
int	nUNKNOWN
double	eP
double	eN
double	ePip
double	ePim
double	ePi0
double	eOther
double	vtx_x
double	vtx_y
double	vtx_z
double	fErecoNue
double	fRecoLepEnNue
double	fRecoHadEnNue
int	fRecoMethodNue
double	fErecoNumu
double	fRecoLepEnNumu
double	fRecoHadEnNumu
int	fRecoMethodNumu
int	fLongestTrackContNumu
int	fTrackMomMethodNumu
double	fMVAResult
double	fMVAResultNue
double	fMVAResultNumu
double	fCVNResultIsAntineutrino
double	fCVNResultNue
double	fCVNResultNumu
double	fCVNResultNutau
double	fCVNResultNC
double	fCVNResult0Protons
double	fCVNResult1Protons
double	fCVNResult2Protons
double	fCVNResultNProtons
double	fCVNResult0Pions
double	fCVNResult1Pions
double	fCVNResult2Pions
double	fCVNResultNPions
double	fCVNResult0Pizeros
double	fCVNResult1Pizeros
double	fCVNResult2Pizeros
double	fCVNResultNPizeros
double	fCVNResult0Neutrons
double	fCVNResult1Neutrons
double	fCVNResult2Neutrons
double	fCVNResultNNeutrons
double	fRegCNNNueE
double	meta_pot
int	meta_run
int	meta_subrun
int	meta_version
systtools::provider_list_t	fSystProviders

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 61 of file CAFMaker_module.cc.

Constructor & Destructor Documentation

dunemva::CAFMaker::CAFMaker ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 145 of file CAFMaker_module.cc.

    : EDAnalyzer(pset)
  {
    this->reconfigure(pset);
  }

dunemva::CAFMaker::~CAFMaker ( )

virtual

Definition at line 151 of file CAFMaker_module.cc.

{}

Member Function Documentation

void dunemva::CAFMaker::analyze ( art::Event const &  evt )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 323 of file CAFMaker_module.cc.

  {
    auto pidin = evt.getHandle<dunemva::MVASelectPID>(fMVASelectLabel);
    auto pidinnue = evt.getHandle<dunemva::MVASelectPID>(fMVASelectNueLabel);
    auto pidinnumu = evt.getHandle<dunemva::MVASelectPID>(fMVASelectNumuLabel);
    art::InputTag itag1(fCVNLabel, "cvnresult");
    auto cvnin = evt.getHandle<std::vector<cvn::Result>>(itag1);
    art::InputTag itag2(fRegCNNLabel, "regcnnresult");
    auto regcnnin = evt.getHandle<std::vector<cnn::RegCNNResult>>(itag2);
    auto ereconuein = evt.getHandle<dune::EnergyRecoOutput>(fEnergyRecoNueLabel);
    auto ereconumuin = evt.getHandle<dune::EnergyRecoOutput>(fEnergyRecoNumuLabel);

    fRun = evt.id().run();
    fSubrun = evt.id().subRun();
    fEvent = evt.id().event();
    meta_run = fRun;
    meta_subrun = fSubrun;

    if( !pidin.failedToGet() ) {
      fMVAResult = pidin->pid;

      //Fill MVA reco stuff
      fErecoNue          = ereconuein->fNuLorentzVector.E();
      fRecoLepEnNue      = ereconuein->fLepLorentzVector.E();
      fRecoHadEnNue      = ereconuein->fHadLorentzVector.E();
      fRecoMethodNue     = ereconuein->recoMethodUsed;
      fErecoNumu         = ereconumuin->fNuLorentzVector.E();
      fRecoLepEnNumu     = ereconumuin->fLepLorentzVector.E();
      fRecoHadEnNumu     = ereconumuin->fHadLorentzVector.E();
      fRecoMethodNumu    = ereconumuin->recoMethodUsed;
      fLongestTrackContNumu  = ereconumuin->longestTrackContained;
      fTrackMomMethodNumu    = ereconumuin->trackMomMethod;
    }

    if( !pidinnue.failedToGet() ) {
      fMVAResultNue = pidinnue->pid;
    }

    if( !pidinnumu.failedToGet() ) {
      fMVAResultNumu = pidinnumu->pid;
    }

    if( !cvnin.failedToGet() ) {
      //using i = cvn::Interaction;
      //if(cvnin->empty() || (*cvnin)[0].fOutput.size() <= i::kNutauOther){
      if(cvnin->empty()){
        fCVNResultIsAntineutrino = fCVNResultNue = fCVNResultNumu = fCVNResultNutau = fCVNResultNC = \
        fCVNResult0Protons = fCVNResult1Protons = fCVNResult2Protons = fCVNResultNProtons = \
        fCVNResult0Pions = fCVNResult1Pions = fCVNResult2Pions = fCVNResultNPions = \
        fCVNResult0Pizeros = fCVNResult1Pizeros = fCVNResult2Pizeros = fCVNResultNPizeros = \
        fCVNResult0Neutrons = fCVNResult1Neutrons = fCVNResult2Neutrons = fCVNResultNNeutrons = -3;
      }
      else{
        //const std::vector<float>& v = (*cvnin)[0].fOutput;
        //fCVNResultNue = v[i::kNueQE] + v[i::kNueRes] + v[i::kNueDIS] + v[i::kNueOther];
        //fCVNResultNumu = v[i::kNumuQE] + v[i::kNumuRes] + v[i::kNumuDIS] + v[i::kNumuOther];
        //fCVNResultNutau = v[i::kNutauQE] + v[i::kNutauRes] + v[i::kNutauDIS] + v[i::kNutauOther]

        fCVNResultIsAntineutrino = (*cvnin)[0].GetIsAntineutrinoProbability();

        fCVNResultNue = (*cvnin)[0].GetNueProbability();
        fCVNResultNumu = (*cvnin)[0].GetNumuProbability();
        fCVNResultNutau = (*cvnin)[0].GetNutauProbability();
        fCVNResultNC = (*cvnin)[0].GetNCProbability();

        fCVNResult0Protons = (*cvnin)[0].Get0protonsProbability();
        fCVNResult1Protons = (*cvnin)[0].Get1protonsProbability();
        fCVNResult2Protons = (*cvnin)[0].Get2protonsProbability();
        fCVNResultNProtons = (*cvnin)[0].GetNprotonsProbability();

        fCVNResult0Pions = (*cvnin)[0].Get0pionsProbability();
        fCVNResult1Pions = (*cvnin)[0].Get1pionsProbability();
        fCVNResult2Pions = (*cvnin)[0].Get2pionsProbability();
        fCVNResultNPions = (*cvnin)[0].GetNpionsProbability();

        fCVNResult0Pizeros = (*cvnin)[0].Get0pizerosProbability();
        fCVNResult1Pizeros = (*cvnin)[0].Get1pizerosProbability();
        fCVNResult2Pizeros = (*cvnin)[0].Get2pizerosProbability();
        fCVNResultNPizeros = (*cvnin)[0].GetNpizerosProbability();

        fCVNResult0Neutrons = (*cvnin)[0].Get0neutronsProbability();
        fCVNResult1Neutrons = (*cvnin)[0].Get1neutronsProbability();
        fCVNResult2Neutrons = (*cvnin)[0].Get2neutronsProbability();
        fCVNResultNNeutrons = (*cvnin)[0].GetNneutronsProbability();
      }
    }

    fRegCNNNueE = -1.;  // initializing
    if(!regcnnin.failedToGet()){
      if (!regcnnin->empty()){
        const std::vector<float>& v = (*regcnnin)[0].fOutput;
        fRegCNNNueE = v[0];
      }
    }

    std::vector< art::Ptr<simb::MCTruth> > truth;
    auto mct = evt.getHandle< std::vector<simb::MCTruth> >("generator");
    if ( mct )
      art::fill_ptr_vector(truth, mct);
    else
      mf::LogWarning("CAFMaker") << "No MCTruth.";

    std::vector< art::Ptr<simb::MCFlux> > flux;
    auto mcf = evt.getHandle< std::vector<simb::MCFlux> >("generator");
    if ( mcf )
      art::fill_ptr_vector(flux, mcf);
    else
      mf::LogWarning("CAFMaker") << "No MCFlux.";
/*
    std::vector< art::Ptr<simb::GTruth> > gtru;
    auto gt = evt.getHandle< std::vector<simb::GTruth> >("generator");
    if ( gt )
      art::fill_ptr_vector(gtru, gt);
    else
      mf::LogWarning("CAFMaker") << "No GTruth.";
*/

    for(size_t i=0; i<truth.size(); i++){

      if(i>1){
        mf::LogWarning("CAFMaker") << "Skipping MC truth index " << i;
        continue;
      }

      fIsFD     = 1; // always FD
      fIsFHC    = 999; // don't know how to get this?
      fIsCC     = !(truth[i]->GetNeutrino().CCNC());  // ccnc is 0=CC 1=NC
      fNuPDG    = truth[i]->GetNeutrino().Nu().PdgCode();
      fNuPDGunosc = flux[i]->fntype;
      fMode     = truth[i]->GetNeutrino().Mode(); //0=QE/El, 1=RES, 2=DIS, 3=Coherent production; this is different than mode in ND
      fEv       = truth[i]->GetNeutrino().Nu().E();
      fQ2       = truth[i]->GetNeutrino().QSqr();
      fW        = truth[i]->GetNeutrino().W();
      fX        = truth[i]->GetNeutrino().X();
      fY        = truth[i]->GetNeutrino().Y();
      fNuMomX   = truth[i]->GetNeutrino().Nu().Momentum().X();
      fNuMomY   = truth[i]->GetNeutrino().Nu().Momentum().Y();
      fNuMomZ   = truth[i]->GetNeutrino().Nu().Momentum().Z();

      vtx_x     = truth[i]->GetNeutrino().Lepton().Vx();
      vtx_y     = truth[i]->GetNeutrino().Lepton().Vy();
      vtx_z     = truth[i]->GetNeutrino().Lepton().Vz();

      //Lepton stuff
      fLepPDG     = truth[i]->GetNeutrino().Lepton().PdgCode();
      fLepMomX    = truth[i]->GetNeutrino().Lepton().Momentum().X();
      fLepMomY    = truth[i]->GetNeutrino().Lepton().Momentum().Y();
      fLepMomZ    = truth[i]->GetNeutrino().Lepton().Momentum().Z();
      fLepE       = truth[i]->GetNeutrino().Lepton().Momentum().T();
      fLepNuAngle = truth[i]->GetNeutrino().Nu().Momentum().Vect().Angle(truth[i]->GetNeutrino().Lepton().Momentum().Vect());

      // TODO
      fWeight = 0.;
      fOscPro = 0.;
      //      fOscPro   = fMVAAlg.OscPro(fCCNC,fBeamPdg,fNuPDG,fEtrue);

      nP        = 0;
      nN        = 0;
      nPip      = 0;
      nPim      = 0;
      nPi0      = 0;
      nKp       = 0;
      nKm       = 0;
      nK0       = 0;
      nEM       = 0;
      nOtherHad = 0;
      nNucleus  = 0;
      nUNKNOWN  = 0;

      eP = 0.;
      eN = 0.;
      ePip = 0.;
      ePim = 0.;
      ePi0 = 0.;
      eOther = 0.;

      for( int p = 0; p < truth[i]->NParticles(); p++ ) {
        if( truth[i]->GetParticle(p).StatusCode() == genie::kIStHadronInTheNucleus ) {

          int pdg = truth[i]->GetParticle(p).PdgCode();
          double ke = truth[i]->GetParticle(p).E() - truth[i]->GetParticle(p).Mass();
          if     ( pdg == genie::kPdgProton ) {
            nP++;
            eP += ke;
          } else if( pdg == genie::kPdgNeutron ) {
            nN++;
            eN += ke;
          } else if( pdg == genie::kPdgPiP ) {
            nPip++;
            ePip += ke;
          } else if( pdg == genie::kPdgPiM ) {
            nPim++;
            ePim += ke;
          } else if( pdg == genie::kPdgPi0 ) {
            nPi0++;
            ePi0 += ke;
          } else if( pdg == genie::kPdgKP ) {
            nKp++;
            eOther += ke;
          } else if( pdg == genie::kPdgKM ) {
            nKm++;
            eOther += ke;
          } else if( pdg == genie::kPdgK0 || pdg == genie::kPdgAntiK0 || pdg == genie::kPdgK0L || pdg == genie::kPdgK0S ) {
            nK0++;
            eOther += ke;
          } else if( pdg == genie::kPdgGamma ) {
            nEM++;
            eOther += ke;
          } else if( genie::pdg::IsHadron(pdg) ) {
            nOtherHad++; // charm mesons, strange and charm baryons, antibaryons, etc.
            eOther += ke;
          } else if( genie::pdg::IsIon(pdg) ) {
            nNucleus++;
          } else {
            nUNKNOWN++;
          }

        }
      }

      // Reweighting variables
      //systtools::ScrubUnityEventResponses(er);

      // struct ParamResponses { 
      //   paramId_t pid;
      //   std::vector<double> responses;
      // }
      // typedef std::vector<ParamResponses> event_unit_response_t;
      // typedef std::vector<event_unit_response_t> EventResponse;

      for( auto &sp : fSystProviders ) {
        std::unique_ptr<systtools::EventResponse> syst_resp = sp->GetEventResponse(evt);
        if( !syst_resp ) {
          std::cout << "[ERROR]: Got nullptr systtools::EventResponse from provider "
                    << sp->GetFullyQualifiedName();
          continue;
        }

        for( systtools::EventResponse::iterator itResp = syst_resp->begin(); itResp != syst_resp->end(); ++itResp ) {
          systtools::event_unit_response_t resp = *itResp;
          for( systtools::event_unit_response_t::iterator it = resp.begin(); it != resp.end(); ++it ) {
            fNwgt[(*it).pid] = (*it).responses.size();
            //fCvWgts[(*it).pid] = (*it).CV_weight;
            for( unsigned int i = 0; i < (*it).responses.size(); ++i ) {
              fWgts[(*it).pid][i] = (*it).responses[i];
            }
          }
        }
      }
    } // loop through MC truth i

    fTree->Fill();
    return;
  }

void dunemva::CAFMaker::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 182 of file CAFMaker_module.cc.

  {

    art::ServiceHandle<art::TFileService> tfs;
    fTree = tfs->make<TTree>("caf", "caf");
    fMetaTree = tfs->make<TTree>("meta", "meta");

    // book-keeping
    fTree->Branch("run",         &fRun,        "run/I");
    fTree->Branch("subrun",      &fSubrun,     "subrun/I");
    fTree->Branch("event",       &fEvent,      "event/I");
    fTree->Branch("isFD",        &fIsFD,       "isFD/I");
    fTree->Branch("isFHC",       &fIsFHC,      "isFHC/I");
    fTree->Branch("isCC",        &fIsCC,       "isCC/I");

    // true interaction quantities
    fTree->Branch("nuPDG",        &fNuPDG,        "nuPDG/I");
    fTree->Branch("nuPDGunosc",   &fNuPDGunosc,   "nuPDGunosc/I");
    fTree->Branch("NuMomX",       &fNuMomX,       "NuMomX/D");
    fTree->Branch("NuMomY",       &fNuMomY,       "NuMomY/D");
    fTree->Branch("NuMomZ",       &fNuMomZ,       "NuMomZ/D");
    fTree->Branch("Ev",           &fEv,           "Ev/D");
    fTree->Branch("mode",         &fMode,         "mode/I");
    fTree->Branch("LepPDG",       &fLepPDG,       "LepPDG/I");
    fTree->Branch("LepMomX",      &fLepMomX,      "LepMomX/D");
    fTree->Branch("LepMomY",      &fLepMomY,      "LepMomY/D");
    fTree->Branch("LepMomZ",      &fLepMomZ,      "LepMomZ/D");
    fTree->Branch("LepE",         &fLepE,         "LepE/D");
    fTree->Branch("LepNuAngle",   &fLepNuAngle,   "LepNuAngle/D");
    fTree->Branch("Q2",           &fQ2,           "Q2/D");
    fTree->Branch("W",            &fW,            "W/D");
    fTree->Branch("X",            &fX,            "X/D");
    fTree->Branch("Y",            &fY,            "Y/D");

    // FS particle counts
    fTree->Branch("nP",        &nP,         "nP/I");
    fTree->Branch("nN",        &nN,         "nN/I");
    fTree->Branch("nipip",     &nPip,       "nipip/I");
    fTree->Branch("nipim",     &nPim,       "nipim/I");
    fTree->Branch("nipi0",     &nPi0,       "nipi0/I");
    fTree->Branch("nikp",      &nKp,        "nikp/I");
    fTree->Branch("nikm",      &nKm,        "nikm/I");
    fTree->Branch("nik0",      &nK0,        "nik0/I");
    fTree->Branch("niem",      &nEM,        "niem/I");
    fTree->Branch("niother",   &nOtherHad,  "niother/I");
    fTree->Branch("nNucleus",  &nNucleus,   "nNucleus/I");
    fTree->Branch("nUNKNOWN",  &nUNKNOWN,   "nUNKNOWN/I");
    fTree->Branch("eP",        &eP,         "eP/D");
    fTree->Branch("eN",        &eN,         "eN/D");
    fTree->Branch("ePip",      &ePip,       "ePip/D");
    fTree->Branch("ePim",      &ePim,       "ePim/D");
    fTree->Branch("ePi0",      &ePi0,       "ePi0/D");
    fTree->Branch("eOther",    &eOther,     "eOther/D");

    // vertex position
    fTree->Branch("vtx_x",   &vtx_x,    "vtx_x/D");
    fTree->Branch("vtx_y",   &vtx_y,    "vtx_y/D");
    fTree->Branch("vtx_z",   &vtx_z,    "vtx_z/D");

    // Reco variables
    fTree->Branch("mvaresult",   &fMVAResult,  "mvaresult/D");
    fTree->Branch("mvanue",      &fMVAResultNue,  "mvanue/D");
    fTree->Branch("mvanumu",     &fMVAResultNumu, "mvanumu/D");

    fTree->Branch("cvnisantineutrino", &fCVNResultIsAntineutrino, "cvnisantineutrino/D");
    fTree->Branch("cvnnue",            &fCVNResultNue,            "cvnnue/D");
    fTree->Branch("cvnnumu",           &fCVNResultNumu,           "cvnnumu/D");
    fTree->Branch("cvnnutau",          &fCVNResultNutau,          "cvnnutau/D");
    fTree->Branch("cvnnc",             &fCVNResultNC,             "cvnnc/D");
    fTree->Branch("cvn0protons",       &fCVNResult0Protons,       "cvn0protons/D");
    fTree->Branch("cvn1protons",       &fCVNResult1Protons,       "cvn1protons/D");
    fTree->Branch("cvn2protons",       &fCVNResult2Protons,       "cvn2protons/D");
    fTree->Branch("cvnNprotons",       &fCVNResultNProtons,       "cvnNprotons/D");
    fTree->Branch("cvn0pions",         &fCVNResult0Pions,         "cvn0pions/D");
    fTree->Branch("cvn1pions",         &fCVNResult1Pions,         "cvn1pions/D");
    fTree->Branch("cvn2pions",         &fCVNResult2Pions,         "cvn2pions/D");
    fTree->Branch("cvnNpions",         &fCVNResultNPions,         "cvnNpions/D");
    fTree->Branch("cvn0pizeros",       &fCVNResult0Pizeros,       "cvn0pizeros/D");
    fTree->Branch("cvn1pizeros",       &fCVNResult1Pizeros,       "cvn1pizeros/D");
    fTree->Branch("cvn2pizeros",       &fCVNResult2Pizeros,       "cvn2pizeros/D");
    fTree->Branch("cvnNpizeros",       &fCVNResultNPizeros,       "cvnNpizeros/D");
    fTree->Branch("cvn0neutrons",      &fCVNResult0Neutrons,      "cvn0neutrons/D");
    fTree->Branch("cvn1neutrons",      &fCVNResult1Neutrons,      "cvn1neutrons/D");
    fTree->Branch("cvn2neutrons",      &fCVNResult2Neutrons,      "cvn2neutrons/D");
    fTree->Branch("cvnNneutrons",      &fCVNResultNNeutrons,      "cvnNneutrons/D");

    fTree->Branch("RegCNNNueE",  &fRegCNNNueE,   "RegCNNNueE/D");
    fTree->Branch("weight",      &fWeight,     "weight/D");
    fTree->Branch("oscpro",      &fOscPro,     "oscpro/F");

    fTree->Branch("Ev_reco_nue",      &fErecoNue,        "Ev_reco_nue/D");
    fTree->Branch("RecoLepEnNue",     &fRecoLepEnNue,    "RecoLepEnNue/D");
    fTree->Branch("RecoHadEnNue",     &fRecoHadEnNue,    "RecoHadEnNue/D");
    fTree->Branch("RecoMethodNue",    &fRecoMethodNue,   "RecoMethodNue/I");
    fTree->Branch("Ev_reco_numu",     &fErecoNumu,       "Ev_reco_numu/D");
    fTree->Branch("RecoLepEnNumu",    &fRecoLepEnNumu,   "RecoLepEnNumu/D");
    fTree->Branch("RecoHadEnNumu",    &fRecoHadEnNumu,   "RecoHadEnNumu/D");
    fTree->Branch("RecoMethodNumu",   &fRecoMethodNumu,  "RecoMethodNumu/I");
    fTree->Branch("LongestTrackContNumu",  &fLongestTrackContNumu, "LongestTrackContNumu/I");
    fTree->Branch("TrackMomMethodNumu",    &fTrackMomMethodNumu,   "TrackMomMethodNumu/I");

    fMetaTree->Branch("pot", &meta_pot, "pot/D");
    fMetaTree->Branch("run", &meta_run, "run/I");
    fMetaTree->Branch("subrun", &meta_subrun, "subrun/I");
    fMetaTree->Branch("version", &meta_version, "version/I");

    // make DUNErw variables
    for( auto &sp : fSystProviders ) {
      systtools::SystMetaData metaData = sp->GetSystMetaData();
      for( systtools::SystMetaData::iterator itMeta = metaData.begin(); itMeta != metaData.end(); ++itMeta ) {      
        systtools::SystParamHeader head = *itMeta;
        std::string name = head.prettyName;
        unsigned int parId = head.systParamId;
        std::cout << "Adding reweight branch " << parId << " for " << name << " with " << head.paramVariations.size() << " shifts" << std::endl;
        fTree->Branch( Form("%s_nshifts", name.c_str()), &fNwgt[parId], Form("%s_nshifts/I", name.c_str()) );
        fTree->Branch( Form("%s_cvwgt", name.c_str()), &fCvWgts[parId], Form("%s_cvwgt/D", name.c_str()) );
        fTree->Branch( Form("wgt_%s", name.c_str()), fWgts[parId], Form("wgt_%s[%s_nshifts]/D", name.c_str(), name.c_str()) );
      }
    }

    // initialize weight variables -- some won't ever be set
    for( int i = 0; i < knShifts; ++i ) {
      fNwgt[i] = 0;
      fCvWgts[i] = 1.;
      for( int j = 0; j < kmaxRwgts; ++j ) {
        fWgts[i][j] = 0.;
      }
    }

    meta_pot = 0.;
    meta_version = 1;
  }

void dunemva::CAFMaker::beginSubRun ( const art::SubRun &  sr )

override

Definition at line 316 of file CAFMaker_module.cc.

  {
    auto pots = sr.getHandle< sumdata::POTSummary >("generator");
    if( pots ) meta_pot += pots->totpot;
  }

void dunemva::CAFMaker::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 584 of file CAFMaker_module.cc.

  {
    fMetaTree->Fill();
  }

void dunemva::CAFMaker::endSubRun ( const art::SubRun &  sr )

override

Definition at line 581 of file CAFMaker_module.cc.

                                             {
  }

void dunemva::CAFMaker::reconfigure

(

fhicl::ParameterSet const &

pset

)

Definition at line 154 of file CAFMaker_module.cc.

  {

    fMVASelectLabel = pset.get<std::string>("MVASelectLabel");
    fMVASelectNueLabel = pset.get<std::string>("MVASelectNueLabel");
    fMVASelectNumuLabel = pset.get<std::string>("MVASelectNumuLabel");
    fCVNLabel = pset.get<std::string>("CVNLabel");
    fRegCNNLabel = pset.get<std::string>("RegCNNLabel");

    fEnergyRecoNueLabel = pset.get<std::string>("EnergyRecoNueLabel");
    fEnergyRecoNumuLabel = pset.get<std::string>("EnergyRecoNumuLabel");

    // Get DUNErw stuff from its fhicl, which should be included on the CAFMaker config file
    //if( !pset.has_key("generated_systematic_provider_configuration") ) {
    //  std::cout << "[ERROR]: Could not find producer key: "
    //               "\"generated_systematic_provider_configuration\". This should "
    //               "contain a list of configured systematic providers generated by "
    //               "GenerateSystProviderConfig." << std::endl;
    //  return;
    //}

    fhicl::ParameterSet syst_provider_config = pset.get<fhicl::ParameterSet>("generated_systematic_provider_configuration");

    fSystProviders = systtools::ConfigureISystProvidersFromParameterHeaders(syst_provider_config);
  }

Member Data Documentation

double dunemva::CAFMaker::eN

private

Definition at line 107 of file CAFMaker_module.cc.

double dunemva::CAFMaker::eOther

private

Definition at line 107 of file CAFMaker_module.cc.

double dunemva::CAFMaker::eP

private

Definition at line 107 of file CAFMaker_module.cc.

double dunemva::CAFMaker::ePi0

private

Definition at line 107 of file CAFMaker_module.cc.

double dunemva::CAFMaker::ePim

private

Definition at line 107 of file CAFMaker_module.cc.

double dunemva::CAFMaker::ePip

private

Definition at line 107 of file CAFMaker_module.cc.

std::string dunemva::CAFMaker::fCVNLabel

private

Definition at line 80 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult0Neutrons

private

Definition at line 132 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult0Pions

private

Definition at line 130 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult0Pizeros

private

Definition at line 131 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult0Protons

private

Definition at line 129 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult1Neutrons

private

Definition at line 132 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult1Pions

private

Definition at line 130 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult1Pizeros

private

Definition at line 131 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult1Protons

private

Definition at line 129 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult2Neutrons

private

Definition at line 132 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult2Pions

private

Definition at line 130 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult2Pizeros

private

Definition at line 131 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResult2Protons

private

Definition at line 129 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultIsAntineutrino

private

Definition at line 127 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultNC

private

Definition at line 128 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultNNeutrons

private

Definition at line 132 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultNPions

private

Definition at line 130 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultNPizeros

private

Definition at line 131 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultNProtons

private

Definition at line 129 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultNue

private

Definition at line 128 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultNumu

private

Definition at line 128 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCVNResultNutau

private

Definition at line 128 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fCvWgts[knShifts]

private

Definition at line 94 of file CAFMaker_module.cc.

std::string dunemva::CAFMaker::fEnergyRecoNueLabel

private

Definition at line 83 of file CAFMaker_module.cc.

std::string dunemva::CAFMaker::fEnergyRecoNumuLabel

private

Definition at line 84 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fErecoNue

private

Definition at line 111 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fErecoNumu

private

Definition at line 115 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fEv

private

Definition at line 104 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fEvent

private

Definition at line 101 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fIsCC

private

Definition at line 103 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fIsFD

private

Definition at line 99 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fIsFHC

private

Definition at line 99 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fLepE

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fLepMomX

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fLepMomY

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fLepMomZ

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fLepNuAngle

private

Definition at line 104 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fLepPDG

private

Definition at line 103 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fLongestTrackContNumu

private

Definition at line 119 of file CAFMaker_module.cc.

TTree* dunemva::CAFMaker::fMetaTree

private

Definition at line 90 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fMode

private

Definition at line 103 of file CAFMaker_module.cc.

std::string dunemva::CAFMaker::fMVAMethod

private

Definition at line 85 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fMVAResult

private

Definition at line 122 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fMVAResultNue

private

Definition at line 123 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fMVAResultNumu

private

Definition at line 124 of file CAFMaker_module.cc.

std::string dunemva::CAFMaker::fMVASelectLabel

private

Definition at line 76 of file CAFMaker_module.cc.

std::string dunemva::CAFMaker::fMVASelectNueLabel

private

Definition at line 77 of file CAFMaker_module.cc.

std::string dunemva::CAFMaker::fMVASelectNumuLabel

private

Definition at line 78 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fNuMomX

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fNuMomY

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fNuMomZ

private

Definition at line 104 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fNuPDG

private

Definition at line 103 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fNuPDGunosc

private

Definition at line 103 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fNwgt[knShifts]

private

Definition at line 93 of file CAFMaker_module.cc.

float dunemva::CAFMaker::fOscPro

private

Definition at line 87 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fQ2

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fRecoHadEnNue

private

Definition at line 113 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fRecoHadEnNumu

private

Definition at line 117 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fRecoLepEnNue

private

Definition at line 112 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fRecoLepEnNumu

private

Definition at line 116 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fRecoMethodNue

private

Definition at line 114 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fRecoMethodNumu

private

Definition at line 118 of file CAFMaker_module.cc.

std::string dunemva::CAFMaker::fRegCNNLabel

private

Definition at line 81 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fRegCNNNueE

private

Definition at line 134 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fRun

private

Definition at line 101 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fSubrun

private

Definition at line 101 of file CAFMaker_module.cc.

systtools::provider_list_t dunemva::CAFMaker::fSystProviders

private

Definition at line 139 of file CAFMaker_module.cc.

int dunemva::CAFMaker::fTrackMomMethodNumu

private

Definition at line 120 of file CAFMaker_module.cc.

TTree* dunemva::CAFMaker::fTree

private

Definition at line 89 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fW

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fWeight

private

Definition at line 88 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fWgts[knShifts][kmaxRwgts]

private

Definition at line 95 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fX

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::fY

private

Definition at line 104 of file CAFMaker_module.cc.

double dunemva::CAFMaker::meta_pot

private

Definition at line 136 of file CAFMaker_module.cc.

int dunemva::CAFMaker::meta_run

private

Definition at line 137 of file CAFMaker_module.cc.

int dunemva::CAFMaker::meta_subrun

private

Definition at line 137 of file CAFMaker_module.cc.

int dunemva::CAFMaker::meta_version

private

Definition at line 137 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nEM

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nK0

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nKm

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nKp

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nN

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nNucleus

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nOtherHad

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nP

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nPi0

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nPim

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nPip

private

Definition at line 106 of file CAFMaker_module.cc.

int dunemva::CAFMaker::nUNKNOWN

private

Definition at line 106 of file CAFMaker_module.cc.

double dunemva::CAFMaker::vtx_x

private

Definition at line 108 of file CAFMaker_module.cc.

double dunemva::CAFMaker::vtx_y

private

Definition at line 108 of file CAFMaker_module.cc.

double dunemva::CAFMaker::vtx_z

private

Definition at line 108 of file CAFMaker_module.cc.

---

The documentation for this class was generated from the following file:

• duneana/duneana/CAFMaker/CAFMaker_module.cc