MucMcHitBuilder.cxx

Go to the documentation of this file.

#include "RawDataCnv/EventManagement/MucMcHitBuilder.h"
#include "Identifier/MucID.h"
#include <fstream>
#include <iostream>
using namespace std;
 
MucMcHitBuilder::MucMcHitBuilder() : Builder() { initialize( Builder::m_confFile ); }
 
void MucMcHitBuilder::unPack( vector<uint32_t>::const_iterator& iter,
                              vector<uint32_t>::const_iterator& eiter, MucTruth_t& mt ) {
  uint32_t helpVal = *( iter++ );
  assert( iter != eiter );
  mt.x = *( iter++ );
  assert( iter != eiter );
  mt.y = *( iter++ );
  assert( iter != eiter );
  mt.z = *( iter++ );
  assert( iter != eiter );
  mt.px = *( iter++ );
  assert( iter != eiter );
  mt.py = *( iter++ );
  assert( iter != eiter );
  mt.pz = *( iter++ );
 
  mt.trackIndex = ( helpVal & m_trackIndexMask ) >> m_trackIndexIndex;
  mt.partId     = ( helpVal & m_partIdMask ) >> m_partIdIndex;
  mt.segId      = ( helpVal & m_segIdMask ) >> m_segIdIndex;
  mt.gapId      = ( helpVal & m_gapIdMask ) >> m_gapIdIndex;
  mt.stripId    = ( helpVal & m_stripIdMask ) >> m_stripIdIndex;
 
  return;
}
 
StatusCode MucMcHitBuilder::pack( MucMcHitCol* mucMcHitCol, WriteRawEvent*& re ) {
  /*
 
    if (mucMcHitCol == NULL) {
      cout << "MucMcHitBuilder::pack cant get MucMcHitCol" << endl;
      return StatusCode::FAILURE;
    }
 
    vector<uint32_t> *mucReMcHitVec = new vector<uint32_t>;
    MucTruth_t m_MucTruth;
 
    MucMcHitCol::const_iterator pMucMcHit = mucMcHitCol->begin();
    for ( ; pMucMcHit != mucMcHitCol->end(); pMucMcHit++) {
      //make the MucTruth data
      makeMucTruth(pMucMcHit, m_MucTruth);
      //pack integers in mucTruth
      uint32_t helpVal = (m_MucTruth.trackIndex<<m_trackIndexIndex) & m_trackIndexMask;
      helpVal |= ((m_MucTruth.partId<<m_partIdIndex) & m_partIdMask);
      helpVal |= ((m_MucTruth.segId<<m_segIdIndex) & m_segIdMask);
      helpVal |= ((m_MucTruth.gapId<<m_gapIdIndex) & m_gapIdMask);
      helpVal |= ((m_MucTruth.stripId<<m_stripIdIndex) & m_stripIdMask);
      //fill the McHit vector
      mucReMcHitVec->push_back(helpVal);
      mucReMcHitVec->push_back(m_MucTruth.x);
      mucReMcHitVec->push_back(m_MucTruth.y);
      mucReMcHitVec->push_back(m_MucTruth.z);
      mucReMcHitVec->push_back(m_MucTruth.px);
      mucReMcHitVec->push_back(m_MucTruth.py);
      mucReMcHitVec->push_back(m_MucTruth.pz);
    }
 
    OfflineEventFormat::SubDetectorHeader sh(OfflineEventFormat::MUCTRUTH);
    SubRawEvent *sub = new SubRawEvent(sh, mucReMcHitVec);
    re->append(sub);
    */
 
  return StatusCode::SUCCESS;
}
 
// initialize re2te tables
 
StatusCode MucMcHitBuilder::initialize( string& initFile ) {
  ifstream f;
 
  // read init file
  f.open( initFile.c_str() );
 
  if ( f.bad() )
  {
    cerr << "Error: could not open file " << initFile << endl;
    return StatusCode::FAILURE;
  }
 
  if ( !Builder::find( f, "##MucTruthConf", initFile ) )
  {
    cerr << "Error: could not find '##MucTruthConf' in file " << initFile << endl;
    return StatusCode::FAILURE;
  }
 
  if ( !Builder::expect( f, "#MucTruthShift", initFile ) ||
       !Builder::expectInt( f, "trackIndex", initFile, m_trackIndexIndex, m_trackIndexMask ) ||
       !Builder::expectInt( f, "partId", initFile, m_partIdIndex, m_partIdMask ) ||
       !Builder::expectInt( f, "segId", initFile, m_segIdIndex, m_segIdMask ) ||
       !Builder::expectInt( f, "gapId", initFile, m_gapIdIndex, m_gapIdMask ) ||
       !Builder::expectInt( f, "stripId", initFile, m_stripIdIndex, m_stripIdMask ) ||
       !Builder::expect( f, "#MucTruthCoeff", initFile ) ||
       !Builder::expectLong( f, "x", initFile, m_xCoeff ) ||
       !Builder::expectLong( f, "y", initFile, m_yCoeff ) ||
       !Builder::expectLong( f, "z", initFile, m_zCoeff ) ||
       !Builder::expectLong( f, "px", initFile, m_pxCoeff ) ||
       !Builder::expectLong( f, "py", initFile, m_pyCoeff ) ||
       !Builder::expectLong( f, "pz", initFile, m_pzCoeff ) )
    return StatusCode::FAILURE;
 
  f.close();
 
  return StatusCode::SUCCESS;
}
 
uint32_t MucMcHitBuilder::getTEID( uint32_t reid ) { return 0; }
 
uint32_t MucMcHitBuilder::getREID( uint32_t teid ) { return 0; }
 
void MucMcHitBuilder::makeMucTruth( MucMcHitCol::const_iterator& pMucMcHit, MucTruth_t& mt ) {
  Identifier ident = ( *pMucMcHit )->identify();
 
  mt.trackIndex = ( *pMucMcHit )->getTrackIndex() % 1000;
  mt.partId     = MucID::part( ident );
  mt.segId      = MucID::seg( ident );
  mt.gapId      = MucID::gap( ident );
  mt.stripId    = MucID::strip( ident );
  mt.x          = int( ( *pMucMcHit )->getPositionX() * m_xCoeff );
  mt.y          = int( ( *pMucMcHit )->getPositionY() * m_yCoeff );
  mt.z          = int( ( *pMucMcHit )->getPositionZ() * m_zCoeff );
  mt.px         = int( ( *pMucMcHit )->getPx() * m_pxCoeff );
  mt.py         = int( ( *pMucMcHit )->getPy() * m_pyCoeff );
  mt.pz         = int( ( *pMucMcHit )->getPz() * m_pzCoeff );
 
  return;
}