NtpSREvent.h

Go to the documentation of this file.

//
// NtpSREvent
//
// NtpSREvent is an ntuple class to hold event data
//
// Based on Roy Lee's CandEventSR/EventSRTTree class

#ifndef NTPSREVENT_H
#define NTPSREVENT_H

#include <iosfwd> // ostream fwd decl'
#include "CandNtupleSR/NtpSRVertex.h"
#include "CandNtupleSR/NtpSRPlane.h"
#include "CandNtupleSR/NtpSRBleach.h"
#include "CandNtupleSR/NtpSRWindow.h"
#include "CandNtupleSR/NtpSRStripPulseHeight.h"

class NtpSREvent;
std::ostream &operator << (std::ostream& os, const NtpSREvent& event);

class NtpSREvent : public TObject {

 public:
  NtpSREvent(): index(0),slc(0),ndigit(0),nstpcnt(0),nstrip(0),stp(0),
       stpph0sigmap(0),stpph0mip(0),stpph0gev(0),stpph1sigmap(0),
       stpph1mip(0),stpph1gev(0),nshower(0),shw(0),ntrack(0),trk(0),
       contained(1) {} // def const'r
  NtpSREvent(Int_t nstripinit,Int_t nshowerinit, Int_t ntrackinit); 
  virtual ~NtpSREvent() { 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 ind);
  void AddShowerAt(Int_t showerindex,Int_t ind);
  void AddTrackAt(Int_t trackindex,Int_t ind);

  void Clear(Option_t* = "") { ClearStrips(); 
                               nshower = 0; if ( shw ) delete [] shw; shw = 0;
                               ntrack = 0; if ( trk ) delete [] trk; trk = 0; }
  void ClearStrips() { nstrip = 0; if ( stp ) delete [] stp; stp = 0; }

  void SetPh(Float_t sigmap,Float_t mip,Float_t gev,
             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 event in event array
  Short_t slc;    // index of associated slice in slice array
  Int_t ndigit;   // number of digits in event
  // 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
  // array (stp) has been deleted.
  Int_t nstpcnt;  // number of strips in event
  Int_t nstrip;   // number of strips in strip array
  Int_t* stp;     //[nstrip] indices to associated strips in strip array
  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)
  Int_t nshower;  // number of showers in event
  Int_t* shw;     //[nshower] indices to associated showers in shower array
  Int_t ntrack;   // number of tracks in event
  Int_t* trk;     //[ntrack] indices to associated tracks in track array
  Int_t primshw;  // index to primary shower in shw array
  // Further description of primshw from Jim's commit notes:
  // Normally, this index is zero (primary is first element of array).  
  // However, if the index is less than 0, no shower is identified as the 
  // primary.  In the current primary shower algorithm, this means that 
  // there are no showers either with energies greater than 2 GeV, or, if 
  // a track exists, with a vertex within 0.5 m of the track vertex.
  Int_t primtrk;  // index to primary track in trk array
  Int_t contained;

  NtpSRStripPulseHeight ph; // summed strip pulse height
  NtpSRPlane    plane;    // range of planes 
  NtpSRVertex   vtx;      // vertex  
  NtpSRVertex   end;      // event end point, set to end of primary track
  NtpSRBleach   bleach;   // values used to cut low ph near detector evts
  NtpSRWindow   win;      // total+unassigned charge info in pln+time window around evt

  ClassDef(NtpSREvent,8)
};

inline void NtpSREvent::SetPh(Float_t sigmap, Float_t mip, Float_t gev,
                               UInt_t istrip, UShort_t iend) {
  if ( (Int_t)istrip >= nstrip ) return;
  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 // NTPSREVENT_H

---