---

TriggerRecon.cxx

Go to the documentation of this file.

//  $Header: /nfs/slac/g/glast/ground/cvs/reconstruction/src/TriggerRecon.cxx,v 1.13 2000/12/11 16:38:45 burnett Exp $
//

#include "reconstruction/TriggerRecon.h"

#include "data/CsIData.h"
#include "data/GlastData.h"

#include "DetectorData.h"

#include "reconstruction/ReconVisitor.h"

#include <algorithm>
using namespace std;


map<idents::ModuleId, int> TriggerRecon::TrigTowers;
const int TriggerRecon::MaxNumTowers = 60;


TriggerRecon::TriggerRecon ()
: Recon()
{
 setTitle("TriggerRecon: ver 2.1");     // hsg: changed version to 2.1 with level 2 addition
#if 0 // THB: remove 

 i_Triage_Time          = data.addElement("Triage_Time");
 i_bits                 = data.addElement("Trig_Bits");
 i_Energy_max_Twr       = data.addElement("Trig_max_Twr_E");
 i_Energy_max_Xtl       = data.addElement("Trig_max_Xtl_E");
 i_Count                = data.addElement("Trig_1x1");
 i_CV_Count             = data.addElement("Trig_cv_Vetos");
 i_Gd_Trig              = data.addElement("Trig_1x1_Gd_Trig");
 i_Level2_Trig          = data.addElement("Trig_Level2_Bits");
 i_Dead_Time            = data.addElement("Dead_Time");
 i_Max_FE_Buff          = data.addElement("Max_FE_Buff");
 i_Total_Evt_Size       = data.addElement("Total_Evt_Size");
 i_Max_Evt_Size         = data.addElement("Max_Evt_Size");
 i_Max_TEM_Thrpt        = data.addElement("Max_TEM_Thrpt");
 i_Max_RO_Time          = data.addElement("Max_RO_Time");
 i_N_Hits               = data.addElement("N_Tracker_Hits");

 i_L2_Layers            = data.addElement("L2_Layers");
 i_L2_Top               = data.addElement("L2_Top");
 i_L2_Mid               = data.addElement("L2_Mid");
 i_L2_Bot               = data.addElement("L2_Bot");
 i_L2_Dsize             = data.addElement("L2_Dsize");
 i_L3_flag              = data.addElement("L3_flag");
 i_L3_vetoFlag          = data.addElement("L3_vetoFlag");
 i_L3_nTracks           = data.addElement("L3_nTracks");
#endif
}

void TriggerRecon::accept (ReconVisitor& rv)
{
   rv.visit(this);
}

void TriggerRecon::reconstruct (const GlastData* glastdata)
{
    if(state() == DONE) return;
#if 0 //THB: put may problems to solve here

    // Get energy info from the CsI data
    const CsIData* csi = glastdata->getCsIData();

    map<ModuleId, float, less<ModuleId> > tower_sum;

    float   max_e_tower=0,
            max_e_xtal=0;

    for(int layer=0; layer < Calorimeter::numberOfLayers(); layer++)
    {
        int nh = csi->nHits(layer);
        for(int i=0; i<nh; i++) {
            float e = csi->energy(layer, i);
            max_e_xtal = max(max_e_xtal, e);
            ModuleId id = csi->moduleId(layer,i);
            tower_sum[id] += e;
                max_e_tower = std::max(max_e_tower, tower_sum[id]);
        }
    }

    int noCV=0,  no1x1=0, no1x1Trig=0, bitsL2=0, bitsL1=0, maxThrpt = 0,
        maxFEBuff = 0, totalEvtSize = 0;
    unsigned    nhits = 0;
    double       maxREvt = 0, maxROTime = 0;

    GlastEvent* evt = glastdata->glastEvent();
    if (!evt) { setState(DONE); return; }
    // get readout data statistics, if available. otherwise we are done here.

    const TEMDataMap&       tem_data = evt->tem_data();
    // get access to stored readout data...

    // statistics
    for( TEMDataMap::const_iterator it = tem_data.begin(); it != tem_data.end(); ++it) {
        if ((*it).second) {
            const TEMEvent& temevt = *(*it).second;
            unsigned        bit = temevt.level1bits();

            noCV +=      ((bit & Level1::b_ACDL) !=0);
            no1x1 +=     ((bit & Level1::b_Track) != 0);
            no1x1Trig += ((bit & Level1::b_Gamma) != 0);
            bitsL1 |= temevt.level1bits();
            const Level2& l2 = *temevt.level2();
            bitsL2 |= l2.flag();
            maxREvt = std::max(maxREvt, temevt.readevt_time() + temevt.sendevt_time());
            maxThrpt = std::max<int>(maxThrpt, temevt.throughput()*16);
            maxFEBuff = std::max<int>(maxFEBuff, temevt.fe_depth());
            nhits = std::max( nhits, temevt.data().nhits());
            
            totalEvtSize += (temevt.dram()) ? temevt.dram()->size()*16 : 0;
            maxROTime = std::max( maxROTime, temevt.readevt_time() + temevt.sendevt_time());
            
        
            if((bit & Level1::b_Track) != 0) {  // This tower triggered  (3 in a row in tracker)
                TrigTowers[(*it).first] = 1;
            }
        }
    }


    data.addData(i_Triage_Time   , evt->key());
    data.addData(i_bits          , bitsL1 & 127); // mask off system bits for now!
    data.addData(i_Energy_max_Twr, max_e_tower);
    data.addData(i_Energy_max_Xtl, max_e_xtal);
    data.addData(i_Count         , no1x1);
    data.addData(i_CV_Count      , noCV);
    data.addData(i_Gd_Trig       , no1x1Trig);
    data.addData(i_Level2_Trig   , bitsL2);
    data.addData(i_Dead_Time     , evt->dead_time());
    data.addData(i_Max_TEM_Thrpt , maxThrpt);
    data.addData(i_Max_FE_Buff   , maxFEBuff);
    data.addData(i_Max_RO_Time   , maxROTime);
    data.addData(i_Total_Evt_Size, totalEvtSize);
    data.addData(i_N_Hits        , nhits);
#if 0
    //Perform Level 3 trigger
    L3T* l3 = L3T::instance();
    l3->wakeup(evt);

    data.addData(i_L3_flag       , l3->get_flag());
    data.addData(i_L3_vetoFlag   , l3->get_vetoFlag());
    data.addData(i_L3_nTracks    , l3->nTracks());
#endif
    TEMDataMap::const_iterator    it13 = tem_data.find(ModuleId(13));
    if (it13 != tem_data.end()) {
    const TEMEvent& t13evt = *(*it13).second;
    //unused    const Level2& t13L2 = *t13evt.level2();
/*
    data.addData(i_L2_Layers     , t13L2.c_layers);
    data.addData(i_L2_Top       , t13L2.c_top);
    data.addData(i_L2_Mid       , t13L2.c_mid);
    data.addData(i_L2_Bot       , t13L2.c_bot); */
    data.addData(i_L2_Dsize     , (t13evt.dram()) ? t13evt.dram()->size() : 0);
    }
    setState(DONE);
    if( mode()==2 && bitsL2!=1 && ((bitsL1&64)==0))
        setState(FAIL);
#endif
    return;
}

void TriggerRecon::clear ()
{
   // not STL (not in SGI version)   TrigTowers.assign(MaxNumTowers, (int)0);
   // hsg: use map for TrigTowers - same functionality.
   TrigTowers.clear();

   Recon::clear();
   data.clear();
}

---