#include <HLTMhtProducer.h>

Inheritance diagram for HLTMhtProducer:

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

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

	~HLTMhtProducer ()

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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
double	etaJet_

edm::InputTag	inputJetTag_

double	minPtJet_

bool	usePt_

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

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

Author: Gheorghe Lungu

Definition at line 19 of file HLTMhtProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 30 of file HLTMhtProducer.cc.

References etaJet_, edm::ParameterSet::getParameter(), inputJetTag_, minPtJet_, and usePt_.

 {
   inputJetTag_ = iConfig.getParameter< edm::InputTag > ("inputJetTag");
   minPtJet_= iConfig.getParameter<double> ("minPtJet");
   etaJet_= iConfig.getParameter<double> ("etaJet");
   usePt_= iConfig.getParameter<bool>("usePt");
 
   //register your products
   produces<reco::METCollection>();
 }

HLTMhtProducer::~HLTMhtProducer ( )

Definition at line 41 of file HLTMhtProducer.cc.

41 {}

Member Function Documentation

void HLTMhtProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 43 of file HLTMhtProducer.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("inputJetTag",edm::InputTag("hltMCJetCorJetIcone5HF07"));
   desc.add<double>("minPtJet",20.0);
   desc.add<double>("etaJet",9999.0);
   desc.add<bool>("usePt",true);
   descriptions.add("hltMhtProducer",desc);
 }

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

virtual

Implements edm::EDProducer.

Definition at line 54 of file HLTMhtProducer.cc.

References abs, funct::cos(), etaJet_, edm::Event::getByLabel(), inputJetTag_, metsig::jet, fwrapper::jets, minPtJet_, edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, query::result, funct::sin(), mathSSE::sqrt(), and usePt_.

 {
   using namespace std;
   using namespace edm;
   using namespace reco;
 
   auto_ptr<reco::METCollection> result (new reco::METCollection); 
 
   math::XYZPoint vtx(0,0,0);
   
   //Handle<CaloJetCollection> recocalojets;
   edm::Handle<edm::View<reco::Jet> > jets;
   iEvent.getByLabel(inputJetTag_,jets);
 
   // look at all candidates,  check cuts and add to result object
   double mhtx=0., mhty=0., mht;
   double jetVar;
   
   if(jets->size() > 0){
     // events with at least one jet
     //for (CaloJetCollection::const_iterator jet = jets->begin(); jet != jets->end(); jet++) {
     for(edm::View<reco::Jet>::const_iterator jet = jets->begin(); jet != jets->end(); jet++ ) {
       jetVar = jet->pt();
       if (!usePt_) jetVar = jet->et();
 
       //---get MHT
       if (jetVar > minPtJet_ && std::abs(jet->eta()) < etaJet_) {
     mhtx -= jetVar*cos(jet->phi());
     mhty -= jetVar*sin(jet->phi());
       }
     }
     mht = sqrt(mhtx*mhtx + mhty*mhty);
 
     math::XYZTLorentzVector mhtVec(mhtx,mhty,0,mht);
     reco::MET mhtobj(mhtVec,vtx);
     result->push_back( mhtobj );
     
   } // events with at least one jet
   
     
   // put object into the Event
   iEvent.put(result);
 
 }