L1RCTRelValAnalyzer.cc

Go to the documentation of this file.

// system include files
#include <memory>

// 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/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/L1CaloTrigger/interface/L1CaloCollections.h"

#include "L1Trigger/RegionalCaloTrigger/interface/L1RCTRelValAnalyzer.h"

using std::string;
using std::cout;
using std::endl;

//
// constructors and destructor
//
L1RCTRelValAnalyzer::L1RCTRelValAnalyzer(const edm::ParameterSet& iConfig) :
  rctEmCandsLabel(iConfig.getParameter<edm::InputTag>("rctEmCandsLabel")),
  rctRegionsLabel(iConfig.getParameter<edm::InputTag>("rctRegionsLabel"))
{
   //now do what ever initialization is needed

  edm::Service<TFileService> fs;
  h_emRank = fs->make<TH1F>( "emRank", "emRank", 64, 0., 64. );
  h_emIeta = fs->make<TH1F>( "emOccupancyIeta", "emOccupancyIeta", 22, 0., 22. );
  h_emIphi = fs->make<TH1F>( "emOccupancyIphi", "emOccupancyIphi", 18, 0., 18. );
  h_emIsoOccIetaIphi = fs->make<TH2F>( "emIsoOccupancy2D", "emIsoOccupancy2D", 22, 0., 22.,
                       18, 0., 18. );
  h_emNonIsoOccIetaIphi = fs->make<TH2F>( "emNonIsoOccupancy2D", "emNonIsoOccupancy2D", 22, 0., 
                      22., 18, 0., 18. );
  
  h_regionSum = fs->make<TH1F>( "regionSum", "regionSum", 100, 0., 100. );
  h_regionSumIetaIphi = fs->make<TH2F>( "regionSumEtWeighted2D", "regionSumEtWeighted2D", 22, 
                      0., 22., 18, 0., 18. );
  h_regionOccIetaIphi = fs->make<TH2F>( "regionOccupancy2D", "regionOccupancy2D", 
                    22, 0., 22., 18, 0., 18. );
}


L1RCTRelValAnalyzer::~L1RCTRelValAnalyzer()
{

   // do anything here that needs to be done at destruction time
   // (e.g. close files, deallocate resources etc.)

}


//
// member functions
//

// ------------ method called to produce the data  ------------
void
L1RCTRelValAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{
   using namespace edm;
#ifdef THIS_IS_AN_EVENT_EXAMPLE
   Handle<ExampleData> pIn;
   iEvent.getByLabel("example",pIn);
#endif

#ifdef THIS_IS_AN_EVENTSETUP_EXAMPLE
   ESHandle<SetupData> pSetup;
   iSetup.get<SetupRecord>().get(pSetup);
#endif

   // as in L1GctTestAnalyzer.cc
   Handle<L1CaloEmCollection> rctEmCands;
   Handle<L1CaloRegionCollection> rctRegions;

   L1CaloEmCollection::const_iterator em;
   L1CaloRegionCollection::const_iterator rgn;

   iEvent.getByLabel(rctEmCandsLabel, rctEmCands);
   iEvent.getByLabel(rctRegionsLabel, rctRegions);

   for (em=rctEmCands->begin(); em!=rctEmCands->end(); em++)
     {
       if ((*em).rank() > 0)
     {
       h_emRank->Fill( (*em).rank() );
       h_emIeta->Fill( (*em).regionId().ieta() );
       h_emIphi->Fill( (*em).regionId().iphi() );
       if ((*em).isolated())
         {
           h_emIsoOccIetaIphi->Fill( (*em).regionId().ieta(),
                    (*em).regionId().iphi() );
         }
       else
         {
           h_emNonIsoOccIetaIphi->Fill( (*em).regionId().ieta(),
                       (*em).regionId().iphi() );
         }
     }
     }
   
   for (rgn=rctRegions->begin(); rgn!=rctRegions->end(); rgn++){
     if ( (*rgn).et() > 0 )
       {
     h_regionSum->Fill( (*rgn).et() );
     h_regionSumIetaIphi->Fill( (*rgn).gctEta(), (*rgn).gctPhi(),
                      (*rgn).et() );
     h_regionOccIetaIphi->Fill( (*rgn).gctEta(), (*rgn).gctPhi() );
       }
   }
}