#include <EcalTrigPrimProducer.h>

Inheritance diagram for EcalTrigPrimProducer:

Public Member Functions
void	beginRun (const edm::Run &run, const edm::EventSetup &es) override

	EcalTrigPrimProducer (const edm::ParameterSet &conf)

void	endRun (const edm::Run &, const edm::EventSetup &) override

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

	~EcalTrigPrimProducer () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

Private Member Functions
unsigned long long	getRecords (edm::EventSetup const &setup)

Private Attributes
std::unique_ptr< EcalTrigPrimFunctionalAlgo >	algo_

bool	barrelOnly_

int	binOfMaximum_

unsigned long long	cacheID_

bool	debug_

bool	famos_

bool	fillBinOfMaximumFromHistory_

bool	tcpFormat_

edm::EDGetTokenT< EBDigiCollection >	tokenEB_

edm::EDGetTokenT< EEDigiCollection >	tokenEE_

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

EcalTrigPrimProducer produces a EcalTrigPrimDigiCollection The barrel code does a detailed simulation The code for the endcap is simulated in a rough way, due to missing strip geometry

Author: Ursula Berthon, Stephanie Baffioni, LLR Palaiseau

Version: 1st Version may 2006; 2nd Version jul 2006

EcalTrigPrimProducer produces a EcalTrigPrimDigiCollection Simulation as close as possible to hardware Main algorithm is EcalTrigPrimFunctionalAlgo which is now templated to take EBdataFrames/EEDataFrames as input

Author: Ursula Berthon, Stephanie Baffioni, Pascal Paganini, LLR Palaiseau

Version: 1st Version may 2006; 2nd Version jul 2006; 3rd Version nov 2006; 4th Version apr 2007 full endcap

Definition at line 35 of file EcalTrigPrimProducer.h.

Constructor & Destructor Documentation

EcalTrigPrimProducer::EcalTrigPrimProducer ( const edm::ParameterSet & conf )

explicit

Definition at line 64 of file EcalTrigPrimProducer.cc.

References tcpFormat_.

     : barrelOnly_(iConfig.getParameter<bool>("BarrelOnly")),
       tcpFormat_(iConfig.getParameter<bool>("TcpOutput")),
       debug_(iConfig.getParameter<bool>("Debug")),
       famos_(iConfig.getParameter<bool>("Famos")),
       tokenEB_(consumes<EBDigiCollection>(
           edm::InputTag(iConfig.getParameter<std::string>("Label"),
                         iConfig.getParameter<std::string>("InstanceEB")))),
       tokenEE_(consumes<EEDigiCollection>(
           edm::InputTag(iConfig.getParameter<std::string>("Label"),
                         iConfig.getParameter<std::string>("InstanceEE")))),
       binOfMaximum_(iConfig.getParameter<int>("binOfMaximum")),
       fillBinOfMaximumFromHistory_(-1 == binOfMaximum_) {
   // register your products
   produces<EcalTrigPrimDigiCollection>();
   if (tcpFormat_)
     produces<EcalTrigPrimDigiCollection>("formatTCP");
 }

EcalTrigPrimProducer::~EcalTrigPrimProducer ( )

override

Definition at line 241 of file EcalTrigPrimProducer.cc.

241 {}

Member Function Documentation

void EcalTrigPrimProducer::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	es
	)

override

Definition at line 132 of file EcalTrigPrimProducer.cc.

References algo_, barrelOnly_, binOfMaximum_, cacheID_, debug_, famos_, fillBinOfMaximumFromHistory_, findBinOfMaximum(), getRecords(), edm::Run::processHistory(), and tcpFormat_.

                                                                 {
   // ProcessHistory is guaranteed to be constant for an entire Run
   binOfMaximum_ = findBinOfMaximum(fillBinOfMaximumFromHistory_, binOfMaximum_,
                                    run.processHistory());
 
   algo_.reset(new EcalTrigPrimFunctionalAlgo(setup, binOfMaximum_, tcpFormat_,
                                              barrelOnly_, debug_, famos_));
 
   // get a first version of the records
   cacheID_ = this->getRecords(setup);
 }

void EcalTrigPrimProducer::endRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

override

Definition at line 145 of file EcalTrigPrimProducer.cc.

References algo_.

                                                               {
   algo_.reset();
 }

void EcalTrigPrimProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 344 of file EcalTrigPrimProducer.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addDefault(), kComment, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                 {
 
   edm::ParameterSetDescription desc;
   desc.add<bool>("BarrelOnly", false);
   desc.add<bool>("TcpOutput", false);
   desc.add<bool>("Debug", false);
   desc.add<bool>("Famos", false);
   desc.add<std::string>("Label", "simEcalUnsuppressedDigis");
   desc.add<std::string>("InstanceEB", "");
   desc.add<std::string>("InstanceEE", "");
   const std::string kComment(
       "A value of -1 will make the module lookup the value of 'binOfMaximum' "
       "from the module 'ecalUnsuppressedDigis' from the process history. "
       "Allowed values are -1 and from 1-10.");
   // The code before the existence of fillDescriptions did something special if
   // 'binOfMaximum' was missing. This replicates that behavior.
   desc.add<int>("binOfMaximum", -1)->setComment(kComment);
   descriptions.addDefault(desc);
 }

unsigned long long EcalTrigPrimProducer::getRecords ( edm::EventSetup const & setup )

private

Definition at line 151 of file EcalTrigPrimProducer.cc.

References algo_, edm::EventSetup::get(), and edm::ESHandle< T >::product().

Referenced by beginRun(), and produce().

                                                            {
   // get Eventsetup records
 
   // for EcalFenixStrip...
   // get parameter records for xtals
   edm::ESHandle<EcalTPGLinearizationConst> theEcalTPGLinearization_handle;
   setup.get<EcalTPGLinearizationConstRcd>().get(theEcalTPGLinearization_handle);
   const EcalTPGLinearizationConst *ecaltpLin =
       theEcalTPGLinearization_handle.product();
   edm::ESHandle<EcalTPGPedestals> theEcalTPGPedestals_handle;
   setup.get<EcalTPGPedestalsRcd>().get(theEcalTPGPedestals_handle);
   const EcalTPGPedestals *ecaltpPed = theEcalTPGPedestals_handle.product();
   edm::ESHandle<EcalTPGCrystalStatus> theEcalTPGCrystalStatus_handle;
   setup.get<EcalTPGCrystalStatusRcd>().get(theEcalTPGCrystalStatus_handle);
   const EcalTPGCrystalStatus *ecaltpgBadX =
       theEcalTPGCrystalStatus_handle.product();
 
   // for strips
   edm::ESHandle<EcalTPGSlidingWindow> theEcalTPGSlidingWindow_handle;
   setup.get<EcalTPGSlidingWindowRcd>().get(theEcalTPGSlidingWindow_handle);
   const EcalTPGSlidingWindow *ecaltpgSlidW =
       theEcalTPGSlidingWindow_handle.product();
   edm::ESHandle<EcalTPGWeightIdMap> theEcalTPGWEightIdMap_handle;
   setup.get<EcalTPGWeightIdMapRcd>().get(theEcalTPGWEightIdMap_handle);
   const EcalTPGWeightIdMap *ecaltpgWeightMap =
       theEcalTPGWEightIdMap_handle.product();
   edm::ESHandle<EcalTPGWeightGroup> theEcalTPGWEightGroup_handle;
   setup.get<EcalTPGWeightGroupRcd>().get(theEcalTPGWEightGroup_handle);
   const EcalTPGWeightGroup *ecaltpgWeightGroup =
       theEcalTPGWEightGroup_handle.product();
   edm::ESHandle<EcalTPGFineGrainStripEE> theEcalTPGFineGrainStripEE_handle;
   setup.get<EcalTPGFineGrainStripEERcd>().get(
       theEcalTPGFineGrainStripEE_handle);
   const EcalTPGFineGrainStripEE *ecaltpgFgStripEE =
       theEcalTPGFineGrainStripEE_handle.product();
   edm::ESHandle<EcalTPGStripStatus> theEcalTPGStripStatus_handle;
   setup.get<EcalTPGStripStatusRcd>().get(theEcalTPGStripStatus_handle);
   const EcalTPGStripStatus *ecaltpgStripStatus =
       theEcalTPGStripStatus_handle.product();
 
   algo_->setPointers(ecaltpLin, ecaltpPed, ecaltpgSlidW, ecaltpgWeightMap,
                      ecaltpgWeightGroup, ecaltpgFgStripEE, ecaltpgBadX,
                      ecaltpgStripStatus);
 
   // .. and for EcalFenixTcp
   // get parameter records for towers
   edm::ESHandle<EcalTPGFineGrainEBGroup> theEcalTPGFineGrainEBGroup_handle;
   setup.get<EcalTPGFineGrainEBGroupRcd>().get(
       theEcalTPGFineGrainEBGroup_handle);
   const EcalTPGFineGrainEBGroup *ecaltpgFgEBGroup =
       theEcalTPGFineGrainEBGroup_handle.product();
 
   edm::ESHandle<EcalTPGLutGroup> theEcalTPGLutGroup_handle;
   setup.get<EcalTPGLutGroupRcd>().get(theEcalTPGLutGroup_handle);
   const EcalTPGLutGroup *ecaltpgLutGroup = theEcalTPGLutGroup_handle.product();
 
   edm::ESHandle<EcalTPGLutIdMap> theEcalTPGLutIdMap_handle;
   setup.get<EcalTPGLutIdMapRcd>().get(theEcalTPGLutIdMap_handle);
   const EcalTPGLutIdMap *ecaltpgLut = theEcalTPGLutIdMap_handle.product();
 
   edm::ESHandle<EcalTPGFineGrainEBIdMap> theEcalTPGFineGrainEBIdMap_handle;
   setup.get<EcalTPGFineGrainEBIdMapRcd>().get(
       theEcalTPGFineGrainEBIdMap_handle);
   const EcalTPGFineGrainEBIdMap *ecaltpgFineGrainEB =
       theEcalTPGFineGrainEBIdMap_handle.product();
 
   edm::ESHandle<EcalTPGFineGrainTowerEE> theEcalTPGFineGrainTowerEE_handle;
   setup.get<EcalTPGFineGrainTowerEERcd>().get(
       theEcalTPGFineGrainTowerEE_handle);
   const EcalTPGFineGrainTowerEE *ecaltpgFineGrainTowerEE =
       theEcalTPGFineGrainTowerEE_handle.product();
 
   edm::ESHandle<EcalTPGTowerStatus> theEcalTPGTowerStatus_handle;
   setup.get<EcalTPGTowerStatusRcd>().get(theEcalTPGTowerStatus_handle);
   const EcalTPGTowerStatus *ecaltpgBadTT =
       theEcalTPGTowerStatus_handle.product();
 
   edm::ESHandle<EcalTPGSpike> theEcalTPGSpike_handle;
   setup.get<EcalTPGSpikeRcd>().get(theEcalTPGSpike_handle);
   const EcalTPGSpike *ecaltpgSpike = theEcalTPGSpike_handle.product();
 
   algo_->setPointers2(ecaltpgFgEBGroup, ecaltpgLutGroup, ecaltpgLut,
                       ecaltpgFineGrainEB, ecaltpgFineGrainTowerEE, ecaltpgBadTT,
                       ecaltpgSpike);
 
   // we will suppose that everything is to be updated if the
   // EcalTPGLinearizationConstRcd has changed
   return setup.get<EcalTPGLinearizationConstRcd>().cacheIdentifier();
 }

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

