TruthCuts.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/Cut.h"

namespace ana
{
  /// Is this a Neutral %Current event?
  extern const Cut kIsNC;
  extern const Cut kIsNue;
  extern const Cut kIsNumu;

  extern const NuTruthCut kIsTau_NT;
  extern const Cut kIsTau;

  extern const NuTruthCut kIsNC_NT;

  //----------------------------------------------------------------------
  /// Helper for defining true CC event cuts
  class CCFlavSel
  {
  public:
    CCFlavSel(int pdg, int pdgorig) : fPdg(pdg), fPdgOrig(pdgorig)
    {
    }

    bool operator()(const caf::SRProxy* sr) const;
  protected:
    int fPdg, fPdgOrig;
  };

  // Finally, the function argument to the Cut constructor can be a "functor"
  // object (one with operator()). This allows similar logic but with different
  // constants to be easily duplicated.

  /// Select CC \f$ \nu_\mu\to\nu_e \f$
  const Cut kIsSig    (CCFlavSel(12, 14));
  /// Select CC \f$ \nu_\mu\to\nu_\mu \f$
  const Cut kIsNumuCC (CCFlavSel(14, 14));
  /// Select CC \f$ \nu_e\to\nu_e \f$
  const Cut kIsBeamNue(CCFlavSel(12, 12));
  /// Select CC \f$ \nu_e\to\nu_\mu \f$
  const Cut kIsNumuApp(CCFlavSel(14, 12));
  /// Select CC \f$ \nu_\mu\to\nu_\tau \f$
  const Cut kIsTauFromMu(CCFlavSel(16, 14));
  /// Select CC \f$ \nu_e\to\nu_\tau \f$
  const Cut kIsTauFromE(CCFlavSel(16, 12));

  extern const NuTruthCut kIsNumuCC_NT;
  extern const NuTruthCut kIsBeamNue_NT;
  extern const NuTruthCut kIsSig_NT;

  /// Is this truly an antineutrino?
  extern const Cut kIsAntiNu;
  extern const Cut kIsNu;

  extern const Cut kHasNu;

  extern const NuTruthCut kIsAntiNu_NT;

  /// Check if MC slice has neutrino information (useful for in-and-out tests)
  extern const Cut kHasNeutrino;

  /// Select truly quasielastic events
  extern const NuTruthCut kIsQE_NT;
  extern const Cut kIsQE;

  extern const NuTruthCut kIsRes_NT;
  extern const Cut kIsRes;

  extern const NuTruthCut kIsDIS_NT;
  extern const Cut kIsDIS;

  extern const NuTruthCut kIsCoh_NT;
  extern const Cut kIsCoh;

  //// Select neutrinos by their parents type
  extern const Cut kIsNumuFromKaon;
  extern const Cut kIsNumuFromPion;

  extern const Cut kIsPion;
  extern const Cut kIsNumuFromMuon;

  extern const Cut kIsKaon;
  extern const Cut kIsMuon;

  extern const Cut kIsNueFromKaon;
  extern const Cut kIsNueFromPion;
  extern const Cut kIsNueFromMuon;

  extern const NuTruthCut kIsDytmanMEC_NT;
  extern const Cut kIsDytmanMEC;

  extern const Cut kIsVtxCont;
  extern const NuTruthCut kIsVtxCont_NT;
}