/data/refman/pasoursint/CMSSW_5_2_7_hltpatch1/src/L1Trigger/L1ExtraFromDigis/src/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
// $Id: L1ExtraTestAnalyzer.cc,v 1.4 2010/02/11 00:12:51 wmtan Exp $
//
//


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

// 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/Common/interface/Handle.h"
#include "DataFormats/L1Trigger/interface/L1EmParticle.h"
#include "DataFormats/L1Trigger/interface/L1JetParticle.h"
#include "DataFormats/L1Trigger/interface/L1MuonParticle.h"
#include "DataFormats/L1Trigger/interface/L1EtMissParticle.h"
#include "DataFormats/L1Trigger/interface/L1ParticleMap.h"
#include "DataFormats/L1Trigger/interface/L1ParticleMapFwd.h"
#include "DataFormats/L1Trigger/interface/L1HFRings.h"
#include "DataFormats/L1Trigger/interface/L1HFRingsFwd.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"

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

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

using namespace std ;

//
// class decleration
//

class L1ExtraTestAnalyzer : public edm::EDAnalyzer {
   public:
      explicit L1ExtraTestAnalyzer(const edm::ParameterSet&);
      ~L1ExtraTestAnalyzer();


      virtual void analyze(const edm::Event&, const edm::EventSetup&);
   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,
            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);