NtpSREvent.cxx

Go to the documentation of this file.

// NtpSREvent
//
// NtpSREvent is an ntuple class to hold the event data
//

#include <iostream>
using namespace std;

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

ClassImp(NtpSREvent)

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

CVSID("$Id: NtpSREvent.cxx,v 1.13 2007/09/11 21:01:19 schubert Exp $");

NtpSREvent::NtpSREvent(Int_t nstripinit, Int_t nshowerinit, Int_t ntrackinit): 
  index(0),slc(0),ndigit(0),nstpcnt(nstripinit),nstrip(nstripinit),stp(0),
  stpph0sigmap(0),stpph0mip(0),stpph0gev(0),stpph1sigmap(0),stpph1mip(0),
  stpph1gev(0),nshower(nshowerinit),shw(0),ntrack(ntrackinit),trk(0),
  primshw(-1),primtrk(-1), contained(1) {

  // Normal constructor

  if ( nstrip ) {
    stp = new Int_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];
    
    // Initialize all indices to -1
    for (Int_t i = 0; i < nstrip; i++ ) { 
      stp[i] = -1;
      stpph0sigmap[i] = -1;
      stpph0mip[i]    = -1;
      stpph0gev[i]    = -1;
      stpph1sigmap[i] = -1;
      stpph1mip[i]    = -1;
      stpph1gev[i]    = -1;
    }
    
  }

  if ( nshower ) {
    shw = new Int_t[nshower];
    // Initialize all indices to -1
    for (Int_t i = 0; i < nshower; i++ ) shw[i] = -1;
  }

  if ( ntrack ) {
    trk = new Int_t[ntrack];
    // Initialize all indices to -1
    for (Int_t i = 0; i < ntrack; i++ ) trk[i] = -1;
  }



}

void NtpSREvent::AddTrackAt(Int_t trackindex, Int_t ind) {
  //
  //  Purpose:  Add trackindex to trk array at position ind.
  //
  //  Arguments: trackindex: index of track in assoc trk TClonesArray. 
  //             ind: index in trk array, cannot exceed bounds
  //                  0->ntrack-1.
  //
  //  Return:  none.
  //
  //  Contact:   S. Kasahara
  // 

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

  trk[ind] = trackindex;
  return;
    
}

void NtpSREvent::AddShowerAt(Int_t showerindex, Int_t ind) {
  //
  //  Purpose:  Add showerindex to shw array at position ind.
  //
  //  Arguments: showerindex: index of shower in assoc shw TClonesArray. 
  //             ind: index in shw array, cannot exceed bounds
  //                  0->nshower-1.
  //
  //  Return:  none.
  //
  //  Contact:   S. Kasahara
  // 

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

  shw[ind] = showerindex;
  return;
    
}

void NtpSREvent::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;
    
}


std::ostream& NtpSREvent::Print(std::ostream& os) const {
  //
  //  Purpose:  Print event data on ostream.
  //

  os << "NtpSREvent::Print" << endl;
  os << "index " << index << " slc " << slc << " ndigit " << ndigit 
     << " nstpcnt " << nstpcnt << " nstrip " << nstrip << " pln " << plane;
  os << "Event has " << nshower << " shower";
  if ( nshower > 1 ) {
    os << "s w/indices ";
    for ( int is = 0; is < nshower-1; is++ ) os << shw[is] << ",";
    os << shw[nshower-1] << " and ";
  }
  else if ( nshower == 1 ) {
    os << " w/index " << shw[nshower-1] << " and ";
  }
  else {
    os << "s and ";
  }
  
  os << ntrack << " track";
  if ( ntrack > 1 ) {
    os << "s w/indices ";
    for ( int it = 0; it < ntrack-1; it++ ) os << trk[it] << ",";
    os << trk[ntrack-1];
  }
  else if ( ntrack == 1 ) {
    os << " w/index " << trk[ntrack-1];
  }
  else {
    os << "s";
  }
  os << endl;
  os << ph << "vtx " << vtx << "end " << end;
  os << bleach;
  os << win;
  return os;

}

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

  Print(std::cout);
  return;

}

---