#include <L1RCTProducer.h>

Inheritance diagram for L1RCTProducer:

Public Member Functions
void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, const edm::EventSetup &context) final

void	beginRun (edm::Run const &r, const edm::EventSetup &c) final

const std::vector< int >	getFedVectorFromOmds (const edm::EventSetup &)

const std::vector< int >	getFedVectorFromRunInfo (const edm::EventSetup &)

	L1RCTProducer (const edm::ParameterSet &ps)

void	printFedVector (const std::vector< int > &)

void	printUpdatedFedMask ()

void	printUpdatedFedMaskVerbose ()

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

void	updateConfiguration (const edm::EventSetup &)

void	updateFedVector (const edm::EventSetup &, bool getFromOmds, int)

	~L1RCTProducer () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Types
enum	crateSection { c_min, ebOddFed = c_min, ebEvenFed, eeFed, hbheFed, hfFed, hfFedUp, c_max = hfFedUp }

Private Attributes
std::vector< int >	bunchCrossings

std::string	conditionsLabel

std::vector< edm::InputTag >	ecalDigis

L1RCTChannelMask *	fedUpdatedMask

bool	getFedsFromOmds

std::vector< edm::InputTag >	hcalDigis

unsigned int	queryDelayInLS

unsigned int	queryIntervalInLS

L1RCT *	rct

L1RCTLookupTables *	rctLookupTables

bool	useEcal

bool	useHcal

Static Private Attributes
static const int	crateFED [18][6]

static const int	maxBarrel = 17

static const int	maxEndcap = 28

static const int	maxHF = 32

static const int	minBarrel = 1

static const int	minEndcap = 17

static const int	minHF = 29

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 44 of file L1RCTProducer.h.

Member Enumeration Documentation

enum L1RCTProducer::crateSection

private

Enumerator
c_min
ebOddFed
ebEvenFed
eeFed
hbheFed
hfFed
hfFedUp
c_max

Definition at line 79 of file L1RCTProducer.h.

                     {
     c_min,
     ebOddFed = c_min,
     ebEvenFed,
     eeFed,
     hbheFed,
     hfFed,
     hfFedUp,
     c_max = hfFedUp
   };

Constructor & Destructor Documentation

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

explicit

Definition at line 25 of file L1RCTProducer.cc.

References edm::EDConsumerBase::consumes(), ecalDigis, and hcalDigis.

     : rctLookupTables(new L1RCTLookupTables), rct(new L1RCT(rctLookupTables)),
       useEcal(conf.getParameter<bool>("useEcal")),
       useHcal(conf.getParameter<bool>("useHcal")),
       ecalDigis(conf.getParameter<std::vector<edm::InputTag>>("ecalDigis")),
       hcalDigis(conf.getParameter<std::vector<edm::InputTag>>("hcalDigis")),
       bunchCrossings(conf.getParameter<std::vector<int>>("BunchCrossings")),
       getFedsFromOmds(conf.getParameter<bool>("getFedsFromOmds")),
       queryDelayInLS(conf.getParameter<unsigned int>("queryDelayInLS")),
       queryIntervalInLS(conf.getParameter<unsigned int>("queryIntervalInLS")),
       conditionsLabel(conf.getParameter<std::string>("conditionsLabel")),
       fedUpdatedMask(nullptr) {
   produces<L1CaloEmCollection>();
   produces<L1CaloRegionCollection>();
 
   for (unsigned int ihc = 0; ihc < hcalDigis.size(); ihc++) {
     consumes<edm::SortedCollection<
         HcalTriggerPrimitiveDigi,
         edm::StrictWeakOrdering<HcalTriggerPrimitiveDigi>>>(hcalDigis[ihc]);
   }
 
   for (unsigned int iec = 0; iec < ecalDigis.size(); iec++) {
     consumes<edm::SortedCollection<
         EcalTriggerPrimitiveDigi,
         edm::StrictWeakOrdering<EcalTriggerPrimitiveDigi>>>(ecalDigis[iec]);
   }
 }

L1RCTProducer::~L1RCTProducer ( )

override

