#include <L1ExtraMEtMixerPlugin.h>

Inheritance diagram for L1ExtraMEtMixerPlugin:

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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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

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)

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.

   : L1ExtraMixerPluginBase(cfg),
     srcMuons_(cfg.getParameter<edm::InputTag>("srcMuons")),
     srcDistanceMapMuPlus_(cfg.getParameter<edm::InputTag>("distanceMapMuPlus")),
     srcDistanceMapMuMinus_(cfg.getParameter<edm::InputTag>("distanceMapMuMinus")),
     sfAbsEtaLt12_(cfg.getParameter<double>("H_Calo_AbsEtaLt12")),
     sfAbsEta12to17_(cfg.getParameter<double>("H_Calo_AbsEta12to17")),
     sfAbsEtaGt17_(cfg.getParameter<double>("H_Calo_AbsEtaGt17"))
 {}

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 assert(), funct::cos(), DENSITY_BRASS, DENSITY_PBWO4, DetId::Ecal, l1extra::L1EtMissParticle::etTotal(), edm::hlt::Exception, edm::Event::getByLabel(), getDeDxForPbWO4(), DetId::Hcal, HcalOuter, 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().

 {
   //std::cout << "<L1ExtraMEtMixerPlugin::produce>:" << std::endl;
   //std::cout << " src1 = " << src1_ << std::endl;
   //std::cout << " src2 = " << src2_ << std::endl;
   //std::cout << " instanceLabel = " << instanceLabel_ << std::endl;
 
   edm::Handle<l1extra::L1EtMissParticleCollection> met1;
   evt.getByLabel(src1_, met1);
   
   edm::Handle<l1extra::L1EtMissParticleCollection> met2;
   evt.getByLabel(src2_, met2);
 
   edm::Handle<reco::CandidateCollection> muons;
   evt.getByLabel(srcMuons_, muons);
 
   edm::Handle<detIdToFloatMap> distanceMapMuPlus;
   evt.getByLabel(srcDistanceMapMuPlus_, distanceMapMuPlus);
   edm::Handle<detIdToFloatMap> distanceMapMuMinus;
   evt.getByLabel(srcDistanceMapMuMinus_, distanceMapMuMinus);
 
   std::auto_ptr<l1extra::L1EtMissParticleCollection> metSum(new l1extra::L1EtMissParticleCollection());
 
   // CV: entries in MET collections refer to different bunch-crossings.
   //     The number of entries in the two MET collections is not necessarily the same
   //    --> match entries by bunch-crossing number
   std::set<int> bx_set;
   for ( l1extra::L1EtMissParticleCollection::const_iterator met1_i = met1->begin();
     met1_i != met1->end(); ++met1_i ) {
     bx_set.insert(met1_i->bx());
   }
   for ( l1extra::L1EtMissParticleCollection::const_iterator met2_i = met2->begin();
     met2_i != met2->end(); ++met2_i ) {
     bx_set.insert(met2_i->bx());
   }
 
   std::vector<int> bx_vector;
   for ( std::set<int>::const_iterator bx = bx_set.begin();
     bx != bx_set.end(); ++bx ) {
     bx_vector.push_back(*bx);
   }
   std::sort(bx_vector.begin(), bx_vector.end());
 
   for ( std::vector<int>::const_iterator bx = bx_vector.begin();
     bx != bx_vector.end(); ++bx ) {
     bool errorFlag = false;
 
     const l1extra::L1EtMissParticle* met1_bx = 0;
     for ( l1extra::L1EtMissParticleCollection::const_iterator met1_i = met1->begin();
       met1_i != met1->end(); ++met1_i ) {
       if ( met1_i->bx() == (*bx) ) {
     if ( met1_bx ) errorFlag = true;
     met1_bx = &(*met1_i);
       }
     }
 
     const l1extra::L1EtMissParticle* met2_bx = 0;
     for ( l1extra::L1EtMissParticleCollection::const_iterator met2_i = met2->begin();
       met2_i != met2->end(); ++met2_i ) {
       if ( met2_i->bx() == (*bx) ) {
     if ( met2_bx ) errorFlag = true;
     met2_bx = &(*met2_i);
       }
     }
 
     if ( errorFlag )
       throw cms::Exception("L1ExtraMEtMixer::produce")
     << " Failed to find unique match of MET objects for BX = " << (*bx) << " !!\n";
     assert(met1_bx || met2_bx);
     
      // CV: check that both MET objects are of the same type
     if ( met1_bx && met2_bx && met1_bx->type() != met2_bx->type() )
       throw cms::Exception("L1ExtraMEtMixer::produce")
     << " Mismatch in type between MET objects stored in collections 'src1' and 'src2' !!\n";
 
     double metSumPx = 0.;
     double metSumPy = 0.;
     double metSumEt = 0.;
     int type = -1;
     if ( met1_bx ) {
       //std::cout << "met1 (BX = " << (*bx) << "): Px = " << met1_bx->px() << ", Py = " << met1_bx->py() 
       //      << " (Et = " << met1_bx->etMiss() << ", Pt = " << met1_bx->pt() << ", phi = " << met1_bx->phi() << ", sumEt = " << met1_bx->etTotal() << ")" << std::endl;
       metSumPx += met1_bx->px();
       metSumPy += met1_bx->py();
       metSumEt += met1_bx->etTotal();
       type = met1_bx->type();
     }
     if ( met2_bx ) {
       //std::cout << "met2 (BX = " << (*bx) << "): Px = " << met2_bx->px() << ", Py = " << met2_bx->py() 
       //      << " (Et = " << met2_bx->etMiss() << ", Pt = " << met2_bx->pt() << ", phi = " << met2_bx->phi() << ", sumEt = " << met2_bx->etTotal() << ")" << std::endl;
       metSumPx += met2_bx->px();
       metSumPy += met2_bx->py();
       metSumEt += met2_bx->etTotal();
       type = met2_bx->type();
     }
 
     // Subtract contribution of muons that were replaced by simulated taus
     if ( (*bx) == 0 ) {
       for ( reco::CandidateCollection::const_iterator muon = muons->begin(); 
         muon != muons->end(); ++muon ) {
     //std::cout << "muon: Pt = " << muon->pt() << ", eta = " << muon->eta() << ", phi = " << muon->phi() 
     //      << " (Px = " << muon->px() << ", Py = " << muon->py() << ")" << std::endl;
 
     double distance_ecal = 0.;
     double distance_hcal = 0.;
     const detIdToFloatMap& distanceMap = ( muon->charge() > 0 ) ? 
       (*distanceMapMuPlus) : (*distanceMapMuMinus);
     for ( detIdToFloatMap::const_iterator dist_iter = distanceMap.begin(); 
           dist_iter != distanceMap.end(); ++dist_iter ) {
       const DetId& id = dist_iter->first;
       if ( id.det() == DetId::Ecal                              ) distance_ecal += dist_iter->second;
       if ( id.det() == DetId::Hcal && id.subdetId() != HcalOuter) distance_hcal += dist_iter->second;
     }
     //std::cout << "distance_ecal = " << distance_ecal << std::endl;
     //std::cout << "distance_hcal = " << distance_hcal << std::endl;
 
     // This is what we would expect from theory:
     double dedx = getDeDxForPbWO4(muon->p());
     //std::cout << "dedx = " << dedx << std::endl;
     double muonEnLoss_theory = dedx*(distance_ecal*DENSITY_PBWO4 + distance_hcal*DENSITY_BRASS);
 
     // We correct the muon energy loss computed from theory
     // by the ratio of average reconstructed L1Extra MET (projected in direction of the muon)
     // in single muon gun Monte Carlo events to the theory expecation.
     // The ratio depends on eta:
     double muonEnLoss = muonEnLoss_theory;
     //std::cout << "muonEnLoss (before eta-correction) = " << muonEnLoss << std::endl;
     if      ( fabs(muon->eta()) < 1.2 ) muonEnLoss *= sfAbsEtaLt12_;
     else if ( fabs(muon->eta()) < 1.7 ) muonEnLoss *= sfAbsEta12to17_;
     else                                muonEnLoss *= sfAbsEtaGt17_;
     //std::cout << "muonEnLoss (after eta-correction) = " << muonEnLoss << std::endl;
 
     double muonEnLoss_transverse = muonEnLoss*sin(muon->theta());
     if ( muonEnLoss_transverse > metSumEt ) muonEnLoss_transverse = metSumEt;
     //std::cout << "muonEnLoss_transverse = " << muonEnLoss_transverse << std::endl;
 
     // The direction is given by the muon direction
     double muonMetCorrPx    = muonEnLoss_transverse*cos(muon->phi());
     double muonMetCorrPy    = muonEnLoss_transverse*sin(muon->phi());
     double muonMetCorrSumEt = muonEnLoss_transverse;
     //std::cout << "muonMetCorr: Px = " << muonMetCorrPx << ", Py = " << muonMetCorrPy << " (sumEt = " << muonMetCorrSumEt << ")" << std::endl;
     metSumPx += muonMetCorrPx;
     metSumPy += muonMetCorrPy;
     metSumEt -= muonMetCorrSumEt;
       }
     }
 
     const double metSumPt = TMath::Sqrt(metSumPx*metSumPx + metSumPy*metSumPy);
     // CV: setting edm::Refs to L1Gct objects not implemented yet
     l1extra::L1EtMissParticle metSum_bx(
       reco::Candidate::LorentzVector(metSumPx, metSumPy, 0., metSumPt),
       (l1extra::L1EtMissParticle::EtMissType)type,
       metSumEt,
       edm::Ref<L1GctEtMissCollection>(),
       edm::Ref<L1GctEtTotalCollection>(),
       edm::Ref<L1GctHtMissCollection>(),
       edm::Ref<L1GctEtHadCollection>(),
       (*bx));
     //std::cout << "metSum (BX = " << (*bx) << "): Px = " << metSum_bx.px() << ", Py = " << metSum_bx.py() 
     //          << " (Et = " << metSum_bx.etMiss() << ", Pt = " << metSum_bx.pt() << ", phi = " << metSum_bx.phi() << ")" << std::endl;
     metSum->push_back(metSum_bx);                          
   }
   
   evt.put(metSum, instanceLabel_);
 }

void L1ExtraMEtMixerPlugin::registerProducts ( edm::EDProducer & producer )

virtual

Implements L1ExtraMixerPluginBase.

Definition at line 26 of file L1ExtraMEtMixerPlugin.cc.

References L1ExtraMixerPluginBase::instanceLabel_.

 {
   producer.produces<l1extra::L1EtMissParticleCollection>(instanceLabel_); 
 }

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

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation