#include <ZmumuEvtSelEffCorrWeightProducer.h>

Inheritance diagram for ZmumuEvtSelEffCorrWeightProducer:

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &)

	ZmumuEvtSelEffCorrWeightProducer (const edm::ParameterSet &)

	~ZmumuEvtSelEffCorrWeightProducer ()

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 Attributes
TH2 *	lutEffCorrEta_

TH2 *	lutEfficiencyPt_

double	maxWeight_

double	minWeight_

edm::InputTag	srcSelectedMuons_

int	verbosity_

TAxis *	xAxisEffCorrEta_

TAxis *	xAxisEfficiencyPt_

TAxis *	yAxisEffCorrEta_

TAxis *	yAxisEfficiencyPt_

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)

Detailed Description

Compute weight factor and uncertainty for correcting Embedded samples for efficiency with which Zmumu events used as input for Embedding production were selected

Authors: Christian Veelken

Version

Revision:: 1.1

Id:: ZmumuEvtSelEffCorrWeightProducer.h,v 1.1 2013/02/21 14:08:41 veelken Exp

Definition at line 32 of file ZmumuEvtSelEffCorrWeightProducer.h.

Constructor & Destructor Documentation

ZmumuEvtSelEffCorrWeightProducer::ZmumuEvtSelEffCorrWeightProducer ( const edm::ParameterSet & cfg )

explicit

Definition at line 17 of file ZmumuEvtSelEffCorrWeightProducer.cc.

References python.multivaluedict::append(), edm::hlt::Exception, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), analyzePatCleaning_cfg::inputFile, writeAntiElectronDiscrMVA_cfi::inputFileName, lutEffCorrEta_, lutEfficiencyPt_, maxWeight_, minWeight_, srcSelectedMuons_, AlCaHLTBitMon_QueryRunRegistry::string, edm::FileInPath::Unknown, verbosity_, xAxisEffCorrEta_, xAxisEfficiencyPt_, yAxisEffCorrEta_, and yAxisEfficiencyPt_.

   : lutEfficiencyPt_(0),
     lutEffCorrEta_(0)
 {
   srcSelectedMuons_ = cfg.getParameter<edm::InputTag>("selectedMuons");
 
   edm::FileInPath inputFileName = cfg.getParameter<edm::FileInPath>("inputFileName");
   if ( inputFileName.location() == edm::FileInPath::Unknown) 
     throw cms::Exception("MuonRadiationCorrWeightProducer") 
       << " Failed to find File = " << inputFileName << " !!\n";
   std::auto_ptr<TFile> inputFile(new TFile(inputFileName.fullPath().data()));
 
   std::string lutEfficiencyPtName = cfg.getParameter<std::string>("lutEfficiencyPt");
   TH2* lutEfficiencyPt = dynamic_cast<TH2*>(inputFile->Get(lutEfficiencyPtName.data()));
   if ( !lutEfficiencyPt ) 
     throw cms::Exception("MuonRadiationCorrWeightProducer") 
       << " Failed to load LUT = " << lutEfficiencyPtName << " from file = " << inputFile->GetName() << " !!\n";
   lutEfficiencyPt_ = (TH2*)lutEfficiencyPt->Clone(std::string(lutEfficiencyPt->GetName()).append("_cloned").data());
   xAxisEfficiencyPt_ = lutEfficiencyPt_->GetXaxis();
   yAxisEfficiencyPt_ = lutEfficiencyPt_->GetYaxis();
 
   std::string lutEffCorrEtaName = cfg.getParameter<std::string>("lutEffCorrEta");
   TH2* lutEffCorrEta = dynamic_cast<TH2*>(inputFile->Get(lutEffCorrEtaName.data()));
   if ( !lutEffCorrEta ) 
     throw cms::Exception("MuonRadiationCorrWeightProducer") 
       << " Failed to load LUT = " << lutEffCorrEtaName << " from file = " << inputFile->GetName() << " !!\n";
   lutEffCorrEta_ = (TH2*)lutEffCorrEta->Clone(std::string(lutEffCorrEta->GetName()).append("_cloned").data());
   xAxisEffCorrEta_ = lutEffCorrEta_->GetXaxis();
   yAxisEffCorrEta_ = lutEffCorrEta_->GetYaxis();
   
   minWeight_ = cfg.getParameter<double>("minWeight"); 
   maxWeight_ = cfg.getParameter<double>("maxWeight");
 
   verbosity_ = ( cfg.exists("verbosity") ) ?
     cfg.getParameter<int>("verbosity") : 0;
 
   produces<double>("weight");
   produces<double>("weightUp");
   produces<double>("weightDown");
 }

ZmumuEvtSelEffCorrWeightProducer::~ZmumuEvtSelEffCorrWeightProducer ( )

Definition at line 58 of file ZmumuEvtSelEffCorrWeightProducer.cc.

References lutEffCorrEta_, and lutEfficiencyPt_.

 {
   delete lutEfficiencyPt_;
   delete lutEffCorrEta_;
 }

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 80 of file ZmumuEvtSelEffCorrWeightProducer.cc.

References gather_cfg::cout, postValidation_cfi::efficiency, reco::Candidate::eta(), getSelMuons(), getTheMuMinus(), getTheMuPlus(), edm::RefToBase< T >::isNonnull(), lutEffCorrEta_, lutEfficiencyPt_, maxWeight_, minWeight_, reco::Candidate::phi(), reco::Candidate::pt(), edm::Event::put(), Clusterizer1DCommons::square(), srcSelectedMuons_, verbosity_, histoStyle::weight, xAxisEffCorrEta_, xAxisEfficiencyPt_, yAxisEffCorrEta_, and yAxisEfficiencyPt_.

 {
   if ( verbosity_ ) {
     std::cout << "<ZmumuEvtSelEffCorrWeightProducer::produce>:" << std::endl;
     std::cout << " srcSelectedMuons = " << srcSelectedMuons_ << std::endl;
   }
 
   double weight = 1.;
   double weightUp = 2.*maxWeight_;
   double weightDown = 0.;
   
   std::vector<reco::CandidateBaseRef> selMuons = getSelMuons(evt, srcSelectedMuons_);
   const reco::CandidateBaseRef muPlus  = getTheMuPlus(selMuons);
   const reco::CandidateBaseRef muMinus = getTheMuMinus(selMuons);
   if ( muPlus.isNonnull() && muMinus.isNonnull() ) {
     if ( verbosity_ ) {
       std::cout << "Mu+: Pt = " << muPlus->pt() << ", eta = " << muPlus->eta() << ", phi = " << muPlus->phi() << std::endl;
       std::cout << "Mu-: Pt = " << muMinus->pt() << ", eta = " << muMinus->eta() << ", phi = " << muMinus->phi() << std::endl;
     }
 
     int binX_efficiencyPt = findBin(xAxisEfficiencyPt_, muPlus->pt());
     int binY_efficiencyPt = findBin(yAxisEfficiencyPt_, muMinus->pt());
     double efficiencyPt = lutEfficiencyPt_->GetBinContent(binX_efficiencyPt, binY_efficiencyPt);
     double efficiencyPtErr = lutEfficiencyPt_->GetBinError(binX_efficiencyPt, binY_efficiencyPt);
 
     int binX_effCorrEta = findBin(xAxisEffCorrEta_, muPlus->eta());
     int binY_effCorrEta = findBin(yAxisEffCorrEta_, muMinus->eta());
     double effCorrEta = lutEffCorrEta_->GetBinContent(binX_effCorrEta, binY_effCorrEta);
     double effCorrEtaErr = lutEffCorrEta_->GetBinError(binX_effCorrEta, binY_effCorrEta);
 
     double efficiency = efficiencyPt*effCorrEta;
     if ( efficiency > 0. ) {
       weight = 1./efficiency;
       if ( weight > maxWeight_ ) weight = maxWeight_;
       if ( weight < minWeight_ ) weight = minWeight_;
       weightUp = weight + TMath::Sqrt(square(efficiencyPtErr*effCorrEta) + square(efficiencyPt*effCorrEtaErr));
       if ( weightUp > (2.*maxWeight_) ) weightUp = 2.*maxWeight_;
       if ( weightUp < weight ) weightUp = weight;
       weightDown = weight - TMath::Sqrt(square(efficiencyPtErr*effCorrEta) + square(efficiencyPt*effCorrEtaErr));
       if ( weightDown > weight ) weightDown = weight;
       if ( weightDown < 0. ) weightDown = 0.;
     } else {
       weight = maxWeight_;
       weightUp = 2.*maxWeight_;
       weightDown = 0.;
     }
   }
   
   if ( verbosity_ ) {
     std::cout << "--> weight = " << weight << " + " << (weightUp - weight) << " - " << (weight - weightDown) << std::endl;
   }
 
   std::auto_ptr<double> weightPtr(new double(weight));
   evt.put(weightPtr, "weight");
   std::auto_ptr<double> weightUpPtr(new double(weightUp));
   evt.put(weightUpPtr, "weightUp");
   std::auto_ptr<double> weightDownPtr(new double(weightDown));
   evt.put(weightDownPtr, "weightDown");
 }