#include <ElectronSqPtTkIsolationProducer.h>

Inheritance diagram for ElectronSqPtTkIsolationProducer:

Public Member Functions
	ElectronSqPtTkIsolationProducer (const edm::ParameterSet &)

virtual void	produce (edm::Event &, const edm::EventSetup &)

	~ElectronSqPtTkIsolationProducer ()

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Attributes
bool	absolut_

edm::ParameterSet	conf_

edm::InputTag	electronProducer_

double	extRadius_

double	intRadius_

double	maxVtxDist_

double	ptMin_

edm::InputTag	trackProducer_

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

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

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

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 14 of file ElectronSqPtTkIsolationProducer.h.

Constructor & Destructor Documentation

ElectronSqPtTkIsolationProducer::ElectronSqPtTkIsolationProducer ( const edm::ParameterSet & config )

explicit

Definition at line 15 of file ElectronSqPtTkIsolationProducer.cc.

References absolut_, conf_, electronProducer_, extRadius_, edm::ParameterSet::getParameter(), intRadius_, maxVtxDist_, ptMin_, and trackProducer_.

                                                                                               : conf_(config)
 {
   // use configuration file to setup input/output collection names
   electronProducer_               = conf_.getParameter<edm::InputTag>("electronProducer");
   
   trackProducer_           = conf_.getParameter<edm::InputTag>("trackProducer");
 
   ptMin_                = conf_.getParameter<double>("ptMin");
   intRadius_            = conf_.getParameter<double>("intRadius");
   extRadius_            = conf_.getParameter<double>("extRadius");
   maxVtxDist_           = conf_.getParameter<double>("maxVtxDist");
 
   absolut_              = conf_.getParameter<bool>("absolut");
 
   //register your products
   produces < reco::CandViewDoubleAssociations>();
 
 }

ElectronSqPtTkIsolationProducer::~ElectronSqPtTkIsolationProducer ( )

Definition at line 34 of file ElectronSqPtTkIsolationProducer.cc.

34 {}

Member Function Documentation

void ElectronSqPtTkIsolationProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements edm::EDProducer.

Definition at line 38 of file ElectronSqPtTkIsolationProducer.cc.

References absolut_, electronProducer_, ExpressReco_HICollisions_FallBack::et, funct::exp(), extRadius_, edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), ElectronSqPtTkIsolation::getPtTracks(), i, intRadius_, maxVtxDist_, edm::Handle< T >::product(), ptMin_, edm::Event::put(), edm::AssociationVector< KeyRefProd, CVal, KeyRef, SizeType, KeyReferenceHelper >::setValue(), funct::sin(), ExpressReco_HICollisions_FallBack::trackCollection, trackProducer_, and testEve_cfg::tracks.

Referenced by python.JSONExport.JsonExport::export(), and python.HTMLExport.HTMLExport::export().

 {
   // Get the  filtered objects
   edm::Handle< reco::GsfElectronCollection> electronHandle;
   iEvent.getByLabel(electronProducer_,electronHandle);
   
   //get the tracks
   edm::Handle<reco::TrackCollection> tracks;
   iEvent.getByLabel(trackProducer_,tracks);
   const reco::TrackCollection* trackCollection = tracks.product();
   
   reco::CandViewDoubleAssociations* isoMap = new reco::CandViewDoubleAssociations(reco::GsfElectronRefProd( electronHandle ));
   
   ElectronSqPtTkIsolation myTkIsolation (extRadius_,intRadius_,ptMin_,maxVtxDist_,trackCollection) ;
   
   for(unsigned int i = 0 ; i < electronHandle->size(); ++i ){
     double isoValue = myTkIsolation.getPtTracks(&(electronHandle->at(i)));
     if(absolut_==true){
       isoMap->setValue(i,isoValue);
     }
     else{
       reco::SuperClusterRef sc = (electronHandle->at(i)).superCluster();
       double et = sc.get()->energy()*sin(2*atan(exp(-sc.get()->eta())));
       isoMap->setValue(i,isoValue/et);
     }
   }
   
   std::auto_ptr<reco::CandViewDoubleAssociations> isolMap(isoMap);
   iEvent.put(isolMap);
 }