CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Types | Private Attributes
L1ExtraMEtMixerPlugin Class Reference

#include <L1ExtraMEtMixerPlugin.h>

Inheritance diagram for L1ExtraMEtMixerPlugin:
L1ExtraMixerPluginBase edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 L1ExtraMEtMixerPlugin (const edm::ParameterSet &)
 
virtual void produce (edm::Event &, const edm::EventSetup &)
 
virtual void registerProducts (edm::EDProducer &)
 
 ~L1ExtraMEtMixerPlugin ()
 
- Public Member Functions inherited from L1ExtraMixerPluginBase
 L1ExtraMixerPluginBase (const edm::ParameterSet &)
 
 ~L1ExtraMixerPluginBase ()
 
- 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
 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
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Types

typedef std::map< uint32_t, float > detIdToFloatMap
 

Private Attributes

double sfAbsEta12to17_
 
double sfAbsEtaGt17_
 
double sfAbsEtaLt12_
 
edm::InputTag srcDistanceMapMuMinus_
 
edm::InputTag srcDistanceMapMuPlus_
 
edm::InputTag srcMuons_
 

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)
 
- Protected Attributes inherited from L1ExtraMixerPluginBase
std::string instanceLabel_
 
edm::InputTag src1_
 
edm::InputTag src2_
 

Detailed Description

Mix L1Extra MET objects, i.e. add MET reconstructed in original Zmumu event and negative vectorial momentum sum of embedded tau decay products

Author
Christian Veelken, LLR
Version
Revision:
1.5
Id:
L1ExtraMEtMixerPlugin.h,v 1.5 2013/01/31 09:07:18 veelken Exp

Definition at line 27 of file L1ExtraMEtMixerPlugin.h.

Member Typedef Documentation

typedef std::map<uint32_t, float> L1ExtraMEtMixerPlugin::detIdToFloatMap
private

Definition at line 37 of file L1ExtraMEtMixerPlugin.h.

Constructor & Destructor Documentation

L1ExtraMEtMixerPlugin::L1ExtraMEtMixerPlugin ( const edm::ParameterSet &  cfg)
explicit

Definition at line 16 of file L1ExtraMEtMixerPlugin.cc.

17  : L1ExtraMixerPluginBase(cfg),
18  srcMuons_(cfg.getParameter<edm::InputTag>("srcMuons")),
19  srcDistanceMapMuPlus_(cfg.getParameter<edm::InputTag>("distanceMapMuPlus")),
20  srcDistanceMapMuMinus_(cfg.getParameter<edm::InputTag>("distanceMapMuMinus")),
21  sfAbsEtaLt12_(cfg.getParameter<double>("H_Calo_AbsEtaLt12")),
22  sfAbsEta12to17_(cfg.getParameter<double>("H_Calo_AbsEta12to17")),
23  sfAbsEtaGt17_(cfg.getParameter<double>("H_Calo_AbsEtaGt17"))
24 {}
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
L1ExtraMixerPluginBase::L1ExtraMixerPluginBase
L1ExtraMixerPluginBase(const edm::ParameterSet &)
Definition: L1ExtraMixerPluginBase.cc:3
L1ExtraMEtMixerPlugin::srcDistanceMapMuMinus_
edm::InputTag srcDistanceMapMuMinus_
Definition: L1ExtraMEtMixerPlugin.h:41
L1ExtraMEtMixerPlugin::srcDistanceMapMuPlus_
edm::InputTag srcDistanceMapMuPlus_
Definition: L1ExtraMEtMixerPlugin.h:40
L1ExtraMEtMixerPlugin::~L1ExtraMEtMixerPlugin ( )
inline

Definition at line 31 of file L1ExtraMEtMixerPlugin.h.

31 {}

Member Function Documentation

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

Implements L1ExtraMixerPluginBase.

Definition at line 31 of file L1ExtraMEtMixerPlugin.cc.

References funct::cos(), DENSITY_BRASS, DENSITY_PBWO4, DetId::Ecal, l1extra::L1EtMissParticle::etTotal(), edm::hlt::Exception, edm::Event::getByLabel(), getDeDxForPbWO4(), DetId::Hcal, HcalOuter, reco::if(), L1ExtraMixerPluginBase::instanceLabel_, metsig::muon, patZpeak::muons, edm::Event::put(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), sfAbsEta12to17_, sfAbsEtaGt17_, sfAbsEtaLt12_, funct::sin(), python.multivaluedict::sort(), L1ExtraMixerPluginBase::src1_, L1ExtraMixerPluginBase::src2_, srcDistanceMapMuMinus_, srcDistanceMapMuPlus_, srcMuons_, and l1extra::L1EtMissParticle::type().