Definition at line 53 of file L1RCTProducer.cc.

References fedUpdatedMask, rct, and rctLookupTables.

                               {
   if (rct != nullptr)
     delete rct;
   if (rctLookupTables != nullptr)
     delete rctLookupTables;
   if (fedUpdatedMask != nullptr)
     delete fedUpdatedMask;
 }

Member Function Documentation

void L1RCTProducer::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		const edm::EventSetup &	context
	)

final

Definition at line 73 of file L1RCTProducer.cc.

References getFedsFromOmds, edm::LuminosityBlockBase::luminosityBlock(), queryDelayInLS, queryIntervalInLS, edm::LuminosityBlockBase::run(), convertSQLiteXML::runNumber, and updateFedVector().

                                                                        {
   // check LS number every LS, if the checkOMDS flag is set AND it's the right
   // LS, update the FED vector from OMDS can pass the flag as the bool??  but
   // only check LS number if flag is true anyhow
   if (getFedsFromOmds) {
     unsigned int nLumi =
         lumiSeg.luminosityBlock(); // doesn't even need the (unsigned int) cast
                                    // because LuminosityBlockNumber_t is already
                                    // an unsigned int
     // LS count starts at 1, want to be able to delay 0 LS's intuitively
     if (((nLumi - 1) == queryDelayInLS) ||
         (queryIntervalInLS > 0 &&
          nLumi % queryIntervalInLS ==
              0)) // to guard against problems if online DQM crashes; every 100
                  // LS is ~20-30 minutes, not too big a load, hopefully not too
                  // long between
     {
       int runNumber = lumiSeg.run();
       //      std::cout << "Lumi section for this FED vector update is " <<
       // nLumi << std::endl;
       updateFedVector(context, true, runNumber); // OMDS
     } else if (queryIntervalInLS <= 0) {
       // don't do interval checking... cout message??
     }
   }
 }

void L1RCTProducer::beginRun	(	edm::Run const &	r,
		const edm::EventSetup &	c
	)

final

Definition at line 62 of file L1RCTProducer.cc.

References edm::RunBase::run(), convertSQLiteXML::runNumber, updateConfiguration(), and updateFedVector().

                                                               {
   //  std::cout << "getFedsFromOmds is " << getFedsFromOmds << std::endl;
 
   updateConfiguration(eventSetup);
 
   int runNumber = run.run();
   updateFedVector(eventSetup, false,
                   runNumber); // RUNINFO ONLY at beginning of run
 }

const std::vector< int > L1RCTProducer::getFedVectorFromOmds ( const edm::EventSetup & eventSetup )

Definition at line 314 of file L1RCTProducer.cc.

References edm::EventSetup::get(), getFedVectorFromRunInfo(), edm::ESHandleBase::isValid(), RunInfo::m_fed_in, edm::ESHandle< T >::product(), and heppy_report::summary.

Referenced by updateFedVector().

                                                                    {
 
   //  std::cout << "Getting FED vector from my specific ES RunInfo object" <<
   //  std::endl;
 
   // get FULL FED vector from RunInfo object specifically created to have OMDS
   // fed vector
   edm::ESHandle<RunInfo> sum;
   eventSetup.get<RunInfoRcd>().get(
       "OmdsFedVector",
       sum); // using label to get my specific instance of RunInfo
   if (sum.isValid()) {
     const RunInfo *summary = sum.product();
     const std::vector<int> fedvector = summary->m_fed_in;
 
     return fedvector;
   } else {
     return getFedVectorFromRunInfo(eventSetup);
   }
 }

const std::vector< int > L1RCTProducer::getFedVectorFromRunInfo ( const edm::EventSetup & eventSetup )

Definition at line 301 of file L1RCTProducer.cc.

References edm::EventSetup::get(), RunInfo::m_fed_in, edm::ESHandle< T >::product(), and heppy_report::summary.

Referenced by getFedVectorFromOmds(), and updateFedVector().

                                                                       {
   //  std::cout << "Getting FED vector from standard RunInfo object" <<
   //  std::endl;
   // get FULL FED vector from RUNINFO
   edm::ESHandle<RunInfo> sum;
   eventSetup.get<RunInfoRcd>().get(sum);
   const RunInfo *summary = sum.product();
   const std::vector<int> fedvector = summary->m_fed_in;
 
   return fedvector;
 }

