GctTimingAnalyzer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    GctTimingAnalyzer
// Class:      GctTimingAnalyzer
// 
//
// Original Author:  Alex Tapper
//         Created:  Mon Apr 21 14:21:06 CEST 2008
//
//

// Include file
#include "L1Trigger/L1GctAnalyzer/interface/GctTimingAnalyzer.h"

GctTimingAnalyzer::GctTimingAnalyzer(const edm::ParameterSet& iConfig):
  m_outputFileName(iConfig.getUntrackedParameter<std::string>("outFile", "gctTiming.txt")),
  m_gctSource(iConfig.getUntrackedParameter<edm::InputTag>("gctSource")),
  m_isoEmSource(iConfig.getUntrackedParameter<edm::InputTag>("isoEmSource")),
  m_nonIsoEmSource(iConfig.getUntrackedParameter<edm::InputTag>("nonIsoEmSource")),
  m_cenJetsSource(iConfig.getUntrackedParameter<edm::InputTag>("cenJetsSource")),
  m_forJetsSource(iConfig.getUntrackedParameter<edm::InputTag>("forJetsSource")),
  m_tauJetsSource(iConfig.getUntrackedParameter<edm::InputTag>("tauJetsSource")),
  m_doInternal(iConfig.getUntrackedParameter<bool>("doInternal")),
  m_doElectrons(iConfig.getUntrackedParameter<bool>("doElectrons")),
  m_doJets(iConfig.getUntrackedParameter<bool>("doJets")),
  m_doHFRings(iConfig.getUntrackedParameter<bool>("doHFRings")),
  m_doESums(iConfig.getUntrackedParameter<bool>("doESums")),
  m_evtNum(0)
{
  m_outputFile.open(m_outputFileName.c_str());
}

GctTimingAnalyzer::~GctTimingAnalyzer()
{
  m_outputFile.close();
}

void GctTimingAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  using namespace edm;
  using namespace std;
  
  // Electrons
  if(m_doElectrons){

    // Isolated EM cands in GCT output
    Handle<L1GctEmCandCollection> isoEm; 
    iEvent.getByLabel(m_isoEmSource,isoEm);    

    for (L1GctEmCandCollection::const_iterator em=isoEm->begin(); em!=isoEm->end(); em++){
      if (em->rank()>0) {
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*em) << std::endl; 
      }
    }
  
    // Non-Isolated EM cands in GCT output
    Handle<L1GctEmCandCollection> nonIsoEm; 
    iEvent.getByLabel(m_nonIsoEmSource,nonIsoEm);    
    
    for (L1GctEmCandCollection::const_iterator em=nonIsoEm->begin(); em!=nonIsoEm->end(); em++){
      if (em->rank()>0) {
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*em) << std::endl; 
      }
    }

    if (m_doInternal){
      // Internal GCT EM cands
      Handle<L1GctInternEmCandCollection> internEm; 
      iEvent.getByLabel(m_gctSource,internEm);    

      if (internEm.isValid()) {
        for (L1GctInternEmCandCollection::const_iterator em=internEm->begin(); em!=internEm->end(); em++){
          if (em->rank()>0) {
            m_outputFile << "BX = " << dec << m_evtNum << " " << (*em) << std::endl; 
          }
        }
      }
    }

    // RCT EM cands
    Handle<L1CaloEmCollection> rctEm; 
    iEvent.getByLabel(m_gctSource,rctEm);    
    
    for (L1CaloEmCollection::const_iterator em=rctEm->begin(); em!=rctEm->end(); em++){
      if (em->rank()>0) {
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*em) << std::endl; 
      }
    }
  }     
  
  // Jets
  if (m_doJets){

    // Central jet cands in GCT output
    Handle<L1GctJetCandCollection> cenJets; 
    iEvent.getByLabel(m_cenJetsSource,cenJets);    

    for (L1GctJetCandCollection::const_iterator cj=cenJets->begin(); cj!=cenJets->end(); cj++){
      if (cj->rank()>0) {
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*cj) << std::endl; 
      }
    }
    
    // Forward jet cands in GCT output
    Handle<L1GctJetCandCollection> forJets; 
    iEvent.getByLabel(m_forJetsSource,forJets);    
    
    for (L1GctJetCandCollection::const_iterator fj=forJets->begin(); fj!=forJets->end(); fj++){
      if (fj->rank()>0) {
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*fj) << std::endl; 
      }
    }
    
    // Tau jet cands in GCT output
    Handle<L1GctJetCandCollection> tauJets; 
    iEvent.getByLabel(m_tauJetsSource,tauJets);    
    
    for (L1GctJetCandCollection::const_iterator tj=tauJets->begin(); tj!=tauJets->end(); tj++){
      if (tj->rank()>0) {
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*tj) << std::endl; 
      }
    }
  
    if (m_doInternal){

      // Internal GCT jet cands
      Handle<L1GctInternJetDataCollection> internJets; 
      iEvent.getByLabel(m_gctSource,internJets);    
    
      if (internJets.isValid()){
        for (L1GctInternJetDataCollection::const_iterator j=internJets->begin(); j!=internJets->end(); j++){
          if ((j->et()>0) || (j->rank()>0)) {
            m_outputFile << "BX = " << dec << m_evtNum << " " << (*j) << std::endl; 
          }
        }
      }
    }
  }
  
  // RCT regions
  Handle<L1CaloRegionCollection> rctRn; 
  iEvent.getByLabel(m_gctSource,rctRn);    
    
  for (L1CaloRegionCollection::const_iterator rn=rctRn->begin(); rn!=rctRn->end(); rn++){
    if (rn->et()>0) {
      m_outputFile << "BX = " << dec << m_evtNum << " " << (*rn) << std::endl; 
    }
  }

  // HF Rings
  if (m_doHFRings){

    // HFRing counts
    Handle<L1GctHFBitCountsCollection> hfBitCnt;
    iEvent.getByLabel(m_gctSource,hfBitCnt);

    for (L1GctHFBitCountsCollection::const_iterator jc=hfBitCnt->begin(); jc!=hfBitCnt->end(); jc++){
      if (jc->bitCount(0) || jc->bitCount(1) || jc->bitCount(2) || jc->bitCount(3)){
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*jc) << std::endl; 
      }
    }

    // HFRing Et sums
    Handle<L1GctHFRingEtSumsCollection> hfEtSums;
    iEvent.getByLabel(m_gctSource,hfEtSums);

    for (L1GctHFRingEtSumsCollection::const_iterator js=hfEtSums->begin(); js!=hfEtSums->end(); js++){
      if (js->etSum(0) || js->etSum(1) || js->etSum(2) || js->etSum(3)){
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*js) << std::endl; 
      }
    }

    if (m_doInternal){
      // Internal HF data
      Handle<L1GctInternHFDataCollection> internHF;
      iEvent.getByLabel(m_gctSource,internHF);
      
      if (internHF.isValid()){
        for (L1GctInternHFDataCollection::const_iterator hf=internHF->begin(); hf!=internHF->end(); hf++){
          if (hf->value(0) || hf->value(1) || hf->value(2) || hf->value(3)){
            m_outputFile << "BX = " << dec << m_evtNum << " " << (*hf) << std::endl; 
          }
        }
      }
    }    
  }

  // HT, MET and ET
  if (m_doESums){

    // MET
    Handle<L1GctEtMissCollection> missEt;
    iEvent.getByLabel(m_gctSource,missEt);

    for (L1GctEtMissCollection::const_iterator met=missEt->begin(); met!=missEt->end(); met++){
      if (met->et()>0){
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*met) << std::endl; 
      }
    }

    // Total Et
    Handle<L1GctEtTotalCollection> totEt;
    iEvent.getByLabel(m_gctSource,totEt);
    
    for (L1GctEtTotalCollection::const_iterator tet=totEt->begin(); tet!=totEt->end(); tet++){
      if (tet->et()>0){
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*tet) << std::endl; 
      }
    }
    
    // Ht
    Handle<L1GctEtHadCollection> hadEt;
    iEvent.getByLabel(m_gctSource,hadEt);

    for (L1GctEtHadCollection::const_iterator ht=hadEt->begin(); ht!=hadEt->end(); ht++){
      if (ht->et()>0){
        m_outputFile << "BX = " << dec << m_evtNum << " " << (*ht) << std::endl; 
      }
    }

    if (m_doInternal){
      // Internal GCT Et sums
      Handle<L1GctInternEtSumCollection> Et;
      iEvent.getByLabel(m_gctSource,Et);
      
      if (Et.isValid()){
        for (L1GctInternEtSumCollection::const_iterator e=Et->begin(); e!=Et->end(); e++){
          if (e->et()>0){
            m_outputFile << "BX = " << dec << m_evtNum << " " << (*e) << std::endl;
          }
        }
      }
    }    
  }

  m_evtNum++;

}