NtpMREvent.cxx

Go to the documentation of this file.

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

#include <iostream>
using namespace std;

#include "MessageService/MsgService.h"
#include "MuonRemoval/NtpMREvent.h"
#include "MuonRemoval/AlgRmMu.h"

ClassImp(NtpMREvent)

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

CVSID("$Id: NtpMREvent.cxx,v 1.8 2008/09/17 02:56:19 tjyang Exp $");

void NtpMREvent::AddStripAt(Int_t stripindex, Int_t mrinfo, Int_t ind) {

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

Bool_t NtpMREvent::StpIsTrueMu(Int_t stripindex)
{
  if(stripindex>=nstrip) return false;
  return (mrstp[stripindex] & RmMuMask::kRMMU_ISMU_MASK);
}

Bool_t NtpMREvent::StpIsTrueShw(Int_t stripindex)
{
  if(stripindex>=nstrip) return false;
  return (mrstp[stripindex] & RmMuMask::kRMMU_ISSHW_MASK);
}

Bool_t NtpMREvent::StpIsScaled(Int_t stripindex)
{
  if(stripindex>=nstrip) return false;
  return (mrstp[stripindex] & RmMuMask::kRMMU_ISSCALED_MASK);
}

Bool_t NtpMREvent::StpIsInRecoTrk(Int_t stripindex)
{
  if(stripindex>=nstrip) return false;
  return (mrstp[stripindex] & RmMuMask::kRMMU_INRECOTRK_MASK);
}

Float_t  NtpMREvent::StpScaleFactor(Int_t stripindex)
{
  if(stripindex>=nstrip) return 0.0;
  return ( float(mrstp[stripindex]>>RmMuMask::kRMMU_NUM_SHIFT) / 
           RmMuMask::kRMMU_SCL_FACT_SCALE );
}

Bool_t NtpMREvent::StpIsRetained(Int_t stripindex)
{
  if(stripindex>=nstrip) return false;
  return (mrstp[stripindex] & RmMuMask::kRMMU_ISRETAIN_MASK);
}

Bool_t NtpMREvent::StpIsMCElec(Int_t stripindex)
{
  if(stripindex>=nstrip) return false;
  return (mrstp[stripindex] & RmMuMask::kRMMU_ISMCELEC_MASK);
}

NtpMREvent::NtpMREvent(): 
  index(0),ndigit(0),nstrip(0),stp(0),mrstp(0),
  orig_event(-1),best_event(-1),best_purity(0),best_complete(0),
  best_purity_phw(0),best_complete_phw(0),
  elec_complete(0),elec_complete_phw(0),
  best_elec_complete_phw(0), best_elec_complete(0), 
  best_elec_event(0), comb_complete(0), comb_complete_phw(0),
  best_comb_complete(0), best_comb_complete_phw(0),
  best_comb_event(0),
  shwvtxx(0),shwvtxy(0),shwvtxz(0),shwendx(0),shwendy(0),shwendz(0),
  shwvtxplane(0),shwendplane(0),shwnplane(0),shwcharge(0),
  vtxx(0),vtxy(0),vtxz(0),vtxdistance(),endx(0),endy(0),endz(0),enddistance(0),
  vtxp(0),endp(0),npln(0),prng(0),pcrv(0),pvdx(0),pvdy(0),pvdz(0),zenith(0),azimuth(0),
  fitp(0),endc(0),pass(0),pmux(0),pmuy(0),
  pmuz(0),mxpl(0),qp(0),mrmpmux(0),mrmpmuy(0),
  mrmpmuz(0),mrmQ2(0),mrmEshw(0)
{
}

NtpMREvent::NtpMREvent(Int_t ninit) : 
  index(0),ndigit(0),nstrip(ninit),stp(0),mrstp(0),
  orig_event(-1),best_event(-1),best_purity(0),best_complete(0),
  best_purity_phw(0),best_complete_phw(0),
  elec_complete(0),elec_complete_phw(0),
  best_elec_complete_phw(0), best_elec_complete(0),
  best_elec_event(0), comb_complete(0), comb_complete_phw(0),
  best_comb_complete(0), best_comb_complete_phw(0),
  best_comb_event(0),
  shwvtxx(0),shwvtxy(0),shwvtxz(0),shwendx(0),shwendy(0),shwendz(0),
  shwvtxplane(0),shwendplane(0),shwnplane(0),shwcharge(0),
  vtxx(0),vtxy(0),vtxz(0),vtxdistance(),endx(0),endy(0),endz(0),enddistance(0),
  vtxp(0),endp(0),npln(0),prng(0),pcrv(0),pvdx(0),pvdy(0),pvdz(0),zenith(0),azimuth(0),
  fitp(0),endc(0),pass(0),pmux(0),pmuy(0),
  pmuz(0),mxpl(0),qp(0),mrmpmux(0),mrmpmuy(0),
  mrmpmuz(0),mrmQ2(0),mrmEshw(0)
{
  if ( nstrip ) {
    stp = new Int_t[nstrip];
    mrstp = new Int_t[nstrip];
    for ( Int_t i = 0; i < nstrip; i++ ) {
      stp[i] = -1;
      mrstp[i] = -1;
    }
  }
}

std::ostream& NtpMREvent::Print(std::ostream& os) const {
  
  os << "NtpMREvent::Print" << endl;
  return os;
  
}

void NtpMREvent::Print(Option_t* /* option */) const {
  
  Print(std::cout);
  return;

}

---