#include <CaloTowerCreatorForTauHLT.h>

Inheritance diagram for CaloTowerCreatorForTauHLT:

Public Member Functions
	CaloTowerCreatorForTauHLT (const edm::ParameterSet &)
	constructor from parameter set More...

	~CaloTowerCreatorForTauHLT () override
	destructor More...

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

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

virtual	~ProducerBase () noexcept(false)

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

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)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &desc)

Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::StreamID sid, edm::Event &evt, const edm::EventSetup &stp) const override
	process one event More...

Private Attributes
const double	mCone
	use only towers in cone mCone around L1 candidate for regional jet reco More...

const double	mEThreshold
	E threshold. More...

const double	mEtThreshold
	imitator of L1 seeds More...

const int	mTauId

const edm::EDGetTokenT< l1extra::L1JetParticleCollection >	mTauTrigger_token
	label of tau trigger type analysis More...

const edm::EDGetTokenT< CaloTowerCollection >	mtowers_token
	label of source collection More...

const int	mVerbose
	verbosity More...

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	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

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

Framework module that produces a collection of calo towers in the region of interest for Tau HLT reconnstruction, depending on tau type trigger: Tau1 - take location of 1st L1 Tau Tau2 - take location of 2nd L1 Tau; if does not exists, take location of 1st Calo Tower ETau - take L1 Tau candidate which is not collinear to HLT (or L1) electron candidate.

Author: A. Nikitenko. IC. based on L. Lista and J. Mans

Definition at line 32 of file CaloTowerCreatorForTauHLT.h.

Constructor & Destructor Documentation

CaloTowerCreatorForTauHLT::CaloTowerCreatorForTauHLT ( const edm::ParameterSet & p )

constructor from parameter set

Definition at line 21 of file CaloTowerCreatorForTauHLT.cc.

   :
   mVerbose (p.getUntrackedParameter<int> ("verbose", 0)),
   mtowers_token (consumes<CaloTowerCollection>(p.getParameter<InputTag>("towers") ) ),
   mCone (p.getParameter<double> ("UseTowersInCone")),
   mTauTrigger_token (consumes<L1JetParticleCollection>(p.getParameter<InputTag>("TauTrigger") ) ),
   mEtThreshold (p.getParameter<double> ("minimumEt")),
   mEThreshold (p.getParameter<double> ("minimumE")),
   mTauId (p.getParameter<int> ("TauId"))
 {
   produces<CaloTowerCollection>();
 }

CaloTowerCreatorForTauHLT::~CaloTowerCreatorForTauHLT ( )

override

destructor

Definition at line 34 of file CaloTowerCreatorForTauHLT.cc.

34 {

35 }

Member Function Documentation

void CaloTowerCreatorForTauHLT::fillDescriptions ( edm::ConfigurationDescriptions & desc )

static

Definition at line 88 of file CaloTowerCreatorForTauHLT.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), edm::ConfigurationDescriptions::setComment(), and edm::ParameterDescriptionNode::setComment().

                                                                                       {
 
   edm::ParameterSetDescription aDesc;
   aDesc.add<edm::InputTag>("TauTrigger",edm::InputTag("l1extraParticles","Tau"))->setComment("L1ExtraJet collection for seeding");
   aDesc.add<int>("TauId",0)->setComment("Item from L1ExtraJet collection used for seeding");
   aDesc.add<edm::InputTag>("towers",edm::InputTag("towerMaker"))->setComment("Input tower collection");
   aDesc.add<double>("UseTowersInCone",0.8)->setComment("Radius of cone around seed");
   aDesc.add<double>("minimumE",0.8)->setComment("Minimum tower energy");
   aDesc.add<double>("minimumEt",0.5)->setComment("Minimum tower ET");
   aDesc.addUntracked<int>("verbose",0)->setComment("Verbosity level; 0=silent");
   desc.setComment("Produce tower collection around L1ExtraJetParticle seed.");
   desc.add("caloTowerMakerHLT",aDesc);
 
 }

void CaloTowerCreatorForTauHLT::produce	(	edm::StreamID	sid,
		edm::Event &	evt,
		const edm::EventSetup &	stp
	)		const

overrideprivate

process one event

Definition at line 37 of file CaloTowerCreatorForTauHLT.cc.

References eleHcalExtractorBlocks_cff::caloTowers, delta, HiRegitMuonDetachedTripletStep_cff::DeltaR, reco::LeafCandidate::energy(), CaloTower::et(), reco::LeafCandidate::eta(), edm::Event::getByToken(), training_settings::idx, mCone, mEThreshold, mEtThreshold, eostools::move(), mTauId, mTauTrigger_token, mtowers_token, mVerbose, AlCaHLTBitMon_ParallelJobs::p, reco::LeafCandidate::phi(), edm::Event::put(), and edm::SortedCollection< T, SORT >::size().

                                                                                                {
   edm::Handle<CaloTowerCollection> caloTowers;
   evt.getByToken( mtowers_token, caloTowers );
 
   // imitate L1 seeds
   edm::Handle<L1JetParticleCollection> jetsgen;
   evt.getByToken( mTauTrigger_token, jetsgen);
 
   std::unique_ptr<CaloTowerCollection> cands( new CaloTowerCollection );
   cands->reserve( caloTowers->size() );
   
   int idTau =0;
   L1JetParticleCollection::const_iterator myL1Jet = jetsgen->begin();
   for(;myL1Jet != jetsgen->end();myL1Jet++)
     {
       if(idTau == mTauId)
     {
       double Sum08 = 0.;
       
       unsigned idx = 0;
       for (; idx < caloTowers->size(); idx++) {
         const CaloTower* cal = &((*caloTowers) [idx]);
         bool isAccepted = false;
         if (mVerbose == 2) {
           edm::LogInfo("JetDebugInfo") << "CaloTowerCreatorForTauHLT::produce-> " << idx << " tower et/eta/phi/e: " 
             << cal->et() << '/' << cal->eta() << '/' << cal->phi() << '/' << cal->energy() << " is...";
         }
         if (cal->et() >= mEtThreshold && cal->energy() >= mEThreshold ) {
           math::PtEtaPhiELorentzVector p( cal->et(), cal->eta(), cal->phi(), cal->energy() );
           double delta  = ROOT::Math::VectorUtil::DeltaR((*myL1Jet).p4().Vect(), p);
           
           if(delta < mCone) {
         isAccepted = true;
         Sum08 += cal->et(); 
         cands->push_back( *cal );
           }
         }
         if (mVerbose == 2){
           if (isAccepted) edm::LogInfo("JetDebugInfo") << "accepted \n";
           else edm::LogInfo("JetDebugInfo") << "rejected \n";
         }
       }
 
     }
       idTau++;
     }
 
   evt.put(std::move(cands));
   
 }