EcalTrigPrimProducer Class Reference

#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	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () 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::EcalTrigPrimProducer ( const edm::ParameterSet &  conf )

explicit

Definition at line 67 of file EcalTrigPrimProducer.cc.

    : 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");
}

References tcpFormat_.

~EcalTrigPrimProducer()

EcalTrigPrimProducer::~EcalTrigPrimProducer ( )

override

Definition at line 221 of file EcalTrigPrimProducer.cc.

{}

Member Function Documentation

beginRun()

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

override

Definition at line 125 of file EcalTrigPrimProducer.cc.

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

References algo_, barrelOnly_, binOfMaximum_, cacheID_, debug_, famos_, fillBinOfMaximumFromHistory_, findBinOfMaximum(), getRecords(), writedatasetfile::run, singleTopDQM_cfi::setup, and tcpFormat_.

endRun()

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

override

Definition at line 135 of file EcalTrigPrimProducer.cc.

{ algo_.reset(); }

References algo_.

fillDescriptions()

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

static

Definition at line 312 of file EcalTrigPrimProducer.cc.

                                                                                      {
  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);
}

References edm::ConfigurationDescriptions::addDefault(), submitPVResolutionJobs::desc, kComment, and AlCaHLTBitMon_QueryRunRegistry::string.

getRecords()

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

private

Definition at line 137 of file EcalTrigPrimProducer.cc.

                                                                              {
  // 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();
}

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

Referenced by beginRun(), and produce().

produce()

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

override

Definition at line 224 of file EcalTrigPrimProducer.cc.

                                                                           {
  // 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");
}

References algo_, Reference_intrackfit_cff::barrel, barrelOnly_, cacheID_, edm::eventsetup::EventSetupRecord::cacheIdentifier(), MillePedeFileConverter_cfg::e, SimL1EmulatorDM_cff::eeDigis, makeMuonMisalignmentScenario::endcap, Exception, edm::EventSetup::get(), getRecords(), mps_fire::i, SummaryClient_cfi::labels, LogDebug, eostools::move(), print(), edm::Handle< T >::product(), edm::DataFrameContainer::size(), tcpFormat_, tokenEB_, and tokenEE_.

Member Data Documentation

algo_

std::unique_ptr<EcalTrigPrimFunctionalAlgo> EcalTrigPrimProducer::algo_

private

Definition at line 48 of file EcalTrigPrimProducer.h.

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

barrelOnly_

bool EcalTrigPrimProducer::barrelOnly_

private

Definition at line 49 of file EcalTrigPrimProducer.h.

Referenced by beginRun(), and produce().

binOfMaximum_

int EcalTrigPrimProducer::binOfMaximum_

private

Definition at line 56 of file EcalTrigPrimProducer.h.

Referenced by beginRun().

cacheID_

unsigned long long EcalTrigPrimProducer::cacheID_

private

Definition at line 61 of file EcalTrigPrimProducer.h.

Referenced by beginRun(), and produce().

debug_

bool EcalTrigPrimProducer::debug_

private

Definition at line 51 of file EcalTrigPrimProducer.h.

Referenced by beginRun().

famos_

bool EcalTrigPrimProducer::famos_

private

Definition at line 52 of file EcalTrigPrimProducer.h.

Referenced by beginRun().

fillBinOfMaximumFromHistory_

bool EcalTrigPrimProducer::fillBinOfMaximumFromHistory_

private

Definition at line 57 of file EcalTrigPrimProducer.h.

Referenced by beginRun().

tcpFormat_

bool EcalTrigPrimProducer::tcpFormat_

private

Definition at line 50 of file EcalTrigPrimProducer.h.

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

tokenEB_

edm::EDGetTokenT<EBDigiCollection> EcalTrigPrimProducer::tokenEB_

private

Definition at line 53 of file EcalTrigPrimProducer.h.

Referenced by produce().

tokenEE_

edm::EDGetTokenT<EEDigiCollection> EcalTrigPrimProducer::tokenEE_

private

Definition at line 54 of file EcalTrigPrimProducer.h.

Referenced by produce().