/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/L1Trigger/L1GctAnalyzer/src/DumpGctDigis.cc

Go to the documentation of this file.

#include "L1Trigger/L1GctAnalyzer/interface/DumpGctDigis.h"

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

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.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/L1CaloTrigger/interface/L1CaloCollections.h"

#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctCollections.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctEtSums.h"

using std::string;
using std::ios;
using std::endl;

//
// constructors and destructor
//
DumpGctDigis::DumpGctDigis(const edm::ParameterSet& iConfig) :
  rawLabel_( iConfig.getUntrackedParameter<edm::InputTag>("rawInput", edm::InputTag("L1GctRawDigis") ) ),
  emuRctLabel_( iConfig.getUntrackedParameter<edm::InputTag>("emuRctInput", edm::InputTag("L1RctEmuDigis") ) ),
  emuGctLabel_( iConfig.getUntrackedParameter<edm::InputTag>("emuGctInput", edm::InputTag("L1GctEmuDigis") ) ),
  outFilename_( iConfig.getUntrackedParameter<string>("outFile", "gctAnalyzer.txt") ),
  doHW_( iConfig.getUntrackedParameter<bool>("doHardware", true) ),
  doEmu_( iConfig.getUntrackedParameter<bool>("doEmulated", true) ),
  doRctEM_( iConfig.getUntrackedParameter<bool>("doRctEm", true) ),
  doEM_( iConfig.getUntrackedParameter<bool>("doEm", true) ),
  doRegions_( iConfig.getUntrackedParameter<bool>("doRegions", false) ),
  doJets_( iConfig.getUntrackedParameter<bool>("doJets", false) ),
  doInternEM_( iConfig.getUntrackedParameter<bool>("doInternEm", true) ),
  doFibres_( iConfig.getUntrackedParameter<bool>("doFibres", false) ),
  doEnergySums_( iConfig.getUntrackedParameter<bool>("doEnergySums", false) ),
  rctEmMinRank_( iConfig.getUntrackedParameter<unsigned>("rctEmMinRank", 0) )
{
  //now do what ever initialization is needed

  outFile_.open(outFilename_.c_str(), ios::out);

}


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

  outFile_.close();

}


//
// member functions
//

// ------------ method called to produce the data  ------------
void
DumpGctDigis::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  using namespace edm;
  
  outFile_ << "Run :" << iEvent.id().run() << "  Event :" << iEvent.id().event() << endl;
  
  // EM
  if (doRctEM_ && doHW_) { doRctEM(iEvent, rawLabel_); }
  if (doRctEM_ && doEmu_) { doRctEM(iEvent, emuRctLabel_); }
  if (doEM_ && doHW_) { doEM(iEvent, rawLabel_); }
  if (doEM_ && doEmu_){ doEM(iEvent, emuGctLabel_); }

  // Jets
  if (doRegions_ && doHW_) { doRegions(iEvent, rawLabel_); }
  if (doRegions_ && doEmu_) { doRegions(iEvent, emuRctLabel_); }
  if (doJets_ && doHW_) { doJets(iEvent, rawLabel_); }
  if (doJets_ && doEmu_) { doJets(iEvent, emuGctLabel_); }

  // Energy Sums
  if (doEnergySums_ && doHW_) { doEnergySums(iEvent, rawLabel_); }
  if (doEnergySums_ && doEmu_) { doEnergySums(iEvent, emuGctLabel_); }

  // debugging
  if (doInternEM_ && doHW_) { doInternEM(iEvent, rawLabel_); }
  if (doFibres_ && doHW_) { doFibres(iEvent, rawLabel_); }

}

void DumpGctDigis::doEM(const edm::Event& iEvent, const edm::InputTag& label) {

  using namespace edm;

  Handle<L1GctEmCandCollection> isoEm;
  Handle<L1GctEmCandCollection> nonIsoEm;

  L1GctEmCandCollection::const_iterator ie;
  L1GctEmCandCollection::const_iterator ne;
  
  iEvent.getByLabel(label.label(),"isoEm",isoEm);
  iEvent.getByLabel(label.label(),"nonIsoEm",nonIsoEm);

  outFile_ << "Iso EM from : " << label.label() << endl;
  for (ie=isoEm->begin(); ie!=isoEm->end(); ie++) {
    outFile_ << (*ie) 
             << " ieta(detID)=" << ie->regionId().ieta()
             << " iphi(detID)=" << ie->regionId().iphi()
             << endl;
  } 
  outFile_ << endl;
  
  outFile_ << "Non-iso EM from : " << label.label() << endl;
  for (ne=nonIsoEm->begin(); ne!=nonIsoEm->end(); ne++) {
    outFile_ << (*ne) 
             << " ieta(detID)=" << ne->regionId().ieta()
             << " iphi(detID)=" << ne->regionId().iphi()
             << endl;
  } 
  outFile_ << endl;

}

void DumpGctDigis::doRctEM(const edm::Event& iEvent, const edm::InputTag& label) {

  using namespace edm;

  Handle<L1CaloEmCollection> em;

  L1CaloEmCollection::const_iterator e;
 
  iEvent.getByLabel(label, em);

  outFile_ << "RCT EM from : " << label.label() << endl;
  for (e=em->begin(); e!=em->end(); e++) {
    if (e->rank() >= rctEmMinRank_) {
      outFile_ << (*e) 
               << " ieta(detID)=" << e->regionId().ieta()
               << " iphi(detID)=" << e->regionId().iphi()
               << endl;
    }
  } 
  outFile_ << endl;
  
}


void DumpGctDigis::doRegions(const edm::Event& iEvent, const edm::InputTag& label) {

  using namespace edm;

  Handle<L1CaloRegionCollection> rgns;

  L1CaloRegionCollection::const_iterator r;
  
  iEvent.getByLabel(label, rgns);

  outFile_ << "Regions from : " << label.label() << endl;
  for (r=rgns->begin(); r!=rgns->end(); r++) {
    outFile_ << (*r) << endl;
  } 
  outFile_ << endl;

}


void DumpGctDigis::doJets(const edm::Event& iEvent, const edm::InputTag& label) {

  using namespace edm;

  Handle<L1GctJetCandCollection> cenJets;
  Handle<L1GctJetCandCollection> forJets;
  Handle<L1GctJetCandCollection> tauJets;
  
  L1GctJetCandCollection::const_iterator cj;
  L1GctJetCandCollection::const_iterator fj;
  L1GctJetCandCollection::const_iterator tj;
  
  const std::string labelStr = label.label();
  
  iEvent.getByLabel(labelStr,"cenJets",cenJets);
  iEvent.getByLabel(labelStr,"forJets",forJets);
  iEvent.getByLabel(labelStr,"tauJets",tauJets);
  
  outFile_ << "Central jets from : " << labelStr << endl;
  for (cj=cenJets->begin(); cj!=cenJets->end(); cj++) {
    outFile_ << (*cj) << endl;
  } 
  outFile_ << endl;
  
  outFile_ << "Forward jets from : " << labelStr << endl;
  for (fj=forJets->begin(); fj!=forJets->end(); fj++) {
    outFile_ << (*fj) << endl;
  } 
  outFile_ << endl;
  
  outFile_ << "Tau jets from : " << labelStr << endl;
  for (tj=tauJets->begin(); tj!=tauJets->end(); tj++) {
    outFile_ << (*tj) << endl;
  }
}


void DumpGctDigis::doInternEM(const edm::Event& iEvent, const edm::InputTag& label) {

  using namespace edm;

  Handle<L1GctInternEmCandCollection> em;

  L1GctInternEmCandCollection::const_iterator e;
  
  iEvent.getByLabel(label, em);

  outFile_ << "Internal EM from : " << label.label() << endl;
  for (e=em->begin(); e!=em->end(); e++) {
    outFile_ << (*e) 
             << " ieta(detID)=" << e->regionId().ieta()
             << " iphi(detID)=" << e->regionId().iphi()
             << endl;
  } 
  outFile_ << endl;
  
}


void DumpGctDigis::doFibres(const edm::Event& iEvent, const edm::InputTag& label) {

  using namespace edm;

  Handle<L1GctFibreCollection> fibres;

  L1GctFibreCollection::const_iterator f;
  
  iEvent.getByLabel(label, fibres);

  outFile_ << "Fibres from : " << label.label() << endl;
  for (f=fibres->begin(); f!=fibres->end(); f++) {
    outFile_ << (*f) << endl;
  } 
  outFile_ << endl;
  
}

void DumpGctDigis::doEnergySums(const edm::Event& iEvent, const edm::InputTag& label)
{
  using namespace edm;
  
  Handle<L1GctEtTotalCollection> etTotal;
  Handle<L1GctEtHadCollection> etHad;
  Handle<L1GctEtMissCollection> etMiss;
  Handle<L1GctHtMissCollection> htMiss;
  
  iEvent.getByLabel(label, etTotal);
  iEvent.getByLabel(label, etHad);
  iEvent.getByLabel(label, etMiss);
  iEvent.getByLabel(label, htMiss);
  
  outFile_ << "Energy sums from: " << label.label() << endl;
  
  L1GctEtTotalCollection::const_iterator et;
  for (et=etTotal->begin(); et!=etTotal->end(); et++){
    outFile_ << *(et) << endl;
  }

  L1GctEtHadCollection::const_iterator ht;
  for (ht=etHad->begin(); ht!=etHad->end(); ht++){
    outFile_ << *(ht) << endl;
  }

  L1GctEtMissCollection::const_iterator met;
  for (met=etMiss->begin(); met!=etMiss->end(); met++){
    outFile_ << *(met) << endl;
  }

  L1GctHtMissCollection::const_iterator mht;
  for (mht=htMiss->begin(); mht!=htMiss->end(); mht++){
    outFile_ << *(mht) << endl;
  }

}