EventInfoByteStreamAuxCnv Class Reference

ByteStream converter for xAOD::EventAuxInfo. More...

#include <EventInfoByteStreamAuxCnv.h>

Inheritance diagram for EventInfoByteStreamAuxCnv:

Collaboration diagram for EventInfoByteStreamAuxCnv:

Public Member Functions
	EventInfoByteStreamAuxCnv (ISvcLocator *svcloc)
virtual	~EventInfoByteStreamAuxCnv () override=default
virtual StatusCode	initialize () override
virtual StatusCode	finalize () override
virtual StatusCode	createObj (IOpaqueAddress *pAddr, DataObject *&pObj) override
	converter method to create object
virtual StatusCode	createRep (DataObject *pObj, IOpaqueAddress *&pAddr) override
	converter method to write object
virtual long	repSvcType () const override
	Storage type and class ID.
void	addTrack (TIDA::Track *t)
const std::vector< TIDA::Track * > &	tracks () const
void	clear ()
void	selectTracks (const Rec::TrackParticleContainer *trigtracks)
void	selectTracks (const Analysis::MuonContainer *muontracks)
bool	msgLvl (const MSG::Level lvl) const
	Test the output level.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.
void	setLevel (MSG::Level lvl)
	Change the current logging level.

Static Public Member Functions
static long	storageType ()
static CLID	classID ()

Protected Member Functions
double	phiCorr (double phi)
void	ipCorr (double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)

Protected Attributes
double	m_beamX
double	m_beamY
double	m_beamZ
std::vector< TIDA::Track * >	m_tracks

Private Member Functions
std::string	ascTime (unsigned int t)
	convert timestamp to ascii time.
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
ServiceHandle< IROBDataProviderSvc >	m_robDataProvider
	RODDataProviderSvc handle.
ServiceHandle< StoreGateSvc >	m_mdSvc
	TDS handle.
ToolHandle< IMCEventInfoByteStreamTool >	m_mcEventInfoTool {"MCEventInfoByteStreamTool"}
	Tool for decoding MC EventInfo from ROB.
bool	m_isSimulation
bool	m_isTestbeam
bool	m_isCalibration
std::string	m_nm
	Message source name.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels).
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging).

Detailed Description

ByteStream converter for xAOD::EventAuxInfo.

Reading: EventInfo is built from the RawEvent header. For MC ByteStream files, MC-specific fields are decoded from a dedicated ROB fragment via MCEventInfoByteStreamTool.

Writing: For MC files, MC-specific EventInfo fields are encoded to a dedicated ROB fragment via MCEventInfoByteStreamTool.

Definition at line 35 of file EventInfoByteStreamAuxCnv.h.

Constructor & Destructor Documentation

EventInfoByteStreamAuxCnv()

EventInfoByteStreamAuxCnv::EventInfoByteStreamAuxCnv

(

ISvcLocator *

svcloc

)

Definition at line 31 of file EventInfoByteStreamAuxCnv.cxx.

  : Converter(storageType(), classID(), svcloc)
  , AthMessaging(svcloc != nullptr ? msgSvc() : nullptr, "EventInfoByteStreamAuxCnv")
  , m_robDataProvider("ROBDataProviderSvc", "EventInfoByteStreamAuxCnv")
  , m_mdSvc("InputMetaDataStore", "EventInfoByteStreamAuxCnv")
  , m_isSimulation(false)
  , m_isTestbeam(false)
  , m_isCalibration(false)
{
}

~EventInfoByteStreamAuxCnv()

virtual EventInfoByteStreamAuxCnv::~EventInfoByteStreamAuxCnv ( )

overridevirtualdefault

Member Function Documentation

addTrack()

void Converter::addTrack ( TIDA::Track * t )

inlineinherited

Definition at line 43 of file Converter.h.

                              {
    m_tracks.push_back(t);
  }

ascTime()

std::string EventInfoByteStreamAuxCnv::ascTime ( unsigned int t )

private

convert timestamp to ascii time.

Definition at line 336 of file EventInfoByteStreamAuxCnv.cxx.

