CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
ElseMETProducer.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: METProducers
4 // Class: ElseMETProducer
5 //
6 //
7 
8 //____________________________________________________________________________||
9 #include "RecoMET/METProducers/interface/ElseMETProducer.h"
10 
11 #include "FWCore/ParameterSet/interface/ParameterSet.h"
12 #include "FWCore/Framework/interface/Event.h"
13 #include "FWCore/Framework/interface/EventSetup.h"
14 #include "FWCore/Framework/interface/ConsumesCollector.h"
15 
16 #include "FWCore/Utilities/interface/InputTag.h"
17 
18 #include "DataFormats/Math/interface/LorentzVector.h"
19 #include "DataFormats/Math/interface/Point3D.h"
20 #include "DataFormats/Common/interface/Handle.h"
21 #include "DataFormats/METReco/interface/METFwd.h"
22 #include "DataFormats/METReco/interface/MET.h"
23 #include "DataFormats/METReco/interface/CommonMETData.h"
24 
25 #include "RecoMET/METAlgorithms/interface/METAlgo.h"
26 
27 #include <string.h>
28 
29 //____________________________________________________________________________||
30 namespace cms
31 {
32 
33 //____________________________________________________________________________||
34  ElseMETProducer::ElseMETProducer(const edm::ParameterSet& iConfig)
35  : inputToken_(consumes<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("src")))
36  , globalThreshold_(iConfig.getParameter<double>("globalThreshold"))
37  {
38  std::string alias = iConfig.exists("alias") ? iConfig.getParameter<std::string>("alias") : "";
39 
40  produces<reco::METCollection>().setBranchAlias(alias);
41  }
42 
43 
44 //____________________________________________________________________________||
45  void ElseMETProducer::produce(edm::Event& event, const edm::EventSetup& setup)
46  {
47  edm::Handle<edm::View<reco::Candidate> > input;
48  event.getByToken(inputToken_, input);
49 
50  METAlgo algo;
51  CommonMETData commonMETdata = algo.run(*input.product(), globalThreshold_);
52 
53  math::XYZTLorentzVector p4(commonMETdata.mex, commonMETdata.mey, 0.0, commonMETdata.met);
54  math::XYZPoint vtx(0,0,0);
55  reco::MET met(commonMETdata.sumet, p4, vtx);
56  std::auto_ptr<reco::METCollection> metcoll;
57  metcoll.reset(new reco::METCollection);
58  metcoll->push_back(met);
59  event.put(metcoll);
60  }
61 
62 //____________________________________________________________________________||
63  DEFINE_FWK_MODULE(ElseMETProducer);
64 }
65 
66 //____________________________________________________________________________||
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
cms::ElseMETProducer::ElseMETProducer
ElseMETProducer(const edm::ParameterSet &)
Definition: ElseMETProducer.cc:34
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:758
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
reco::METCollection
std::vector< reco::MET > METCollection
collection of MET objects
Definition: METCollection.h:23
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
input
static std::string const input
Definition: EdmProvDump.cc:44
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
METAlgo::run
CommonMETData run(const edm::View< reco::Candidate > &candidates, double globalThreshold=0.0)
Definition: METAlgo.cc:16
reco::MET
Definition: MET.h:42
p4
double p4[4]
Definition: TauolaWrapper.h:92
cms::ElseMETProducer::produce
virtual void produce(edm::Event &, const edm::EventSetup &) override
Definition: ElseMETProducer.cc:45
CommonMETData
Structure containing data common to all types of MET.
Definition: CommonMETData.h:12
event
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
edm::Handle::product
T const * product() const
Definition: Handle.h:81
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
caloMETBenchmarkGeneric_cfi.met
string met
Definition: caloMETBenchmarkGeneric_cfi.py:3
cms::ElseMETProducer::inputToken_
edm::EDGetTokenT< edm::View< reco::Candidate > > inputToken_
Definition: ElseMETProducer.h:43
HcalObjRepresent::setup
void setup(std::vector< TH2F > &depth, std::string name, std::string units="")
Definition: HcalObjRepresent.h:628
Generated for CMSSW Reference Manual by doxygen 1.8.5