#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::ProducerBase
ProducesCollector	producesCollector ()

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

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

L1RCTProducer::eeFed

Definition: L1RCTProducer.h:78

L1RCTProducer::hbheFed

Definition: L1RCTProducer.h:78

L1RCTProducer::c_max

Definition: L1RCTProducer.h:78

L1RCTProducer::ebOddFed

Definition: L1RCTProducer.h:78

L1RCTProducer::hfFedUp

Definition: L1RCTProducer.h:78

L1RCTProducer::c_min

Definition: L1RCTProducer.h:78

L1RCTProducer::hfFed

Definition: L1RCTProducer.h:78

L1RCTProducer::ebEvenFed

Definition: L1RCTProducer.h:78

Constructor & Destructor Documentation

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

explicit

Definition at line 33 of file L1RCTProducer.cc.

References 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 60 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 79 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 69 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 309 of file L1RCTProducer.cc.

References edm::EventSetup::get(), getFedVectorFromRunInfo(), edm::ESHandleBase::isValid(), RunInfo::m_fed_in, edm::ESHandle< T >::product(), and edmLumisInFiles::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 297 of file L1RCTProducer.cc.

References edm::EventSetup::get(), RunInfo::m_fed_in, edm::ESHandle< T >::product(), and edmLumisInFiles::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 384 of file L1RCTProducer.cc.

References filterCSVwithJSON::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 391 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 400 of file L1RCTProducer.cc.

References gather_cfg::cout, L1RCTChannelMask::ecalMask, fedUpdatedMask, L1RCTChannelMask::hcalMask, L1RCTChannelMask::hfMask, mps_fire::i, dqmiolumiharvest::j, and dqmdumpme::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 328 of file L1RCTProducer.cc.

References bunchCrossings, L1RCT::digiInput(), bsc_activity_cfg::ecal, ecalDigis, Exception, edm::Event::getByLabel(), L1RCT::getIsolatedEGObjects(), L1RCT::getNonisolatedEGObjects(), L1RCT::getRegions(), patCandidatesForDimuonsSequences_cff::hcal, hcalDigis, mps_fire::i, dqmiolumiharvest::j, 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 104 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 139 of file L1RCTProducer.cc.

References c_max, c_min, crateFED, fwrapper::cs, ebEvenFed, ebOddFed, bsc_activity_cfg::ecal, L1RCTChannelMask::ecalMask, eeFed, distMuonTCMETValueMapProducer_cff::eta_max, egammaIdentification::eta_min, fedUpdatedMask, spr::find(), edm::EventSetup::get(), getFedVectorFromOmds(), getFedVectorFromRunInfo(), hbheFed, L1RCTChannelMask::hcalMask, hfFed, hfFedUp, L1RCTChannelMask::hfMask, mps_fire::i, LEDCalibrationChannels::ieta, dqmiolumiharvest::j, dqmdumpme::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 69 of file L1RCTProducer.h.

Referenced by produce().

std::string L1RCTProducer::conditionsLabel

private

Definition at line 73 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}}