Inheritance diagram for FSRWeightProducer:

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

	~FSRWeightProducer ()

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)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
double	alphaRatio (double)

virtual void	beginJob () override

virtual void	endJob () override

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

Private Attributes
edm::EDGetTokenT< reco::GenParticleCollection >	genToken_

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

Definition at line 20 of file FSRWeightProducer.cc.

Constructor & Destructor Documentation

FSRWeightProducer::FSRWeightProducer ( const edm::ParameterSet & pset )

explicit

Definition at line 36 of file FSRWeightProducer.cc.

References genToken_, and edm::ParameterSet::getUntrackedParameter().

                                                                 {
       genToken_ = consumes<reco::GenParticleCollection>(pset.getUntrackedParameter<edm::InputTag> ("GenTag", edm::InputTag("genParticles")));
 
       produces<double>();
 }

FSRWeightProducer::~FSRWeightProducer ( )

Definition at line 43 of file FSRWeightProducer.cc.

43 {}

Member Function Documentation

double FSRWeightProducer::alphaRatio ( double pt )

private

Definition at line 125 of file FSRWeightProducer.cc.

References patCaloMETCorrections_cff::A, TtFullHadDaughter::B, patCaloMETCorrections_cff::C, DEFINE_FWK_MODULE, cmsBatch::log, and M_PI.

Referenced by produce().

                                               {
 
       double pigaga = 0.;
 
       // Leptonic contribution (just one loop, precise at < 0.3% level)
       const double alphapi = 1/137.036/M_PI;
       const double mass_e = 0.0005;
       const double mass_mu = 0.106;
       const double mass_tau = 1.777;
       const double mass_Z = 91.2;
       if (pt>mass_e) pigaga += alphapi * (2*log(pt/mass_e)/3.-5./9.);
       if (pt>mass_mu) pigaga += alphapi * (2*log(pt/mass_mu)/3.-5./9.);
       if (pt>mass_tau) pigaga += alphapi * (2*log(pt/mass_tau)/3.-5./9.);
 
       // Hadronic vaccum contribution
       // Using simple effective parametrization from Physics Letters B 513 (2001) 46.
       // Top contribution neglected
       double A = 0.;
       double B = 0.;
       double C = 0.;
       if (pt<0.7) {
             A = 0.0; B = 0.0023092; C = 3.9925370;
       } else if (pt<2.0) {
             A = 0.0; B = 0.0022333; C = 4.2191779;
       } else if (pt<4.0) {
             A = 0.0; B = 0.0024402; C = 3.2496684;
       } else if (pt<10.0) {
             A = 0.0; B = 0.0027340; C = 2.0995092;
       } else if (pt<mass_Z) {
             A = 0.0010485; B = 0.0029431; C = 1.0;
       } else if (pt<10000.) {
             A = 0.0012234; B = 0.0029237; C = 1.0;
       } else {
             A = 0.0016894; B = 0.0028984; C = 1.0;
       }
       pigaga += A + B*log(1.+C*pt*pt);
 
       // Done
       return 1./(1.-pigaga);
 }

void FSRWeightProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 46 of file FSRWeightProducer.cc.

46 {}

void FSRWeightProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 49 of file FSRWeightProducer.cc.

49 {}

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

overrideprivatevirtual

Definition at line 52 of file FSRWeightProducer.cc.

References funct::abs(), alphaRatio(), funct::cos(), reco::CompositeRefCandidateT< D >::daughter(), reco::CompositeRefCandidateT< D >::daughterRef(), delta, reco::Candidate::energy(), reco::LeafCandidate::energy(), GenHFHadronMatcher_cfi::genParticles, genToken_, edm::Event::getByToken(), mps_fire::i, edm::EventBase::isRealData(), softMuonTagInfos_cfi::leptonId, reco::Candidate::mass(), reco::LeafCandidate::mass(), reco::CompositeRefCandidateT< D >::mother(), eostools::move(), reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::CompositeRefCandidateT< D >::numberOfMothers(), reco::Candidate::pdgId(), reco::LeafCandidate::pdgId(), reco::Candidate::phi(), reco::LeafCandidate::phi(), muons2muons_cfi::photon, funct::pow(), edm::Event::put(), funct::sin(), mathSSE::sqrt(), reco::LeafCandidate::status(), reco::Candidate::theta(), reco::LeafCandidate::theta(), and mps_merge::weight.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                                       {
 
       if (iEvent.isRealData()) return;
 
       edm::Handle<reco::GenParticleCollection> genParticles;
       iEvent.getByToken(genToken_, genParticles);
 
       std::unique_ptr<double> weight (new double);
 
       // Set a default weight to start with
       (*weight) = 1.;
 
       unsigned int gensize = genParticles->size();
       for (unsigned int i = 0; i<gensize; ++i) {
             const reco::GenParticle& lepton = (*genParticles)[i];
             if (lepton.status()!=3) continue;
             int leptonId = lepton.pdgId();
             if (abs(leptonId)!=11 && abs(leptonId)!=13 && abs(leptonId)!=15) continue;
             if (lepton.numberOfMothers()!=1) continue;
             const reco::Candidate * boson = lepton.mother();
             int bosonId = abs(boson->pdgId());
             if (bosonId!=23  && bosonId!=24) continue;
             double bosonMass = boson->mass();
             double leptonMass = lepton.mass();
             double leptonEnergy = lepton.energy();
             double cosLeptonTheta = cos(lepton.theta());
             double sinLeptonTheta = sin(lepton.theta());
             double leptonPhi = lepton.phi();
 
             int trueKey = i;
             if (lepton.numberOfDaughters()==0) {
                   continue;
             } else if (lepton.numberOfDaughters()==1) {
                   int otherleptonKey = lepton.daughterRef(0).key();
                   const reco::GenParticle& otherlepton = (*genParticles)[otherleptonKey];
                   if (otherlepton.pdgId()!=leptonId) continue;
                   if (otherlepton.numberOfDaughters()<=1) continue;
                   trueKey = otherleptonKey;
             }
 
             const reco::GenParticle& trueLepton = (*genParticles)[trueKey];
             unsigned int nDaughters = trueLepton.numberOfDaughters();
 
             for (unsigned int j = 0; j<nDaughters; ++j) {
                   const reco::Candidate * photon = trueLepton.daughter(j);
                   if (photon->pdgId()!=22) continue;
                   double photonEnergy = photon->energy();
                   double cosPhotonTheta = cos(photon->theta());
                   double sinPhotonTheta = sin(photon->theta());
                   double photonPhi = photon->phi();
                   double costheta = sinLeptonTheta*sinPhotonTheta*cos(leptonPhi-photonPhi)
                                   + cosLeptonTheta*cosPhotonTheta;
                   // Missing O(alpha) terms in soft-collinear approach
                   // Only for W, from hep-ph/0303260
                   if (bosonId==24) {
                         double betaLepton = sqrt(1-pow(leptonMass/leptonEnergy,2));
                         double delta = - 8*photonEnergy *(1-betaLepton*costheta)
                           / pow(bosonMass,3)
                           / (1-pow(leptonMass/bosonMass,2))
                           / (4-pow(leptonMass/bosonMass,2))
                           * leptonEnergy * (pow(leptonMass,2)/bosonMass+2*photonEnergy);
                         (*weight) *= (1 + delta);
                   }
                   // Missing NLO QED orders in QED parton shower approach
                   // Change coupling scale from 0 to kT to estimate this effect
                   (*weight) *= alphaRatio(photonEnergy*sqrt(1-pow(costheta,2)));
             }
       }
 
 
       iEvent.put(std::move(weight));
 }

Member Data Documentation

edm::EDGetTokenT<reco::GenParticleCollection> FSRWeightProducer::genToken_

private

Definition at line 31 of file FSRWeightProducer.cc.

Referenced by FSRWeightProducer(), and produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation