In This Package:

SimHeaderCnvHits.cc

Go to the documentation of this file.

#include "SimHeaderCnv.h"
#include "PerSimEvent/PerSimHitHeader.h"
#include "Event/SimHitHeader.h"

#include "PerSimEvent/PerSimPmtHit.h"
#include "PerSimEvent/PerSimRpcHit.h"
#include "Event/SimPmtHit.h"
#include "Event/SimRpcHit.h"

#include <iostream>

using namespace DayaBay;
using namespace std;


DayaBay::SimHitCollection* SimHeaderCnv::convert(vector<SimTrack*>& timap,
                                                 const PerSimHitCollection& phc)
{
    const PerSimHitCollection::phit_container& phits = phc.hits;

    size_t nphits = phits.size();

    SimHitCollection* hc = new SimHitCollection;
    hc->setDetector(Detector(phc.detId));

    if (!nphits) return hc; // empty

    bool is_pmt = (0 != dynamic_cast<const PerSimPmtHit*>(phits[0]));

    SimHitCollection::hit_container hits;

    if (is_pmt) {
        for (size_t ihit=0; ihit < nphits; ++ihit) {
            const PerSimPmtHit* phit = dynamic_cast<const PerSimPmtHit*>(phits[ihit]);
            SimPmtHit* pmt_hit = new SimPmtHit;
            pmt_hit->setHc(hc);
            pmt_hit->setDir(phit->dir);
            pmt_hit->setPol(phit->pol);
            pmt_hit->setWavelength(phit->wavelength);
            pmt_hit->setType(phit->type);
            hits.push_back(pmt_hit);
        }
    }
    else {
        for (size_t ihit=0; ihit < nphits; ++ihit) {
            const PerSimRpcHit* phit = dynamic_cast<const PerSimRpcHit*>(phits[ihit]);
            SimRpcHit* rpc_hit = new SimRpcHit;
            rpc_hit->setHc(hc);
            rpc_hit->setEnergyDep(phit->energy);
            hits.push_back(rpc_hit);
        }
    }

    // Now copy the common part
    for (size_t ihit=0; ihit < nphits; ++ihit) {
        const PerSimHit *phit = phits[ihit];

        if (0 == (0x0000ff00&phit->sensDetId) && (0x000000ff&phit->sensDetId) > 6) {
            cerr << "SimHeaderCnvHits: warning got bogus PMTid: " << (void*)phit->sensDetId << endl;
        }

        SimHit *hit = hits[ihit];
        hit->setSensDetId(phit->sensDetId);
        hit->setHitTime(phit->time);
        hit->setLocalPos(phit->pos);
        hit->setWeight(phit->weight);

        int index = phit->ancestor.index;
        if (phit->ancestor.index < 0) {
            // emit error to log
            continue;
        }
        int count = phit->ancestor.count;
        hit->setAncestor(SimTrackReference(timap[index],count));
    }


    hc->setCollection(hits);
    return hc;
}

PerSimHitCollection* SimHeaderCnv::convert(map<const SimTrack*,int>& trackMap,
                                           const DayaBay::SimHitCollection& hc)
{
    MsgStream log(msgSvc(), "SimHeaderCnvHits::convert(t2p)");

    const SimHitCollection::hit_container &hits = hc.collection();

    PerSimHitCollection* phc = new PerSimHitCollection;
    phc->detId = hc.detector().siteDetPackedData();

    size_t nhits = hits.size();

    if (!nhits) return phc; // empty

    bool is_pmt = (0 != dynamic_cast<const SimPmtHit*>(hits[0]));

    PerSimHitCollection::phit_container& phits = phc->hits;

    // First create the output hit of specific type and fill specific data
    if (is_pmt) {
        for (size_t ihit=0; ihit < nhits; ++ihit) {
            const SimPmtHit *hit = dynamic_cast<const SimPmtHit*>(hits[ihit]);
            PerSimPmtHit* pmt_hit = new PerSimPmtHit(hit->dir(),
                                                     hit->pol(),
                                                     hit->wavelength(),
                                                     hit->type());
            phits.push_back(pmt_hit);
        }
    }
    else {
        for (size_t ihit=0; ihit < nhits; ++ihit) {
            const SimRpcHit *hit = dynamic_cast<const SimRpcHit*>(hits[ihit]);
            PerSimRpcHit* rpc_hit = new PerSimRpcHit(hit->energyDep());
            phits.push_back(rpc_hit);
        }
    }

    // Now copy the common part
    int nfailed = 0;
    for (size_t ihit=0; ihit < nhits; ++ihit) {
        const SimHit *hit = hits[ihit];
        PerSimHit *phit = phits[ihit];
        phit->sensDetId = hit->sensDetId();
        phit->time = hit->hitTime();
        phit->pos = hit->localPos();
        phit->weight = hit->weight();

        const SimTrack* ancestor = hit->ancestor().track();
        if (!ancestor) {
            ++nfailed;
            continue;
        }

        map<const SimTrack*,int>::const_iterator ait = trackMap.find(ancestor);
        if (ait == trackMap.end()) {
            log << MSG::ERROR
                << "No ancestor track index found for #" 
                << ancestor->trackId()
                << endreq;
            continue;
        }
        int count = hit->ancestor().indirection();
        phit->ancestor = PerSimIndirection(ait->second,count);
    }
    if (nfailed) {
        log << MSG::DEBUG
            << "Failed to find ancestor for hit " << nfailed << " times"  << endreq;
    }
        
    return phc;
}

---