ANtpMCManipulator.cxx

Go to the documentation of this file.

//$Id: ANtpMCManipulator.cxx,v 1.3 2005/11/17 22:15:19 brebel Exp $
//
//ANtpMCManipulator.cxx
//
//B. Rebel 10/2004

#include "AnalysisNtuples/Module/ANtpMCManipulator.h"
#include "MessageService/MsgService.h"

#include <cassert>
#include <algorithm>

ClassImp(ANtpMCManipulator)

CVSID("$Id: ANtpMCManipulator.cxx,v 1.3 2005/11/17 22:15:19 brebel Exp $");

//-------------------------------------------------------------------
ANtpMCManipulator::ANtpMCManipulator() :
  fMCArray(0),
  fStdArray(0),
  fTHEventArray(0),
  fTHTrackArray(0),
  fTHShowerArray(0),
  fTHStripArray(0)
{
    
  MSG("ANtpMCManipulator", Msg::kDebug) << "ANtpMCManipulator::Constructor" << endl;
  
}

//-------------------------------------------------------------------
void ANtpMCManipulator::Initialize(TClonesArray *mcArray, TClonesArray *stdArray,
                                   TClonesArray  *thEventArray, TClonesArray *thTrackArray, 
                                   TClonesArray *thShowerArray, TClonesArray *thStripArray)
{
  fMCArray = mcArray;
  fStdArray = stdArray;
  fTHEventArray = thEventArray;
  fTHTrackArray = thTrackArray;
  fTHShowerArray = thShowerArray;
  fTHStripArray = thStripArray;

  MSG("ANtpMCManipulator", Msg::kDebug) << "ANtpMCManipulator::Constructor" << endl;    
}


//----------------------------------------------------------------------
ANtpMCManipulator::~ANtpMCManipulator()
{ 
  MSG("ANtpMCManipulator", Msg::kDebug) << "ANtpMCManipulator::Destructor" << endl; 
}

//------------------------------------------------------------------------
NtpMCTruth *ANtpMCManipulator::GetNtpMCTruth(Int_t index) 
{
  NtpMCTruth *ntpMCTruth = 0;

  if(!fMCArray) return ntpMCTruth;

  if(index < 0 && fMCArray->GetEntries()>0) index = 0;
  if(index < fMCArray->GetEntries()+1)
    ntpMCTruth = dynamic_cast<NtpMCTruth *>(fMCArray->At(index));
  
  return ntpMCTruth;
}

//------------------------------------------------------------------------
NtpMCStdHep *ANtpMCManipulator::GetNtpMCStdHep(Int_t index) 
{
  NtpMCStdHep *ntpMCStdHep = 0;

  if(!fStdArray) return ntpMCStdHep;

  if(index < 0 && fStdArray->GetEntries()>0) index = 0;
  if(index < fStdArray->GetEntries()+1)
    ntpMCStdHep = dynamic_cast<NtpMCStdHep *>(fStdArray->At(index));
  
  return ntpMCStdHep;
}

//------------------------------------------------------------------------
NtpTHEvent *ANtpMCManipulator::GetNtpTHEvent(Int_t index) 
{
  NtpTHEvent *ntpTHEvent = 0;

  if(!fTHEventArray) return ntpTHEvent;

  if(index < 0 && fTHEventArray->GetEntries()>0) index = 0;
  if(index < fTHEventArray->GetLast()+1)
    ntpTHEvent = dynamic_cast<NtpTHEvent *>(fTHEventArray->At(index));
  
  return ntpTHEvent;
}

//------------------------------------------------------------------------
NtpTHTrack *ANtpMCManipulator::GetNtpTHTrack(Int_t index) 
{
  NtpTHTrack *ntpTHTrack = 0;

  if(!fTHTrackArray) return ntpTHTrack;

  if(index < 0 && fTHTrackArray->GetEntries()>0) index = 0;
  if(index < fTHTrackArray->GetLast()+1)
    ntpTHTrack = dynamic_cast<NtpTHTrack *>(fTHTrackArray->At(index));
  
  return ntpTHTrack;
}

//------------------------------------------------------------------------
NtpTHShower *ANtpMCManipulator::GetNtpTHShower(Int_t index) 
{
  NtpTHShower *ntpTHShower = 0;

  if(!fTHShowerArray) return ntpTHShower;

  if(index < 0 && fTHShowerArray->GetEntries()>0) index = 0;
  if(index < fTHShowerArray->GetLast()+1)
    ntpTHShower = dynamic_cast<NtpTHShower *>(fTHShowerArray->At(index));
  
  return ntpTHShower;
}

//------------------------------------------------------------------------
NtpTHStrip *ANtpMCManipulator::GetNtpTHStrip(Int_t index) 
{
  NtpTHStrip *ntpTHStrip = 0;

  if(!fTHStripArray) return ntpTHStrip;

  if(index < 0 && fTHStripArray->GetEntries()>0) index = 0;
  if(index < fTHStripArray->GetLast()+1)
    ntpTHStrip = dynamic_cast<NtpTHStrip *>(fTHStripArray->At(index));
  
  return ntpTHStrip;
}

---