#include <HcalTrigPrimDigiProducer.h>

Inheritance diagram for HcalTrigPrimDigiProducer:

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

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

virtual	~HcalTrigPrimDigiProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
bool	HFEMB_

std::vector< edm::InputTag >	inputLabel_
	input tags for HCAL digis More...

edm::InputTag	inputTagFEDRaw_
	input tag for FEDRawDataCollection More...

edm::ParameterSet	LongShortCut_

double	LongShortOffset_

double	LongShortSlope_

double	MinLongEnergy_

double	MinShortEnergy_

bool	runFrontEndFormatError_

bool	runZS_

HcalTriggerPrimitiveAlgo	theAlgo_

edm::EDGetTokenT < HBHEDigiCollection >	tok_hbhe_

edm::EDGetTokenT < HFDigiCollection >	tok_hf_

edm::EDGetTokenT < FEDRawDataCollection >	tok_raw_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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::EDProducer
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 12 of file HcalTrigPrimDigiProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 28 of file HcalTrigPrimDigiProducer.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HFEMB_, inputLabel_, inputTagFEDRaw_, LongShortCut_, LongShortOffset_, LongShortSlope_, MinLongEnergy_, MinShortEnergy_, runFrontEndFormatError_, HcalTriggerPrimitiveAlgo::setPeakFinderAlgorithm(), theAlgo_, tok_hbhe_, tok_hf_, and tok_raw_.

 : 
   theAlgo_(ps.getParameter<bool>("peakFilter"),
       ps.getParameter<std::vector<double> >("weights"),
       ps.getParameter<int>("latency"),
       ps.getParameter<uint32_t>("FG_threshold"),
           ps.getParameter<uint32_t>("ZS_threshold"),
       ps.getParameter<int>("numberOfSamples"),
       ps.getParameter<int>("numberOfPresamples"),
       ps.getParameter<int>("numberOfSamplesHF"),
       ps.getParameter<int>("numberOfPresamplesHF"),
       ps.getParameter<uint32_t>("MinSignalThreshold"),
       ps.getParameter<uint32_t>("PMTNoiseThreshold")
    ),
   inputLabel_(ps.getParameter<std::vector<edm::InputTag> >("inputLabel")),
   inputTagFEDRaw_(ps.getParameter<edm::InputTag> ("InputTagFEDRaw")),
   runZS_(ps.getParameter<bool>("RunZS")),
   runFrontEndFormatError_(ps.getParameter<bool>("FrontEndFormatError"))
 {
     HFEMB_ = false;
     if(ps.exists("LSConfig"))
     {
         LongShortCut_ = ps.getUntrackedParameter<edm::ParameterSet>("LSConfig");
         HFEMB_ = LongShortCut_.getParameter<bool>("HcalFeatureHFEMBit");
         MinLongEnergy_ = LongShortCut_.getParameter<double>("Min_Long_Energy"); //minimum long energy
         MinShortEnergy_ = LongShortCut_.getParameter<double>("Min_Short_Energy"); //minimum short energy
         LongShortSlope_ = LongShortCut_.getParameter<double>("Long_vrs_Short_Slope"); //slope of the line that cuts are based on
         LongShortOffset_ = LongShortCut_.getParameter<double>("Long_Short_Offset"); //offset of line
     }
   // register for data access
   if (runFrontEndFormatError_) {
     tok_raw_ = consumes<FEDRawDataCollection>(inputTagFEDRaw_);
   }
   tok_hbhe_ = consumes<HBHEDigiCollection>(inputLabel_[0]);
   tok_hf_ = consumes<HFDigiCollection>(inputLabel_[1]);
 
    produces<HcalTrigPrimDigiCollection>();
    theAlgo_.setPeakFinderAlgorithm(ps.getParameter<int>("PeakFinderAlgorithm"));
 }

virtual HcalTrigPrimDigiProducer::~HcalTrigPrimDigiProducer ( )

inlinevirtual

Definition at line 17 of file HcalTrigPrimDigiProducer.h.

17 {}

Member Function Documentation

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

virtual

Produces the EDM products,

Implements edm::EDProducer.

Definition at line 69 of file HcalTrigPrimDigiProducer.cc.

References edm::EventSetup::get(), edm::Event::getByToken(), HFEMB_, inputLabel_, inputTagFEDRaw_, edm::HandleBase::isValid(), LongShortOffset_, LongShortSlope_, MinLongEnergy_, MinShortEnergy_, edm::Handle< T >::product(), edm::ESHandle< class >::product(), edm::Event::put(), query::result, HcalTriggerPrimitiveAlgo::run(), HcalTriggerPrimitiveAlgo::runFEFormatError(), runFrontEndFormatError_, HcalTriggerPrimitiveAlgo::runZS(), runZS_, theAlgo_, tok_hbhe_, tok_hf_, and tok_raw_.

                                                                                         {
 
   // Step A: get the conditions, for the decoding
   edm::ESHandle<HcalTPGCoder> inputCoder;
   eventSetup.get<HcalTPGRecord>().get(inputCoder);
 
   edm::ESHandle<CaloTPGTranscoder> outTranscoder;
   eventSetup.get<CaloTPGRecord>().get(outTranscoder);
 
   edm::ESHandle<HcalLutMetadata> lutMetadata;
   eventSetup.get<HcalLutMetadataRcd>().get(lutMetadata);
   float rctlsb = lutMetadata->getRctLsb();
 
   edm::ESHandle<HcalTrigTowerGeometry> pG;
   eventSetup.get<CaloGeometryRecord>().get(pG);
   
   // Step B: Create empty output
   std::auto_ptr<HcalTrigPrimDigiCollection> result(new HcalTrigPrimDigiCollection());
 
   edm::Handle<HBHEDigiCollection> hbheDigis;
   edm::Handle<HFDigiCollection>   hfDigis;
 
   iEvent.getByToken(tok_hbhe_,hbheDigis);
   iEvent.getByToken(tok_hf_,hfDigis);
 
   // protect here against missing input collections
   // there is no protection in HcalTriggerPrimitiveAlgo
 
   if (!hbheDigis.isValid()) {
       edm::LogInfo("HcalTrigPrimDigiProducer")
               << "\nWarning: HBHEDigiCollection with input tag "
               << inputLabel_[0]
               << "\nrequested in configuration, but not found in the event."
               << "\nQuit returning empty product." << std::endl;
 
       // put empty HcalTrigPrimDigiCollection in the event
       iEvent.put(result);
 
       return;
   }
 
   if (!hfDigis.isValid()) {
       edm::LogInfo("HcalTrigPrimDigiProducer")
               << "\nWarning: HFDigiCollection with input tag "
               << inputLabel_[1]
               << "\nrequested in configuration, but not found in the event."
               << "\nQuit returning empty product." << std::endl;
 
       // put empty HcalTrigPrimDigiCollection in the event
       iEvent.put(result);
 
       return;
   }
 
 
     edm::ESHandle < HcalDbService > pSetup;
     eventSetup.get<HcalDbRecord> ().get(pSetup);
 
     HcalFeatureBit* hfembit = 0;
 
     if(HFEMB_)
     {
         hfembit = new HcalFeatureHFEMBit(MinShortEnergy_, MinLongEnergy_, LongShortSlope_, LongShortOffset_, *pSetup); //inputs values that cut will be based on
         theAlgo_.run(inputCoder.product(), outTranscoder->getHcalCompressor().get(),
                 *hbheDigis, *hfDigis, *result, &(*pG), rctlsb, hfembit);
 
     }// Step C: Invoke the algorithm, passing in inputs and getting back outputs.
     else
     {
         theAlgo_.run(inputCoder.product(), outTranscoder->getHcalCompressor().get(),
                 *hbheDigis, *hfDigis, *result, &(*pG), rctlsb);
     }
 
   // Step C: Invoke the algorithm, passing in inputs and getting back outputs.
  
 
   // Step C.1: Run FE Format Error / ZS for real data.
   if (runFrontEndFormatError_) {
 
        
         const HcalElectronicsMap *emap = pSetup->getHcalMapping();
 
         edm::Handle < FEDRawDataCollection > fedHandle;
         iEvent.getByToken(tok_raw_, fedHandle);
 
         if (fedHandle.isValid() && emap != 0) {
             theAlgo_.runFEFormatError(fedHandle.product(), emap, *result);
         } else {
             edm::LogInfo("HcalTrigPrimDigiProducer")
                     << "\nWarning: FEDRawDataCollection with input tag "
                     << inputTagFEDRaw_
                     << "\nrequested in configuration, but not found in the event."
                     << "\nQuit returning empty product." << std::endl;
 
             // produce empty HcalTrigPrimDigiCollection and put it in the event
             std::auto_ptr < HcalTrigPrimDigiCollection > emptyResult(
                     new HcalTrigPrimDigiCollection());
 
             iEvent.put(emptyResult);
 
             return;
         }
 
   }
 
   if (runZS_) theAlgo_.runZS(*result);
 
   //  edm::LogInfo("HcalTrigPrimDigiProducer") << "HcalTrigPrims: " << result->size();
 
   // Step D: Put outputs into event
   iEvent.put(result);
 }