void L1RCTProducer::printFedVector ( const std::vector< int > & fedVector )

Definition at line 396 of file L1RCTProducer.cc.

References popcon2dropbox::copy(), and gather_cfg::cout.

                                                                   {
   std::cout << "Contents of given fedVector: ";
   std::copy(fedVector.begin(), fedVector.end(),
             std::ostream_iterator<int>(std::cout, ", "));
   std::cout << std::endl;
 }

void L1RCTProducer::printUpdatedFedMask ( )

Definition at line 404 of file L1RCTProducer.cc.

References gather_cfg::cout, fedUpdatedMask, and L1RCTChannelMask::print().

                                         {
   if (fedUpdatedMask != nullptr) {
     fedUpdatedMask->print(std::cout);
   } else {
     std::cout
         << "Trying to print contents of fedUpdatedMask, but it doesn't exist!"
         << std::endl;
   }
 }

void L1RCTProducer::printUpdatedFedMaskVerbose ( )

Definition at line 415 of file L1RCTProducer.cc.

References gather_cfg::cout, L1RCTChannelMask::ecalMask, fedUpdatedMask, L1RCTChannelMask::hcalMask, L1RCTChannelMask::hfMask, mps_fire::i, and gen::k.

                                                {
   if (fedUpdatedMask != nullptr) {
     // print contents of fedvector
     std::cout << "Contents of fedUpdatedMask: ";
     //       std::copy(fedUpdatedMask.begin(), fedUpdatedMask.end(),
     //       std::ostream_iterator<int>(std::cout, ", "));
     std::cout << "--> ECAL mask: " << std::endl;
     for (int i = 0; i < 18; i++) {
       for (int j = 0; j < 2; j++) {
         for (int k = 0; k < 28; k++) {
           std::cout << fedUpdatedMask->ecalMask[i][j][k] << ", ";
         }
       }
     }
     std::cout << "--> HCAL mask: " << std::endl;
     for (int i = 0; i < 18; i++) {
       for (int j = 0; j < 2; j++) {
         for (int k = 0; k < 28; k++) {
           std::cout << fedUpdatedMask->hcalMask[i][j][k] << ", ";
         }
       }
     }
     std::cout << "--> HF mask: " << std::endl;
     for (int i = 0; i < 18; i++) {
       for (int j = 0; j < 2; j++) {
         for (int k = 0; k < 4; k++) {
           std::cout << fedUpdatedMask->hfMask[i][j][k] << ", ";
         }
       }
     }
 
     std::cout << std::endl;
   } else {
     // print error message
     std::cout
         << "Trying to print contents of fedUpdatedMask, but it doesn't exist!"
         << std::endl;
   }
 }

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

final

Definition at line 335 of file L1RCTProducer.cc.