32 {
33  //std::cout << "<L1ExtraMEtMixerPlugin::produce>:" << std::endl;
34  //std::cout << " src1 = " << src1_ << std::endl;
35  //std::cout << " src2 = " << src2_ << std::endl;
36  //std::cout << " instanceLabel = " << instanceLabel_ << std::endl;
37 
38  edm::Handle<l1extra::L1EtMissParticleCollection> met1;
39  evt.getByLabel(src1_, met1);
40 
41  edm::Handle<l1extra::L1EtMissParticleCollection> met2;
42  evt.getByLabel(src2_, met2);
43 
44  edm::Handle<reco::CandidateCollection> muons;
45  evt.getByLabel(srcMuons_, muons);
46 
47  edm::Handle<detIdToFloatMap> distanceMapMuPlus;
48  evt.getByLabel(srcDistanceMapMuPlus_, distanceMapMuPlus);
49  edm::Handle<detIdToFloatMap> distanceMapMuMinus;
50  evt.getByLabel(srcDistanceMapMuMinus_, distanceMapMuMinus);
51 
52  std::auto_ptr<l1extra::L1EtMissParticleCollection> metSum(new l1extra::L1EtMissParticleCollection());
53 
54  // CV: entries in MET collections refer to different bunch-crossings.
55  // The number of entries in the two MET collections is not necessarily the same
56  // --> match entries by bunch-crossing number
57  std::set<int> bx_set;
58  for ( l1extra::L1EtMissParticleCollection::const_iterator met1_i = met1->begin();
59  met1_i != met1->end(); ++met1_i ) {
60  bx_set.insert(met1_i->bx());
61  }
62  for ( l1extra::L1EtMissParticleCollection::const_iterator met2_i = met2->begin();
63  met2_i != met2->end(); ++met2_i ) {
64  bx_set.insert(met2_i->bx());
65  }
66 
67  std::vector<int> bx_vector;
68  for ( std::set<int>::const_iterator bx = bx_set.begin();
69  bx != bx_set.end(); ++bx ) {
70  bx_vector.push_back(*bx);
71  }
72  std::sort(bx_vector.begin(), bx_vector.end());
73 
74  for ( std::vector<int>::const_iterator bx = bx_vector.begin();
75  bx != bx_vector.end(); ++bx ) {
76  bool errorFlag = false;
77 
78  const l1extra::L1EtMissParticle* met1_bx = 0;
79  for ( l1extra::L1EtMissParticleCollection::const_iterator met1_i = met1->begin();
80  met1_i != met1->end(); ++met1_i ) {
81  if ( met1_i->bx() == (*bx) ) {
82  if ( met1_bx ) errorFlag = true;
83  met1_bx = &(*met1_i);
84  }
85  }
86 
87  const l1extra::L1EtMissParticle* met2_bx = 0;
88  for ( l1extra::L1EtMissParticleCollection::const_iterator met2_i = met2->begin();
89  met2_i != met2->end(); ++met2_i ) {
90  if ( met2_i->bx() == (*bx) ) {
91  if ( met2_bx ) errorFlag = true;
92  met2_bx = &(*met2_i);
93  }
94  }
95 
96  if ( errorFlag )
97  throw cms::Exception("L1ExtraMEtMixer::produce")
98  << " Failed to find unique match of MET objects for BX = " << (*bx) << " !!\n";
99  assert(met1_bx || met2_bx);
100 
101  // CV: check that both MET objects are of the same type
102  if ( met1_bx && met2_bx && met1_bx->type() != met2_bx->type() )
103  throw cms::Exception("L1ExtraMEtMixer::produce")
104  << " Mismatch in type between MET objects stored in collections 'src1' and 'src2' !!\n";
105 
106  double metSumPx = 0.;
107  double metSumPy = 0.;
108  double metSumEt = 0.;
109  int type = -1;
110  if ( met1_bx ) {
111  //std::cout << "met1 (BX = " << (*bx) << "): Px = " << met1_bx->px() << ", Py = " << met1_bx->py()
112  // << " (Et = " << met1_bx->etMiss() << ", Pt = " << met1_bx->pt() << ", phi = " << met1_bx->phi() << ", sumEt = " << met1_bx->etTotal() << ")" << std::endl;
113  metSumPx += met1_bx->px();
114  metSumPy += met1_bx->py();
115  metSumEt += met1_bx->etTotal();
116  type = met1_bx->type();
117  }
118  if ( met2_bx ) {
119  //std::cout << "met2 (BX = " << (*bx) << "): Px = " << met2_bx->px() << ", Py = " << met2_bx->py()
120  // << " (Et = " << met2_bx->etMiss() << ", Pt = " << met2_bx->pt() << ", phi = " << met2_bx->phi() << ", sumEt = " << met2_bx->etTotal() << ")" << std::endl;
121  metSumPx += met2_bx->px();
122  metSumPy += met2_bx->py();
123  metSumEt += met2_bx->etTotal();
124  type = met2_bx->type();
125  }
126 
127  // Subtract contribution of muons that were replaced by simulated taus
128  if ( (*bx) == 0 ) {
129  for ( reco::CandidateCollection::const_iterator muon = muons->begin();
130  muon != muons->end(); ++muon ) {
131  //std::cout << "muon: Pt = " << muon->pt() << ", eta = " << muon->eta() << ", phi = " << muon->phi()
132  // << " (Px = " << muon->px() << ", Py = " << muon->py() << ")" << std::endl;
133 
134  double distance_ecal = 0.;
135  double distance_hcal = 0.;
136  const detIdToFloatMap& distanceMap = ( muon->charge() > 0 ) ?
137  (*distanceMapMuPlus) : (*distanceMapMuMinus);
138  for ( detIdToFloatMap::const_iterator dist_iter = distanceMap.begin();
139  dist_iter != distanceMap.end(); ++dist_iter ) {
140  const DetId& id = dist_iter->first;
141  if ( id.det() == DetId::Ecal ) distance_ecal += dist_iter->second;
142  if ( id.det() == DetId::Hcal && id.subdetId() != HcalOuter) distance_hcal += dist_iter->second;
143  }
144  //std::cout << "distance_ecal = " << distance_ecal << std::endl;
145  //std::cout << "distance_hcal = " << distance_hcal << std::endl;
146 
147  // This is what we would expect from theory:
148  double dedx = getDeDxForPbWO4(muon->p());
149  //std::cout << "dedx = " << dedx << std::endl;
150  double muonEnLoss_theory = dedx*(distance_ecal*DENSITY_PBWO4 + distance_hcal*DENSITY_BRASS);
151 
152  // We correct the muon energy loss computed from theory
153  // by the ratio of average reconstructed L1Extra MET (projected in direction of the muon)
154  // in single muon gun Monte Carlo events to the theory expecation.
155  // The ratio depends on eta:
156  double muonEnLoss = muonEnLoss_theory;
157  //std::cout << "muonEnLoss (before eta-correction) = " << muonEnLoss << std::endl;
158  if ( fabs(muon->eta()) < 1.2 ) muonEnLoss *= sfAbsEtaLt12_;
159  else if ( fabs(muon->eta()) < 1.7 ) muonEnLoss *= sfAbsEta12to17_;
160  else muonEnLoss *= sfAbsEtaGt17_;
161  //std::cout << "muonEnLoss (after eta-correction) = " << muonEnLoss << std::endl;
162 
163  double muonEnLoss_transverse = muonEnLoss*sin(muon->theta());
164  if ( muonEnLoss_transverse > metSumEt ) muonEnLoss_transverse = metSumEt;
165  //std::cout << "muonEnLoss_transverse = " << muonEnLoss_transverse << std::endl;
166 
167  // The direction is given by the muon direction
168  double muonMetCorrPx = muonEnLoss_transverse*cos(muon->phi());
169  double muonMetCorrPy = muonEnLoss_transverse*sin(muon->phi());
170  double muonMetCorrSumEt = muonEnLoss_transverse;
171  //std::cout << "muonMetCorr: Px = " << muonMetCorrPx << ", Py = " << muonMetCorrPy << " (sumEt = " << muonMetCorrSumEt << ")" << std::endl;
172  metSumPx += muonMetCorrPx;
173  metSumPy += muonMetCorrPy;
174  metSumEt -= muonMetCorrSumEt;
175  }
176  }
177 
178  const double metSumPt = TMath::Sqrt(metSumPx*metSumPx + metSumPy*metSumPy);
179  // CV: setting edm::Refs to L1Gct objects not implemented yet
180  l1extra::L1EtMissParticle metSum_bx(
181  reco::Candidate::LorentzVector(metSumPx, metSumPy, 0., metSumPt),
182  (l1extra::L1EtMissParticle::EtMissType)type,
183  metSumEt,
184  edm::Ref<L1GctEtMissCollection>(),
185  edm::Ref<L1GctEtTotalCollection>(),
186  edm::Ref<L1GctHtMissCollection>(),
187  edm::Ref<L1GctEtHadCollection>(),
188  (*bx));
189  //std::cout << "metSum (BX = " << (*bx) << "): Px = " << metSum_bx.px() << ", Py = " << metSum_bx.py()
190  // << " (Et = " << metSum_bx.etMiss() << ", Pt = " << metSum_bx.pt() << ", phi = " << metSum_bx.phi() << ")" << std::endl;
191  metSum->push_back(metSum_bx);
192  }
193 
194  evt.put(metSum, instanceLabel_);
195 }
DENSITY_PBWO4
const double DENSITY_PBWO4
Definition: embeddingAuxFunctions.h:54
type
type
Definition: HCALResponse.h:21
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
l1extra::L1EtMissParticle::type
EtMissType type() const
Definition: L1EtMissParticle.h:68
reco::LeafCandidate::py
virtual double py() const GCC11_FINAL
y coordinate of momentum vector
Definition: LeafCandidate.h:149
DENSITY_BRASS
const double DENSITY_BRASS
Definition: embeddingAuxFunctions.h:55
L1ExtraMEtMixerPlugin::detIdToFloatMap
std::map< uint32_t, float > detIdToFloatMap
Definition: L1ExtraMEtMixerPlugin.h:37
edm::Event::put
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:116
getDeDxForPbWO4
double getDeDxForPbWO4(double)
Definition: embeddingAuxFunctions.cc:320
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
reco::LeafCandidate::px
virtual double px() const GCC11_FINAL
x coordinate of momentum vector
Definition: LeafCandidate.h:147
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:390
DetId
Definition: DetId.h:18
reco::Candidate::LorentzVector
math::XYZTLorentzVector LorentzVector
Lorentz vector.
Definition: Candidate.h:41
patZpeak.muons
tuple muons
Definition: patZpeak.py:38
reco::if
if(dp >Float(M_PI)) dp-
l1extra::L1EtMissParticle::etTotal
const double & etTotal() const
Definition: L1EtMissParticle.h:75
L1ExtraMEtMixerPlugin::srcDistanceMapMuMinus_
edm::InputTag srcDistanceMapMuMinus_
Definition: L1ExtraMEtMixerPlugin.h:41
l1extra::L1EtMissParticleCollection
std::vector< L1EtMissParticle > L1EtMissParticleCollection
Definition: L1EtMissParticleFwd.h:33
L1ExtraMEtMixerPlugin::srcDistanceMapMuPlus_
edm::InputTag srcDistanceMapMuPlus_
Definition: L1ExtraMEtMixerPlugin.h:40
void L1ExtraMEtMixerPlugin::registerProducts ( edm::EDProducer &  producer)
virtual

