BmnHit.cxx

Go to the documentation of this file.

 
#include "BmnHit.h"
 
#include "FairRootManager.h"
 
#include <TBranch.h>
#include <TChain.h>
#include <TClonesArray.h>
#include <map>
#include <string>
 
BmnHit::BmnHit()
    : FairHit()
{
    fUsing = kFALSE;
    fFlag = kTRUE;
    fIndex = -1;
    fDetId = kREF;
    fStation = -1;
    fModule = -1;
    fType = -1;
    fLength = -1.0;
    fResX = -1000.;
    fResY = -1000.;
    fCovXY = -1000.;
    fPseudo = false;
}
 
BmnHit::BmnHit(Int_t detID, TVector3 pos, TVector3 dpos, Int_t index)
    : FairHit(detID, pos, dpos, index)
{
    fUsing = kFALSE;
    fFlag = kTRUE;
    fIndex = -1;
    fDetId = kREF;
    fStation = -1;
    fModule = -1;
    fType = -1;
    fLength = -1.0;
    fResX = -1000.;
    fResY = -1000.;
    fCovXY = -1000.;
    fPseudo = false;
}
 
BmnHit::~BmnHit() {}
 
// Strip cluster of the hit in the given strip layer, taken from the cluster branch of this layer
StripCluster* BmnHit::GetCluster(const char* cluster_branch_name, StripLayerType layer)
{
    // a pseudo hit is reconstructed from a lower layer cluster only
    if (layer == UpperStripLayer && fPseudo)
        return nullptr;
    Int_t cluster_index = (layer == LowerStripLayer) ? fLowerClusterIndex : fUpperClusterIndex;
 
    // each branch is looked up only once for the current FairRootManager, because every GetObject call
    // for a branch which is not in memory takes a new slot of the internal array of FairRootManager
    static thread_local FairRootManager* cached_manager = nullptr;
    static thread_local map<string, TClonesArray*> cluster_arrays;
 
    FairRootManager* ioman = FairRootManager::Instance();
    if (ioman != cached_manager) {
        cluster_arrays.clear();
        cached_manager = ioman;
    }
 
    auto it = cluster_arrays.find(cluster_branch_name);
    if (it == cluster_arrays.end()) {
        TClonesArray* clusters = dynamic_cast<TClonesArray*>(ioman->GetObject(cluster_branch_name));
        // a branch of the input file activated only now does not contain the current event yet, so it is read here
        TChain* input_chain = ioman->GetInChain();
        if (clusters && input_chain && input_chain->GetTree()) {
            TBranch* branch = input_chain->GetTree()->GetBranch(cluster_branch_name);
            if (branch)
                branch->GetEntry(input_chain->GetTree()->GetReadEntry());
        }
        it = cluster_arrays.emplace(cluster_branch_name, clusters).first;
    }
 
    TClonesArray* clusters = it->second;
    if (!clusters || cluster_index < 0 || cluster_index >= clusters->GetEntriesFast())
        return nullptr;
    return static_cast<StripCluster*>(clusters->UncheckedAt(cluster_index));
}