SRCosRej.h

Go to the documentation of this file.

#ifndef SRCOSREJ_H
#define SRCOSREJ_H


namespace caf
{
    /// Output from Cosmic Rejection (CosRej) module
    class SRCosRej{
    public:
      SRCosRej();
      ~SRCosRej();

      void setDefault();

      /// \name Track Projections
      //@
      float   mindist;          ///< minimum projected distance to edge of any Kalman track with > 15 hits
      int     mincell;          ///< minimum projected cells to edge of any Kalman track with > 15 hits
      //@}

      /// \name Other Track Information
      //@
      float   anglekal;         ///< cos of angle of best ReMId Kalman track
      float   anglecos;         ///< cos of angle of Cosmic Track
      int     nkal3d;           ///< number of 3D tracks as determined by Kalman tracker
      float   costhetatrue;     ///< cosine of angle between cosmic trk dir and true dir of most contributing particle
      float   kalthetatrue;     ///< cosine of angle between kalman trk dir and true dir of most contributing particle
      int     pdgbest;          ///< pdg code of most contributing true particle to best kalman trk
      //@}


      /// \name PIDs
      //@
      float   numucontpid2020;    ///< cosmic rejection PID for contained events; 2020 Analysis
      float   numucontpid2019;    ///< cosmic rejection PID for contained events; 2019 Analysis
      float   numuunconttunedpid; ///< tuned cosmic rejection PID for uncontained events - Jose's;
      //@}

      /// \name Scattering Variables
      //@
      float   scatt;            ///< sum of all all angular deviation (looping over kaltrk traj pts) / trk len
      float   fscattmax;        ///< maximum scattering angle (Fernanda)
      float   fscattsum;        ///< sum of scattering angles (Fernanda)
      float   fscattext;        ///< scattering variable (Fernanda)
      float   fscattsig;        ///< sigma of scattering variable distribution (Fernanda)
      //@}

      /// \name Energy and Activity
      //@
      float   eratio;           ///< ratio of best kalman track GeV / slice GeV
      float   hadE;             ///< GeV sum of non (Kalman best ReMId) track E
      //@}

      //      track timing fits - disabled until fast timing is implemented
      float   kfitspeed;        ///< fit inverse speed (ns/cm) of Hough timing fit (TimingFit module) for best ReMId Kalman track
      float   cfitspeed;        ///< fit inverse speed (ns/cm) of Hough timing fit (TimingFit module) for cosmic track
      float   kdirscore;        ///< chisq score difference between +c and -c assumptions in Hough timing fit (TimingFit module) for best ReMId Kalman track
      float   cdirscore;        ///< chisq score difference between +c and -c assumptions in Hough timing fit (TimingFit module) for cosmic track
      float   kscorediff;       ///< chisq score difference between best +c/-c assumption and free fit in Hough timing fit (TimingFit module) for best ReMId Kalman track
      float   cscorediff;       ///< chisq score difference between best +c/-c assumption and free fit in Hough timing fit (TimingFit module) for Cosmic track
      float   kalslope;         ///< slope of track timing fit (kalmantrack)
      float   kalchisq;         ///< chisq value of slope fit (kalmantrack)
      float   kalchidiff;       ///< chisq difference of track timing fits (kalmantrack)
      float   cosslope;         ///< slope of track timing fit (cosmictrack)
      float   coschisq;         ///< chisq value of slope fit (cosmictrack)
      float   coschidiff;       ///< chisq difference of track timing fits (cosmictrack)

    };

} // end namespace

#endif // SRCOSREJ_H
//////////////////////////////////////////////////////////////////////////////