#include <PreshowerAnalyzer.h>

Inheritance diagram for PreshowerAnalyzer:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

virtual void	endJob ()

	PreshowerAnalyzer (const edm::ParameterSet &)

	~PreshowerAnalyzer ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Attributes
double	calib_planeX_

double	calib_planeY_

float	EmaxDE_

float	EmaxSC_

float	EminDE_

float	EminSC_

double	gamma_

TH1F *	h1_esDeltaE

TH1F *	h1_esE_x

TH1F *	h1_esE_y

TH1F *	h1_esEta_x

TH1F *	h1_esEta_y

TH1F *	h1_esNhits_x

TH1F *	h1_esNhits_y

TH1F *	h1_esPhi_x

TH1F *	h1_esPhi_y

TH1F *	h1_islandEESCEnergy1

TH1F *	h1_islandEESCEnergy2

TH1F *	h1_nclu_x

TH1F *	h1_nclu_y

std::string	islandEndcapSuperClusterCollection1_

std::string	islandEndcapSuperClusterCollection2_

std::string	islandEndcapSuperClusterProducer1_

std::string	islandEndcapSuperClusterProducer2_

double	mip_

int	nBinDE_

int	nBinSC_

int	nEvt_

std::string	outputFile_

std::string	preshClusterCollectionX_

std::string	preshClusterCollectionY_

std::string	preshClusterProducer_

TFile *	rootFile_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Description: Analyzer to fetch collection of objects from event and make simple plots

Implementation: \

Author: : Shahram Rahatlou, May 2006

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

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

Definition at line 35 of file PreshowerAnalyzer.h.

Constructor & Destructor Documentation

PreshowerAnalyzer::PreshowerAnalyzer ( const edm::ParameterSet & ps )

explicit

Definition at line 30 of file PreshowerAnalyzer.cc.

References calib_planeX_, calib_planeY_, EmaxDE_, EmaxSC_, EminDE_, EminSC_, gamma_, edm::ParameterSet::getParameter(), islandEndcapSuperClusterCollection1_, islandEndcapSuperClusterCollection2_, islandEndcapSuperClusterProducer1_, islandEndcapSuperClusterProducer2_, mip_, nBinDE_, nBinSC_, nEvt_, outputFile_, preshClusterCollectionX_, preshClusterCollectionY_, preshClusterProducer_, and rootFile_.

 {
  
   EminDE_ = ps.getParameter<double>("EminDE");
   EmaxDE_ = ps.getParameter<double>("EmaxDE");
   nBinDE_ = ps.getParameter<int>("nBinDE");
 
   EminSC_ = ps.getParameter<double>("EminSC");
   EmaxSC_ = ps.getParameter<double>("EmaxSC");
   nBinSC_ = ps.getParameter<int>("nBinSC");
 
   preshClusterCollectionX_ = ps.getParameter<std::string>("preshClusterCollectionX");
   preshClusterCollectionY_ = ps.getParameter<std::string>("preshClusterCollectionY");
   preshClusterProducer_   = ps.getParameter<std::string>("preshClusterProducer");
 
   islandEndcapSuperClusterCollection1_ = ps.getParameter<std::string>("islandEndcapSuperClusterCollection1");
   islandEndcapSuperClusterProducer1_   = ps.getParameter<std::string>("islandEndcapSuperClusterProducer1");
 
   islandEndcapSuperClusterCollection2_ = ps.getParameter<std::string>("islandEndcapSuperClusterCollection2");
   islandEndcapSuperClusterProducer2_   = ps.getParameter<std::string>("islandEndcapSuperClusterProducer2");
 
   outputFile_   = ps.getParameter<std::string>("outputFile");
   rootFile_ = TFile::Open(outputFile_.c_str(),"RECREATE"); // open output file to store histograms
 
   // calibration parameters:
   calib_planeX_ = ps.getParameter<double>("preshCalibPlaneX");
   calib_planeY_ = ps.getParameter<double>("preshCalibPlaneY");
   gamma_        = ps.getParameter<double>("preshCalibGamma");
   mip_          = ps.getParameter<double>("preshCalibMIP");
 
   nEvt_ = 0; 
 }

PreshowerAnalyzer::~PreshowerAnalyzer ( )

Definition at line 65 of file PreshowerAnalyzer.cc.