References bunchCrossings, L1RCT::digiInput(), digitizers_cfi::ecal, ecalDigis, Exception, edm::Event::getByLabel(), L1RCT::getIsolatedEGObjects(), L1RCT::getNonisolatedEGObjects(), L1RCT::getRegions(), digitizers_cfi::hcal, hcalDigis, mps_fire::i, eostools::move(), L1RCT::processEvent(), rct, simplePhotonAnalyzer_cfi::sample, useEcal, and useHcal.

                                                              {
 
   std::unique_ptr<L1CaloEmCollection> rctEmCands(new L1CaloEmCollection);
   std::unique_ptr<L1CaloRegionCollection> rctRegions(
       new L1CaloRegionCollection);
 
   if (!(ecalDigis.size() == hcalDigis.size() &&
         hcalDigis.size() == bunchCrossings.size()))
     throw cms::Exception("BadInput")
         << "From what I see the number of your your ECAL input digi "
            "collections.\n"
         << "is different from the size of your HCAL digi input collections\n"
         << "or the size of your BX factor collection"
         << "They must be the same to correspond to the same Bxs\n"
         << "It does not matter if one of them is empty\n";
 
   // loop through and process each bx
   for (unsigned short sample = 0; sample < bunchCrossings.size(); sample++) {
     edm::Handle<EcalTrigPrimDigiCollection> ecal;
     edm::Handle<HcalTrigPrimDigiCollection> hcal;
 
     EcalTrigPrimDigiCollection ecalIn;
     HcalTrigPrimDigiCollection hcalIn;
 
     if (useHcal && event.getByLabel(hcalDigis[sample], hcal))
       hcalIn = *hcal;
 
     if (useEcal && event.getByLabel(ecalDigis[sample], ecal))
       ecalIn = *ecal;
 
     rct->digiInput(ecalIn, hcalIn);
     rct->processEvent();
 
     // Stuff to create
     for (int j = 0; j < 18; j++) {
       L1CaloEmCollection isolatedEGObjects = rct->getIsolatedEGObjects(j);
       L1CaloEmCollection nonisolatedEGObjects = rct->getNonisolatedEGObjects(j);
       for (int i = 0; i < 4; i++) {
         isolatedEGObjects.at(i).setBx(bunchCrossings[sample]);
         nonisolatedEGObjects.at(i).setBx(bunchCrossings[sample]);
         rctEmCands->push_back(isolatedEGObjects.at(i));
         rctEmCands->push_back(nonisolatedEGObjects.at(i));
       }
     }
 
     for (int i = 0; i < 18; i++) {
       std::vector<L1CaloRegion> regions = rct->getRegions(i);
       for (int j = 0; j < 22; j++) {
         regions.at(j).setBx(bunchCrossings[sample]);
         rctRegions->push_back(regions.at(j));
       }
     }
   }
 
   // putting stuff back into event
   event.put(std::move(rctEmCands));
   event.put(std::move(rctRegions));
 }

void L1RCTProducer::updateConfiguration ( const edm::EventSetup & eventSetup )

Definition at line 101 of file L1RCTProducer.cc.

References conditionsLabel, MillePedeFileConverter_cfg::e, edm::EventSetup::get(), h, edm::ESHandle< T >::product(), alignCSCRings::r, rctLookupTables, alignCSCRings::s, L1RCTLookupTables::setEcalScale(), L1RCTLookupTables::setHcalScale(), L1RCTLookupTables::setL1CaloEtScale(), and L1RCTLookupTables::setRCTParameters().

Referenced by beginRun().

                                                                        {
   // Refresh configuration information every event
   // Hopefully, this does not take too much time
   // There should be a call back function in future to
   // handle changes in configuration
   // parameters to configure RCT (thresholds, etc)
   edm::ESHandle<L1RCTParameters> rctParameters;
   eventSetup.get<L1RCTParametersRcd>().get(conditionsLabel, rctParameters);
   const L1RCTParameters *r = rctParameters.product();
 
   // SCALES
 
   // energy scale to convert eGamma output
   edm::ESHandle<L1CaloEtScale> emScale;
   eventSetup.get<L1EmEtScaleRcd>().get(conditionsLabel, emScale);
   const L1CaloEtScale *s = emScale.product();
 
   // get energy scale to convert input from ECAL
   edm::ESHandle<L1CaloEcalScale> ecalScale;
   eventSetup.get<L1CaloEcalScaleRcd>().get(conditionsLabel, ecalScale);
   const L1CaloEcalScale *e = ecalScale.product();
 
   // get energy scale to convert input from HCAL
   edm::ESHandle<L1CaloHcalScale> hcalScale;
   eventSetup.get<L1CaloHcalScaleRcd>().get(conditionsLabel, hcalScale);
   const L1CaloHcalScale *h = hcalScale.product();
 
   // set scales
   rctLookupTables->setEcalScale(e);
   rctLookupTables->setHcalScale(h);
 
   rctLookupTables->setRCTParameters(r);
   rctLookupTables->setL1CaloEtScale(s);
 }

void L1RCTProducer::updateFedVector	(	const edm::EventSetup &	eventSetup,
		bool	getFromOmds,
		int	runNumber
	)