override

Definition at line 244 of file EcalTrigPrimProducer.cc.

References algo_, Reference_intrackfit_cff::barrel, barrelOnly_, cacheID_, edm::eventsetup::EventSetupRecord::cacheIdentifier(), makeMuonMisalignmentScenario::endcap, Exception, edm::EventSetup::get(), edm::Event::getByToken(), getRecords(), mps_fire::i, edm::EventBase::id(), tablePrinter::labels, LogDebug, edm::ProductLabels::module, eostools::move(), edm::print(), edm::Handle< T >::product(), edm::ProductLabels::productInstance, edm::Event::put(), edm::DataFrameContainer::size(), tcpFormat_, tokenEB_, and tokenEE_.

                                                                 {
 
   // update constants if necessary
   if (iSetup.get<EcalTPGLinearizationConstRcd>().cacheIdentifier() != cacheID_)
     cacheID_ = this->getRecords(iSetup);
 
   // get input collections
 
   edm::Handle<EBDigiCollection> ebDigis;
   edm::Handle<EEDigiCollection> eeDigis;
   bool barrel = true;
   bool endcap = true;
   if (barrelOnly_)
     endcap = false;
 
   if (!e.getByToken(tokenEB_, ebDigis)) {
     barrel = false;
     edm::EDConsumerBase::Labels labels;
     labelsForToken(tokenEB_, labels);
     edm::LogWarning("EcalTPG")
         << " Couldnt find Barrel dataframes with producer " << labels.module
         << " and label " << labels.productInstance << "!!!";
   }
   if (!barrelOnly_) {
     if (!e.getByToken(tokenEE_, eeDigis)) {
       endcap = false;
       edm::EDConsumerBase::Labels labels;
       labelsForToken(tokenEE_, labels);
       edm::LogWarning("EcalTPG")
           << " Couldnt find Endcap dataframes with producer " << labels.module
           << " and label " << labels.productInstance << "!!!";
     }
   }
   if (!barrel && !endcap) {
     edm::EDConsumerBase::Labels labels;
     labelsForToken(tokenEB_, labels);
     throw cms::Exception(" ProductNotFound")
         << "No EBDataFrames(EEDataFrames) with producer " << labels.module
         << " and label " << labels.productInstance << " found in input!!\n";
   }
 
   if (!barrelOnly_)
     LogDebug("EcalTPG") << " =================> Treating event  " << e.id()
                         << ", Number of EBDataFrames "
                         << ebDigis.product()->size()
                         << ", Number of EEDataFrames "
                         << eeDigis.product()->size();
   else
     LogDebug("EcalTPG") << " =================> Treating event  " << e.id()
                         << ", Number of EBDataFrames "
                         << ebDigis.product()->size();
 
   auto pOut = std::make_unique<EcalTrigPrimDigiCollection>();
   auto pOutTcp = std::make_unique<EcalTrigPrimDigiCollection>();
 
   // invoke algorithm
 
   const EBDigiCollection *ebdc = nullptr;
   const EEDigiCollection *eedc = nullptr;
   if (barrel) {
     ebdc = ebDigis.product();
     algo_->run(iSetup, ebdc, *pOut, *pOutTcp);
   }
 
   if (endcap) {
     eedc = eeDigis.product();
     algo_->run(iSetup, eedc, *pOut, *pOutTcp);
   }
 
   edm::LogInfo("produce") << "For Barrel + Endcap, " << pOut->size()
                           << " TP  Digis were produced";
 
   //  debug prints if TP >0
 
   for (unsigned int i = 0; i < pOut->size(); ++i) {
     bool print = false;
     for (int isam = 0; isam < (*pOut)[i].size(); ++isam) {
       if ((*pOut)[i][isam].raw())
         print = true;
     }
     if (print)
       LogDebug("EcalTPG") << " For tower  " << (((*pOut)[i])).id() << ", TP is "
                           << (*pOut)[i];
   }
   if (barrelOnly_)
     LogDebug("EcalTPG") << "\n =================> For Barrel , " << pOut->size()
                         << " TP  Digis were produced (including zero ones)";
   else
     LogDebug("EcalTPG") << "\n =================> For Barrel + Endcap, "
                         << pOut->size()
                         << " TP  Digis were produced (including zero ones)";
 
   // put result into the Event
 
   e.put(std::move(pOut));
   if (tcpFormat_)
     e.put(std::move(pOutTcp), "formatTCP");
 }

