L1TStage2Layer1Producer Class Reference

Inheritance diagram for L1TStage2Layer1Producer:

Public Member Functions
	L1TStage2Layer1Producer (const edm::ParameterSet &ps)
	~L1TStage2Layer1Producer ()
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 (const std::string &iProcessName, std::vector< const char * > &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 ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	beginJob () override
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) override
virtual void	endJob () override
virtual void	endRun (edm::Run const &, edm::EventSetup const &) override
virtual void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
int	bxFirst_
int	bxLast_
std::vector< edm::EDGetToken >	ecalToken_
Stage2Layer1FirmwareFactory	factory_
unsigned	fwv_
std::vector< edm::EDGetToken >	hcalToken_
int	ietaMax_
int	ietaMin_
int	iphiMax_
int	iphiMin_
CaloParamsHelper *	params_
unsigned long long	paramsCacheId_
boost::shared_ptr < Stage2PreProcessor >	processor_
int	verbosity_

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
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 62 of file L1TStage2Layer1Producer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 103 of file L1TStage2Layer1Producer.cc.

References bxFirst_, bxLast_, ecalToken_, fwv_, edm::ParameterSet::getParameter(), hcalToken_, and params_.

                                                                          :
  verbosity_(ps.getParameter<int>("verbosity")),
  bxFirst_(ps.getParameter<int>("bxFirst")),
  bxLast_(ps.getParameter<int>("bxLast")),
  ietaMin_(-32),
  ietaMax_(32),
  iphiMin_(1),
  iphiMax_(72),
  ecalToken_(bxLast_+1-bxFirst_),
  hcalToken_(bxLast_+1-bxFirst_),
  paramsCacheId_(0),
  params_(0)
{

  // register what you produce
  produces<CaloTowerBxCollection> ();

  // register what you consume and keep token for later access:
  for (int ibx=0; ibx<bxLast_+1-bxFirst_; ibx++) {
    ecalToken_[ibx] = consumes<EcalTrigPrimDigiCollection>(ps.getParameter<edm::InputTag>("ecalToken"));
    hcalToken_[ibx] = consumes<HcalTrigPrimDigiCollection>(ps.getParameter<edm::InputTag>("hcalToken"));
  }

  // placeholder for the parameters
  params_ = new CaloParamsHelper;

  // set firmware version from python config for now
  fwv_ = ps.getParameter<int>("firmware");

}

L1TStage2Layer1Producer::~L1TStage2Layer1Producer

(

)

Definition at line 134 of file L1TStage2Layer1Producer.cc.

References params_.

                                                  {

  delete params_;

}

Member Function Documentation

void L1TStage2Layer1Producer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 273 of file L1TStage2Layer1Producer.cc.

{
}