{
   struct tm t;
   t.tm_sec   = tstamp;
   t.tm_min   = 0;
   t.tm_hour  = 0;
   t.tm_mday  = 0;
   t.tm_mon   = 0;
   t.tm_year  = 70;
   t.tm_wday  = 00;
   t.tm_yday  = 00;
   t.tm_isdst = 0;
   time_t ct = mktime(&t);
   struct tm t2;
   gmtime_r(&ct, &t2);
   char buf[32];
   asctime_r (&t2, buf);
   return std::string (buf);
}

classID()

CLID EventInfoByteStreamAuxCnv::classID ( )

static

Definition at line 42 of file EventInfoByteStreamAuxCnv.cxx.

{
  return ClassID_traits<xAOD::EventAuxInfo>::ID();
}

clear()

void Converter::clear ( )

inlineinherited

Definition at line 51 of file Converter.h.

{m_tracks.clear();} 

createObj()

StatusCode EventInfoByteStreamAuxCnv::createObj ( IOpaqueAddress * pAddr, DataObject *& pObj )

overridevirtual

converter method to create object

Definition at line 123 of file EventInfoByteStreamAuxCnv.cxx.

{
  ByteStreamAddress* pRE_Addr{};
  ATH_CHECK( (pRE_Addr = dynamic_cast<ByteStreamAddress*>(pAddr)) != nullptr );
 
  const EventContext& ctx = pRE_Addr->getEventContext();
  ATH_MSG_DEBUG("Creating Objects");
 
  // get RawEvent
  const RawEvent* re = m_robDataProvider->getEvent(ctx);
  if (!re) {
    ATH_MSG_ERROR("Can not get RawEvent ");
    return StatusCode::FAILURE;
  }
 
  // Run Number
  int runNumber = re->run_no();
 
  // Event Number
  uint64_t eventNumber;
  if (re->version() < 0x03010000) {
    eventNumber=re->lvl1_id();
  } else {
    eventNumber=re->global_id();
  }
  
  // Time Stamp
  uint32_t bc_time_sec = re->bc_time_seconds();
  uint32_t bc_time_ns  = re->bc_time_nanoseconds();
  // bc_time_ns should be lt 1e9.
  if (bc_time_ns > 1000000000) {
    if (runNumber < 20920) {
      //  In M4 we saw otherwise, because sec/ns was swapped in raw data.
      uint32_t temp = bc_time_ns;
      bc_time_ns = bc_time_sec;
      bc_time_sec= temp;
      ATH_MSG_DEBUG("bc_time  second/nanosecond swapped, sec/ns = " << bc_time_sec << " " << bc_time_ns);
    } 
    else { 
      // For later runs, the nanosecond clock sometimes is not reset, making it overrun 1e9. Round it off to 1e9
      ATH_MSG_WARNING("bc_time nanosecond number larger than 1e9, it is " << bc_time_ns << ", reset it to 1 sec");
      bc_time_ns = 1000000000;
    }
  }
 
  // luminosity block number
  uint16_t lumiBlock = re->lumi_block();
 
  // bunch crossing identifier
  uint16_t bcID = re->bc_id();
 
  unsigned int detMask0 = 0xFFFFFFFF, detMask1 = 0xFFFFFFFF, detMask2 = 0xFFFFFFFF, detMask3 = 0xFFFFFFFF;
  // Get ByteStream Metadata from Input MetaData Store
  const ByteStreamMetadataContainer* metadatacont{nullptr};
  StatusCode status = m_mdSvc->retrieve(metadatacont, "ByteStreamMetadata");
  if (!status.isSuccess()) {
    ATH_MSG_WARNING("Unable to retrieve Input MetaData for ByteStream");
  } 
  else {
    const ByteStreamMetadata* metadata = *(metadatacont->begin());
    uint64_t detectorMask = metadata->getDetectorMask();
    detMask0 = (unsigned int)(detectorMask & 0x00000000FFFFFFFF);
    detMask1 = (unsigned int)(detectorMask >> 32);
    uint64_t detectorMask2 = metadata->getDetectorMask2();
    detMask2 = (unsigned int)(detectorMask2 & 0x00000000FFFFFFFF);
    detMask3 = (unsigned int)(detectorMask2 >> 32);
  }
 
 
  xAOD::EventInfo evtInfo;
  xAOD::EventAuxInfo* pEvtInfoAux = new xAOD::EventAuxInfo();
  evtInfo.setStore(pEvtInfoAux);
 
  evtInfo.setRunNumber(runNumber);
  evtInfo.setEventNumber(eventNumber);
  evtInfo.setLumiBlock(lumiBlock);
  evtInfo.setTimeStamp(bc_time_sec);
  evtInfo.setTimeStampNSOffset(bc_time_ns);
  evtInfo.setBCID(bcID);
  evtInfo.setDetectorMask(detMask0,detMask1);
  evtInfo.setDetectorMaskExt(detMask2,detMask3);
 
  // Set default MC values (will be overwritten if MC ROB is present)
  evtInfo.setMCChannelNumber(0);
  evtInfo.setMCEventNumber(0);
  evtInfo.setMCEventWeights(std::vector<float>(1,1));
 
  // Set Event Type
  uint32_t eventTypeBitmask{0};
  if (m_isSimulation) {
    eventTypeBitmask |= xAOD::EventInfo::IS_SIMULATION;
  }
  if (m_isTestbeam) {
    eventTypeBitmask |= xAOD::EventInfo::IS_TESTBEAM;
  }
  if (m_isCalibration) {
    eventTypeBitmask |= xAOD::EventInfo::IS_CALIBRATION;
  }
  evtInfo.setEventTypeBitmask(eventTypeBitmask);
 
  // Try to decode MC EventInfo from dedicated ROB fragment (for MC ByteStream files)
  bool mcEventInfoDecoded = false;
  if (m_isSimulation && m_mcEventInfoTool.isEnabled()) {
    // Request the MC EventInfo ROB
    std::vector<uint32_t> mcRobIds = m_mcEventInfoTool->robIds();
    std::vector<const OFFLINE_FRAGMENTS_NAMESPACE::ROBFragment*> mcRobFragments;
    m_robDataProvider->getROBData(ctx, mcRobIds, mcRobFragments, "EventInfoByteStreamAuxCnv");
 
    if (!mcRobFragments.empty()) {
      ATH_MSG_DEBUG("Found MC EventInfo ROB fragment, decoding...");
      StatusCode sc = m_mcEventInfoTool->convertFromBS(mcRobFragments[0], evtInfo);
      if (sc.isFailure()) {
        ATH_MSG_WARNING("Failed to decode MC EventInfo from ROB fragment, using default values");
      } else {
        ATH_MSG_DEBUG("Successfully decoded MC EventInfo from ROB fragment");
        mcEventInfoDecoded = true;
      }
    } else {
      ATH_MSG_DEBUG("No MC EventInfo ROB fragment found, using default values");
    }
  }
 
  // Trigger Info
  const OFFLINE_FRAGMENTS_NAMESPACE::DataType* buffer;
  // status element
  re->status(buffer);
  uint32_t statusElement = *buffer;
 
  // extended LVL1ID
  uint32_t extendedLevel1ID = re->lvl1_id();
 
  // LVL1 trigger type
  uint32_t level1TriggerType = re->lvl1_trigger_type();
 
  // stream tag
  std::vector<xAOD::EventInfo::StreamTag> streamTags;
  std::vector<eformat::helper::StreamTag> onl_streamTags;
  re->stream_tag(buffer);
  eformat::helper::decode(re->nstream_tag(), buffer, onl_streamTags);
  for (const eformat::helper::StreamTag& onl_streamTag : onl_streamTags) {
    std::set<uint32_t> tmp_off_dets = std::set<uint32_t>();
    if (!onl_streamTag.dets.empty()) {
      std::set<eformat::SubDetector> tmp_onl_dets = onl_streamTag.dets;
      for (const eformat::SubDetector& subdet : tmp_onl_dets) {
        tmp_off_dets.insert((uint32_t) subdet);
      }
    }
    streamTags.push_back(xAOD::EventInfo::StreamTag(onl_streamTag.name
                            , onl_streamTag.type
                            , onl_streamTag.obeys_lumiblock
                            , onl_streamTag.robs
                            , tmp_off_dets)
             );
  }
 
  evtInfo.setStatusElement(statusElement);
  // Only set extendedLevel1ID from RawEvent header if MC EventInfo wasn't decoded
  // (MC EventInfo ROB contains the original extendedLevel1ID value)
  if (!mcEventInfoDecoded) {
    evtInfo.setExtendedLevel1ID(extendedLevel1ID);
  }
  evtInfo.setLevel1TriggerType(level1TriggerType);
  evtInfo.setStreamTags(streamTags);
 
  // record EventInfo
  evtInfo.setEventFlags(xAOD::EventInfo::Core, m_robDataProvider->getEventStatus(ctx));
  pObj = SG::asStorable(pEvtInfoAux);
 
  ATH_MSG_DEBUG(" New xAOD::EventAuxInfo made, run/event= " << runNumber
        << " " << eventNumber
        << " Time stamp  = " << ascTime(bc_time_sec) 
        );
 
  return StatusCode::SUCCESS;
}

createRep()

StatusCode EventInfoByteStreamAuxCnv::createRep ( DataObject * pObj, IOpaqueAddress *& pAddr )

overridevirtual

converter method to write object

Definition at line 299 of file EventInfoByteStreamAuxCnv.cxx.

{
  ATH_MSG_DEBUG("createRep for xAOD::EventAuxInfo");
 
  // For MC ByteStream files, create the MC EventInfo ROB fragment
  if (m_isSimulation && m_mcEventInfoTool.isEnabled()) {
    ATH_MSG_DEBUG("Creating MC EventInfo ROB fragment for simulation");
 
    // Get RawEventWrite via ByteStreamCnvSvc
    SmartIF<IByteStreamEventAccess> byteStreamCnvSvc(service("ByteStreamCnvSvc"));
    if (!byteStreamCnvSvc.isValid()) {
      ATH_MSG_ERROR("Failed to retrieve ByteStreamCnvSvc");
      return StatusCode::FAILURE;
    }
    RawEventWrite* re = byteStreamCnvSvc->getRawEvent();
    if (!re) {
      ATH_MSG_ERROR("Failed to get RawEventWrite");
      return StatusCode::FAILURE;
    }
 
    // Create MC EventInfo ROB fragment
    std::vector<OFFLINE_FRAGMENTS_NAMESPACE_WRITE::ROBFragment*> mcRobs;
    ATH_CHECK(m_mcEventInfoTool->convertToBS(mcRobs, Gaudi::Hive::currentContext()));
 
    // Append ROB fragments to the event
    for (auto* rob : mcRobs) {
      // Set LVL1 ID and trigger type from the full event
      rob->rod_lvl1_id(re->lvl1_id());
      rob->rod_lvl1_type(re->lvl1_trigger_type());
      re->append(rob);
      ATH_MSG_DEBUG("Added MC EventInfo ROB fragment to output event");
    }
  }
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode EventInfoByteStreamAuxCnv::finalize ( )

overridevirtual

Definition at line 112 of file EventInfoByteStreamAuxCnv.cxx.

{
  ATH_MSG_DEBUG("Finalize");
  
  StatusCode sc = Converter::finalize();
  if (sc.isFailure()) {
    ATH_MSG_WARNING("Converter::finalize() failed");
  }
  return sc;
}

initialize()

StatusCode EventInfoByteStreamAuxCnv::initialize ( )

overridevirtual

Definition at line 52 of file EventInfoByteStreamAuxCnv.cxx.

{
  ATH_MSG_DEBUG("Initialize");
 
  CHECK(Converter::initialize());
 
  CHECK(m_robDataProvider.retrieve());
  CHECK(m_mdSvc.retrieve());
 
  SmartIF<IProperty> byteStreamCnvSvc(service("ByteStreamCnvSvc"));
  CHECK( byteStreamCnvSvc.isValid() );
 
  SimpleProperty<bool> propIsSimulation("IsSimulation", m_isSimulation);
  StatusCode sc = byteStreamCnvSvc->getProperty(&propIsSimulation);
  if (sc.isSuccess()) {
    m_isSimulation = propIsSimulation.value();
    ATH_MSG_INFO("IsSimulation : " << m_isSimulation);
  } 
  else {
    ATH_MSG_ERROR("Cannot get IsSimulation");
    return sc;
  }
 
  SimpleProperty<bool> propIsTestbeam("IsTestbeam", m_isTestbeam);
  sc = byteStreamCnvSvc->getProperty(&propIsTestbeam);
  if (sc.isSuccess()) {
    m_isTestbeam = propIsTestbeam.value();
    ATH_MSG_INFO("IsTestbeam : " << m_isTestbeam);
  } 
  else {
    ATH_MSG_ERROR("Cannot get IsTestbeam");
    return sc;
  }
 
  SimpleProperty<bool> propIsCalibration("IsCalibration", m_isCalibration);
  sc = byteStreamCnvSvc->getProperty(&propIsCalibration);
  if (sc.isSuccess()) {
    m_isCalibration = propIsCalibration.value();
    ATH_MSG_INFO("IsCalibration : " << m_isCalibration);
  }
  else {
    ATH_MSG_ERROR("Cannot get IsCalibration");
    return sc;
  }
 
  // Retrieve MC EventInfo decoder tool (optional - only needed for MC BS files)
  if (m_isSimulation) {
    if (m_mcEventInfoTool.retrieve().isFailure()) {
      ATH_MSG_WARNING("Failed to retrieve MCEventInfoByteStreamTool - MC EventInfo will use default values");
      m_mcEventInfoTool.disable();
    } else {
      ATH_MSG_DEBUG("MCEventInfoByteStreamTool retrieved successfully");
    }
  } else {
    m_mcEventInfoTool.disable();
  }
 
  return StatusCode::SUCCESS;
}

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

ipCorr()

void Converter::ipCorr ( double d0, double z0, double & d0c, double & z0c, double phi0, double eta, double pt )

inlineprotectedinherited

Definition at line 200 of file Converter.h.

                                                                                                  {
    
    double sn  = sin(phi0);
    double cs  = cos(phi0);
    double sd0 = (d0 != 0 ? d0/fabs(d0) : 1);
    double spt = (pt != 0 ? pt/fabs(pt) : 1);
    
    if (fabs(pt) >= 1*CLHEP::TeV) {
 
      d0c = d0 + m_beamX*sn - m_beamY*cs;
 
    } else {
      
      double rc = fabs(pt)*15.0/(9.0*1.042);
      
      double xc = (fabs(d0)-spt*sd0*rc)*cos(phi0+M_PI/2*sd0) - m_beamX;
      double yc = (fabs(d0)-spt*sd0*rc)*sin(phi0+M_PI/2*sd0) - m_beamY;
      
      double newphi;
      double xd01,yd01,xd02,yd02;
      
      if (xc == 0) {
    xd01 = 0; yd01 = rc+yc;
    xd02 = 0; yd02 = -rc+yc; 
      } else {
    xd01 = xc+yc/xc*yc+sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd01 = yc/xc*xd01;
    xd02 = xc+yc/xc*yc-sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd02 = yc/xc*xd02;
      }
      
      double r1 = sqrt(xd01*xd01+yd01*yd01);
      double r2 = sqrt(xd02*xd02+yd02*yd02);
      
      double phiV;
      
      if (r1 < r2)
    phiV = atan2(yd01,xd01); 
      else
    phiV = atan2(yd02,xd02); 
      
      double phi1 = phiCorr(phiV+M_PI/2);
      double phi2 = phiCorr(phiV-M_PI/2);
      
      if (fabs(phiCorr(phi1-phi0))<=fabs(phiCorr(phi2-phi0)))
    newphi = phi1;
      else
    newphi = phi2;
      
      d0c = fabs(sqrt(xc*xc+yc*yc)-rc)*sin(phiV-newphi);
      
      double theta=2*atan2(exp(-eta),1);
      double theta_save=theta;
      theta = theta - (1+spt)/2*M_PI;
      if (theta>0) theta = theta_save;
      
      double deltaz0= -spt*rc/tan(theta)*phiCorr(phiCorr(newphi)-phiCorr(phi0));
      z0c = z0 + deltaz0;
      
    }
  }

msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 167 of file AthMessaging.h.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 182 of file AthMessaging.h.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  // If user did not set explicit message level we have to initialize
  // the messaging and retrieve the default via the MessageSvc.
  if (m_lvl==MSG::NIL && !m_initialized.test_and_set()) initMessaging();
 
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

phiCorr()

double Converter::phiCorr ( double phi )

inlineprotectedinherited

Definition at line 191 of file Converter.h.

                            {
    
    if (phi < -TMath::Pi()) phi += 2*TMath::Pi();
    if (phi >  TMath::Pi()) phi -= 2*TMath::Pi();
    
    return phi;
  }

repSvcType()

virtual long EventInfoByteStreamAuxCnv::repSvcType ( ) const

inlineoverridevirtual

Storage type and class ID.

Definition at line 50 of file EventInfoByteStreamAuxCnv.h.

{ return i_repSvcType(); }

selectTracks() [1/2]

void Converter::selectTracks ( const Analysis::MuonContainer * muontracks )

inlineinherited

Definition at line 120 of file Converter.h.

                                                               {
    
    static int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
    
    Analysis::MuonContainer::const_iterator muonitr = muontracks->begin();
    Analysis::MuonContainer::const_iterator muonend = muontracks->end();
    
    while(muonitr!=muonend) { 
      
      bool hasInDet = (*muonitr)->hasInDetTrackParticle();
      const Rec::TrackParticle *trackitr=NULL; 
      if(hasInDet) trackitr=(*muonitr)->inDetTrackParticle();
      muonitr++;
      if(!hasInDet) continue;
      
      const Trk::MeasuredPerigee* measPer = trackitr->measuredPerigee();
      
      CLHEP::HepVector perigeeParams = measPer->parameters();    
      double pT    = measPer->pT(); 
      double eta   = measPer->eta();
      double phi   = perigeeParams[Trk::phi0];
      double z0    = perigeeParams[Trk::z0];
      double d0    = perigeeParams[Trk::d0];
      //double theta = perigeeParams[Trk::theta];
      
      double deta = 1;
      double dphi = 1;
      double dz0  = 1;
      double dd0  = 1;
      double dpT  = 1;    
      
      // Check number of hits
      // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
      // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
      // stereo clusters making a spacepoint are two "hits"
      const Trk::TrackSummary *summary = trackitr->trackSummary();
      int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits); 
      int nPixelHits  = 2*summary->get(Trk::numberOfPixelHits);  
      int nSctHits    = summary->get(Trk::numberOfSCTHits); 
      int nStrawHits  = summary->get(Trk::numberOfTRTHits);
      int nTrHits     = summary->get(Trk::numberOfTRTHighThresholdHits);
      
      int nSiHits     = nPixelHits + nSctHits;
      
      const Trk::FitQuality *quality   = trackitr->fitQuality();
      double chi2 = quality->chiSquared();
      double dof  = quality->quality->numberDoF();
 
      unsigned bitmap = 0;
      
      for ( int ih=0 ; ih<20 ; ih++ ) {
    if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );    
      }
      
      // Shift d0 and z0 according to beam position
      ipCorr(d0, z0, d0, z0, phi, eta, pT);
      
      // Create and save Track  
      TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
                                 deta, dphi, dz0, dd0, dpT,
                                 nBlayerHits, nPixelHits, nSctHits, nSiHits,
                                 nStrawHits, nTrHits, bitmap, 0,
                                 -1) ;
      
      addTrack(t);
    }
    
  }

