#include <AlCaDiJetsProducer.h>

Inheritance diagram for cms::AlCaDiJetsProducer:

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

virtual void	beginJob ()

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

	~AlCaDiJetsProducer ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
bool	allowMissingInputs_

std::vector< edm::InputTag >	ecalLabels_

edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT < HFRecHitCollection >	tok_hf_

edm::EDGetTokenT < HORecHitCollection >	tok_ho_

edm::EDGetTokenT < reco::CaloJetCollection >	tok_jets_

std::vector< edm::EDGetTokenT < EcalRecHitCollection > >	toks_ecal_

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

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)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

Definition at line 42 of file AlCaDiJetsProducer.h.

Constructor & Destructor Documentation

cms::AlCaDiJetsProducer::AlCaDiJetsProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 14 of file AlCaDiJetsProducer.cc.

References edm::ParameterSet::getParameter(), and i.

 {
    tok_jets_ = consumes<CaloJetCollection>(iConfig.getParameter<edm::InputTag>("jetsInput"));
    ecalLabels_=iConfig.getParameter<std::vector<edm::InputTag> >("ecalInputs");
    tok_hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
    tok_ho_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInput"));
    tok_hf_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInput")); 
    allowMissingInputs_ = true;
 
   // fill ecal tokens from input labels
    const unsigned nLabels = ecalLabels_.size();
    for ( unsigned i=0; i != nLabels; i++ ) 
      toks_ecal_.push_back(consumes<EcalRecHitCollection>(ecalLabels_[i]));
 
 //register your products
    produces<CaloJetCollection>("DiJetsBackToBackCollection");
    produces<EcalRecHitCollection>("DiJetsEcalRecHitCollection");
    produces<HBHERecHitCollection>("DiJetsHBHERecHitCollection");
    produces<HORecHitCollection>("DiJetsHORecHitCollection");
    produces<HFRecHitCollection>("DiJetsHFRecHitCollection");
 
 }

cms::AlCaDiJetsProducer::~AlCaDiJetsProducer ( )

Definition at line 40 of file AlCaDiJetsProducer.cc.

 {
  
 }

Member Function Documentation

void cms::AlCaDiJetsProducer::beginJob ( void )

virtual

Reimplemented from edm::EDProducer.

Definition at line 36 of file AlCaDiJetsProducer.cc.

37 {

38 }

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

virtual

Implements edm::EDProducer.

Definition at line 49 of file AlCaDiJetsProducer.cc.