Definition at line 136 of file L1RCTProducer.cc.

References c_max, c_min, crateFED, fwrapper::cs, ebEvenFed, ebOddFed, digitizers_cfi::ecal, L1RCTChannelMask::ecalMask, eeFed, MuonTCMETValueMapProducer_cff::eta_max, MuonME0RecHits_cfi::eta_min, fedUpdatedMask, spr::find(), edm::EventSetup::get(), getFedVectorFromOmds(), getFedVectorFromRunInfo(), hbheFed, L1RCTChannelMask::hcalMask, hfFed, hfFedUp, L1RCTChannelMask::hfMask, mps_fire::i, gen::k, maxBarrel, maxEndcap, maxHF, minBarrel, minEndcap, minHF, edm::ESHandle< T >::product(), rctLookupTables, L1RCTLookupTables::setChannelMask(), and L1RCTLookupTables::setNoisyChannelMask().

Referenced by beginLuminosityBlock(), and beginRun().

 {
   // list of RCT channels to mask
   edm::ESHandle<L1RCTChannelMask> channelMask;
   eventSetup.get<L1RCTChannelMaskRcd>().get(channelMask);
   const L1RCTChannelMask *cEs = channelMask.product();
 
   // list of Noisy RCT channels to mask
   edm::ESHandle<L1RCTNoisyChannelMask> hotChannelMask;
   eventSetup.get<L1RCTNoisyChannelMaskRcd>().get(hotChannelMask);
   const L1RCTNoisyChannelMask *cEsNoise = hotChannelMask.product();
   rctLookupTables->setNoisyChannelMask(cEsNoise);
 
   // Update the channel mask according to the FED VECTOR
   // This is the beginning of run. We delete the old
   // create the new and set it in the LUTs
 
   if (fedUpdatedMask != nullptr)
     delete fedUpdatedMask;
 
   fedUpdatedMask = new L1RCTChannelMask();
   // copy a constant object
   for (int i = 0; i < 18; i++) {
     for (int j = 0; j < 2; j++) {
       for (int k = 0; k < 28; k++) {
         fedUpdatedMask->ecalMask[i][j][k] = cEs->ecalMask[i][j][k];
         fedUpdatedMask->hcalMask[i][j][k] = cEs->hcalMask[i][j][k];
       }
       for (int k = 0; k < 4; k++) {
         fedUpdatedMask->hfMask[i][j][k] = cEs->hfMask[i][j][k];
       }
     }
   }
 
   //   // adding fed mask into channel mask
 
   const std::vector<int> Feds = getFromOmds
                                     ? getFedVectorFromOmds(eventSetup)
                                     : getFedVectorFromRunInfo(eventSetup);
   // so can create/initialize/assign const quantity in one line accounting for
   // if statement wikipedia says this is exactly what it's for:
   // http://en.wikipedia.org/wiki/%3F:#C.2B.2B
 
   //   std::cout << "Contents of ";
   //   std::cout << (getFromOmds ? "OMDS RunInfo" : "standard RunInfo");
   //   std::cout << " FED vector" << std::endl;
   //   printFedVector(Feds);
 
   bool useUpgradedHF = false;
 
   std::vector<int> caloFeds; // pare down the feds to the interesting ones
   // is this unneccesary?
   // Mike B : This will decrease the find speed so better do it
   for (std::vector<int>::const_iterator cf = Feds.begin(); cf != Feds.end();
        ++cf) {
     int fedNum = *cf;
     if ((fedNum > 600 && fedNum < 724) || fedNum == 1118 || fedNum == 1120 ||
         fedNum == 1122)
       caloFeds.push_back(fedNum);
 
     if (fedNum == 1118 || fedNum == 1120 || fedNum == 1122)
       useUpgradedHF = true;
   }
 
   for (int cr = 0; cr < 18; ++cr) {
 
     for (crateSection cs = c_min; cs <= c_max; cs = crateSection(cs + 1)) {
       bool fedFound = false;
 
       // Try to find the FED
       std::vector<int>::iterator fv =
           std::find(caloFeds.begin(), caloFeds.end(), crateFED[cr][cs]);
       if (fv != caloFeds.end())
         fedFound = true;
 
       if (!fedFound) {
         int eta_min = 0;
         int eta_max = 0;
         bool phi_even[2] = {false}; //, phi_odd = false;
         bool ecal = false;
 
         switch (cs) {
         case ebEvenFed:
           eta_min = minBarrel;
           eta_max = maxBarrel;
           phi_even[0] = true;
           ecal = true;
           break;
 
         case ebOddFed:
           eta_min = minBarrel;
           eta_max = maxBarrel;
           phi_even[1] = true;
           ecal = true;
           break;
 
         case eeFed:
           eta_min = minEndcap;
           eta_max = maxEndcap;
           phi_even[0] = true;
           phi_even[1] = true;
           ecal = true;
           break;
 
         case hbheFed:
           eta_min = minBarrel;
           eta_max = maxEndcap;
           phi_even[0] = true;
           phi_even[1] = true;
           ecal = false;
           break;
 
         case hfFed:
           if (useUpgradedHF)
             break;
 
           eta_min = minHF;
           eta_max = maxHF;
 
           phi_even[0] = true;
           phi_even[1] = true;
           ecal = false;
           break;
 
         case hfFedUp:
           if (!useUpgradedHF)
             break;
 
           eta_min = minHF;
           eta_max = maxHF;
 
           phi_even[0] = true;
           phi_even[1] = true;
           ecal = false;
           break;
 
         default:
           break;
         }
         for (int ieta = eta_min; ieta <= eta_max; ++ieta) {
           if (ieta <= 28) // barrel and endcap
             for (int even = 0; even <= 1; even++) {
               if (phi_even[even]) {
                 if (ecal)
                   fedUpdatedMask->ecalMask[cr][even][ieta - 1] = true;
                 else
                   fedUpdatedMask->hcalMask[cr][even][ieta - 1] = true;
               }
             }
           else
             for (int even = 0; even <= 1; even++)
               if (phi_even[even])
                 fedUpdatedMask->hfMask[cr][even][ieta - 29] = true;
         }
       }
     }
   }
 
   rctLookupTables->setChannelMask(fedUpdatedMask);
 }

