RawDataEmcMcHitCnv.hpp

Go to the documentation of this file.

#pragma once
 
#include "RawDataBaseCnv.h"
 
#include "Identifier/EmcID.h"
#include "McTruth/EmcMcHit.h"
#include "RawDataCnv/ClassID_temp.h"
#include "RawDataCnv/EventManagement/EmcMcHitBuilder.h"
#include "RawDataCnv/EventManagement/McTruthDataStructure.h"
 
extern const CLID& CLID_EmcMcHitCol;
 
class RawDataEmcMcHitCnv : public RawDataBaseCnv {
public:
  RawDataEmcMcHitCnv( ISvcLocator* svc ) : RawDataBaseCnv( svc, classID() ) {}
 
  static const CLID& classID() { return CLID_EmcMcHitCol; }
 
  StatusCode createObj( IOpaqueAddress* pAddr, DataObject*& pObj ) {
 
    // This converter will create an empty EmcMcHitCol on the TDS
    EmcMcHitCol* emcMcHitCol = new EmcMcHitCol;
    pObj                     = emcMcHitCol;
 
    RAWEVENT* evt = m_inputSvc->currentEvent();
    if ( evt == NULL )
    {
      // MsgStream log(msgSvc(), "RawDataEmcMcHitCnv");
      // log << MSG::ERROR << "RawDataCnv has no event!" << endmsg;
      return StatusCode::FAILURE;
    }
 
    EmcMcHit*  emcMcHit;
    EmcTruth_t m_EmcTruth;
 
    assert( ( evt->getEmcTruth().size() % 8 ) == 0 );
 
    std::vector<uint32_t>::const_iterator iter  = evt->getEmcTruth().begin();
    std::vector<uint32_t>::const_iterator eiter = evt->getEmcTruth().end();
 
    for ( int emcMcHitId = 0; iter != eiter; emcMcHitId++ )
    {
      // retrieve the McTruth data
      m_emcMcHitBuilder.unPack( iter, eiter, m_EmcTruth );
      // construct the identifier
      Identifier ident =
          EmcID::crystal_id( m_EmcTruth.partId, m_EmcTruth.numTheta, m_EmcTruth.numPhi );
      // construct the EmcMcHit
      emcMcHit = new EmcMcHit(
          ident, m_EmcTruth.trackIndex, m_EmcTruth.x / ( m_emcMcHitBuilder.m_xCoeff * 1.0 ),
          m_EmcTruth.y / ( m_emcMcHitBuilder.m_yCoeff * 1.0 ),
          m_EmcTruth.z / ( m_emcMcHitBuilder.m_zCoeff * 1.0 ),
          m_EmcTruth.px / ( m_emcMcHitBuilder.m_pxCoeff * 1.0 ),
          m_EmcTruth.py / ( m_emcMcHitBuilder.m_pyCoeff * 1.0 ),
          m_EmcTruth.pz / ( m_emcMcHitBuilder.m_pzCoeff * 1.0 ),
          m_EmcTruth.totalEdep / ( m_emcMcHitBuilder.m_totalEdepCoeff * 1.0 ) );
      // And add the stuff to the container
      emcMcHitCol->push_back( emcMcHit );
    }
 
    return StatusCode::SUCCESS;
  }
 
  StatusCode createRep( DataObject* pObj, IOpaqueAddress*& pAddr ) {
    // convert PixelRaw in the container into ByteStream
    // MsgStream log(msgSvc(), "RawDataEmcMcHitCnv");
 
    WriteRawEvent*& re = m_cnvSvc->getRawEvent();
 
    if ( re == 0 )
    {
      // log << " get RawEvent failed !" << endmsg;
      return StatusCode::FAILURE;
    }
 
    SmartDataPtr<EmcMcHitCol> mcHitCol( dataProvider(), EventModel::MC::EmcMcHitCol );
 
    if ( mcHitCol == 0 )
    {
      // log << "no EmcMcHitCol found" << endmsg;
      return StatusCode::FAILURE;
    }
 
    StatusCode sc = m_emcMcHitBuilder.pack( mcHitCol, re );
 
    return sc;
  }
 
private:
  /// EmcMcHitBuilder for packing or unpacking data
  EmcMcHitBuilder m_emcMcHitBuilder;
};