#include <PhysicsTools/TagProbeMassProducer/src/TagProbeMassProducer.cc>

Inheritance diagram for TagProbeMassProducer:

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

	~TagProbeMassProducer () 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
void	beginJob () override

void	endJob () override

bool	isPassingProbe (const unsigned int iprobe) const

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

Private Attributes
double	delRMaxCut_

double	delRMinCut_

double	massMaxCut_

double	massMinCut_

edm::InputTag	passingProbeCollection_

edm::Handle< reco::CandidateView >	passingProbes

edm::InputTag	probeCollection_

edm::Handle< reco::CandidateView >	probes

bool	requireOS_

edm::InputTag	tagCollection_

edm::Handle< reco::CandidateView >	tags

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

Description: <one line="" class="" summary>="">

Usage: <usage>

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 34 of file TagProbeMassProducer.h.

Constructor & Destructor Documentation

TagProbeMassProducer::TagProbeMassProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 37 of file TagProbeMassProducer.cc.

References delRMaxCut_, delRMinCut_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), massMaxCut_, massMinCut_, passingProbeCollection_, probeCollection_, requireOS_, and tagCollection_.

 {
    tagCollection_   = iConfig.getParameter<edm::InputTag>("TagCollection");
    probeCollection_ = iConfig.getParameter<edm::InputTag>("ProbeCollection");
    passingProbeCollection_ = iConfig.getParameter<edm::InputTag>("PassingProbeCollection");
 
    massMinCut_      = iConfig.getUntrackedParameter<double>("MassMinCut",50.0);
    massMaxCut_      = iConfig.getUntrackedParameter<double>("MassMaxCut",120.0);
    delRMinCut_      = iConfig.getUntrackedParameter<double>("DelRMinCut",0.0);
    delRMaxCut_      = iConfig.getUntrackedParameter<double>("DelRMaxCut",10000.0);
 
    requireOS_       = iConfig.getUntrackedParameter<bool>("RequireOS",true);
 
    produces<std::vector<float> >("TPmass");
 }

TagProbeMassProducer::~TagProbeMassProducer ( )

override

Definition at line 54 of file TagProbeMassProducer.cc.

 {
  
 }

Member Function Documentation

void TagProbeMassProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 134 of file TagProbeMassProducer.cc.

135 {

136 }

void TagProbeMassProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 140 of file TagProbeMassProducer.cc.

140 {

141 }

bool TagProbeMassProducer::isPassingProbe ( const unsigned int iprobe ) const

private

Definition at line 144 of file TagProbeMassProducer.cc.

References DEFINE_FWK_MODULE, passingProbes, probes, edm::View< T >::refAt(), and edm::View< T >::size().

Referenced by produce().

                                                                           {
 
   if (iProbe > probes->size()) return false;
 
   edm::RefToBase<reco::Candidate> probeRef = probes->refAt(iProbe);
   edm::RefToBase<reco::Candidate> passingProbeRef;
 
   unsigned int numPassingProbes = passingProbes->size();
 
   for (unsigned int iPassProbe = 0; iPassProbe < numPassingProbes; ++iPassProbe) {
     passingProbeRef = passingProbes->refAt(iPassProbe);
     if (passingProbeRef == probeRef) {
       return true;      
     }
   }
   return false;
 }

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

overrideprivate

Definition at line 66 of file TagProbeMassProducer.cc.

References edm::RefToBaseVector< T >::begin(), delRMaxCut_, delRMinCut_, edm::RefToBaseVector< T >::end(), edm::Event::getByLabel(), mps_fire::i, reco::tau::disc::InvariantMass(), isPassingProbe(), edm::HandleBase::isValid(), massMaxCut_, massMinCut_, eostools::move(), passingProbeCollection_, passingProbes, probeCollection_, probes, edm::RefToBaseVector< T >::push_back(), edm::Event::put(), edm::View< T >::refAt(), requireOS_, Validation_hcalonly_cfi::sign, edm::View< T >::size(), GlobalPosition_Frontier_DevDB_cff::tag, tagCollection_, and tags.

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

 {
 
    // We need the output Muon association collection to fill
    std::unique_ptr<std::vector<float> > TPmass( new std::vector<float>);
 
    if ( !iEvent.getByLabel( tagCollection_, tags ) ) {
       edm::LogWarning("TagProbe") << "Could not extract tag muons with input tag "
                  << tagCollection_;
    }
 
    if ( !iEvent.getByLabel( probeCollection_, probes ) ) {
       edm::LogWarning("TagProbe") << "Could not extract probe muons with input tag "
                  << probeCollection_;
    }
 
    if ( !iEvent.getByLabel( passingProbeCollection_, passingProbes ) ) {
       edm::LogWarning("TagProbe") << "Could not extract passing probe muons with input tag "
                  << passingProbeCollection_;
    }
 
    // Loop over Tag and associate with Probes
    if( tags.isValid() && probes.isValid() )
    {
       edm::RefToBaseVector<reco::Candidate> vtags;
       for (size_t i = 0; i < tags->size(); ++i) {
          vtags.push_back(tags->refAt(i));
       }
       edm::RefToBaseVector<reco::Candidate> vprobes;
       for (size_t i = 0; i < probes->size(); ++i) {
          vprobes.push_back(probes->refAt(i));
       }
 
       int itag = 0;
       edm::RefToBaseVector<reco::Candidate>::const_iterator tag = vtags.begin();
       for( ; tag != vtags.end(); ++tag, ++itag ) 
       {  
      int iprobe = 0;
      edm::RefToBaseVector<reco::Candidate>::const_iterator probe = vprobes.begin();
      for( ; probe != vprobes.end(); ++probe, ++iprobe ) 
      {
         // Tag-Probe invariant mass cut
         double invMass = ROOT::Math::VectorUtil::InvariantMass((*tag)->p4(), (*probe)->p4());
         if( invMass < massMinCut_ ) continue;
         if( invMass > massMaxCut_ ) continue;
 
         // Tag-Probe deltaR cut
         double delR = reco::deltaR<double>((*tag)->eta(),(*tag)->phi(),(*probe)->eta(),(*probe)->phi());
         if( delR < delRMinCut_ ) continue;
         if( delR > delRMaxCut_ ) continue;
 
         // Tag-Probe opposite sign
         int sign = (*tag)->charge() * (*probe)->charge();
         if( requireOS_ && sign > 0 ) continue;
 
         bool isPassing = isPassingProbe (iprobe);
 
         if (isPassing) TPmass->push_back(invMass);
      }   
       }
    }
 
    // Finally put the tag probe collection in the event
    iEvent.put(std::move(TPmass),"TPmass");
 }