FSRWeightProducer Class Reference

Inheritance diagram for FSRWeightProducer:

Public Member Functions
	FSRWeightProducer (const edm::ParameterSet &)
	~FSRWeightProducer () override
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
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 ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
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
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
double	alphaRatio (double)
void	beginJob () override
void	endJob () override
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
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
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)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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 35 of file FSRWeightProducer.cc.

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

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

  produces<double>();
}

FSRWeightProducer::~FSRWeightProducer ( )

override

Definition at line 43 of file FSRWeightProducer.cc.

{}

Member Function Documentation

double FSRWeightProducer::alphaRatio ( double  pt )

private

Definition at line 128 of file FSRWeightProducer.cc.

References MaterialEffects_cfi::A, TtFullHadDaughter::B, gen::C, DEFINE_FWK_MODULE, dqm-mbProfile::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.

{}

void FSRWeightProducer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 49 of file FSRWeightProducer.cc.

{}

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

overrideprivate

Definition at line 52 of file FSRWeightProducer.cc.

References funct::abs(), alphaRatio(), funct::cos(), reco::CompositeRefCandidateT< D >::daughter(), reco::CompositeRefCandidateT< D >::daughterRef(), dumpMFGeometry_cfg::delta, reco::Candidate::energy(), reco::LeafCandidate::energy(), genParticles2HepMC_cfi::genParticles, genToken_, edm::Event::getByToken(), mps_fire::i, edm::EventBase::isRealData(), dqmiolumiharvest::j, HLT_2018_cff::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, MaterialEffects_cfi::photonEnergy, HLTObjectMonitor_cfi::photonPhi, 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().