L1ExtraTestAnalyzer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    L1ExtraTestAnalyzer
// Class:      L1ExtraTestAnalyzer
//
//
// Original Author:  Werner Sun
//         Created:  Fri Jul 28 14:22:31 EDT 2006
//
//

// system include files
#include <memory>
#include <string>

// user include files
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
#include "DataFormats/L1Trigger/interface/L1EmParticle.h"
#include "DataFormats/L1Trigger/interface/L1EtMissParticle.h"
#include "DataFormats/L1Trigger/interface/L1HFRings.h"
#include "DataFormats/L1Trigger/interface/L1HFRingsFwd.h"
#include "DataFormats/L1Trigger/interface/L1JetParticle.h"
#include "DataFormats/L1Trigger/interface/L1MuonParticle.h"
#include "DataFormats/L1Trigger/interface/L1ParticleMap.h"
#include "DataFormats/L1Trigger/interface/L1ParticleMapFwd.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "TFile.h"
#include "TH1.h"

using namespace std;

//
// class decleration
//

class L1ExtraTestAnalyzer : public edm::one::EDAnalyzer<> {
public:
  explicit L1ExtraTestAnalyzer(const edm::ParameterSet &);
  ~L1ExtraTestAnalyzer() override;

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

private:
  // ----------member data ---------------------------

  edm::InputTag isoEmSource_;
  edm::InputTag nonIsoEmSource_;
  edm::InputTag cenJetSource_;
  edm::InputTag forJetSource_;
  edm::InputTag tauJetSource_;
  edm::InputTag muonSource_;
  edm::InputTag etMissSource_;
  edm::InputTag htMissSource_;
  edm::InputTag hfRingsSource_;
  edm::InputTag gtReadoutSource_;
  edm::InputTag particleMapSource_;

  TFile file_;
  TH1F hist_;
};

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
L1ExtraTestAnalyzer::L1ExtraTestAnalyzer(const edm::ParameterSet &iConfig)
    : isoEmSource_(iConfig.getParameter<edm::InputTag>("isolatedEmSource")),
      nonIsoEmSource_(iConfig.getParameter<edm::InputTag>("nonIsolatedEmSource")),
      cenJetSource_(iConfig.getParameter<edm::InputTag>("centralJetSource")),
      forJetSource_(iConfig.getParameter<edm::InputTag>("forwardJetSource")),
      tauJetSource_(iConfig.getParameter<edm::InputTag>("tauJetSource")),
      muonSource_(iConfig.getParameter<edm::InputTag>("muonSource")),
      etMissSource_(iConfig.getParameter<edm::InputTag>("etMissSource")),
      htMissSource_(iConfig.getParameter<edm::InputTag>("htMissSource")),
      hfRingsSource_(iConfig.getParameter<edm::InputTag>("hfRingsSource")),
      gtReadoutSource_(iConfig.getParameter<edm::InputTag>("gtReadoutSource")),
      particleMapSource_(iConfig.getParameter<edm::InputTag>("particleMapSource")),
      file_("l1extra.root", "RECREATE"),
      hist_("triggers",
            "Triggers",
            2 * l1extra::L1ParticleMap::kNumOfL1TriggerTypes + 1,
            -0.75,
            static_cast<double>(l1extra::L1ParticleMap::kNumOfL1TriggerTypes) + 0.5 - 0.75) {
  // now do what ever initialization is needed
}

L1ExtraTestAnalyzer::~L1ExtraTestAnalyzer() {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)

  file_.cd();
  hist_.Write();
}

//
// member functions
//

