TrgTrackToTrack.cpp

Go to the documentation of this file.

#include "TrgTrackToTrack.h"
#include "TrgConverter/TrgTrackStreamer.h"
#include "Event/TrackStreamer.h"

#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "CLHEP/Units/SystemOfUnits.h"

//-----------------------------------------------------------------------------
// Implementation file for class : TrgTrackToTrack
//
// 13/01/2005 : Jose A. Hernando
//-----------------------------------------------------------------------------

// Declaration of the Algorithm Factory
static const  AlgFactory<TrgTrackToTrack>           s_trgTrackToTrackFact ;
const        IAlgFactory& TrgTrackToTrackFactory  = s_trgTrackToTrackFact ;


TrgTrackToTrack::TrgTrackToTrack( const std::string& name,
                                  ISvcLocator* pSvcLocator)
  : GaudiAlgorithm ( name , pSvcLocator )
  , m_trgContainer()
{
  declareProperty( "TrgTracksName" ,
                   m_trgTracksName  = TrgTrackLocation::Long );

  declareProperty( "TracksLocation" ,
                   m_tracksLocation = "/Event/Rec/Tr/NewTracks" );

}

TrgTrackToTrack::~TrgTrackToTrack() {};

StatusCode TrgTrackToTrack::initialize() {

  // initialize the gaudi-algorithm
  StatusCode sc = GaudiAlgorithm::initialize();
  if ( sc.isFailure() ) return sc;
  
  // get the dataProvider
  TrgDataProvider* trgData = tool<TrgDataProvider>( "TrgDataProvider" );
  
  m_trgContainer = trgData->tracks(m_trgTracksName);

  release(trgData);

  info() << " -- parameters of TrgTrackToTrack -- " << endreq;
  info() << " TrgTrack input  container " << m_trgTracksName << endreq;
  info() << "    Track output contaner "  << m_tracksLocation << endreq;
  info() << " ----------------------------------- " << endreq;

  return StatusCode::SUCCESS;
}

StatusCode TrgTrackToTrack::execute() {

  StatusCode sc = StatusCode::SUCCESS;

  Tracks* tracks = new Tracks();

  for (TrgTrackIterator it = m_trgContainer->begin();
       it != m_trgContainer->end(); it++) {
    const TrgTrack& gtrack = *(*it);
    Track* track = new Track();
    TrgTrackConverter::toTrack(gtrack,*track);
    tracks->insert(track);
    testLHCbIDs(gtrack);
    debug() << str::trgTrack(gtrack) << endreq;
    debug() << str::track(*track) << endreq;
  }

  debug() << " number of TrgTracks " << m_trgContainer->size() << endreq;
  debug() << " number converted Tracks " << tracks->size() << endreq;

  // locate the tracks in the TES
  sc = put( tracks, m_tracksLocation);

  return sc;
  
}

StatusCode TrgTrackToTrack::finalize() {
  return StatusCode::SUCCESS;
}


void TrgTrackToTrack::testLHCbIDs(const TrgTrack& gtrack) {

  std::vector<LHCbID> ids = TrgTrackConverter::lhcbIDs(gtrack);
  int i = 0;

  const std::vector< std::pair<VeloChannelID,int> >& vref 
    = gtrack.veloReference();
  for (std::vector< std::pair<VeloChannelID,int> >::const_iterator
         it = vref.begin(); it != vref.end(); it++) {
    const std::pair<VeloChannelID,int>& velo = *it;
    LHCbID id = ids[i];
    debug() << " detector: velo " << velo.first << " size " << velo.second
            << endreq;
    debug() << " channelID " << velo.first.channelID() << endreq;
    debug() << " LHCbID " << id.detectorType() << ", "
            << id.channelID() <<"," << id.spareBits()<< endreq;
    // id.setSpareBits((unsigned) velo.second);
    i++;
  }
  
  const std::vector< std::pair<ITChannelID,int> >& itref 
    = gtrack.itReference();
  for (std::vector< std::pair<ITChannelID,int> >::const_iterator
         it = itref.begin(); it != itref.end(); it++) {
    const std::pair<ITChannelID,int>& si = *it;
    LHCbID id = ids[i];
    debug() << " detector: it " << si.first << " size " << si.second
            << endreq;
    debug() << " channelID " << si.first.channelID() << endreq;
    debug() << " LHCbID " << id.detectorType() << ", "
            << id.channelID() <<"," << id.spareBits()<< endreq;
    i++;
  }
  
  
  const std::vector< std::pair<OTChannelID,int> >& otref 
    = gtrack.otReference();
  for (std::vector< std::pair<OTChannelID,int> >::const_iterator
         it = otref.begin(); it != otref.end(); it++) {
    const std::pair<OTChannelID,int>& ot = *it;
    LHCbID id = ids[i];
    debug() << " detector: ot " << ot.first << " size " << ot.second
            << endreq;
    debug() << " channelID " << ot.first.channelID()  << endreq;
    debug() << " LHCbID " << id.detectorType() << ", "
            << id.channelID() <<"," << id.spareBits()<< endreq;
    // id.setSpareBits((unsigned) ot.second);
    i++;
  }
}

---