NtpSRTrack.h

Go to the documentation of this file.

//
// NtpSRTrack
//
// NtpSRTrack is an ntuple class to hold 3-d track data
//
// Based on Roy Lee's CandEventSR/TrackSRTTree class.

#ifndef NTPSRTRACK_H
#define NTPSRTRACK_H

#include <iosfwd> // ostream fwd decl'

#include "CandNtupleSR/NtpSRTrackPlane.h"
#include "CandNtupleSR/NtpSRStripPulseHeight.h"
#include "CandNtupleSR/NtpSRVertex.h"
#include "CandNtupleSR/NtpSRFiducial.h"
#include "CandNtupleSR/NtpSRTrackTime.h"
#include "CandNtupleSR/NtpSRMomentum.h"
#include "CandNtupleSR/NtpSRFitTrack.h"
#include "CandNtupleSR/NtpSRCosmicRay.h"

class NtpSRTrack;
std::ostream &operator << (std::ostream& os, const NtpSRTrack& track);

class NtpSRTrack : public TObject {

 public:
  NtpSRTrack() : index(0),slc(0),ndigit(0),nstpcnt(0),nstrip(0),stp(0),
                 stpfit(0),stpu(0),
                 stpv(0),stpx(0),stpy(0),stpz(0),stpds(0),stpfitchi2(0),
                 stpfitprechi2(0),stpfitqp(0),stpph0sigmap(0),stpph0mip(0),
                 stpph0gev(0),stpph1sigmap(0),stpph1mip(0),stpph1gev(0),
                 stpattn0c0(0),stpattn1c0(0),stpt0(0),stpt1(0),stptcal0t0(0),
                 stptcal1t0(0),ds(0),range(0),cputime(0), contained(0) {} // def const'r
  NtpSRTrack(Int_t nstripinit); // normal const'r
  NtpSRTrack(const NtpSRTrack& o); // copy ctor

  virtual ~NtpSRTrack() { this -> Clear(); }

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

  // State changing methods
  void Clear(Option_t* = "") { ClearStrips(); }
  void ClearStrips();
  
  void AddStripAt(Int_t stripindex,Int_t ind);
  void SetPh(Float_t sigmap,Float_t mip,Float_t gev,
                                        UInt_t istrip,UShort_t iend);
  void SetTime(Double_t time, UInt_t istrip, UShort_t iend);
  void SetCalT0(Double_t calt0, UInt_t istrip, UShort_t iend);
  void SetAttnC0(Double_t attnc0, UInt_t istrip, UShort_t iend);
  

 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 track in track array
  Short_t slc;          // index of associated slice in slice array
  Int_t  ndigit;        // number of digits in track
  // nstpcnt will represent the number of strips in track even in "slim"
  // ntuple output file for which nstrip is set to 0 and the strip 
  // arrays (stp,stpfit, etc.) have been deleted. 
  Int_t  nstpcnt;       // number of strips in track
  Int_t  nstrip;        // number of strips in strip track arrays
  Int_t  *stp;          //[nstrip] indices of assoc. strips in strip array 
  Byte_t *stpfit;       //[nstrip] 1 if hit strip was used in final track fit
  // Track positional information at plane associated with each strip
  Float_t *stpu;        //[nstrip] U-position of track at strip plane (m) 
  Float_t *stpv;        //[nstrip] V-position of track at strip plane (m)
  Float_t *stpx;        //[nstrip] Z-position of track at strip plane (m)
  Float_t *stpy;        //[nstrip] Y-position of track at strip plane (m)
  Float_t *stpz;        //[nstrip] Z-position of track at strip plane (m)
  Float_t *stpds;       //[nstrip] track path length from track end (m)
  Float_t *stpfitchi2;   //[nstrip] final chi2 at strip plane
  Float_t *stpfitprechi2;//[nstrip] pre chi2 from Kalman filter at strip plane
  Float_t *stpfitqp;     //[nstrip] charge over momentum at strip plane (GeV)
  Float_t *stpph0sigmap;//[nstrip] ph corrected for fiber attenuation (east)
  Float_t *stpph0mip;   //[nstrip] ph in MIPs      (east)
  Float_t *stpph0gev;   //[nstrip] ph in GeV       (east)
  Float_t *stpph1sigmap;//[nstrip] ph corrected for fiber attenuation (west)
  Float_t *stpph1mip;   //[nstrip] ph in MIPs      (west)
  Float_t *stpph1gev;   //[nstrip] ph in GeV       (west)
  Float_t *stpattn0c0;  //[nstrip] C0 from mapper results (east)
  Float_t *stpattn1c0;  //[nstrip] C0 from mapper results (west)
  Double_t *stpt0; //[nstrip] 3-d time corr. for propagation delays (sec) (E)
  Double_t *stpt1; //[nstrip] 3-d time corr. for propagation delays (sec) (W)
  Double_t *stptcal0t0; //[nstrip] T0 offset (sec) (east)
  Double_t *stptcal1t0; //[nstrip] T0 offset (sec) (west)
  Float_t ds;    // total path length of track from vertex to end (m)
  Float_t range; // amount of material traversed from vertex to end (g/cm**2)
  Float_t cputime; // CPU time to create track list which contains this track
  Int_t contained; // 0 = not contained, 1 = contained
  NtpSRStripPulseHeight ph; // summed strip pulse height
  NtpSRTrackPlane  plane;   // range of planes, timing used to determine start
  NtpSRVertex  vtx;      // vertex coordinates
  NtpSRVertex  end;      // coordinates at end plane
  NtpSRVertex  lin;      // coordinates from linear fit
  NtpSRFiducial fidvtx;  // containment information for vertex
  NtpSRFiducial fidend;  // containment information for track end
  NtpSRFiducial fidall;  // containment information for distance of closest
                         // approach over all 3D points on the track
  NtpSRTrackTime time;   // track timing data
  NtpSRMomentum momentum; // track momentum data
  NtpSRFitTrack fit;     // data from track fit
  NtpSRCosmicRay cr;     // cosmic ray data

  ClassDef(NtpSRTrack,8)
};

inline void NtpSRTrack::SetTime(Double_t time, UInt_t istrip, UShort_t iend) {
  if      ( iend == 0 ) stpt0[istrip] = time;
  else if ( iend == 1 ) stpt1[istrip] = time;
  return;
}
inline void NtpSRTrack::SetCalT0(Double_t calt0, UInt_t istrip,UShort_t iend) {
  if      ( iend == 0 ) stptcal0t0[istrip] = calt0;
  else if ( iend == 1 ) stptcal1t0[istrip] = calt0;
  return;
}
inline void NtpSRTrack::SetAttnC0(Double_t attnc0,UInt_t istrip,UShort_t iend){
  if      ( iend == 0 ) stpattn0c0[istrip] = attnc0;
  else if ( iend == 1 ) stpattn1c0[istrip] = attnc0;
  return;
}
inline void NtpSRTrack::SetPh(Float_t sigmap, Float_t mip, Float_t gev, 
                              UInt_t istrip, UShort_t iend) {
  if ( iend == 0 ) {
    stpph0sigmap[istrip] = sigmap;
    stpph0mip[istrip] = mip;
    stpph0gev[istrip] = gev;
  }
  else if ( iend == 1 ) {
    stpph1sigmap[istrip] = sigmap;
    stpph1mip[istrip] = mip;
    stpph1gev[istrip] = gev;
  }

  return;
}


#endif // NTPSRTRACK_H

---