MCint.cxx

Go to the documentation of this file.

#include "MCint.h"
#include "Mint.h"
#include "Jint.h"

#include <JobControl/JobC.h>
#include <MinosObjectMap/MomNavigator.h>
#include <Plex/PlexStripEndId.h>
#include <Plex/PlexHandle.h>
#include <UgliGeometry/UgliGeomHandle.h>
#include <CandDigit/CandDigitListHandle.h>
#include <CandDigit/CandDigitHandle.h>
#include <Record/SimSnarlRecord.h>
#include <RerootExodus/RerootExodus.h>
#include <REROOT_Classes/REROOT_FLSHit.h>
#include <REROOT_Classes/REROOT_FLSDigit.h>
#include <REROOT_Classes/REROOT_StdHep.h>

#include <Digitization/DigiScintHit.h>

#include <DataUtil/CDL2STL.h>
using namespace DataUtil;

#include <TClonesArray.h>
#include <TParticle.h>

#include <map>
#include <vector>
using namespace std;

MCint::MCint(Mint& mint)
    : fMint(mint)
    , fDirty(false)
{
}

MCint::~MCint()
{
}

void MCint::Update()
{
    fDirty = true;
}
void MCint::Clear()
{
    fScintHits.clear();
}

template <class T>
static void clones_array_to_vector(const TClonesArray& ca, vector<T*>& v)
{
    unsigned int siz = ca.GetLast()+1;
    v.clear();
    v.reserve(siz);
    for (unsigned int ind=0; ind<siz; ++ind) {
        T* t = dynamic_cast<T*>(ca.UncheckedAt(ind));
        if (!t) continue;
        v.push_back(t);
    }
    if (siz != v.size()) {
        cerr << "Warning: not all elements of TClonesArray were of the right type\n";
    }
}
void MCint::GetData()
{
    if (!fDirty) return;
    this->Clear();
    
    const MomNavigator& mom = fMint.GetJobC().Mom;
    SimSnarlRecord* ssr = dynamic_cast<SimSnarlRecord*>(mom.GetFragment("SimSnarlRecord"));
    if (!ssr) {
        cerr << "No SimSnarlRecord in this data\n";
        return;
    }

    const TClonesArray* arr = dynamic_cast<const TClonesArray*>
        (ssr->FindComponent("TClonesArray","DigiScintHits"));
    if (!arr) {
        cerr << "No DigiScintHit TClonesArray in this data\n";
        return;
    }

    clones_array_to_vector(*arr,fScintHits);
    fDirty = false;
}

const vector<DigiScintHit*>& MCint::GetScintHits() 
{
    this->GetData();
    return fScintHits;
}

---