BPFVars.h

Go to the documentation of this file.

#pragma once

// Definitions for specific BPF variables.

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

#include "CAFAna/Vars/CVNProngVars.h"

namespace ana
{
  const Cut kBPFValidMuon = (kCVNMuonIdx >= 0);

  //=============================//
  // BPF single-variables:
  //=============================//

  /// \brief:   Best muon PID for the slice
  ///
  /// Event-level muon PID from the highest scoring BPF muon track, using CVN-prong as a PID.
  ///
  /// Suggested range:  complete coverage from 0 -> 1
  ///
  // const Var kBPFBestMuonPID - Depreciated Var, please use kCVNBestMuonScore in CVNProngVars.h



  /// \brief:   Number of 3D Prongs
  ///
  /// Reconstructed number of hits (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonNpng;

  /// \brief:   Number of hits in the best muon track
  ///
  /// Reconstructed number of hits (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonNhits;

  /// \brief:   Reco Start X-position for the best muon track
  ///
  /// Reconstructed Start X-position (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonStartx;

  /// \brief:   Reco Start Y-position for the best muon track
  ///
  /// Reconstructed Start Y-position (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonStarty;

  /// \brief:   Reco Start Z-position for the best muon track
  ///
  /// Reconstructed Start Z-position (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonStartz;

  /// \brief:   Reco End X-position for the best muon track
  ///
  /// Reconstructed End X-position (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonStopx;

  /// \brief:   Reco End Y-position for the best muon track
  ///
  /// Reconstructed End Y-position (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonStopy;

  /// \brief:   Reco End Z-position for the best muon track
  ///
  /// Reconstructed End Z-position (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonStopz;

  /// \brief:   Reco Start X-direction for the best muon track
  ///
  /// Reconstructed Start X-direction (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonDirx;

  /// \brief:   Reco Start Y-direction for the best muon track
  ///
  /// Reconstructed Start Y-direction (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonDiry;

  /// \brief:   Cosine of Reco Start Y-direction for the best muon track
  ///
  /// Cosine of the reconstructed Start Y-direction (under the muon assumption)
  /// for the BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonCosDirY;

  /// \brief:   Reco Start Z-direction for the best muon track
  ///
  /// Reconstructed Start Z-direction (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonDirz;

  /// \brief:   Reco X-momentum for the best muon track
  ///
  /// Reconstructed X-momentum (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonPx;

  /// \brief:   Reco Y-momentum for the best muon track
  ///
  /// Reconstructed Y-momentum (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonPy;

  /// \brief:   Reco Z-momentum for the best muon track
  ///
  /// Reconstructed Z-momentum (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonPz;

  /// \brief:   Cal energy for the best muon track
  ///
  /// Calorimetric energy (under the muon assumption) for the
  /// BPF track with the highest muon PID.
  extern const Var kBPFBestMuonCalE;

  /// \brief:   Track energy for the best muon track
  ///
  /// BPF track energy (under the muon assumption) for the
  /// BPF track with the highest muon PID
  extern const Var kBPFBestMuonTrkE;

  /// \brief:   Reco length for the best muon track
  ///
  /// Reconstructed track length (under the muon assumption) for the
  /// BPF track with the highest muon PID.
  extern const Var kBPFBestMuonTrkLen;

  /// \brief:   Number of missing planes for best muon track
  /// 
  /// Number of missing planes (under the muon assumption) for the
  /// BPF track with the highest muon PID
  extern const Var kBPFBestMuonTrkNPlaneGap;

  /// \brief:   Number of planes for best muon track
  /// 
  /// Number of planes (under the muon assumption) for the
  /// BPF track with the highest muon PID
  extern const Var kBPFBestMuonTrkNPlanes;


  /// \brief:   Pt/P for Best Muon track
  ///
  extern const Var kBPFBestMuonPtP;

  /// \brief:   Cos(NuMI beam dir) w.r.t. Best BPF Muon track
  ///
  extern const Var kBPFBestMuonCosNumi;

  /// \brief:   EXAMPLE VAR:  longest track length under the muon assumption
  ///
  /// Reconstructed track length (under the muon assumption) for the
  /// longest BPF track in the slice. This varibale is provided
  /// partially as an example of how to loop over all prongs & tracks.
  extern const Var kBPFLongestMuonTrkLen;

  /// \brief:   chi2T for the best muon track
  ///
  /// Input variable "chi2T" (total track fit chi^2) to the muon PID from the best muon track.
  ///
  /// Suggested range:  0 -> 20 (although values can occur outside of this range)
  extern const Var kBPFBestMuonChi2T;

  /// \brief:   dEdXLL for the best muon track
  ///
  /// Input variable "dEdXLL" (dE/dx log-likelihood) to the muon PID from the best muon track.
  ///
  /// Suggested range:  -10 -> 0 (although values can occur outside of this range)
  extern const Var kBPFBestMuondEdXLL;

  /// \brief:   hitRatio for the best muon track
  ///
  /// Input variable "hitRatio" to the muon PID from the best muon track.
  ///
  /// Suggested range:  complete coverage from 0 -> 1
  extern const Var kBPFBestMuonHitRatio;



  /// \brief:   Position at the MuCat transition for the best muon track
  ///
  /// Interpolated position at the MuCat transition (under the muon assumption) for the
  /// BPF track with the highest CVN-prong muon PID.
  extern const Var kBPFBestMuonPosAtTrans;

  /// \brief:   Projected number of cells from endpoint forward to detector edge in ND for the best muon track
  ///
  extern const Var kBPFBestMuonFwdCellND;

  /// \brief:   Projected number of cells from endpoint forward to detector edge for the best muon track
  ///
  extern const Var kBPFBestMuonFwdCell;

  /// \brief:   Projected number of cells from endpoint backward to detector edge in ND for the best muon track
  ///
  extern const Var kBPFBestMuonBakCellND;

  /// \brief:   Projected number of cells from endpoint backward to detector edge for the best muon track
  ///
  extern const Var kBPFBestMuonBakCell;


  //=============================//
  // BPF multi-variables:
  //=============================//



  /// \brief:   EXAMPLE MULTIVAR:  reco E for all "pion-like" tracks
  ///
  /// This multivar is provided as an example of how to plot a value for a
  /// set of tracks that match user defined criteria. The example is to plot
  /// the BPF reconstructed energy under the pion assumption for all "pion-like"
  /// tracks. "pion-like" is definied simply as all tracks < 2 m long and NOT
  /// identified as the muon by the BPF muon PID.
  extern const MultiVar kBPFPionRecoE;

  /// \brief:   track length for all BPF tracks under the muon assumption
  ///
  /// The length of every track in the slice made using the muon assumption.
  extern const MultiVar kBPFAllTrkLenMuon;

  /// \brief:   track length for all BPF tracks under the pion assumption
  ///
  /// The length of every track in the slice made using the pion assumption.
  extern const MultiVar kBPFAllTrkLenPion;

  /// \brief:   track length for all BPF tracks under the proton assumption
  ///
  /// The length of every track in the slice made using the proton assumption.
  extern const MultiVar kBPFAllTrkLenProton;

} // namespace