Member Data Documentation

std::unique_ptr<EcalTrigPrimFunctionalAlgo> EcalTrigPrimProducer::algo_

private

Definition at line 48 of file EcalTrigPrimProducer.h.

Referenced by beginRun(), endRun(), getRecords(), and produce().

bool EcalTrigPrimProducer::barrelOnly_

private

Definition at line 49 of file EcalTrigPrimProducer.h.

Referenced by beginRun(), and produce().

int EcalTrigPrimProducer::binOfMaximum_

private

Definition at line 56 of file EcalTrigPrimProducer.h.

Referenced by beginRun().

unsigned long long EcalTrigPrimProducer::cacheID_

private

Definition at line 61 of file EcalTrigPrimProducer.h.

Referenced by beginRun(), and produce().

bool EcalTrigPrimProducer::debug_

private

Definition at line 51 of file EcalTrigPrimProducer.h.

Referenced by beginRun().

bool EcalTrigPrimProducer::famos_

private

Definition at line 52 of file EcalTrigPrimProducer.h.

Referenced by beginRun().

bool EcalTrigPrimProducer::fillBinOfMaximumFromHistory_

private

Definition at line 57 of file EcalTrigPrimProducer.h.

Referenced by beginRun().

bool EcalTrigPrimProducer::tcpFormat_

private

Definition at line 50 of file EcalTrigPrimProducer.h.

Referenced by beginRun(), EcalTrigPrimProducer(), and produce().

edm::EDGetTokenT<EBDigiCollection> EcalTrigPrimProducer::tokenEB_

private

Definition at line 53 of file EcalTrigPrimProducer.h.

Referenced by produce().

edm::EDGetTokenT<EEDigiCollection> EcalTrigPrimProducer::tokenEE_

private

Definition at line 54 of file EcalTrigPrimProducer.h.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation