ANtpSnarlManipulator.cxx

Go to the documentation of this file.

//$Id: ANtpSnarlManipulator.cxx,v 1.2 2005/05/04 16:09:58 minoscvs Exp $
//
//ANtpSnarlManipulator.cxx
//
//B. Rebel 10/2004

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

#include <cassert>
#include <algorithm>

ClassImp(ANtpSnarlManipulator)

CVSID("$Id: ANtpSnarlManipulator.cxx,v 1.2 2005/05/04 16:09:58 minoscvs Exp $");

//-------------------------------------------------------------------
ANtpSnarlManipulator::ANtpSnarlManipulator() :
  fEventArray(0),
  fTrackArray(0),
  fShowerArray(0),
  fStripArray(0),
  fUseLengthTrack(0),
  fUsePHShower(0),
  fUseNPlanesShower(0),
  fUsePHTrack(0),
  fUseNPlanesTrack(0)
{
    
  MSG("ANtpSnarlManipulator", Msg::kDebug) << "ANtpSnarlManipulator::Constructor" << endl;
  
}

//-------------------------------------------------------------------
void ANtpSnarlManipulator::Initialize(TClonesArray  *eventArray, TClonesArray *trackArray, 
                                      TClonesArray *showerArray, TClonesArray *stripArray)
{
  fEventArray = eventArray;
  fTrackArray = trackArray;
  fShowerArray = showerArray;
  fStripArray = stripArray;

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


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

//------------------------------------------------------------------------
NtpSREvent *ANtpSnarlManipulator::GetEvent(Int_t index) 
{
  NtpSREvent *ntpEvent = 0;
      
  if(index < fEventArray->GetLast()+1)
    ntpEvent = dynamic_cast<NtpSREvent *>(fEventArray->At(index));
  
  return ntpEvent;
}

//------------------------------------------------------------------------
NtpSRTrack *ANtpSnarlManipulator::GetTrack(Int_t index) 
{
  NtpSRTrack *ntpTrack = 0;
      
  if(index < fTrackArray->GetLast()+1)
    ntpTrack = dynamic_cast<NtpSRTrack *>(fTrackArray->At(index));
  
  return ntpTrack;
}

//------------------------------------------------------------------------
NtpSRShower *ANtpSnarlManipulator::GetShower(Int_t index) 
{
  NtpSRShower *ntpShower = 0;

  if(index < fShowerArray->GetLast()+1)
    ntpShower = dynamic_cast<NtpSRShower *>(fShowerArray->At(index));
  
  return ntpShower;
}

//-------------------------------------------------------------------------
NtpSRTrack *ANtpSnarlManipulator::GetPrimaryTrack() 
{
  NtpSRTrack *primaryTrack = 0;
  NtpSRTrack *ntpTrack = 0;

  for(Int_t i = 0; i < fTrackArray->GetLast()+1; ++i){
    
    ntpTrack = dynamic_cast<NtpSRTrack *>(fTrackArray->At(i));
    
    //if you dont already have a candidate, make this track it
    if(!primaryTrack) primaryTrack = ntpTrack;
    
    //check to see if this track is better than the best so far.  the fUseXXX
    //variables are set to 1. or 0. depending on if they should be used to 
    //determine which is the primary track
    if(ntpTrack->ds > primaryTrack->ds*fUseLengthTrack
       && ntpTrack->plane.n > primaryTrack->plane.n*fUseNPlanesTrack
       && ntpTrack->ph.sigmap > primaryTrack->ph.sigmap*fUsePHTrack)
      primaryTrack = ntpTrack;
    
  }//end loop over tracks in the event
  
  return primaryTrack;
  
}

//-------------------------------------------------------------------------
NtpSRShower *ANtpSnarlManipulator::GetPrimaryShower() 
{
  NtpSRShower *ntpShower = 0;
  NtpSRShower *primaryShower = 0;
  
  for(Int_t i = 0; i < fShowerArray->GetLast()+1; ++i){
    
    ntpShower = dynamic_cast<NtpSRShower *>(fShowerArray->At(i));
    
    //if you dont already have a candidate, make this Shower it
    if(!primaryShower) primaryShower = ntpShower;
    
    //check to see if this Shower is better than the best so far.  the fUseXXX
    //variables are set to 1. or 0. depending on if they should be used to 
    //determine which is the primary Shower
    if(ntpShower->plane.n > primaryShower->plane.n*fUseNPlanesShower
       && ntpShower->ph.sigmap > primaryShower->ph.sigmap*fUsePHShower)
      primaryShower = ntpShower;
    
  }//end loop over Showers in the event
  
  return primaryShower;
  
}

//------------------------------------------------------------------------
NtpSRStrip *ANtpSnarlManipulator::GetStrip(Int_t index) 
{
  NtpSRStrip *ntpStrip = 0;

  if(index < fStripArray->GetLast()+1)
    ntpStrip = dynamic_cast<NtpSRStrip *>(fStripArray->At(index));
  
  return ntpStrip;
}

//--------------------------------------------------------------------------
void ANtpSnarlManipulator::SetPrimaryTrackCriteria(Int_t useNPlanes, Int_t useLength,
                                                   Int_t usePH)
{
  fUseLengthTrack = useLength;
  fUseNPlanesTrack = useNPlanes;
  fUsePHTrack = usePH;
  
  return;
}

//--------------------------------------------------------------------------
void ANtpSnarlManipulator::SetPrimaryShowerCriteria(Int_t useNPlanes, Int_t usePH)
{
  fUseNPlanesShower = useNPlanes;
  fUsePHShower = usePH;
  
  return;
}

---