Member Data Documentation

std::vector<int> L1RCTProducer::bunchCrossings

private

Definition at line 70 of file L1RCTProducer.h.

Referenced by produce().

std::string L1RCTProducer::conditionsLabel

private

Definition at line 74 of file L1RCTProducer.h.

Referenced by updateConfiguration().

const int L1RCTProducer::crateFED

staticprivate

Initial value:

= {
    {613, 614, 603, 702, 718, 1118}, {611, 612, 602, 700, 718, 1118},
    {627, 610, 601, 716, 722, 1122}, {625, 626, 609, 714, 722, 1122},
    {623, 624, 608, 712, 722, 1122}, {621, 622, 607, 710, 720, 1120},
    {619, 620, 606, 708, 720, 1120}, {617, 618, 605, 706, 720, 1120},
    {615, 616, 604, 704, 718, 1118}, {631, 632, 648, 703, 719, 1118},
    {629, 630, 647, 701, 719, 1118}, {645, 628, 646, 717, 723, 1122},
    {643, 644, 654, 715, 723, 1122}, {641, 642, 653, 713, 723, 1122},
    {639, 640, 652, 711, 721, 1120}, {637, 638, 651, 709, 721, 1120},
    {635, 636, 650, 707, 721, 1120}, {633, 634, 649, 705, 719, 1118}}

Definition at line 90 of file L1RCTProducer.h.

Referenced by updateFedVector().

std::vector<edm::InputTag> L1RCTProducer::ecalDigis

private

Definition at line 68 of file L1RCTProducer.h.

Referenced by L1RCTProducer(), and produce().

L1RCTChannelMask* L1RCTProducer::fedUpdatedMask

private

Definition at line 77 of file L1RCTProducer.h.

Referenced by printUpdatedFedMask(), printUpdatedFedMaskVerbose(), updateFedVector(), and ~L1RCTProducer().

bool L1RCTProducer::getFedsFromOmds

private