Trigger.cpp

Go to the documentation of this file.

#include "cafe/Trigger.hpp"

#include <algorithm>
#include "cafe/Config.hpp"
#include "cafe/Event.hpp"
#include "tmb_tree/TMBTrigger.hpp"
#include "tmb_tree/TMBTrigBits.hpp"

using namespace std;

namespace cafe {

    Trigger::Trigger(const char *name)
        : Processor(name),
          _vars("_trgname"),
          _numEvents(0),
          _numPassed(0),
          _stat(0)
    {
        Config config(name);
        vector<string> temp = config.getVString("Triggers", " ,");
        _Tracking = config.get("Accounting", false);

        for(int i = 0; i < temp.size(); i++) {
            _triggerNames[temp.at(i)] = 0;
        }

        _numEvents = 0;
        _numPassed = 0;
        if(_triggerNames.size() == 0) {
            warn() << "Trigger[" << name << "] No triggers specified ? All events will fail" 
                   << endl;
        }

        out() << "Trigger[" << name << "] Triggers = ";
        for(map<string, long>::const_iterator it = _triggerNames.begin(); 
            it != _triggerNames.end(); ++it) {
            std::string tmpname = it->first;
            out() << it->first << ' ';
        }
        out() << endl;
        _useTrigBits=config.get("UseTrigBits",false);
        if(_useTrigBits) out()<<"Will use TMBTrigBits instead of TMBTriggers"<<endl;
        _required=config.get("Required",true);
        if(!_required) out()<<"Will veto on trigger list"<<endl;
    }
        
    bool Trigger::processEvent(Event& event)
    {
        ++_numEvents;

        _stat = 0;
        event.get("StatPointer", _stat);
        bool trig_found = !_required;
        Collection<TMBTrigger> triggers = event.getTriggers(_vars);
        if(triggers.size()>0 && !_useTrigBits){

        //
        // This is proportional to the number of triggers in the event
        // and about log(numUserTriggers).
        //


          for(Collection<TMBTrigger>::const_iterator it = triggers.begin(); 
              it != triggers.end();
            ++it) {
            if(_triggerNames.find((*it).getTrgName()) != _triggerNames.end()) {
                trig_found = _required;

                if(_Tracking) {

                    // Find the trigger in the list and increment the counter
                        
                    std::string tmpname = (*it).getTrgName();
                    ++_triggerNames[tmpname];
                }
            }
          }

        }else{
          const TMBTrigBits* trigbits = event.getTrigBits();
          if(!trigbits) return false;
          for(map<string,long>::const_iterator it = _triggerNames.begin();
              it != _triggerNames.end(); ++it){
            std::string tmpname = it->first;
            if(trigbits->triggerFired(tmpname)){
              trig_found=_required;
              if(_Tracking)  ++_triggerNames[tmpname];
            }
          }
        }

        if(trig_found) {
            ++_numPassed;
            if(_stat != 0) {
                _stat->EventSelected(name()); 
            }
            return true;
        }

        return false;
    }


    void Trigger::finish() {

        if(_numEvents > 0) {
            out() << "============================================\n";
            out() << "  Trigger[" << name() << "] SUMMARY\n\n";
            out() << "  " << _numEvents << " events processed.\n";
            out() << "  " << _numPassed << " events passed, pass rate = "
                  << double(_numPassed) / double(_numEvents) * 100.
                  << "%\n\n";
            if(_Tracking) {
                out() << "    Individual Trigger Efficencies\n\n";
             
                for(map<string, long>::const_iterator it = _triggerNames.begin(); 
                    it != _triggerNames.end(); ++it) {
                    std::string tmpname = it->first;
                    out() << "    " << tmpname << ": " 
                          << (double(it->second) / _numEvents) * 100.0 
                          << "% (" << it->second << " events)\n";
                }
                out() << '\n';
            }

            out() << "============================================" << endl;
        }
    }
}
ClassImp(cafe::Trigger)

---