// ------------ method called to produce the data  ------------
void L1ExtraTestAnalyzer::analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) {
  using namespace edm;
  using namespace l1extra;

  static int iev = 0;
  cout << "EVENT " << ++iev << endl;

  cout << "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~" << endl;

  // Isolated EM particles
  Handle<L1EmParticleCollection> isoEmColl;
  iEvent.getByLabel(isoEmSource_, isoEmColl);
  cout << "Number of isolated EM " << isoEmColl->size() << endl;

  for (L1EmParticleCollection::const_iterator emItr = isoEmColl->begin(); emItr != isoEmColl->end(); ++emItr) {
    cout << "  p4 (" << emItr->px() << ", " << emItr->py() << ", " << emItr->pz() << ", " << emItr->energy() << ") et "
         << emItr->et() << " eta " << emItr->eta() << " phi " << emItr->phi() << endl;
  }

  // Non-isolated EM particles
  Handle<L1EmParticleCollection> nonIsoEmColl;
  iEvent.getByLabel(nonIsoEmSource_, nonIsoEmColl);
  cout << "Number of non-isolated EM " << nonIsoEmColl->size() << endl;

  for (L1EmParticleCollection::const_iterator emItr = nonIsoEmColl->begin(); emItr != nonIsoEmColl->end(); ++emItr) {
    cout << "  p4 (" << emItr->px() << ", " << emItr->py() << ", " << emItr->pz() << ", " << emItr->energy() << ") et "
         << emItr->et() << " eta " << emItr->eta() << " phi " << emItr->phi() << endl;
  }

  // Jet particles
  Handle<L1JetParticleCollection> cenJetColl;
  iEvent.getByLabel(cenJetSource_, cenJetColl);
  cout << "Number of central jets " << cenJetColl->size() << endl;

  for (L1JetParticleCollection::const_iterator jetItr = cenJetColl->begin(); jetItr != cenJetColl->end(); ++jetItr) {
    cout << "  p4 (" << jetItr->px() << ", " << jetItr->py() << ", " << jetItr->pz() << ", " << jetItr->energy()
         << ") et " << jetItr->et() << " eta " << jetItr->eta() << " phi " << jetItr->phi() << endl;
  }

  Handle<L1JetParticleCollection> forJetColl;
  iEvent.getByLabel(forJetSource_, forJetColl);
  cout << "Number of forward jets " << forJetColl->size() << endl;

  for (L1JetParticleCollection::const_iterator jetItr = forJetColl->begin(); jetItr != forJetColl->end(); ++jetItr) {
    cout << "  p4 (" << jetItr->px() << ", " << jetItr->py() << ", " << jetItr->pz() << ", " << jetItr->energy()
         << ") et " << jetItr->et() << " eta " << jetItr->eta() << " phi " << jetItr->phi() << endl;
  }

  Handle<L1JetParticleCollection> tauColl;
  iEvent.getByLabel(tauJetSource_, tauColl);
  cout << "Number of tau jets " << tauColl->size() << endl;

  for (L1JetParticleCollection::const_iterator tauItr = tauColl->begin(); tauItr != tauColl->end(); ++tauItr) {
    cout << "  p4 (" << tauItr->px() << ", " << tauItr->py() << ", " << tauItr->pz() << ", " << tauItr->energy()
         << ") et " << tauItr->et() << " eta " << tauItr->eta() << " phi " << tauItr->phi() << endl;
  }

  // Muon particles
  Handle<L1MuonParticleCollection> muColl;
  iEvent.getByLabel(muonSource_, muColl);
  cout << "Number of muons " << muColl->size() << endl;

  for (L1MuonParticleCollection::const_iterator muItr = muColl->begin(); muItr != muColl->end(); ++muItr) {
    cout << "  q " << muItr->charge() << " p4 (" << muItr->px() << ", " << muItr->py() << ", " << muItr->pz() << ", "
         << muItr->energy() << ") et " << muItr->et() << " eta " << muItr->eta() << endl
         << "    phi " << muItr->phi() << " iso " << muItr->isIsolated() << " mip " << muItr->isMip() << " fwd "
         << muItr->isForward() << " rpc " << muItr->isRPC() << endl;
  }

  // MET
  Handle<L1EtMissParticleCollection> etMissColl;
  iEvent.getByLabel(etMissSource_, etMissColl);
  cout << "MET Coll (" << etMissColl->begin()->px() << ", " << etMissColl->begin()->py() << ", "
       << etMissColl->begin()->pz() << ", " << etMissColl->begin()->energy() << ") phi " << etMissColl->begin()->phi()
       << " EtTot " << etMissColl->begin()->etTotal() << endl;

  // MHT
  Handle<L1EtMissParticleCollection> htMissColl;
  iEvent.getByLabel(htMissSource_, htMissColl);
  cout << "MHT Coll (" << htMissColl->begin()->px() << ", " << htMissColl->begin()->py() << ", "
       << htMissColl->begin()->pz() << ", " << htMissColl->begin()->energy() << ") phi " << htMissColl->begin()->phi()
       << " HtTot " << htMissColl->begin()->etTotal() << endl;

  // HF Rings
  Handle<L1HFRingsCollection> hfRingsColl;
  iEvent.getByLabel(hfRingsSource_, hfRingsColl);
  cout << "HF Rings:" << endl;
  for (int i = 0; i < L1HFRings::kNumRings; ++i) {
    cout << "  " << i << ": et sum = " << hfRingsColl->begin()->hfEtSum((L1HFRings::HFRingLabels)i)
         << ", bit count = " << hfRingsColl->begin()->hfBitCount((L1HFRings::HFRingLabels)i) << endl;
  }
  cout << endl;

  //    // L1GlobalTriggerReadoutRecord
  //    Handle< L1GlobalTriggerReadoutRecord > gtRecord ;
  //    iEvent.getByLabel( gtReadoutSource_, gtRecord ) ;
  //    cout << "Global trigger decision " << gtRecord->decision() << endl ;

  cout << endl;
}

// define this as a plug-in
DEFINE_FWK_MODULE(L1ExtraTestAnalyzer);