References PV3DBase< T, PVType, FrameType >::eta(), eta(), edm::EventSetup::get(), edm::Event::getByToken(), CaloGeometry::getPosition(), i, fwrapper::jets, phi, PV3DBase< T, PVType, FrameType >::phi(), reco::LeafCandidate::phi(), pi, edm::ESHandle< class >::product(), edm::Event::put(), query::result, and mathSSE::sqrt().

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

 {
 
   double pi = 4.*atan(1.);
 
    std::auto_ptr<CaloJetCollection> result (new CaloJetCollection); //Corrected jets
    std::auto_ptr<EcalRecHitCollection> miniDiJetsEcalRecHitCollection(new EcalRecHitCollection);
    std::auto_ptr<HBHERecHitCollection> miniDiJetsHBHERecHitCollection(new HBHERecHitCollection);
    std::auto_ptr<HORecHitCollection> miniDiJetsHORecHitCollection(new HORecHitCollection);
    std::auto_ptr<HFRecHitCollection> miniDiJetsHFRecHitCollection(new HFRecHitCollection);
 
    edm::ESHandle<CaloGeometry> pG;
    iSetup.get<CaloGeometryRecord>().get(pG);
    const CaloGeometry* geo = pG.product();
 
   // Jets Collections 
 
    vector<CaloJet> jetv; 
 
    CaloJet fJet1, fJet2, fJet3;
    edm::Handle<CaloJetCollection> jets;                
    iEvent.getByToken(tok_jets_, jets);    
    int iflag_select = 0; 
    if(jets->size()>1){
     fJet1 = (*jets)[0];
     fJet2 = (*jets)[1];
     double dphi = fabs(fJet1.phi() - fJet2.phi());  
     if(dphi > pi){dphi = 2*pi - dphi;}
     double degreedphi = dphi*180./pi; 
     if(fabs(degreedphi-180)<30.){iflag_select = 1;}
    }
    if(iflag_select == 1){
      result->push_back(fJet1);
      result->push_back(fJet2);
      jetv.push_back(fJet1); 
      jetv.push_back(fJet2); 
      if(jets->size()>2){
      fJet3 = (*jets)[2]; 
      result->push_back(fJet3);
      jetv.push_back(fJet3);
      }     
    } else {
      iEvent.put( result, "DiJetsBackToBackCollection");
      iEvent.put( miniDiJetsEcalRecHitCollection,"DiJetsEcalRecHitCollection");
      iEvent.put( miniDiJetsHBHERecHitCollection, "DiJetsHBHERecHitCollection");
      iEvent.put( miniDiJetsHORecHitCollection, "DiJetsHORecHitCollection");
      iEvent.put( miniDiJetsHFRecHitCollection, "DiJetsHFRecHitCollection");
      return;
    }  
   
   // Ecal Collections 
 
    std::vector<edm::EDGetTokenT<EcalRecHitCollection> >::const_iterator i;
    for (i=toks_ecal_.begin(); i!=toks_ecal_.end(); i++) {
    edm::Handle<EcalRecHitCollection> ec;
    iEvent.getByToken(*i,ec);
    for(EcalRecHitCollection::const_iterator ecItr = (*ec).begin();
                                                  ecItr != (*ec).end(); ++ecItr)
         {
 // EcalBarrel = 1, EcalEndcap = 2
           GlobalPoint pos = geo->getPosition(ecItr->detid());
           double phihit = pos.phi();
           double etahit = pos.eta();
           int iflag_select = 0;  
           for(unsigned int i=0; i<jetv.size(); i++){
             double deta = fabs(etahit - jetv[i].eta());
             double dphi = fabs(phihit - jetv[i].phi());
             if(dphi > pi) dphi = 2*pi - dphi;
             double dr = sqrt(deta*deta+dphi*dphi);
             if(dr < 1.4){iflag_select = 1;}
           }  
           if(iflag_select==1){miniDiJetsEcalRecHitCollection->push_back(*ecItr);}
        }
 
    }
 
   // HB & HE Collections 
    
    edm::Handle<HBHERecHitCollection> hbhe;
    iEvent.getByToken(tok_hbhe_,hbhe);
    for(HBHERecHitCollection::const_iterator hbheItr=hbhe->begin(); 
                                                   hbheItr!=hbhe->end(); hbheItr++)
         {
           GlobalPoint pos = geo->getPosition(hbheItr->detid());
           double phihit = pos.phi();
           double etahit = pos.eta();
           int iflag_select = 0;  
           for(unsigned int i=0; i<jetv.size(); i++){
             double deta = fabs(etahit - jetv[i].eta());
             double dphi = fabs(phihit - jetv[i].phi());
             if(dphi > pi) dphi = 2*pi - dphi;
             double dr = sqrt(deta*deta+dphi*dphi);
             if(dr < 1.4){iflag_select = 1;}
           }  
           if(iflag_select==1){miniDiJetsHBHERecHitCollection->push_back(*hbheItr);}
         }
   
 
 // HO Collections 
 
 
    edm::Handle<HORecHitCollection> ho;
    iEvent.getByToken(tok_ho_,ho);
    for(HORecHitCollection::const_iterator hoItr=ho->begin(); 
                                                 hoItr!=ho->end(); hoItr++)
      {
           GlobalPoint pos = geo->getPosition(hoItr->detid());
           double phihit = pos.phi();
           double etahit = pos.eta();
           int iflag_select = 0;  
           for(unsigned int i=0; i<jetv.size(); i++){
             double deta = fabs(etahit - jetv[i].eta());
             double dphi = fabs(phihit - jetv[i].phi());
             if(dphi > pi) dphi = 2*pi - dphi;
             double dr = sqrt(deta*deta+dphi*dphi);
             if(dr < 1.4){iflag_select = 1;}
           }  
           if(iflag_select==1){miniDiJetsHORecHitCollection->push_back(*hoItr);}
         }
   
   // HF Collection
 
  
    edm::Handle<HFRecHitCollection> hf;
    iEvent.getByToken(tok_hf_,hf);
    for(HFRecHitCollection::const_iterator hfItr=hf->begin(); 
                                                 hfItr!=hf->end(); hfItr++)
        {
           GlobalPoint pos = geo->getPosition(hfItr->detid());
           double phihit = pos.phi();
           double etahit = pos.eta();
           int iflag_select = 0;  
           for(unsigned int i=0; i<jetv.size(); i++){
             double deta = fabs(etahit - jetv[i].eta());
             double dphi = fabs(phihit - jetv[i].phi());
             if(dphi > pi) dphi = 2*pi - dphi;
             double dr = sqrt(deta*deta+dphi*dphi);
             if(dr < 1.4){iflag_select = 1;}
           }  
           if(iflag_select==1){miniDiJetsHFRecHitCollection->push_back(*hfItr);}
        }
  
 
   //Put selected information in the event
 
    iEvent.put( result, "DiJetsBackToBackCollection");
    iEvent.put( miniDiJetsEcalRecHitCollection,"DiJetsEcalRecHitCollection");
    iEvent.put( miniDiJetsHBHERecHitCollection, "DiJetsHBHERecHitCollection");
    iEvent.put( miniDiJetsHORecHitCollection, "DiJetsHORecHitCollection");
    iEvent.put( miniDiJetsHFRecHitCollection, "DiJetsHFRecHitCollection");
 }

Member Data Documentation

bool cms::AlCaDiJetsProducer::allowMissingInputs_

private

Definition at line 61 of file AlCaDiJetsProducer.h.

std::vector<edm::InputTag> cms::AlCaDiJetsProducer::ecalLabels_

private

Definition at line 52 of file AlCaDiJetsProducer.h.

edm::EDGetTokenT<HBHERecHitCollection> cms::AlCaDiJetsProducer::tok_hbhe_

private

Definition at line 55 of file AlCaDiJetsProducer.h.

edm::EDGetTokenT<HFRecHitCollection> cms::AlCaDiJetsProducer::tok_hf_

private

Definition at line 57 of file AlCaDiJetsProducer.h.

edm::EDGetTokenT<HORecHitCollection> cms::AlCaDiJetsProducer::tok_ho_

private

Definition at line 56 of file AlCaDiJetsProducer.h.

edm::EDGetTokenT<reco::CaloJetCollection> cms::AlCaDiJetsProducer::tok_jets_

private

Definition at line 54 of file AlCaDiJetsProducer.h.

std::vector<edm::EDGetTokenT<EcalRecHitCollection> > cms::AlCaDiJetsProducer::toks_ecal_

private

Definition at line 59 of file AlCaDiJetsProducer.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation