#include <HcalRealisticZS.h>

Inheritance diagram for HcalRealisticZS:

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

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

	~HcalRealisticZS () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

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_, tok_ho_, and hcalDigisRealistic_cfi::useInstanceLabels.

                                                            :
   inputLabel_(conf.getParameter<std::string>("digiLabel")) {
 
   bool markAndPass=conf.getParameter<bool>("markAndPass");
   bool useInstanceLabels=conf.getParameter<bool>("useInstanceLabels");
 
   // 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_,
                                  useInstanceLabels ? "HFQIE10DigiCollection" : "" ));
   tok_hbheQIE11_ = consumes<QIE11DigiCollection>(edm::InputTag(inputLabel_,
                                    useInstanceLabels ? "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 ( )

override

Definition at line 98 of file HcalRealisticZS.cc.

References algo_.

                                   {
   algo_->clearDbService();
 }

Member Function Documentation

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

override

Definition at line 102 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");
 
 }