L1TBasicDemo.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// L1TBasicDemo:  demonstrate basic access of L1T objects
// 

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/L1Trigger/interface/EGamma.h"
#include "DataFormats/L1Trigger/interface/Tau.h"
#include "DataFormats/L1Trigger/interface/Jet.h"
#include "DataFormats/L1Trigger/interface/Muon.h"
#include "DataFormats/L1Trigger/interface/EtSum.h"
#include "DataFormats/L1Trigger/interface/EtSumHelper.h"

// leaving out the following namespaces, so namespaces are explicit in the demo code:
// using namespace l1t;
// using namespace edm;

class L1TBasicDemo : public edm::EDAnalyzer {
public:
  explicit L1TBasicDemo(const edm::ParameterSet&);
  ~L1TBasicDemo() override;
  
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
  
private:

  void analyze(edm::Event const&, edm::EventSetup const&) override;

  // EDM tokens:
  edm::EDGetTokenT<l1t::EGammaBxCollection> egToken_;
  edm::EDGetTokenT<l1t::TauBxCollection>    tauToken_;
  edm::EDGetTokenT<l1t::JetBxCollection>    jetToken_;
  edm::EDGetTokenT<l1t::EtSumBxCollection>  sumToken_;
  edm::EDGetTokenT<l1t::MuonBxCollection>   muonToken_;

  int trigger_bx_only;

};

L1TBasicDemo::L1TBasicDemo(const edm::ParameterSet& iConfig) { 
  egToken_   = consumes<l1t::EGammaBxCollection> (iConfig.getParameter<edm::InputTag>("EgTag"));
  tauToken_  = consumes<l1t::TauBxCollection>    (iConfig.getParameter<edm::InputTag>("TauTag"));
  jetToken_  = consumes<l1t::JetBxCollection>    (iConfig.getParameter<edm::InputTag>("JetTag"));
  sumToken_  = consumes<l1t::EtSumBxCollection>  (iConfig.getParameter<edm::InputTag>("SumTag"));
  muonToken_ = consumes<l1t::MuonBxCollection>   (iConfig.getParameter<edm::InputTag>("MuonTag"));
  trigger_bx_only = iConfig.getParameter<bool>("UseTriggerBxOnly");
}


L1TBasicDemo::~L1TBasicDemo(){}

void 
L1TBasicDemo::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup)
{
  cout << "INFO:  dumping EGamma BX collection:\n";
  edm::Handle<l1t::EGammaBxCollection> eg;    
  iEvent.getByToken(egToken_, eg);
  if (eg.isValid()){ 
    for (int ibx = eg->getFirstBX(); ibx <= eg->getLastBX(); ++ibx) {
      if (trigger_bx_only && (ibx != 0)) continue;  
      for (auto it=eg->begin(ibx); it!=eg->end(ibx); it++){      
    if (it->et() == 0) continue; // if you don't care about L1T candidates with zero ET.
    cout << "bx:  " << ibx << "  et:  "  << it->et() << "  eta:  "  << it->eta() << "  phi:  "  << it->phi() << "\n";
      }
    }
  } else {
    edm::LogWarning("MissingProduct") << "L1Upgrade e-gamma bx collection not found." << std::endl;
  }

  cout << "INFO:  dumping Tau BX collection:\n";
  edm::Handle<l1t::TauBxCollection> tau;    
  iEvent.getByToken(tauToken_, tau);
  if (tau.isValid()){ 
    for (int ibx = tau->getFirstBX(); ibx <= tau->getLastBX(); ++ibx) {
      if (trigger_bx_only && (ibx != 0)) continue;  
      for (auto it=tau->begin(ibx); it!=tau->end(ibx); it++){      
    if (it->et() == 0) continue; // if you don't care about L1T candidates with zero ET.
    cout << "bx:  " << ibx << "  et:  "  << it->et() << "  eta:  "  << it->eta() << "  phi:  "  << it->phi() << "\n";
      }
    }
  } else {
    edm::LogWarning("MissingProduct") << "L1Upgrade tau bx collection not found." << std::endl;
  }

  cout << "INFO:  dumping Jet BX collection:\n";
  edm::Handle<l1t::JetBxCollection> jet;    
  iEvent.getByToken(jetToken_, jet);
  if (jet.isValid()){ 
    for (int ibx = jet->getFirstBX(); ibx <= jet->getLastBX(); ++ibx) {
      if (trigger_bx_only && (ibx != 0)) continue;  
      for (auto it=jet->begin(ibx); it!=jet->end(ibx); it++){      
    if (it->et() == 0) continue; // if you don't care about L1T candidates with zero ET.
    cout << "bx:  " << ibx << "  et:  "  << it->et() << "  eta:  "  << it->eta() << "  phi:  "  << it->phi() << "\n";
      }
    }
  } else {
    edm::LogWarning("MissingProduct") << "L1Upgrade jet bx collection not found." << std::endl;
  }

  cout << "INFO:  dumping EtSum BX collection:\n";
  edm::Handle<l1t::EtSumBxCollection> sum;    
  iEvent.getByToken(sumToken_, sum);
  if (sum.isValid()){ 
    l1t::EtSumHelper hsum(sum);
    cout << "met:     " << hsum.MissingEt() << "\n";
    cout << "met phi: " << hsum.MissingEtPhi() << "\n";
    cout << "mht:     " << hsum.MissingHt() << "\n";
    cout << "mht phi: " << hsum.MissingHtPhi() << "\n";
    cout << "sum et:  " << hsum.TotalEt() << "\n";
    cout << "sum ht:  " << hsum.TotalHt() << "\n";
  } else {
    edm::LogWarning("MissingProduct") << "L1Upgrade sum bx collection not found." << std::endl;
  }

  cout << "INFO:  dumping Muon BX collection:\n";
  edm::Handle<l1t::MuonBxCollection> muon;    
  iEvent.getByToken(muonToken_, muon);
  if (muon.isValid()){ 
    for (int ibx = muon->getFirstBX(); ibx <= muon->getLastBX(); ++ibx) {
      if (trigger_bx_only && (ibx != 0)) continue;  
      for (auto it=muon->begin(ibx); it!=muon->end(ibx); it++){      
    if (it->et() == 0) continue; // if you don't care about L1T candidates with zero ET.
    cout << "bx:  " << ibx << "  et:  "  << it->et() << "  eta:  "  << it->eta() << "  phi:  "  << it->phi() << "\n";
      }
    }
  } else {
    edm::LogWarning("MissingProduct") << "L1Upgrade muon bx collection not found." << std::endl;
  }
}






void
L1TBasicDemo::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}


DEFINE_FWK_MODULE(L1TBasicDemo);