Vars.h

Go to the documentation of this file.

#pragma once

// Specific variables. When this file gets too crowded it should be split up

// See Cuts.h for a walkthrough of how to define this kind of object

#include "CAFAna/Core/Var.h"
#include "CAFAna/Core/Cut.h"

#include "StandardRecord/Proxy/FwdDeclare.h"

class TH1D;

namespace ana
{
  extern const Cut kIsRHC;

  // Return the run number
  extern const Var kRun;

  // Return the subrun number
  extern const Var kSubrun;

  // Return the cycle number (MC only) 
  extern const Var kCycle;

  // Return the event number
  extern const Var kEvt;

  // Return the slice number
  extern const Var kSlc;

  /// See docdb 23597
  double CalibrationBugCorrectionFactor(const caf::SRHeaderProxy& hdr);

  /// \f$ \nu_e \f$ PID
  extern const Var kLEM;

  /// \f$ \nu_e \f$ PID
  extern const Var kLID;

  /// \f$ \nu_e \f$ PID
  extern const Var kRVP;

  /// \f$ \nu_e \f$ PID
  extern const Var kCVNe;

  extern const Var kCVNe_looseptp;
  extern const Var kCVNe_oldpresel;

  /// 2018 nue PID
  extern const Var kCVNSSe;

  /// \f$ \nu_\mu \f$ PID
  extern const Var kCVNm;

  extern const Var kCVNm_looseptp;
  extern const Var kCVNm_oldpresel;

  /// \f$ \nu_\tau \f$ PID
  extern const Var kCVNt;

  /// \f$ \nu_nc \f$ PID
  extern const Var kCVNnc;

  /// \f$ e- elastic scattering PID - I don't know what the e elastic number is for the final state labels.
  //  Leaving this commented out for now...
  // const Var kCVNeel([](const caf::SRProxy* sr){ return sr->sel.cvn.output[???]; });
  /// \f$ cosmic PID - use this for the final state era of >= prod4 files.
  extern const Var kCVNcos_looseptp;
  extern const Var kCVNcos_oldpresel;

  /// The result of CosVeto
  extern const Var kCosVeto;

  /// Summed calorimetric energy of all hits in slice, uncorrected
  /** As of revision 9885 (2014-06-29), the CalE scale correction is applied in
      CAFMaker, not CAFAna.  For cafs produced with tags after June 29, 2014,
      use this variable instead. */
  extern const Var kCaloE;

  extern const Var kMeanTime;

  /// Number of SlcME's assoc with slice
  extern const Var kNMichels;

  extern const Var kNHit;

  extern const Var kLongestProng;

  /// \f$ \nu_\mu \f$ PID
  extern const Var kRemID;

  /// True neutrino energy
  extern const NuTruthVar kTrueE_NT;
  extern const Var kTrueE;

  /// Neutrino interaction mode
  extern const Var kMode;

  ///True neutrino ancestor (off-target) pT
  extern const Var kTruetpT;

  ///True neutrino ancestor (off-target) pz
  extern const Var kTruetpz;
  
  //True parton Z - energy of most energetic hadron divided by the energy of the hadronic system. 
  extern const Var kPartonZ;

   class GetWtMinervaND
   {
   public:
     GetWtMinervaND(){ }
     double operator() (const caf::SRProxy* sr) const;
   protected:
     void InitializeMinervaND() const;
     mutable TH1D* fRatiosm[4] = {0};
   };

  class GetWtMippNA49NOvAND
   {
   public:
     GetWtMippNA49NOvAND(){ }
     double operator() (const caf::SRProxy* sr) const;
   protected:
     void InitializeMippNA49NOvAND() const;
     mutable TH1D* fRatiosn[4] = {0};
   };
  class GetWtMippNA49NOvAFD
   {
   public:
     GetWtMippNA49NOvAFD(){ }
     double operator() (const caf::SRProxy* sr) const;
   protected:
     void InitializeMippNA49NOvAFD() const;
     mutable TH1D* fRatiosf[4] = {0};
   };

  const Var kNDMinervaRwtVar = Var(GetWtMinervaND());//Minerva wt only for ND

 const Var kNDMippNA49NOvARwtVar = Var(GetWtMippNA49NOvAND());//MiPP+NA49 FROM NOVA wts only for ND
 const Var kFDMippNA49NOvARwtVar = Var(GetWtMippNA49NOvAFD());//Mipp+Na49 From NOvA wts only for FD
 const Var kMippNA49NOvARwtVar = kNDMippNA49NOvARwtVar * kFDMippNA49NOvARwtVar;

} // namespace