TrkTrackFakeReader.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2019 CERN for the benefit of the ATLAS collaboration
*/

 
// Silicon trackers includes
#include "TrkTrackFakeReader.h"
 
//InDetRawData includes
#include "InDetRawData/PixelRDO_Container.h"
#include "InDetRawData/SCT_RDO_Container.h"
#include "InDetRawData/TRT_RDO_Container.h"
#include "InDetRawData/PixelRDO_Collection.h"
#include "InDetRawData/SCT_RDO_Collection.h"
#include "InDetRawData/TRT_RDO_Collection.h"
 
//  #include "InDetRawData/InDetRawDataContainer.h"
//  #include "InDetRawData/InDetRawDataCollection.h"
//  #include "InDetRawData/Pixel1RawData.h"
//  #include "InDetRawData/InDetRawDataCLASS_DEF.h"
//  #include "InDetRawData/PixelRDORawData.h"
//  #include "InDetRawUtils/PixelHid2RESrcID.h"
 
// Id includes
#include "InDetIdentifier/PixelID.h"
#include "InDetIdentifier/SCT_ID.h"
#include "InDetIdentifier/TRT_ID.h"
//#include "InDetDetDescr/InDet_DetDescrManager.h"
 
// test includes
#include "TrkTrack/Track.h"
#include "TrkTrack/TrackCollection.h"
#include "InDetPrepRawData/SiClusterCollection.h"
#include "InDetPrepRawData/SiClusterContainer.h"
#include "InDetPrepRawData/SiCluster.h"
 
#include <map>
#include "TrkTrackFakeWriter.h"
 
#include "StoreGate/ReadCondHandle.h"
 
// Constructor with parameters:
TrkTrackFakeReader::TrkTrackFakeReader(const std::string &name, 
ISvcLocator *pSvcLocator) :
  AthAlgorithm(name,pSvcLocator)
  {}
 
// Initialize method:
StatusCode TrkTrackFakeReader::initialize()
{
  ATH_MSG_INFO( "TrkTrackFakeReader::initialize()"  );
  ATH_CHECK( m_pixelDetEleCollKey.initialize() );
  return StatusCode::SUCCESS;
}
 
// Execute method:
StatusCode TrkTrackFakeReader::execute(const EventContext& ctx) 
{
  ATH_MSG_DEBUG( "TrkTrackFakeReader::execute()"  );
 
  using namespace Trk;
  using namespace InDet;
 
    //Get tracks
  const TrackCollection* tracks = nullptr;
  ATH_CHECK( evtStore()->retrieve(tracks,"Tracks") );
 
  SG::ReadCondHandle<InDetDD::SiDetectorElementCollection> pixelDetEleHandle(m_pixelDetEleCollKey, ctx);
  const InDetDD::SiDetectorElementCollection* elements = *pixelDetEleHandle;
  if (not pixelDetEleHandle.isValid() or elements==nullptr) {
    ATH_MSG_WARNING(m_pixelDetEleCollKey.fullKey() << " is not available.");
  }
 
  for (TrackCollection::const_iterator it = tracks->begin(); it!=tracks->end(); ++it)
  {
      //compare tracks
    Trk::Track* track = FakeTrackBuilder::buildTrack(elements);
      // compareTracks(track,*it);
 
    delete track;
  }
  return StatusCode::SUCCESS;
}
 
 
// Finalize method:
StatusCode TrkTrackFakeReader::finalize() 
{
  ATH_MSG_INFO( "TrkTrackFakeReader::finalize()"  );
  return StatusCode::SUCCESS;
}
 
void TrkTrackFakeReader::compareTracks(const Trk::Track* ref, const Trk::Track* readTrk) {
      // if (**it!=*track) {
      //   log<<MSG::ERROR<<"Tracks differ!"<<endmsg;
      //   
      //   log<<MSG::INFO<<"Reference track: "<<(*track)<<endmsg;
      //   log<<MSG::INFO<<"Track from ESD: "<<(**it)<<endmsg;
      // }
  // Get the messaging service, print where you are
  MsgStream log(msgSvc(), name());
  if ( !(ref->fitQuality () && readTrk->fitQuality() 
    && ref->fitQuality()->chiSquared()==readTrk->fitQuality ()->chiSquared()
  && ref->fitQuality()->numberDoF()==readTrk->fitQuality ()->numberDoF())) {
    ATH_MSG_WARNING( "Track fitqualities differ"  );
  } 
}