References rootFile_.

 {
   delete rootFile_;
 }

Member Function Documentation

void PreshowerAnalyzer::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

virtual

Implements edm::EDAnalyzer.

Definition at line 96 of file PreshowerAnalyzer.cc.

References calib_planeX_, calib_planeY_, gamma_, edm::Event::getByLabel(), h1_esDeltaE, h1_esE_x, h1_esE_y, h1_esEta_x, h1_esEta_y, h1_esNhits_x, h1_esNhits_y, h1_esPhi_x, h1_esPhi_y, h1_islandEESCEnergy1, h1_islandEESCEnergy2, h1_nclu_x, h1_nclu_y, islandEndcapSuperClusterCollection1_, islandEndcapSuperClusterCollection2_, islandEndcapSuperClusterProducer1_, islandEndcapSuperClusterProducer2_, mip_, nEvt_, preshClusterCollectionX_, preshClusterCollectionY_, and preshClusterProducer_.

                                                                          {
 //========================================================================
 
   using namespace edm; // needed for all fwk related classe
 
   //std::cout << "\n .......  Event # " << nEvt_+1 << " is analyzing ....... " << std::endl << std::endl;
 
   // Get island super clusters
   Handle<reco::SuperClusterCollection> pIslandEndcapSuperClusters1;
   evt.getByLabel(islandEndcapSuperClusterProducer1_, islandEndcapSuperClusterCollection1_, pIslandEndcapSuperClusters1);
   const reco::SuperClusterCollection* islandEndcapSuperClusters1 = pIslandEndcapSuperClusters1.product();
   //std::cout << "\n islandEndcapSuperClusters1->size() = " << islandEndcapSuperClusters1->size() << std::endl;
 
   // loop over the super clusters and fill the histogram
   for(reco::SuperClusterCollection::const_iterator aClus = islandEndcapSuperClusters1->begin();
                                                     aClus != islandEndcapSuperClusters1->end(); aClus++) {
     h1_islandEESCEnergy1->Fill( aClus->energy() );
   }
 
   // Get island super clusters
   Handle<reco::SuperClusterCollection> pIslandEndcapSuperClusters2;
   evt.getByLabel(islandEndcapSuperClusterProducer2_, islandEndcapSuperClusterCollection2_, pIslandEndcapSuperClusters2);
   const reco::SuperClusterCollection* islandEndcapSuperClusters2 = pIslandEndcapSuperClusters2.product();
   //std::cout << "\n islandEndcapSuperClusters2->size() = " << islandEndcapSuperClusters2->size() << std::endl;
 
   // loop over the super clusters and fill the histogram
   for(reco::SuperClusterCollection::const_iterator aClus = islandEndcapSuperClusters2->begin();
                                                     aClus != islandEndcapSuperClusters2->end(); aClus++) {
     h1_islandEESCEnergy2->Fill( aClus->energy() );
   }
 
 
   // Get ES clusters in X plane
   Handle<reco::PreshowerClusterCollection> pPreshowerClustersX;
   evt.getByLabel(preshClusterProducer_, preshClusterCollectionX_, pPreshowerClustersX);
   const reco::PreshowerClusterCollection *clustersX = pPreshowerClustersX.product();
   h1_nclu_x->Fill( clustersX->size() );
   //std::cout << "\n pPreshowerClustersX->size() = " << clustersX->size() << std::endl;
 
   Handle<reco::PreshowerClusterCollection> pPreshowerClustersY;
   evt.getByLabel(preshClusterProducer_, preshClusterCollectionY_, pPreshowerClustersY);
   const reco::PreshowerClusterCollection *clustersY = pPreshowerClustersY.product();
   h1_nclu_y->Fill( clustersY->size() );
   //std::cout << "\n pPreshowerClustersY->size() = " << clustersY->size() << std::endl;
 
 
   // loop over the ES clusters and fill the histogram
   float e1 = 0;
   for(reco::PreshowerClusterCollection::const_iterator esClus = clustersX->begin();
                                                        esClus !=clustersX->end(); esClus++) {
       e1 += esClus->energy();
       h1_esE_x->Fill( esClus->energy() );  
       h1_esEta_x->Fill( esClus->eta() );
       h1_esPhi_x->Fill( esClus->phi() );
       h1_esNhits_x->Fill( esClus->nhits() );     
   }
 
   float e2 = 0;
   for(reco::PreshowerClusterCollection::const_iterator esClus = clustersY->begin();
                                                        esClus !=clustersY->end(); esClus++) {
       e2 += esClus->energy();
       h1_esE_y->Fill( esClus->energy() );  
       h1_esEta_y->Fill( esClus->eta() );
       h1_esPhi_y->Fill( esClus->phi() );
       h1_esNhits_y->Fill( esClus->nhits() );     
   }
   
   float deltaE = 0;       
   if(e1+e2 > 1.0e-10) {
       // GeV to #MIPs
       e1 = e1 / mip_;
       e2 = e2 / mip_;
       deltaE = gamma_*(calib_planeX_*e1+calib_planeY_*e2);       
    }
 
   h1_esDeltaE->Fill(deltaE);
 
    nEvt_++;
 
 }

