NtpSRShower.cxx

Go to the documentation of this file.

// NtpSRShower
//
// NtpSRShower is an ntuple class to hold the shower data
//

#include <iostream>
using namespace std;

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

ClassImp(NtpSRShower)

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

CVSID("$Id: NtpSRShower.cxx,v 1.15 2006/08/12 11:14:27 cbs Exp $");

//_____________________________________________________________________________
void NtpSRShower::AddStripAt(Int_t stripindex, Int_t ind) {
  //
  //  Purpose:  Add stripindex to fStripInd array at position ind.
  //
  //  Arguments: stripindex: index of strip in assoc fStrips TClonesArray. 
  //             ind: index in fStripInd array, cannot exceed bounds
  //                  0->fNStrip-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 NtpSRShower::AddClusterAt(Int_t clusterindex, Int_t ind) {
  //
  //  Purpose:  Add clusterindex to clu array at position ind.
  //
  //  Arguments: clusterindex: index of cluster in assoc cluster TClonesArray. 
  //             ind: index in cluster array, cannot exceed bounds
  //                  0->ncluster-1. 
  // 
  //  Return:  none.
  //
  //  Contact:   S. Kasahara
  //

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

  clu[ind] = clusterindex;
  return;

}

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

  if ( nstrip ) {
    stp           = new Int_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]; 
    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] = -1;
      stpu[i] = -999999;
      stpv[i] = -999999;
      stpx[i] = -999999;
      stpy[i] = -999999;
      stpz[i] = -999999;
      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;
    }
  }
  else {
    MSG("NtpSR",Msg::kWarning) << "NtpSRShower called with nstrip = 0" << endl;
  }

  if ( ncluster ) {
    clu = new Int_t[ncluster];
    for ( Int_t i = 0; i < ncluster; i++ ) clu[i] = -1;
  }
  else {
    MSG("NtpSR",Msg::kWarning) << "NtpSRShower called with ncluster = 0" 
                               << endl;
  }
}

//_____________________________________________________________________________
std::ostream& NtpSRShower::Print(std::ostream& os) const {
  //
  //  Purpose:  Print strip data on ostream.
  //

  os << "NtpSRShower::Print" << endl;
  os << "index " << index << " slc " << slc << " ndigit " << ndigit
     << " nstpcnt " << nstpcnt << " nstrip " << nstrip  
     << " ncluster(all,u,v)(" 
     << ncluster << "," << nUcluster << "," << nVcluster << ") "
     << "\npln" << plane << ph << vtx << sss;

  return os;

}

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

  Print(std::cout);
  return;

}

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

  nstrip = 0;
  if ( stp ) delete [] stp; stp = 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 ( 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;
  
}

  

---