NtpSRShower.h

Go to the documentation of this file.

//
// NtpSRShower
//
// NtpSRShower is an ntuple class to hold shower data
//
// Based Roy Lee's CandEventSR/ShowerSRTTree class

#ifndef NTPSRSHOWER_H
#define NTPSRSHOWER_H

#include <iosfwd> // ostream fwd decl'

#include "CandNtupleSR/NtpSRVertex.h"
#include "CandNtupleSR/NtpSRPlane.h"
#include "CandNtupleSR/NtpSRShowerPulseHeight.h"
#include "CandNtupleSR/NtpSRSubShowerSummary.h"

class NtpSRShower;
std::ostream &operator << (std::ostream& os, const NtpSRShower& shower);

class NtpSRShower : public TObject {

 public:
  NtpSRShower(): index(0),slc(0),ndigit(0),nstpcnt(0),nstrip(0),stp(0),
                        stpu(0),stpv(0),stpx(0),stpy(0),stpz(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),
                        nUcluster(0),nVcluster(0),ncluster(0),clu(0),contained(0) {} // def const'r
  NtpSRShower(Int_t nstripinit,Int_t nclusterinit=0); // normal const'r

  virtual ~NtpSRShower() { 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 AddClusterAt(Int_t clusterindex, Int_t ind);
  void Clear(Option_t* = "") { ClearStrips(); ClearClusters(); }
  void ClearStrips();
  void ClearClusters() { ncluster = 0; if ( clu ) delete [] clu; clu = 0; }

  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);

 private:
  // Avoid copy and assignment ctor by making them private and w/no imp
  NtpSRShower(const NtpSRShower& that);
  NtpSRShower& operator=(const NtpSRShower& that);
  
 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 shower in shower array
  Short_t slc;    // index of associated slice in slice array
  Int_t ndigit;   // number of digits in shower
  // 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 shower
  Int_t nstrip;   // number of strips in strip array
  Int_t* stp;     //[nstrip]  indices to associated strips in strip array
  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 *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)
  Short_t nUcluster; // number of U clusters in shower
  Short_t nVcluster; // number of V clusters in shower
  Int_t  ncluster; // number of clusters in shower
  Int_t* clu;     //[ncluster]  indices to associated clusters in cluster array
  Int_t contained;
  NtpSRStripPulseHeight ph; // summed strip pulse height 
  NtpSRShowerPulseHeight shwph; // shower energy members 
  NtpSRPlane  plane;        // range of planes 
  NtpSRVertex vtx;          // vertex - not filled? 
  NtpSRSubShowerSummary sss; //subshower summary block

  ClassDef(NtpSRShower,12)
};

inline void NtpSRShower::SetTime(Double_t time, UInt_t istrip, UShort_t iend) {
  if ( (Int_t)istrip >= nstrip ) return; 
  if      ( iend == 0 ) stpt0[istrip] = time;
  else if ( iend == 1 ) stpt1[istrip] = time;
  return;
}
inline void NtpSRShower::SetCalT0(Double_t calt0, UInt_t istrip,UShort_t iend) {
  if ( (Int_t)istrip >= nstrip ) return; 
  if      ( iend == 0 ) stptcal0t0[istrip] = calt0;
  else if ( iend == 1 ) stptcal1t0[istrip] = calt0;
  return;
}
inline void NtpSRShower::SetAttnC0(Double_t attnc0,UInt_t istrip,UShort_t iend){
  if ( (Int_t)istrip >= nstrip ) return; 
  if      ( iend == 0 ) stpattn0c0[istrip] = attnc0;
  else if ( iend == 1 ) stpattn1c0[istrip] = attnc0;
  return;
}
inline void NtpSRShower::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 // NTPSRSHOWER_H

---