ChopEvaluation.h

Go to the documentation of this file.

#ifndef CHOPEVALUATION_H
#define CHOPEVALUATION_H


#include "CandDigit/CandDigitHandle.h"
#include "DataUtil/Truthifier.h"
#include <TObject.h>
#include <vector>

static const int kNProf = 15;

class profile_element{ 
 public:                                        
  int i;
  float p;
};

class ChopEvaluation //: public TObject
{
 public:
  int   fTdcBeg; // TDC time start
  int   fTdcEnd; // TDC time stop
  int   fPlaneBeg; // start plane
  int   fPlaneEnd; // end plane

  double fSigcor;
  double fSigcor_cal; // Calorimeter
  double fSigcor_calt; // Calorimeter
  double fSigcor_fid; // Calorimeter, "fiducial" volume only.
  double fMeanTime; // charge-weighted
  int    fBigTdc; // Biggest bucket in t-bins
  double fBigTdcE; // Size of biggest bucket
  int    fBigPlane; // Biggest bucket in t-bins
  double fBigPlaneE; // Size of biggest bucket
  int    fPlanes; // Total planes hit
  int    fStrips; // total strips hit
  int    fCalStrips; // total strips hit
  int    fDigits; // Total digits
  
  int    fBestNeutrino1;
  int    fBestNeutrino2;

  profile_element fProfile[kNProf];
  
  // work elements
  std::vector<double> fTdcProfile;      
  std::vector<double> fPlaneProfile;    
  std::vector<double> fNuMatch;         
  std::vector<double> fNuMatch_cal;     
  std::vector<double> fNuMatch_fid;     
  std::vector<double> fNuMatch_calt;    
  std::vector<double> fNuMatch_caltfid; 
  std::vector<double> fNuMatch_strips;  
  std::vector<double> fNuMatch_calstrips;  


  ChopEvaluation();
  void Reset();

  void Evaluate(
                const std::vector<CandDigitHandle>& inList, // digits to evaluate
                const Truthifier& truth,                    // truth.
                const std::vector<int>&    inNeutrinos,     // indecies of neutrinos in truth info
                const std::vector<double>& inNeutrinoTZero, // time of first interaction of each nu
                double inTCut
                );

  //ClassDef(ChopEvaluation,1);
};



class ChopNeutrinoInfo
{
 public:
  int inu;
  float vtxt;
  float vtxtdc;
  float vtxx;
  float vtxy;
  float vtxz;

  float enu;   // neutrino energy
  float emu; // muon
  float eem; // em
  float eshw; // shower
  float evis;  // shower+muon+em
  int   iaction;
  float X;
  float Y;
  float tzero; // time of first hit

  ChopNeutrinoInfo() {
    inu=iaction=-1;
    vtxt=vtxtdc=vtxx=vtxy=vtxz=999999;
    enu=emu=eem=eshw=evis=X=Y=tzero=0;
  }
  //ClassDef(ChopNeutrinoInfo,1);
};

class ChopNeutrinoLeaf : public ChopEvaluation, public ChopNeutrinoInfo
{
 public:
    // Truth:
  int   snarl;
  float viscal; // of neutrino
  float viscalt;
  float visfid;
  float visfidt;
  float sigcorr;
  float strips;
  float calstrips;


  int   chop1;
  float completet1;
  float complete1;
  float complete_tot1;
  float strips1;
  float calstrips1;
  float stripsmatch1;
  float calstripsmatch1;
  float chop1_viscal;
  float chop1_viscalt;
  float chop1_visfid;
  float chop1_sigcorr;
  int   chop1_nu;
  int   chop1_tdc;
  int   chop1_planebeg;
  int   chop1_planeend;


  int   chop2;
  float completet2;
  float complete2;
  float complete_tot2;
  float strips2;
  float calstrips2;
  float stripsmatch2;
  float calstripsmatch2;
  float chop2_viscal;
  float chop2_viscalt;
  float chop2_visfid;
  float chop2_sigcorr;
  int   chop2_nu;
  int   chop2_tdc;
  int   chop2_planebeg;
  int   chop2_planeend;



  ChopNeutrinoLeaf() { Reset(); };
  void Reset() {
    chop1=chop2=-1;
    complete1=complete2=0;
    completet1=completet2=0;
    strips1 = strips2 = 0;
    stripsmatch1 = stripsmatch2 = 0;
    chop1_viscal = chop2_viscal=0;
    chop1_visfid = chop2_visfid = 0;
    chop1_nu = chop2_nu = 0;
    chop1_tdc = chop2_tdc = 0;
  }

};




class ChopTreeLeaf : public ChopEvaluation, public ChopNeutrinoInfo
{
 public:
  int snarl;
  int chop;

  // Evalution:
  float complete1; // Completeness, fraction, of biggest nu event.
  float viscal1; // of neutrino
  float viscalt1;
  float visfidt1;
  float strips1;

  float complete2; // Completeness, fraction, of biggest nu event.
  float viscal2; // of neutrino
  float viscalt2;
  float visfidt2;
  float strips2;

  ChopTreeLeaf() { Reset(); };
  void Reset() {
    snarl = 0;
  }
};


#endif

---