void PreshowerAnalyzer::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 72 of file PreshowerAnalyzer.cc.

References EmaxDE_, EmaxSC_, EminDE_, EminSC_, h1_esDeltaE, h1_esE_x, h1_esE_y, h1_esEta_x, h1_esEta_y, h1_esNhits_x, h1_esNhits_y, h1_esPhi_x, h1_esPhi_y, h1_islandEESCEnergy1, h1_islandEESCEnergy2, h1_nclu_x, h1_nclu_y, nBinDE_, nBinSC_, and rootFile_.

                                  {
 //========================================================================
 
   rootFile_->cd();
 
   h1_esE_x = new TH1F("esE_x"," ES cluster Energy in  X-plane",20, 0, 0.03);
   h1_esE_y = new TH1F("esE_y"," ES cluster Energy in  Y-plane",20, 0, 0.03);
   h1_esEta_x = new TH1F("esEta_x"," ES cluster Eta in X-plane",12, 1.5, 2.7);
   h1_esEta_y = new TH1F("esEta_y"," ES cluster Eta in Y-plane",12, 1.5, 2.7);
   h1_esPhi_x = new TH1F("esPhi_x"," ES cluster Phi in X-plane",20, 0, 6.28);
   h1_esPhi_y = new TH1F("esPhi_y"," ES cluster Phi in Y-plane",20, 0, 6.28);
   h1_esNhits_x = new TH1F("esNhits_x"," ES cluster Nhits in  X-plane",10, 0, 10);
   h1_esNhits_y = new TH1F("esNhits_y"," ES cluster Nhits in  Y-plane",10, 0, 10);
   h1_esDeltaE = new TH1F("esDeltaE"," DeltaE", nBinDE_, EminDE_, EmaxDE_); 
   h1_nclu_x = new TH1F("esNclu_x"," number of ES clusters (for one SC) in X-plane",20, 0, 80);
   h1_nclu_y = new TH1F("esNclu_y"," number of ES clusters (for one SC) in Y-plane",20, 0, 80);
 
   h1_islandEESCEnergy1 = new TH1F("islandEESCEnergy1","Energy of super clusters with island algo - endcap1",nBinSC_,EminSC_,EmaxSC_);
   h1_islandEESCEnergy2 = new TH1F("islandEESCEnergy2","Energy of super clusters with island algo - endcap2",nBinSC_,EminSC_,EmaxSC_);
 }

void PreshowerAnalyzer::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 178 of file PreshowerAnalyzer.cc.

References h1_esDeltaE, h1_esE_x, h1_esE_y, h1_esEta_x, h1_esEta_y, h1_esNhits_x, h1_esNhits_y, h1_esPhi_x, h1_esPhi_y, h1_islandEESCEnergy1, h1_islandEESCEnergy2, h1_nclu_x, h1_nclu_y, and rootFile_.

                                {
 //========================================================================
 
    rootFile_->cd();
 
    h1_esE_x->Write();     
    h1_esE_y->Write();
    h1_esEta_x->Write();
    h1_esEta_y->Write();
    h1_esPhi_x->Write();
    h1_esPhi_y->Write();
    h1_esNhits_x->Write();
    h1_esNhits_y->Write();          
    h1_esDeltaE->Write();
    h1_nclu_x->Write();
    h1_nclu_y->Write();
 
    h1_islandEESCEnergy1->Write();
    h1_islandEESCEnergy2->Write();
 
    rootFile_->Close();
 }