void L1TStage2Layer1Producer::beginRun ( edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 284 of file L1TStage2Layer1Producer.cc.

References l1t::Stage2Layer1FirmwareFactory::create(), factory_, fwv_, edm::EventSetup::get(), LogDebug, params_, paramsCacheId_, processor_, edm::ESHandle< class >::product(), and l1t::CaloParamsHelper::~CaloParamsHelper().

{

  // update parameters and algorithms at run start, if they have changed
  // update params first because the firmware factory relies on pointer to params

  // parameters

  unsigned long long id = iSetup.get<L1TCaloParamsRcd>().cacheIdentifier();

  if (id != paramsCacheId_) {

    paramsCacheId_ = id;

    edm::ESHandle<CaloParams> paramsHandle;
    iSetup.get<L1TCaloParamsRcd>().get(paramsHandle);

    // replace our local copy of the parameters with a new one using placement new
    params_->~CaloParamsHelper();
    params_ = new (params_) CaloParamsHelper(*paramsHandle.product());

    LogDebug("L1TDebug") << *params_ << std::endl;

    if (! params_){
      edm::LogError("l1t|caloStage2") << "Could not retrieve params from Event Setup" << std::endl;
    }

  }

  // firmware

  if ( !processor_ ) { // in future, also check if the firmware cache ID has changed !

    //     m_fwv = ; // get new firmware version in future

    // Set the current algorithm version based on DB pars from database:
    processor_ = factory_.create(fwv_, params_);

    LogDebug("L1TDebug") << "Processor object : " << (processor_?1:0) << std::endl;

    if (! processor_) {
      // we complain here once per run
      edm::LogError("l1t|caloStage2") << "Layer 1 firmware could not be configured.\n";
    }

  }


}

void L1TStage2Layer1Producer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 279 of file L1TStage2Layer1Producer.cc.

                                {
}

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

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 337 of file L1TStage2Layer1Producer.cc.

{
}

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

static

Definition at line 361 of file L1TStage2Layer1Producer.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                    {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

void L1TStage2Layer1Producer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 142 of file L1TStage2Layer1Producer.cc.

References funct::abs(), bxFirst_, bxLast_, ecalToken_, edm::EventSetup::get(), edm::Event::getByToken(), hcalToken_, ietaMax_, ietaMin_, iphiMax_, iphiMin_, LogDebug, params_, processor_, edm::Event::put(), l1t::CaloParamsHelper::towerLsbE(), and l1t::CaloParamsHelper::towerLsbH().

{

  LogDebug("l1t|stage 2") << "L1TStage2Layer1Producer::produce function called..." << std::endl;

  // do event setup
  // get RCT input scale objects
  edm::ESHandle<L1CaloEcalScale> ecalScale;
  iSetup.get<L1CaloEcalScaleRcd>().get(ecalScale);
  //  const L1CaloEcalScale* e = ecalScale.product();

  edm::ESHandle<L1CaloHcalScale> hcalScale;
  iSetup.get<L1CaloHcalScaleRcd>().get(hcalScale);
  //  const L1CaloHcalScale* h = hcalScale.product();

  LogDebug("L1TDebug") << "First BX=" << bxFirst_ << ", last BX=" << bxLast_ << ", LSB(E)=" << params_->towerLsbE() << ", LSB(H)=" << params_->towerLsbH() << std::endl;


  // output collection
  std::auto_ptr<CaloTowerBxCollection> towersColl (new CaloTowerBxCollection);


  // loop over crossings
  for (int bx = bxFirst_; bx < bxLast_+1; ++bx) {

    int ibx = bx-bxFirst_;

    edm::Handle<EcalTrigPrimDigiCollection> ecalTPs;
    edm::Handle<HcalTrigPrimDigiCollection> hcalTPs;

    iEvent.getByToken(hcalToken_[ibx], hcalTPs);
    iEvent.getByToken(ecalToken_[ibx], ecalTPs);

    // create input and output tower vectors for this BX
    std::auto_ptr< std::vector<CaloTower> > localInTowers (new std::vector<CaloTower>(CaloTools::caloTowerHashMax()));
    std::auto_ptr< std::vector<CaloTower> > localOutTowers (new std::vector<CaloTower>()); //this is later filled to the same size as localInTowers

    // loop over ECAL TPs
    EcalTrigPrimDigiCollection::const_iterator ecalItr;
    int nEcal=0;
    for (ecalItr=ecalTPs->begin(); ecalItr!=ecalTPs->end(); ++ecalItr, ++nEcal) {

      int ieta = ecalItr->id().ieta();
      int iphi = ecalItr->id().iphi();

      int ietIn = ecalItr->compressedEt();
      bool ifg  = ecalItr->fineGrain();

      // decompress
      double et = ecalScale->et( ietIn, abs(ieta), (ieta>0) );
      int ietOut = floor( et / params_->towerLsbE() );
      //      int ietOutMask = (int) pow(2,params_->towerNBitsE())-1;

      if (ietIn>0)
    LogDebug("L1TDebug") << " ECAL TP : " << ieta << ", " << iphi << ", " << ietIn << ", " << et << ", " << ietOut << std::endl;

      int itow = CaloTools::caloTowerHash(ieta, iphi);
      localInTowers->at(itow).setHwEtEm(ietOut);// & ietOutMask);
      localInTowers->at(itow).setHwQual( localInTowers->at(itow).hwQual() | (ifg ? 0x4 : 0x0) );

    }

    // loop over HCAL TPs
    HcalTrigPrimDigiCollection::const_iterator hcalItr;
    int nHcal=0;
    for (hcalItr=hcalTPs->begin(); hcalItr!=hcalTPs->end(); ++hcalItr, ++nHcal) {

      int ieta = hcalItr->id().ieta();
      int iphi = hcalItr->id().iphi();

      int ietIn = hcalItr->SOI_compressedEt();
      //int ifg = hcalItr->SOI_fineGrain();

      // decompress
      double et = hcalScale->et( ietIn, abs(ieta), (ieta>0) );
      int ietOut = floor( et / params_->towerLsbH() );
      //      int ietOutMask = (int) pow(2,params_->towerNBitsH() )-1;

      if (ietIn>0)
    LogDebug("L1TDebug") << " HCAL TP : " << ieta << ", " << iphi << ", " << ietIn << ", " << et << ", " << ietOut << std::endl;

      int itow = CaloTools::caloTowerHash(ieta, iphi);
      localInTowers->at(itow).setHwEtHad(ietOut);// & ietOutMask);
      //      towers.at(itow).setHwFGHad(ifg);

    }

    // now calculate remaining tower quantities
    for (int ieta=ietaMin_; ieta<ietaMax_+1; ieta++) {

      for (int iphi=iphiMin_; iphi<iphiMax_+1; iphi++) {

    if(!CaloTools::isValidIEtaIPhi(ieta,iphi)) continue;

    int itow = CaloTools::caloTowerHash(ieta, iphi);

    // get ECAL/HCAL raw numbers
    int ietEcal = localInTowers->at(itow).hwEtEm();
    int ietHcal = localInTowers->at(itow).hwEtHad();

    //  const LorentzVector& p4;
    int iet = ietEcal + ietHcal;   // this is nonsense, temp solution!

    //LogDebug("L1TDebug") << " Tower : " << ieta << ", " << iphi << ", " << iet << ", " << ietEcal << ", " << ietHcal << std::endl;

    localInTowers->at(itow).setHwPt(iet);
    localInTowers->at(itow).setHwEta(ieta);
    localInTowers->at(itow).setHwPhi(iphi);

      }
    }

    // do the decompression
    processor_->processEvent(*localInTowers, *localOutTowers);

    // copy towers to output collection
    for(std::vector<CaloTower>::const_iterator tower = localOutTowers->begin();
    tower != localOutTowers->end();
    ++tower)
      towersColl->push_back(ibx, *tower);

    LogDebug("L1TDebug") << "BX=" << ibx << ", N(Tower in)=" << localInTowers->size() << ", N(Tower out)=" << localOutTowers->size() << std::endl;

  }

  iEvent.put(towersColl);

}

Member Data Documentation

int L1TStage2Layer1Producer::bxFirst_

private

Definition at line 84 of file L1TStage2Layer1Producer.cc.

Referenced by L1TStage2Layer1Producer(), and produce().

int L1TStage2Layer1Producer::bxLast_

private

Definition at line 84 of file L1TStage2Layer1Producer.cc.

Referenced by L1TStage2Layer1Producer(), and produce().

std::vector<edm::EDGetToken> L1TStage2Layer1Producer::ecalToken_

private

Definition at line 87 of file L1TStage2Layer1Producer.cc.

Referenced by L1TStage2Layer1Producer(), and produce().

Stage2Layer1FirmwareFactory L1TStage2Layer1Producer::factory_

private

Definition at line 96 of file L1TStage2Layer1Producer.cc.

Referenced by beginRun().

unsigned L1TStage2Layer1Producer::fwv_

private

Definition at line 92 of file L1TStage2Layer1Producer.cc.

Referenced by beginRun(), and L1TStage2Layer1Producer().

std::vector<edm::EDGetToken> L1TStage2Layer1Producer::hcalToken_

private

Definition at line 88 of file L1TStage2Layer1Producer.cc.

Referenced by L1TStage2Layer1Producer(), and produce().

int L1TStage2Layer1Producer::ietaMax_

private

Definition at line 85 of file L1TStage2Layer1Producer.cc.

Referenced by produce().

int L1TStage2Layer1Producer::ietaMin_

private

Definition at line 85 of file L1TStage2Layer1Producer.cc.

Referenced by produce().

int L1TStage2Layer1Producer::iphiMax_

private

Definition at line 85 of file L1TStage2Layer1Producer.cc.

Referenced by produce().

int L1TStage2Layer1Producer::iphiMin_

private

Definition at line 85 of file L1TStage2Layer1Producer.cc.

Referenced by produce().

CaloParamsHelper* L1TStage2Layer1Producer::params_

private

Definition at line 93 of file L1TStage2Layer1Producer.cc.

Referenced by beginRun(), L1TStage2Layer1Producer(), produce(), and ~L1TStage2Layer1Producer().

unsigned long long L1TStage2Layer1Producer::paramsCacheId_

private

Definition at line 91 of file L1TStage2Layer1Producer.cc.

Referenced by beginRun().

boost::shared_ptr<Stage2PreProcessor> L1TStage2Layer1Producer::processor_

private

Definition at line 97 of file L1TStage2Layer1Producer.cc.

Referenced by beginRun(), and produce().

int L1TStage2Layer1Producer::verbosity_

private

Definition at line 82 of file L1TStage2Layer1Producer.cc.