#include <L1Trigger/RegionalCaloTrigger/plugins/L1RCTRelValAnalyzer.cc>

Inheritance diagram for L1RCTRelValAnalyzer:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
	L1RCTRelValAnalyzer (const edm::ParameterSet &)
	~L1RCTRelValAnalyzer ()
Private Attributes
TH1F *	h_emIeta
TH1F *	h_emIphi
TH2F *	h_emIsoOccIetaIphi
TH2F *	h_emNonIsoOccIetaIphi
TH1F *	h_emRank
TH2F *	h_regionOccIetaIphi
TH1F *	h_regionSum
TH2F *	h_regionSumIetaIphi
edm::InputTag	rctEmCandsLabel
edm::InputTag	rctRegionsLabel

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 45 of file L1RCTRelValAnalyzer.h.

Constructor & Destructor Documentation

L1RCTRelValAnalyzer::L1RCTRelValAnalyzer ( const edm::ParameterSet & iConfig ) [explicit]

Definition at line 25 of file L1RCTRelValAnalyzer.cc.

References h_emIeta, h_emIphi, h_emIsoOccIetaIphi, h_emNonIsoOccIetaIphi, h_emRank, h_regionOccIetaIphi, h_regionSum, and h_regionSumIetaIphi.

                                                                       :
  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 ( )

Definition at line 48 of file L1RCTRelValAnalyzer.cc.

{

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

}

Member Function Documentation

void L1RCTRelValAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		`[virtual]`

Implements edm::EDAnalyzer.

Definition at line 63 of file L1RCTRelValAnalyzer.cc.

References edm::EventSetup::get(), edm::Event::getByLabel(), h_emIeta, h_emIphi, h_emIsoOccIetaIphi, h_emNonIsoOccIetaIphi, h_emRank, h_regionOccIetaIphi, h_regionSum, h_regionSumIetaIphi, rctEmCandsLabel, and rctRegionsLabel.

{
   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() );
       }
   }
}