#include <RecoBTag/JetCrystalsAssociator/src/JetCrystalsAssociator.cc>

Inheritance diagram for JetCrystalsAssociator:

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

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

	~JetCrystalsAssociator ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
std::unique_ptr< JetCrystalsAssociationCollection >	associate (const edm::Handle< CaloJetCollection > &jets, const edm::OrphanHandle< EMLorentzVectorCollection > &myLorentzRecHits) const

Private Attributes
double	m_deltaRCut

edm::InputTag	m_EBRecHits

edm::InputTag	m_EERecHits

edm::InputTag	m_jetsSrc

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

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

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

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

Definition at line 67 of file JetCrystalsAssociator.cc.

Constructor & Destructor Documentation

JetCrystalsAssociator::JetCrystalsAssociator ( const edm::ParameterSet & iConfig )

explicit

Definition at line 89 of file JetCrystalsAssociator.cc.

References edm::ParameterSet::getParameter().

 {
   produces<reco::JetCrystalsAssociationCollection>();
   produces<reco::EMLorentzVectorCollection>();
 
   m_EBRecHits = iConfig.getParameter<edm::InputTag>("EBRecHits");
   m_EERecHits = iConfig.getParameter<edm::InputTag>("EERecHits");
   m_jetsSrc   = iConfig.getParameter<edm::InputTag>("jets");
   m_deltaRCut = iConfig.getParameter<double>("coneSize");
 }

JetCrystalsAssociator::~JetCrystalsAssociator ( )

Definition at line 101 of file JetCrystalsAssociator.cc.

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

Member Function Documentation

std::unique_ptr< JetCrystalsAssociationCollection > JetCrystalsAssociator::associate	(	const edm::Handle< CaloJetCollection > &	jets,
		const edm::OrphanHandle< EMLorentzVectorCollection > &	myLorentzRecHits
	)		const

private

Definition at line 209 of file JetCrystalsAssociator.cc.

References DEFINE_FWK_MODULE, delta, HiRegitMuonDetachedTripletStep_cff::DeltaR, p4, and lumiQTWidget::t.

 {
   // we know we will save an element per input jet
   auto outputCollection = std::make_unique<JetCrystalsAssociationCollection>(jets->size());
 
   //loop on jets and associate
   for (size_t j = 0; j < jets->size(); j++) {
     (*outputCollection)[j].first = edm::RefToBase<Jet>(CaloJetRef(jets, j));
     for (size_t t = 0; t < myLorentzRecHits->size(); t++) {
       double delta = ROOT::Math::VectorUtil::DeltaR((*jets)[j].p4().Vect(), (*myLorentzRecHits)[t]);
       if (delta < m_deltaRCut)
         (*outputCollection)[j].second.push_back( EMLorentzVectorRef(myLorentzRecHits, t) );
     }
   }  
   return outputCollection;
 }

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

overridevirtual

Definition at line 116 of file JetCrystalsAssociator.cc.