Implements L1ExtraMixerPluginBase.

Definition at line 26 of file L1ExtraMEtMixerPlugin.cc.

References L1ExtraMixerPluginBase::instanceLabel_.

27 {
28  producer.produces<l1extra::L1EtMissParticleCollection>(instanceLabel_);
29 }
l1extra::L1EtMissParticleCollection
std::vector< L1EtMissParticle > L1EtMissParticleCollection
Definition: L1EtMissParticleFwd.h:33

Member Data Documentation

double L1ExtraMEtMixerPlugin::sfAbsEta12to17_
private

Definition at line 44 of file L1ExtraMEtMixerPlugin.h.

Referenced by produce().

double L1ExtraMEtMixerPlugin::sfAbsEtaGt17_
private

Definition at line 45 of file L1ExtraMEtMixerPlugin.h.

Referenced by produce().

double L1ExtraMEtMixerPlugin::sfAbsEtaLt12_
private

Definition at line 43 of file L1ExtraMEtMixerPlugin.h.

Referenced by produce().

edm::InputTag L1ExtraMEtMixerPlugin::srcDistanceMapMuMinus_
private

Definition at line 41 of file L1ExtraMEtMixerPlugin.h.

Referenced by produce().

edm::InputTag L1ExtraMEtMixerPlugin::srcDistanceMapMuPlus_
private

Definition at line 40 of file L1ExtraMEtMixerPlugin.h.

Referenced by produce().

edm::InputTag L1ExtraMEtMixerPlugin::srcMuons_
private

Definition at line 39 of file L1ExtraMEtMixerPlugin.h.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.5