NtpSRTrack.cxx

Go to the documentation of this file.

// NtpSRTrack
//
// NtpSRTrack is an ntuple class to hold 3-d track data
//

#include <iostream>
using namespace std;

#include "MessageService/MsgService.h"
#include "CandNtupleSR/NtpSRTrack.h"

ClassImp(NtpSRTrack)

std::ostream& operator << (std::ostream& os, const NtpSRTrack& ntp) 
                                     { return ntp.Print(os); }

CVSID("$Id: NtpSRTrack.cxx,v 1.18 2006/08/03 18:14:19 musser Exp $");

NtpSRTrack::NtpSRTrack(Int_t nstripinit) : index(0),slc(0),ndigit(0),
                 nstpcnt(nstripinit),nstrip(nstripinit),
                 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) {
  //
  // Purpose: Normal constructor
  //

  if ( nstrip ) {
    stp           = new Int_t[nstrip]; 
    stpfit        = new Byte_t[nstrip]; 
    stpu          = new Float_t[nstrip]; 
    stpv          = new Float_t[nstrip]; 
    stpx          = new Float_t[nstrip]; 
    stpy          = new Float_t[nstrip]; 
    stpz          = new Float_t[nstrip]; 
    stpds         = new Float_t[nstrip]; 
    stpfitchi2    = new Float_t[nstrip]; 
    stpfitprechi2 = new Float_t[nstrip]; 
    stpfitqp      = new Float_t[nstrip]; 
    stpph0sigmap  = new Float_t[nstrip]; 
    stpph0mip     = new Float_t[nstrip]; 
    stpph0gev     = new Float_t[nstrip]; 
    stpph1sigmap  = new Float_t[nstrip]; 
    stpph1mip     = new Float_t[nstrip]; 
    stpph1gev     = new Float_t[nstrip]; 
    stpattn0c0    = new Float_t[nstrip]; 
    stpattn1c0    = new Float_t[nstrip]; 
    stpt0         = new Double_t[nstrip]; 
    stpt1         = new Double_t[nstrip]; 
    stptcal0t0    = new Double_t[nstrip]; 
    stptcal1t0    = new Double_t[nstrip]; 
  }
  else {
    MSG("NtpSR",Msg::kWarning) 
      << "NtpSRTrack constructor called with nstrip = 0" << endl;
  }
  
  for (Int_t i = 0; i < nstrip; i++ ) {
       stp[i] = -1;
       stpfit[i] = 1;
       stpu[i] = -999999;
       stpv[i] = -999999;
       stpx[i] = -999999;
       stpy[i] = -999999;
       stpz[i] = -999999;
       stpds[i] = -999999;
       stpfitchi2[i] = 0;
       stpfitprechi2[i] = 0;
       stpfitqp[i] = 0;
       stpph0sigmap[i] = 0.;
       stpph0mip[i] = 0.;
       stpph0gev[i] = 0.;
       stpph1sigmap[i] = 0.;
       stpph1mip[i] = 0.;
       stpph1gev[i] = 0.;
       stpattn0c0[i] = -999999;
       stpattn1c0[i] = -999999;
       stpt0[i] = -999999;
       stpt1[i] = -999999;
       stptcal0t0[i] = -999999;
       stptcal1t0[i] = -999999;
  }

}

NtpSRTrack::NtpSRTrack(const NtpSRTrack& o) // copy ctor
  : TObject(o),
    index(o.index),
    slc(o.slc),
    ndigit(o.ndigit),
    nstpcnt(o.nstpcnt),
    nstrip(o.nstrip),
    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) 
{
  //
  // Purpose: Copy constructor
  //
  if ( nstrip ) {
    stp           = new Int_t[nstrip]; 
    stpfit        = new Byte_t[nstrip]; 
    stpu          = new Float_t[nstrip]; 
    stpv          = new Float_t[nstrip]; 
    stpx          = new Float_t[nstrip]; 
    stpy          = new Float_t[nstrip]; 
    stpz          = new Float_t[nstrip]; 
    stpds         = new Float_t[nstrip]; 
    stpfitchi2    = new Float_t[nstrip]; 
    stpfitprechi2 = new Float_t[nstrip]; 
    stpfitqp      = new Float_t[nstrip]; 
    stpph0sigmap  = new Float_t[nstrip]; 
    stpph0mip     = new Float_t[nstrip]; 
    stpph0gev     = new Float_t[nstrip]; 
    stpph1sigmap  = new Float_t[nstrip]; 
    stpph1mip     = new Float_t[nstrip]; 
    stpph1gev     = new Float_t[nstrip]; 
    stpattn0c0    = new Float_t[nstrip]; 
    stpattn1c0    = new Float_t[nstrip]; 
    stpt0         = new Double_t[nstrip]; 
    stpt1         = new Double_t[nstrip]; 
    stptcal0t0    = new Double_t[nstrip]; 
    stptcal1t0    = new Double_t[nstrip]; 
  }

  for (Int_t i = 0; i < nstrip; i++ ) {
       stp[i] = o.stp[i];
       stpfit[i] = o.stpfit[i];
       stpu[i] = o.stpu[i];
       stpv[i] = o.stpv[i];
       stpx[i] = o.stpx[i];
       stpy[i] = o.stpy[i];
       stpz[i] = o.stpz[i];
       stpds[i] = o.stpds[i];
       stpfitchi2[i] = o.stpfitchi2[i];
       stpfitprechi2[i] = o.stpfitprechi2[i];
       stpfitqp[i] = o.stpfitqp[i];
       stpph0sigmap[i] = o.stpph0sigmap[i];
       stpph0mip[i]    = o.stpph0mip[i];   
       stpph0gev[i]    = o.stpph0gev[i]; 
       stpph1sigmap[i] = o.stpph1sigmap[i]; 
       stpph1mip[i]    = o.stpph1mip[i];
       stpph1gev[i]    = o.stpph1gev[i];
       stpattn0c0[i]   = o.stpattn0c0[i];
       stpattn1c0[i]   = o.stpattn1c0[i];
       stpt0[i]        = o.stpt0[i]; 
       stpt1[i]        = o.stpt1[i]; 
       stptcal0t0[i]   = o.stptcal0t0[i];
       stptcal1t0[i]   = o.stptcal1t0[i];
  }
  ds       = o.ds;
  range    = o.range;
  cputime  = o.cputime;
  ph       = o.ph;
  plane    = o.plane;
  vtx      = o.vtx;             
  end      = o.end;             
  lin      = o.lin;             
  fidvtx   = o.fidvtx;   
  fidend   = o.fidend;   
  fidall   = o.fidall;   
  time     = o.time;     
  momentum = o.momentum; 
  fit      = o.fit;
  cr       = o.cr;
  
}

void NtpSRTrack::AddStripAt(Int_t stripindex, Int_t ind) {
  //
  //  Purpose:  Add stripindex to stp array at position ind.
  //
  //  Arguments: stripindex: index of strip in assoc stp TClonesArray. 
  //             ind: index in stp array, cannot exceed bounds
  //                  0->nstrip-1.
  //
  //  Return:  none.
  //
  //  Contact:   S. Kasahara
  // 

  if ( ind < 0 || ind >= nstrip ) {
    MSG("NtpSR",Msg::kWarning)
      << "Attempt to add strip at index " << ind
      << " outside of array size " << nstrip << " ignored." << endl;
    return;  
  }

  stp[ind] = stripindex;
  return;
    
}

void NtpSRTrack::ClearStrips() {
  //
  //  Purpose:  Clear allocated strip memory. 
  //

  nstrip = 0;
  if ( stp ) delete [] stp; stp = 0;
  if ( stpfit ) delete [] stpfit; stpfit = 0;
  if ( stpu ) delete [] stpu; stpu = 0;
  if ( stpv ) delete [] stpv; stpv = 0;
  if ( stpx ) delete [] stpx; stpx = 0;
  if ( stpy ) delete [] stpy; stpy = 0;
  if ( stpz ) delete [] stpz; stpz = 0;
  if ( stpds ) delete [] stpds; stpds = 0;
  if ( stpfitchi2 ) delete [] stpfitchi2; stpfitchi2 = 0;
  if ( stpfitprechi2 ) delete [] stpfitprechi2; stpfitprechi2 = 0;
  if ( stpfitqp ) delete [] stpfitqp; stpfitqp = 0;
  if ( stpph0sigmap ) delete [] stpph0sigmap; stpph0sigmap = 0;
  if ( stpph0mip ) delete [] stpph0mip; stpph0mip = 0;
  if ( stpph0gev ) delete [] stpph0gev; stpph0gev = 0;
  if ( stpph1sigmap ) delete [] stpph1sigmap; stpph1sigmap = 0;
  if ( stpph1mip ) delete [] stpph1mip; stpph1mip = 0;
  if ( stpph1gev ) delete [] stpph1gev; stpph1gev = 0;
  if ( stpattn0c0 ) delete [] stpattn0c0; stpattn0c0 = 0;
  if ( stpattn1c0 ) delete [] stpattn1c0; stpattn1c0 = 0;
  if ( stpt0 ) delete [] stpt0; stpt0 = 0;
  if ( stpt1 ) delete [] stpt1; stpt1 = 0;
  if ( stptcal0t0 ) delete [] stptcal0t0; stptcal0t0 = 0;
  if ( stptcal1t0 ) delete [] stptcal1t0; stptcal1t0 = 0;
  
  return;
    
}

std::ostream& NtpSRTrack::Print(std::ostream& os) const {
  //
  //  Purpose:  Print track data on ostream.
  //

  os << "NtpSRTrack::Print" << endl;
  os << "index " << index << " slc " << slc << " ndigit " << ndigit
     << " nstpcnt " << nstpcnt << " nstrip " << nstrip  << " ds(m) " 
     << ds << " range(g/cm^2) " 
     << range << " cpu(sec) " << cputime << endl;
  os << plane << ph;
  os << "vtx " << vtx << "end " << end << "lin " << lin;
  os << "fidvtx " << fidvtx << "fidend " << fidend << "fidall " << fidall;
  os << time << momentum << fit << cr << endl;
  
  return os;

}

void NtpSRTrack::Print(Option_t* /* option */) const {
  //
  //  Purpose:  Print track in form supported by TObject::Print.
  //
  //  Arguments: option (not used)
  //

  Print(std::cout);
  return;

}

---