Inheritance diagram for ISRGammaWeightProducer:

Public Member Functions
void	beginJob () override

void	endJob () override

	ISRGammaWeightProducer (const edm::ParameterSet &)

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

	~ISRGammaWeightProducer () 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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

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

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 Attributes
edm::EDGetTokenT< reco::GenParticleCollection >	genToken_

std::vector< double >	isrBinEdges_

std::vector< double >	ptWeights_

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
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 20 of file ISRGammaWeightProducer.cc.

Constructor & Destructor Documentation

◆ ISRGammaWeightProducer()

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

explicit

Definition at line 36 of file ISRGammaWeightProducer.cc.

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

References genToken_, and muonDTDigis_cfi::pset.

◆ ~ISRGammaWeightProducer()

ISRGammaWeightProducer::~ISRGammaWeightProducer ( )

override

Definition at line 44 of file ISRGammaWeightProducer.cc.

44 {}

Member Function Documentation

◆ beginJob()

void ISRGammaWeightProducer::beginJob ( void )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 47 of file ISRGammaWeightProducer.cc.

47 {}

◆ endJob()

void ISRGammaWeightProducer::endJob ( void )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 50 of file ISRGammaWeightProducer.cc.

50 {}

◆ produce()

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

overridevirtual

Implements edm::EDProducer.

Definition at line 53 of file ISRGammaWeightProducer.cc.

                                                                            {
   if (iEvent.isRealData())
     return;
  
   edm::Handle<reco::GenParticleCollection> genParticles;
   iEvent.getByToken(genToken_, genParticles);
   unsigned int gensize = genParticles->size();
  
   std::unique_ptr<double> weight(new double);
  
   // Set a default weight to start with
   (*weight) = 1.;
  
   // Find the boson at the hard scattering level
   const reco::GenParticle* boson = nullptr;
   int parton1Key = -1;
   int parton2Key = -1;
   for (unsigned int i = 0; i < gensize; ++i) {
     const reco::GenParticle& part = (*genParticles)[i];
     int status = abs(part.status());
     if (status != 3)
       continue;
     if (part.numberOfMothers() != 2)
       continue;
     int partId = abs(part.pdgId());
     if (status == 3 && (partId == 23 || abs(partId) == 24)) {
       boson = &(*genParticles)[i];
       parton1Key = part.motherRef(0).key();
       parton2Key = part.motherRef(1).key();
       break;
     }
   }
  
   // Consider only photons near the hard-scattering process
   const reco::GenParticle* photon = nullptr;
   if (boson) {
     for (unsigned int i = 0; i < gensize; ++i) {
       photon = nullptr;
       const reco::GenParticle& part = (*genParticles)[i];
       int status = abs(part.status());
       if (status != 1)
         continue;
       int partId = abs(part.pdgId());
       if (partId != 22)
         continue;
       if (part.numberOfMothers() != 1)
         continue;
       int keyM = part.motherRef(0).key();
       const reco::GenParticle* mother = &(*genParticles)[keyM];
       if (mother->status() != 3)
         continue;
       int mId = mother->pdgId();
       if (abs(mId) > 6 && mId != 2212)
         continue;
       for (unsigned int j = 0; j < mother->numberOfDaughters(); ++j) {
         int keyD = mother->daughterRef(j).key();
         if (keyD == parton1Key || keyD == parton2Key) {
           photon = &part;
           break;
         }
       }
       if (photon)
         break;
     }
   }
  
   if (boson && photon) {
     math::XYZTLorentzVector smom = boson->p4() + photon->p4();
     double s = smom.M2();
     double sqrts = smom.M();
  
     // Go to CM using the boost direction of the boson+photon system
     ROOT::Math::Boost cmboost(smom.BoostToCM());
     math::XYZTLorentzVector photonCM(cmboost(photon->p4()));
     double pcostheta = (smom.x() * photonCM.x() + smom.y() * photonCM.y() + smom.z() * photonCM.z()) / smom.P();
  
     // Determine kinematic invariants
     double t = -sqrts * (photonCM.t() - pcostheta);
     double MV = boson->mass();
     double u = MV * MV - s - t;
     (*weight) = 1. - 2 * t * u / (s * s + MV * MV * MV * MV);
     //printf(">>>>>>>>> s %f t %f u %f, MV %f, weight = %f\n", s, t, u, MV, (*weight));
   }
  
   iEvent.put(std::move(weight));
 }

References funct::abs(), reco::CompositeRefCandidateT< D >::daughterRef(), genParticles2HepMC_cfi::genParticles, genToken_, mps_fire::i, iEvent, dqmiolumiharvest::j, reco::LeafCandidate::mass(), reco::CompositeRefCandidateT< D >::motherRef(), eostools::move(), reco::CompositeRefCandidateT< D >::numberOfDaughters(), reco::LeafCandidate::p4(), reco::LeafCandidate::pdgId(), muons2muons_cfi::photon, alignCSCRings::s, mps_update::status, reco::LeafCandidate::status(), OrderedSet::t, and mps_merge::weight.

Member Data Documentation

◆ genToken_

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

private

Definition at line 30 of file ISRGammaWeightProducer.cc.

Referenced by ISRGammaWeightProducer(), and produce().

◆ isrBinEdges_

std::vector<double> ISRGammaWeightProducer::isrBinEdges_

private

Definition at line 31 of file ISRGammaWeightProducer.cc.

◆ ptWeights_

std::vector<double> ISRGammaWeightProducer::ptWeights_

private

Definition at line 32 of file ISRGammaWeightProducer.cc.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ISRGammaWeightProducer()

◆ ~ISRGammaWeightProducer()

Member Function Documentation

◆ beginJob()

◆ endJob()

◆ produce()

Member Data Documentation

◆ genToken_

◆ isrBinEdges_

◆ ptWeights_