/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Event/RawDataCnv/RawDataCnv-00-04-35/src/EventManagement/EmcMcHitBuilder.cxx

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include "Identifier/EmcID.h"
#include "RawDataCnv/EventManagement/EmcMcHitBuilder.h"
using namespace std;

EmcMcHitBuilder::EmcMcHitBuilder()
   :Builder()
{
  initialize(Builder::m_confFile);
}

void EmcMcHitBuilder::unPack(vector<uint32_t>::const_iterator& iter,
               vector<uint32_t>::const_iterator& eiter, EmcTruth_t& et)
{
  uint32_t helpVal = *(iter++); assert(iter != eiter);
  et.x       = *(iter++); assert(iter != eiter);
  et.y       = *(iter++); assert(iter != eiter);
  et.z       = *(iter++); assert(iter != eiter);
  et.px         = *(iter++); assert(iter != eiter);
  et.py         = *(iter++); assert(iter != eiter);
  et.pz         = *(iter++); assert(iter != eiter);
  et.totalEdep  = *(iter++);

  et.trackIndex = (helpVal&m_trackIndexMask) >> m_trackIndexIndex;
  et.partId     = (helpVal&m_partIdMask) >> m_partIdIndex;
  et.numTheta   = (helpVal&m_numThetaMask) >> m_numThetaIndex;
  et.numPhi     = (helpVal&m_numPhiMask) >> m_numPhiIndex;

  return;
}


StatusCode EmcMcHitBuilder::pack(EmcMcHitCol* emcMcHitCol, WriteRawEvent*& re)
{
  /*
  if (emcMcHitCol == NULL) {
    cout << "EmcMcHitBuilder::pack cant get EmcMcHitCol" << endl;
    return StatusCode::FAILURE;
  }

  vector<uint32_t> *emcReMcHitVec = new vector<uint32_t>;
  EmcTruth_t m_EmcTruth;

  EmcMcHitCol::const_iterator pEmcMcHit = emcMcHitCol->begin();
  for ( ; pEmcMcHit != emcMcHitCol->end(); pEmcMcHit++) {
    // Make the EmcTruth Data
    makeEmcTruth(pEmcMcHit, m_EmcTruth);
    // pack
    uint32_t helpVal = (m_EmcTruth.trackIndex<<m_trackIndexIndex) & m_trackIndexMask;
    helpVal |= ((m_EmcTruth.partId<<m_partIdIndex) & m_partIdMask);
    helpVal |= ((m_EmcTruth.numTheta<<m_numThetaIndex) & m_numThetaMask);
    helpVal |= ((m_EmcTruth.numPhi<<m_numPhiIndex) & m_numPhiMask);
    // fill the McHit vector
    emcReMcHitVec->push_back(helpVal);
    emcReMcHitVec->push_back(m_EmcTruth.x);
    emcReMcHitVec->push_back(m_EmcTruth.y);
    emcReMcHitVec->push_back(m_EmcTruth.z);
    emcReMcHitVec->push_back(m_EmcTruth.px);
    emcReMcHitVec->push_back(m_EmcTruth.py);
    emcReMcHitVec->push_back(m_EmcTruth.pz);
    emcReMcHitVec->push_back(m_EmcTruth.totalEdep);
  }

  OfflineEventFormat::SubDetectorHeader sh(OfflineEventFormat::EMCTRUTH);
  SubRawEvent *sub = new SubRawEvent(sh, emcReMcHitVec);
  re->append(sub);
  */

  return StatusCode::SUCCESS;
}

// initialize re2te tables

StatusCode  EmcMcHitBuilder::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, "##EmcTruthConf", initFile)) {
    cerr << "Error: could not find '##EmcTruthConf' in file " << initFile << endl;
    return StatusCode::FAILURE;
  }

  if ( !Builder::expect( f, "#EmcTruthShift", initFile) ||
       !Builder::expectInt( f, "trackIndex", initFile, m_trackIndexIndex, m_trackIndexMask) ||
       !Builder::expectInt( f, "partId", initFile, m_partIdIndex, m_partIdMask) ||
       !Builder::expectInt( f, "numTheta", initFile, m_numThetaIndex, m_numThetaMask) ||
       !Builder::expectInt( f, "numPhi", initFile, m_numPhiIndex, m_numPhiMask) ||
       !Builder::expect( f, "#EmcTruthCoeff", 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) ||
       !Builder::expectLong( f, "totalEdep", initFile, m_totalEdepCoeff)  )
    return StatusCode::FAILURE;

  f.close();

  return StatusCode::SUCCESS;
}


uint32_t EmcMcHitBuilder::getTEID(uint32_t reid)
{
  return 0;
}

 
uint32_t EmcMcHitBuilder::getREID(uint32_t teid)
{
  return 0;
}

void EmcMcHitBuilder::makeEmcTruth
(EmcMcHitCol::const_iterator& pEmcMcHit, EmcTruth_t& et)
{
  Identifier ident = (*pEmcMcHit)->identify();

  et.trackIndex = (*pEmcMcHit)->getTrackIndex();
  et.partId     = EmcID::barrel_ec(ident);
  et.numTheta   = EmcID::theta_module(ident);
  et.numPhi     = EmcID::phi_module(ident);
  et.x          = int((*pEmcMcHit)->getPositionX() * m_xCoeff);
  et.y          = int((*pEmcMcHit)->getPositionY() * m_yCoeff);
  et.z          = int((*pEmcMcHit)->getPositionZ() * m_zCoeff);
  et.px         = int((*pEmcMcHit)->getPx() * m_pxCoeff);
  et.py         = int((*pEmcMcHit)->getPy() * m_pyCoeff);
  et.pz         = int((*pEmcMcHit)->getPz() * m_pzCoeff);
  et.totalEdep  = int((*pEmcMcHit)->getDepositEnergy() * m_totalEdepCoeff);

  return;
}

---