SRNumuEnergy.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// \file    SRNumuEnergy.h
// \version $Id: SRNumuEnergy.h,v 1.1 2012-11-19 19:45:16 rocco Exp $
// \author  $Author: rocco $
// \date    $Date: 2012-11-19 19:45:16 $
////////////////////////////////////////////////////////////////////////
#ifndef SRNUMUENERGY_H
#define SRNUMUENERGY_H

#include "StandardRecord/SRHadClust.h"
#include "StandardRecord/SRTrueNumuEnergy.h"
#include "StandardRecord/SRBPFEnergy.h"

namespace caf
{
  /// Numu energy estimator output
  class SRNumuEnergy
    {
    public:
      SRNumuEnergy();
      virtual ~SRNumuEnergy();

      float      E;              ///< Neutrino energy, set to match trkccE [GeV]
      float      calccE;         ///< Calorimetric charged current neutrino energy [GeV]
      float      trkqeE;         ///< Track length quasielastic neutrino energy [GeV]
      float      trknonqeE;      ///< Track length non-quasielastic neutrino energy [GeV]
      float      trkccE;         ///< Track length cc neutrino energy [GeV]
      float      shiftedtrkccE;  ///< Track length cc neutrino energy [GeV] (with hadE = hadE*0.77)
      float      angleE;         ///< Quasielastic angle formula neutrino energy [GeV]
      float      angleerror;     ///< 1-sigma error from using quasielastic angle formula for neutrino energy
      float      recomuonE;      ///< Reconstructed muon energy for all neutrino energy estimators [GeV]
      float      ucrecomuonE1trk; ///< Reconstructed muon energy for uncontained single track events at the FarDet from TMVA [GeV]
      float      ucrecomuonE;     ///< Reconstructed muon energy for uncontained non-single track events at the FarDet from TMVA [GeV]
      float      recotrkcchadE;  ///< Reconstructed hadronic energy for track cc neutrino energy estimator [GeV]
      float      hadcalE;        ///< Hadronic calorimetric energy NOT on the muon track[GeV]
      float      hadtrkE;        ///< Hadronic calorimetric energy on the muon track[GeV]
      float      ndtrklenact;    ///< Near detector -- muon track length in active region [cm]
      float      ndtrklencat;    ///< Near detector -- muon track length in muon catcher  [cm]
      float      ndtrkcalactE;   ///< Near detector -- muon calorimetric energy in active region [GeV]
      float      ndtrkcaltranE;  ///< Near detector -- muon calorimetric energy in transition plane [GeV]
      float      ndtrkcalcatE;   ///< Near detector -- muon calorimetric energy in muon catcher [GeV]
      float      ndhadcalactE;   ///< Near detector -- hadronic calorimetric energy NOT on the muon track in active region [GeV]
      float      ndhadcaltranE;  ///< Near detector -- hadronic calorimetric energy NOT on the muon track in transition plane [GeV]
      float      ndhadcalcatE;   ///< Near detector -- hadronic calorimetric energy NOT on the muon track in muon catcher [GeV]
      float      ndhadtrkactE;   ///< Near detector -- hadronic calorimetric energy on the muon track in active region [GeV]
      float      ndhadtrktranE;  ///< Near detector -- hadronic calorimetric energy on the muon track in transition plane [GeV]
      float      ndhadtrkcatE;   ///< Near detector -- hadronic calorimetric energy on the muon track in muon catcher [GeV]
      float      ndtrktranx;     ///< Near detector -- if muon track crosses transition plane, the x location [cm]
      float      ndtrktrany;     ///< Near detector -- if muon track crosses transition plane, the y location [cm]
      SRHadClust hadclust;       ///< Cluster information for hits in slice but not on muon track (track with higest ReMId value)
      SRTrueNumuEnergy mc;       ///< True information used for retuning numu energy fits
      SRBPFEnergy      bpfenergy;///< BreakPointFitter energy estimator

      /* LSTM Energies */
      float      lstmmuon;      ///< Muon Energy predicted by LSTM EE [GeV]
      float      lstmnu;        ///< NuMu Energy predicted by LSTM EE [GeV]

      float      lstmmuon_opphorn; ///< lstmmuon trained with opposite horn current
      float      lstmnu_opphorn;   ///< lstmnu trained with opposite horn current
      
      float      regcvnhadE;    ///< Hadronic Energy predicted by Regression CNN [GeV] 

      float      regcvnhadE_opphorn; ///< regcvnhadE trained with opposite horn current

      virtual void setDefault();
      void setLSTMDefault(bool opphorn = false);
    };

} // end namespace

#endif // SRNUMUENERGY_H
//////////////////////////////////////////////////////////////////////////////