selectTracks() [2/2]

void Converter::selectTracks ( const Rec::TrackParticleContainer * trigtracks )

inlineinherited

Definition at line 54 of file Converter.h.

                                                                   { 
    
    Rec::TrackParticleContainer::const_iterator  trackitr = trigtracks->begin();
    Rec::TrackParticleContainer::const_iterator  trackend = trigtracks->end();
    
    static int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
 
    while ( trackitr!=trackend ) { 
        
      const Trk::MeasuredPerigee* measPer = (*trackitr)->measuredPerigee();
 
      CLHEP::HepVector perigeeParams = measPer->parameters();    
      double pT    = measPer->pT(); 
      double eta   = measPer->eta();
      double phi   = perigeeParams[Trk::phi0];
      double z0    = perigeeParams[Trk::z0];
      double d0    = perigeeParams[Trk::d0];
      //double theta = perigeeParams[Trk::theta];
    
      double deta = 1;
      double dphi = 1;
      double dz0  = 1;
      double dd0  = 1;
      double dpT  = 1;    
 
      // Check number of hits
      // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
      // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
      // stereo clusters making a spacepoint are two "hits"
      const Trk::TrackSummary *summary = (*trackitr)->trackSummary();
      int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits); 
      int nPixelHits  = 2*summary->get(Trk::numberOfPixelHits);  
      int nSctHits    = summary->get(Trk::numberOfSCTHits); 
      int nStrawHits  = summary->get(Trk::numberOfTRTHits);
      int nTrHits     = summary->get(Trk::numberOfTRTHighThresholdHits);
      int nSiHits     = nPixelHits + nSctHits;
      bool expectBL   = false; // Not stored for Rec::TrackParticle
 
      const Trk::FitQuality *quality   = (*trackitr)->fitQuality();
      double chi2 = quality->chiSquared();
      double dof  = quality->numberDoF();
 
      unsigned bitmap = 0;
      
      for ( int ih=0 ; ih<20 ; ih++ ) {
    if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );    
      }
 
      // Shift d0 and z0 according to beam position
      ipCorr(d0, z0, d0, z0, phi, eta, pT);
 
      // Create and save Track  
      TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof, 
                                 deta, dphi, dz0, dd0, dpT,
                                 nBlayerHits, nPixelHits, nSctHits, nSiHits,
                                 nStrawHits, nTrHits, bitmap, 0,
                                 -1,
                                 expectBL) ;
  
      addTrack(t);
      trackitr++;
    }
    
  }

setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

storageType()

long EventInfoByteStreamAuxCnv::storageType ( )

static

Definition at line 47 of file EventInfoByteStreamAuxCnv.cxx.

{
  return ByteStreamAddress::storageType();
}

tracks()

const std::vector< TIDA::Track * > & Converter::tracks ( ) const

inlineinherited

Definition at line 48 of file Converter.h.

{return m_tracks;}

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging).

Definition at line 141 of file AthMessaging.h.

m_beamX

double Converter::m_beamX

protectedinherited

Definition at line 263 of file Converter.h.

m_beamY

double Converter::m_beamY

protectedinherited

Definition at line 263 of file Converter.h.

m_beamZ

double Converter::m_beamZ

protectedinherited

Definition at line 263 of file Converter.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_isCalibration

bool EventInfoByteStreamAuxCnv::m_isCalibration

private

Definition at line 63 of file EventInfoByteStreamAuxCnv.h.

m_isSimulation

bool EventInfoByteStreamAuxCnv::m_isSimulation

private

Definition at line 61 of file EventInfoByteStreamAuxCnv.h.

m_isTestbeam

bool EventInfoByteStreamAuxCnv::m_isTestbeam

private

Definition at line 62 of file EventInfoByteStreamAuxCnv.h.

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

m_mcEventInfoTool

ToolHandle<IMCEventInfoByteStreamTool> EventInfoByteStreamAuxCnv::m_mcEventInfoTool {"MCEventInfoByteStreamTool"}

private

Tool for decoding MC EventInfo from ROB.

Definition at line 58 of file EventInfoByteStreamAuxCnv.h.

{"MCEventInfoByteStreamTool"}; 

m_mdSvc

ServiceHandle<StoreGateSvc> EventInfoByteStreamAuxCnv::m_mdSvc

private

TDS handle.

Definition at line 57 of file EventInfoByteStreamAuxCnv.h.

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels).

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_robDataProvider

ServiceHandle<IROBDataProviderSvc> EventInfoByteStreamAuxCnv::m_robDataProvider

private

RODDataProviderSvc handle.

Definition at line 56 of file EventInfoByteStreamAuxCnv.h.

m_tracks

std::vector<TIDA::Track*> Converter::m_tracks

protectedinherited

Definition at line 264 of file Converter.h.

---

The documentation for this class was generated from the following files:

• EventInfoByteStreamAuxCnv.h
• EventInfoByteStreamAuxCnv.cxx