#include <HcalRealisticZS.h>

Inheritance diagram for HcalRealisticZS:

Public Member Functions
	HcalRealisticZS (const edm::ParameterSet &ps)

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

virtual	~HcalRealisticZS ()

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

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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 Attributes
std::unique_ptr< HcalZSAlgoRealistic >	algo_

std::string	inputLabel_

edm::EDGetTokenT< HBHEDigiCollection >	tok_hbhe_

edm::EDGetTokenT< QIE11DigiCollection >	tok_hbheQIE11_

edm::EDGetTokenT< HFDigiCollection >	tok_hf_

edm::EDGetTokenT< QIE10DigiCollection >	tok_hfQIE10_

edm::EDGetTokenT< HODigiCollection >	tok_ho_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	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::stream::EDProducerBase
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 HcalRealisticZS.h.

Constructor & Destructor Documentation

HcalRealisticZS::HcalRealisticZS ( const edm::ParameterSet & ps )

explicit

Definition at line 12 of file HcalRealisticZS.cc.

References algo_, edm::ParameterSet::getParameter(), inputLabel_, hcalDigis_cfi::markAndPass, tmp, tok_hbhe_, tok_hbheQIE11_, tok_hf_, tok_hfQIE10_, and tok_ho_.

                                                            :
   inputLabel_(conf.getParameter<std::string>("digiLabel")) {
 
   bool markAndPass=conf.getParameter<bool>("markAndPass");
 
   // register for data access
   tok_hbhe_ = consumes<HBHEDigiCollection>(edm::InputTag(inputLabel_));
   tok_ho_ = consumes<HODigiCollection>(edm::InputTag(inputLabel_));
   tok_hf_ = consumes<HFDigiCollection>(edm::InputTag(inputLabel_));
   tok_hfQIE10_ = consumes<QIE10DigiCollection>(edm::InputTag(inputLabel_, "HFQIE10DigiCollection"));
   tok_hbheQIE11_ = consumes<QIE11DigiCollection>(edm::InputTag(inputLabel_, "HBHEQIE11DigiCollection"));
 
 
   std::vector<int> tmp = conf.getParameter<std::vector<int> >("HBregion");
 
   if(tmp[0]<0 || tmp[0]>9 || tmp[1]<0 || tmp[1]>9 || tmp[0]>tmp[1]) {
     edm::LogError("HcalZeroSuppression") << "ZS(HB) region error: " 
                      << tmp[0] << ":" <<tmp[1];
     tmp[0]=0; tmp[1]=9;
   }
 
   std::pair<int,int> HBsearchTS (tmp[0],tmp[1]);
   
   tmp = conf.getParameter<std::vector<int> >("HEregion");
   if(tmp[0]<0 || tmp[0]>9 || tmp[1]<0 || tmp[1]>9 || tmp[0]>tmp[1]) {
     edm::LogError("HcalZeroSuppression") << "ZS(HE) region error: " 
                      << tmp[0] << ":" <<tmp[1];
     tmp[0]=0; tmp[1]=9;
   }
   std::pair<int,int> HEsearchTS (tmp[0],tmp[1]);
   
   tmp = conf.getParameter<std::vector<int> >("HOregion");
   if(tmp[0]<0 || tmp[0]>9 || tmp[1]<0 || tmp[1]>9 || tmp[0]>tmp[1]) {
     edm::LogError("HcalZeroSuppression") << "ZS(HO) region error: " 
                      << tmp[0] << ":" <<tmp[1];
     tmp[0]=0; tmp[1]=9;
   }
   std::pair<int,int> HOsearchTS (tmp[0],tmp[1]);
   
   tmp = conf.getParameter<std::vector<int> >("HFregion");
   if(tmp[0]<0 || tmp[0]>9 || tmp[1]<0 || tmp[1]>9 || tmp[0]>tmp[1]) {
     edm::LogError("HcalZeroSuppression") << "ZS(HF) region error: " 
                      << tmp[0] << ":" <<tmp[1];
     tmp[0]=0; tmp[1]=9;
   }
   std::pair<int,int> HFsearchTS (tmp[0],tmp[1]);
       
   
     //this constructor will be called if useConfigZSvalues is set to 1 in
     //HcalZeroSuppressionProducers/python/hcalDigisRealistic_cfi.py
     //which means that channel-by-channel ZS thresholds from DB will NOT be used
     if ( conf.getParameter<int>("useConfigZSvalues") ) {
 
       algo_.reset(
      new HcalZSAlgoRealistic (markAndPass,
                   conf.getParameter<int>("HBlevel"),
                   conf.getParameter<int>("HElevel"),
                   conf.getParameter<int>("HOlevel"),
                   conf.getParameter<int>("HFlevel"), 
                   HBsearchTS,
                   HEsearchTS,
                   HOsearchTS,
                   HFsearchTS
                   ));
       
     } else {
 
       algo_.reset(
      new HcalZSAlgoRealistic(markAndPass,                 
                  HBsearchTS,
                  HEsearchTS,
                  HOsearchTS,
                  HFsearchTS));    
     }
 
     produces<HBHEDigiCollection>();
     produces<HODigiCollection>();
     produces<HFDigiCollection>();
     produces<QIE10DigiCollection>("HFQIE10DigiCollection");
     produces<QIE11DigiCollection>("HBHEQIE11DigiCollection");
    
 }

HcalRealisticZS::~HcalRealisticZS ( )

virtual

Definition at line 95 of file HcalRealisticZS.cc.

References algo_.

                                   {
   algo_->clearDbService();
 }

Member Function Documentation

void HcalRealisticZS::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

overridevirtual

Definition at line 99 of file HcalRealisticZS.cc.

References algo_, edm::EventSetup::get(), edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, photonIsolationHIProducer_cfi::hf, photonIsolationHIProducer_cfi::ho, eostools::move(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), edm::Event::put(), HcalDataFrameContainer< Digi >::samples(), edm::SortedCollection< T, SORT >::size(), edm::DataFrameContainer::size(), tok_hbhe_, tok_hbheQIE11_, tok_hf_, tok_hfQIE10_, and tok_ho_.

 {
  
   edm::Handle<HBHEDigiCollection> hbhe;    
   edm::Handle<HODigiCollection> ho;    
   edm::Handle<HFDigiCollection> hf;
   edm::Handle<QIE10DigiCollection> hfQIE10;
   edm::Handle<QIE11DigiCollection> hbheQIE11;
 
   edm::ESHandle<HcalDbService> conditions;
   eventSetup.get<HcalDbRecord>().get(conditions);
   algo_->setDbService(conditions.product());
 
   e.getByToken(tok_hbhe_,hbhe);
   
   // create empty output
   std::unique_ptr<HBHEDigiCollection> zs_hbhe(new HBHEDigiCollection);
   
   e.getByToken(tok_ho_,ho);
   
   // create empty output
   std::unique_ptr<HODigiCollection> zs_ho(new HODigiCollection);
   
   e.getByToken(tok_hf_,hf);
   
   // create empty output
   std::unique_ptr<HFDigiCollection> zs_hf(new HFDigiCollection);
   
   e.getByToken(tok_hfQIE10_,hfQIE10);
   e.getByToken(tok_hbheQIE11_,hbheQIE11);
   
   // create empty output
   std::unique_ptr<HBHEDigiCollection> zs_hbheUpgrade(new HBHEDigiCollection);
   std::unique_ptr<HFDigiCollection> zs_hfUpgrade(new HFDigiCollection);
   std::unique_ptr<QIE10DigiCollection> zs_hfQIE10(new QIE10DigiCollection(hfQIE10->samples()));
   std::unique_ptr<QIE11DigiCollection> zs_hbheQIE11(new QIE11DigiCollection(hbheQIE11->samples()));
   
   //run the algorithm
   algo_->suppress(*(hbhe.product()),*zs_hbhe);
   algo_->suppress(*(ho.product()),*zs_ho);
   algo_->suppress(*(hf.product()),*zs_hf);
   algo_->suppress(*(hfQIE10.product()),*zs_hfQIE10);
   algo_->suppress(*(hbheQIE11.product()),*zs_hbheQIE11);
 
   
   edm::LogInfo("HcalZeroSuppression") << "Suppression (HBHE) input " << hbhe->size() << " digis, output " << zs_hbhe->size() << " digis" 
                       <<  " (HO) input " << ho->size() << " digis, output " << zs_ho->size() << " digis"
                       <<  " (HF) input " << hf->size() << " digis, output " << zs_hf->size() << " digis"
                       <<  " (HFQIE10) input " << hfQIE10->size() << " digis, output " << zs_hfQIE10->size() << " digis"
                       <<  " (HBHEQIE11) input " << hbheQIE11->size() << " digis, output " << zs_hbheQIE11->size() << " digis";
   
 
     // return result
     e.put(std::move(zs_hbhe));
     e.put(std::move(zs_ho));
     e.put(std::move(zs_hf));
     e.put(std::move(zs_hfQIE10),"HFQIE10DigiCollection");
     e.put(std::move(zs_hbheQIE11),"HBHEQIE11DigiCollection");
 
 }