---

XANATriggerInfo.cpp

Go to the documentation of this file.

//-------------------------------------------------
// from RunL1GlobalTrigger
//-------------------------------------------------

#include "Utilities/Configuration/interface/Architecture.h"
#include "Utilities/Notification/interface/PackageInitializer.h"
#include "Utilities/Notification/interface/Observer.h" 
#include "Utilities/Notification/interface/Singleton.h"

#include "CARF/G3Event/interface/G3EventProxy.h"

#include "Trigger/L1GlobalTrigger/interface/L1GlobalTrigger.h"
#include "Trigger/L1GlobalTrigger/interface/L1GlobalTriggerSetup.h"

#include "HLT/HLTSteering/interface/HighLevelTriggerSetup.h"
#include "HLT/HLTSteering/interface/HighLevelTrigger.h"
#include "HLT/HLTElements/interface/HighLevelTriggerOrNode.h"
#include "HLT/HLTElements/interface/HighLevelTriggerAndNode.h"
#include "HLT/HLTElements/interface/HighLevelTriggerRandomElement.h"
#include "HLT/HLTElements/interface/HLGL2L25EleTrigger.h"

#include <XANADOO/XANATriggerInfo/interface/XANATriggerInfo.h>
#include <XANADOO/XANAUtilities/interface/XANAToolkit.h>

#include <Utilities/UI/interface/SimpleConfigurable.h> 

#include <TBits.h>
#include <iostream>
#include <vector>

class UserAnalysis : private Observer<G3EventProxy*> {
  
public:
  
  UserAnalysis() { 

    init();

    eventsAnalysed  = 0;
    runsAnalysed    = 0;
    lastrun         = 0;
    eventsTriggered = 0;
  }
  
  virtual ~UserAnalysis() {
    
    xana::out(1) << endl;
    xana::out(1) << endl;
    xana::out(1) << "**** L1GlobalTrigger report "  << endl;
    xana::out(1) << " Number of runs analysed:    " << runsAnalysed    << endl;
    xana::out(1) << " Number of events analysed:  " << eventsAnalysed  << endl;
    xana::out(1) << " Number of triggered events: " << eventsTriggered << endl;
    xana::out(1) << endl;
    xana::out(1) << endl;
  }
  
  void analysis(G3EventProxy* ev) {
    
    //
    // some statistics: count events and runs processed
    //
    eventsAnalysed++;
    if (ev->simSignal()->id().runNumber() != lastrun) {
      lastrun = (unsigned int) ev->simSignal()->id().runNumber();
      runsAnalysed++;
    }
    
    xana::out(1) << endl;
    xana::out(1) << "**** L1GlobalTrigger processing event ";
    xana::out(1) << ev->simSignal()->id().runNumber() << ":"
         << ev->simSignal()->id().eventInRun();
    xana::out(1) << " ****" << endl;
    
    //
    // Set up the L1 Global Trigger 
    //    
    L1GlobalTriggerSetup* setup = Singleton<L1GlobalTriggerSetup>::instance();
    L1GlobalTrigger* trig = setup->GlobalTrigger();
    
    //
    // Run and print L1 Global Trigger result (L1 Accept word)
    // 
    trig->Result().print();      
    
    //
    // Print L1GlobalTrigger input (L1 Trigger Objects)
    // 
    trig->Result().printL1Objects();
    
    //
    // Retrieve one single object info 
    //
    // (second muon)
    //
    xana::out(1) << endl;
    xana::out(1) << "**** RunL1GlobalTrigger:: second muon: " << endl;
    trig->muonCand( 1 ).print();
    xana::out(1) << endl;

    //
    // Retrieve one single object info 
    //
    // (first isolated electron)
    //
    xana::out(1) << endl;
    xana::out(1) << "**** RunL1GlobalTrigger:: first isolated electron: " << endl;
    trig->isolatedElectronCand( 0 ).print();
    xana::out(1) << endl;

    //
    // Count number of triggered events
    //
    if( trig->Result().getGlobalDecision() ) eventsTriggered++;
   
    bool decision = trig->Result().getGlobalDecision();
    const BitArray<128> response = trig->Result().getDecision();
    xana::out(1) << "L1Trigger global decision : " << decision << endl;
    xana::out(1) << "L1Trigger response : " ;
    TBits bits(128);
    for (int i = 0; i < 128; i++) {
      xana::out(1) << (response.element(i) ? '1' : '0');
      bits.SetBitNumber(i, response.element(i));        
    }   
    xana::out(1) << endl;
    SimpleConfigurable<string> l1Menu("lumi2x1033","L1Globaltrigger:TriggerMenu");
    XANATriggerInfo *trigger = new XANATriggerInfo(decision, bits, l1Menu.value());
    xana::out(1) << "===========================================================" << endl;
    
  }
  
private :
  
  void upDate(G3EventProxy * ev) { if (ev != 0) analysis(ev); }
  
private :
  
  unsigned int eventsAnalysed;  // count events that have been analysed
  unsigned int eventsTriggered; // count evennts passed trigger requirement
  unsigned int runsAnalysed;    // count runs
  unsigned int lastrun;
  
};


PKBuilder<UserAnalysis> user("L1GlobalTrigger");     

---