References CaloTower::constituent(), CaloTower::constituentsSize(), DetId::det(), EB, egHLT::errCodes::EBRecHits, DetId::Ecal, EcalBarrel, EcalEndcap, EE, egHLT::errCodes::EERecHits, edm::SortedCollection< T, SORT >::end(), stringResolutionProvider_cfi::et, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), edm::SortedCollection< T, SORT >::find(), geometry, edm::EventSetup::get(), edm::Event::getByLabel(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), fwrapper::jets, M_PI, eostools::move(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::phi(), edm::Event::put(), funct::sin(), DetId::subdetId(), lumiQTWidget::t, PV3DBase< T, PVType, FrameType >::theta(), and theta().

 {
   using namespace edm;
   using namespace reco;
   using namespace std;
 
   // geometry initialization
   ESHandle<CaloGeometry> geometry;
   iSetup.get<CaloGeometryRecord>().get(geometry);
   
   const CaloSubdetectorGeometry* EB = geometry->getSubdetectorGeometry(DetId::Ecal,EcalBarrel);
    const CaloSubdetectorGeometry* EE = geometry->getSubdetectorGeometry(DetId::Ecal,EcalEndcap);
    // end 
    
    Handle<CaloJetCollection> jets;
    iEvent.getByLabel(m_jetsSrc, jets);
    
    // get calo towers collection
    //   Handle<CaloTowerCollection> caloTowers; 
    //   iEvent.getByLabel(m_towersSrc, caloTowers);
 
    // calculation of ECAL isolation
    Handle<EBRecHitCollection> EBRecHits;
    Handle<EERecHitCollection> EERecHits;
    iEvent.getByLabel( m_EBRecHits, EBRecHits );
    iEvent.getByLabel( m_EERecHits, EERecHits );
    
    auto jetRecHits = std::make_unique<EMLorentzVectorCollection>();
    //loop on jets and associate
    for (size_t t = 0; t < jets->size(); t++)
     {
       const std::vector<CaloTowerPtr>  myTowers=(*jets)[t].getCaloConstituents();
       //      cout <<"Jet id "<<t<<endl;
       //      cout <<"Tower size "<<myTowers.size()<<endl;
       for (unsigned int iTower = 0; iTower < myTowers.size(); iTower++)
     {
       CaloTowerPtr theTower = myTowers[iTower];
       size_t numRecHits = theTower->constituentsSize();
     // access CaloRecHits
     for (size_t j = 0; j < numRecHits; j++) {
       DetId RecHitDetID=theTower->constituent(j);
       DetId::Detector DetNum=RecHitDetID.det();
       if( DetNum == DetId::Ecal ){
         int EcalNum =  RecHitDetID.subdetId();
         if( EcalNum == 1 ){
           EBDetId EcalID = RecHitDetID;
           EBRecHitCollection::const_iterator theRecHit=EBRecHits->find(EcalID);
           if(theRecHit != EBRecHits->end()){
         DetId id = theRecHit->detid();
         const CaloCellGeometry* this_cell = EB->getGeometry(id);
         if (this_cell) {
           GlobalPoint posi = this_cell->getPosition();
           double energy = theRecHit->energy();
           double eta = posi.eta();
           double phi = posi.phi();
           double theta = posi.theta();
           if(theta > M_PI) theta = 2 * M_PI- theta;
           double et = energy * sin(theta);
           // cout <<"Et "<<et<<endl;
           EMLorentzVector p(et, eta, phi, energy);
           jetRecHits->push_back(p);
         }
           }
         } else if ( EcalNum == 2 ) {
           EEDetId EcalID = RecHitDetID;
           EERecHitCollection::const_iterator theRecHit=EERecHits->find(EcalID);     
           if(theRecHit != EBRecHits->end()){
         DetId id = theRecHit->detid();
         const CaloCellGeometry* this_cell = EE->getGeometry(id);
         if (this_cell) {
           GlobalPoint posi = this_cell->getPosition();
           double energy = theRecHit->energy();
           double eta = posi.eta();
           double phi = posi.phi();
           double theta = posi.theta();
           if (theta > M_PI) theta = 2 * M_PI - theta;
           double et = energy * sin(theta);
           // cout <<"Et "<<et<<endl;
           EMLorentzVector p(et, eta, phi, energy);
           jetRecHits->push_back(p);
         }
           }
         }
       }
     }
       }
     }
 
   edm::OrphanHandle <reco::EMLorentzVectorCollection> myRecHits = iEvent.put(std::move(jetRecHits));
 
   iEvent.put(associate(jets,myRecHits));
 }

Member Data Documentation

double JetCrystalsAssociator::m_deltaRCut

private

Definition at line 86 of file JetCrystalsAssociator.cc.

edm::InputTag JetCrystalsAssociator::m_EBRecHits

private

Definition at line 83 of file JetCrystalsAssociator.cc.

edm::InputTag JetCrystalsAssociator::m_EERecHits

private

Definition at line 84 of file JetCrystalsAssociator.cc.

edm::InputTag JetCrystalsAssociator::m_jetsSrc

private

Definition at line 82 of file JetCrystalsAssociator.cc.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation