NtpMREvent.h

Go to the documentation of this file.

//
// NtpMREvent
//
// NtpMREvent is an ntuple class to hold MR track data
//

#ifndef NTPMREVENT_H
#define NTPMREVENT_H

#include <iosfwd> // ostream fwd decl'
#include "TObject.h"

class NtpMREvent;
std::ostream &operator << (std::ostream& os, const NtpMREvent& mrevt);

class NtpMREvent : public TObject {

 public:
  NtpMREvent(); //def const'r
  NtpMREvent(Int_t ninit); // normal const'r
  virtual ~NtpMREvent() { this -> Clear(); }

  // State testing methods
  virtual std::ostream& Print(std::ostream& os) const;
  virtual void Print(Option_t* option = "") const; 

  // State changing methods
  void AddStripAt(Int_t stripindex, Int_t mrinfo, Int_t ind);
  void Clear(Option_t* = "") { 
    if(stp) delete [] stp; stp=0; 
    if(mrstp) delete [] mrstp; mrstp=0;
  }

  Bool_t  StpIsTrueMu(Int_t stripindex);    // true muon hit
  Bool_t  StpIsTrueShw(Int_t stripindex);   // true shower hit
  Bool_t  StpIsScaled(Int_t stripindex);    // Scaled by MR procedure
  Bool_t  StpIsInRecoTrk(Int_t stripindex); // Is in reco'd track
  Float_t StpScaleFactor(Int_t stripindex); // scaled by...
  Bool_t  StpIsRetained(Int_t stripindex);  // retained from original event
  Bool_t  StpIsMCElec(Int_t stripindex);    // MC electron added in

 public:
  // Ntuple is treated like a C-struct with public data members and
  // rule-breaking field data members not prefaced by "f" and all
  // lowercase, by popular demand.
  
  UShort_t index; // index of this mr track
  Int_t ndigit;   // number of kept digits in original event
  Int_t nstrip;   // number of kept strips in original event
  Int_t* stp;       //[nstrip]  indices to associated strips in strip array
  Int_t* mrstp;     //[nstrip]  information about strips retained
  Int_t orig_event;  // original event index track was removed from
  Int_t best_event;   // index of best matched event in new event array 
  Float_t best_purity;   // purity of kept event fragment in best_event 
  Float_t best_complete; // completeness of kept event fragment in best_event
  Float_t best_purity_phw;   // PH purity of kept event fragment in best_event 
  Float_t best_complete_phw; // PH completeness of kept event fragment in best_event

  Float_t elec_complete;     // completeness of electron only
  Float_t elec_complete_phw; // PH completeness of electron only

  Float_t best_elec_complete_phw;  
  Float_t best_elec_complete; // sometimes the remnant and elec get split
  Float_t best_elec_event;    //   this tracks where the electron wants to be

  Float_t comb_complete;
  Float_t comb_complete_phw;
  Float_t best_comb_complete;
  Float_t best_comb_complete_phw;
  Float_t best_comb_event;

  // added by J. Ling
  // details of the shower
  Float_t shwvtxx;  // shower vertex x
  Float_t shwvtxy;  // shower vertex y
  Float_t shwvtxz;  // shower vertex z
  Float_t shwendx;  // shower end x
  Float_t shwendy;  // shower end y
  Float_t shwendz;  // shower end z
  Int_t   shwvtxplane;  // shower vertex plane
  Int_t   shwendplane;  // shower end plane
  Int_t   shwnplane;    // shower plane length
  Float_t   shwcharge;    // shower charge, mip
  
  //details of removed track:
  Float_t vtxx;  //track vertex x
  Float_t vtxy;  //track vertex y
  Float_t vtxz;  //track vertex z
  Float_t vtxdistance;  // track vetext distance to the edge
  Float_t endx;  //track end x
  Float_t endy;  //track end y
  Float_t endz;  //track end z
  Float_t enddistance;  // track end distance to the edge
  Int_t vtxp;  //track vertex plane
  Int_t endp;  //track end plane
  Int_t   npln;  //track end - begin plane + 1
  Float_t prng;  //track mom from range
  Float_t pcrv;  //track mom from curvature
  Float_t pvdx;  //track vertex dir cos x
  Float_t pvdy;  //track vertex dir cos y
  Float_t pvdz;  //track vertex dir cos z
  Float_t zenith;  // Zenith angle  
  Float_t azimuth; // azimuth angle
  Int_t   fitp;  //track fit pass
  Int_t   endc;  //track contained?
  Int_t   pass;  //is track contained or, if not, fit is good
  Float_t pmux;  //track x momentum (using appropriate momentum)
  Float_t pmuy;  //track y momentum (using appropriate momentum)
  Float_t pmuz;  //track z momentum (using appropriate momentum)
  Int_t   mxpl;  //max track plane (beyond which all hits are killed by MR)
  Double_t qp;   //q over pcurve
  //variables for MRM study
  Float_t mrmpmux;   //true muon x momentum
  Float_t mrmpmuy;   //true muon y momentum
  Float_t mrmpmuz;   //true muon z momentum
  Float_t mrmQ2;     //converted Q2
  Float_t mrmEshw;   //converted Eshw

  ClassDef(NtpMREvent,6)
};

#endif // NTPMREVENT_H

---