#include <L1RCTSaveInput.h>

Inheritance diagram for L1RCTSaveInput:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

	L1RCTSaveInput (const edm::ParameterSet &)

	~L1RCTSaveInput () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
edm::ESGetToken< L1RCTChannelMask, L1RCTChannelMaskRcd >	channelMaskToken_

edm::EDGetTokenT< EcalTrigPrimDigiCollection >	ecalDigisToken_

edm::ESGetToken< L1CaloEcalScale, L1CaloEcalScaleRcd >	ecalScaleToken_

edm::ESGetToken< L1CaloEtScale, L1EmEtScaleRcd >	emScaleToken_

std::string	fileName

edm::EDGetTokenT< HcalTrigPrimDigiCollection >	hcalDigisToken_

edm::ESGetToken< L1CaloHcalScale, L1CaloHcalScaleRcd >	hcalScaleToken_

std::ofstream	ofs

L1RCT *	rct

L1RCTLookupTables *	rctLookupTables

edm::ESGetToken< L1RCTParameters, L1RCTParametersRcd >	rctParametersToken_

EcalTPGScale::Tokens	tokens_

edm::ESGetToken< CaloTPGTranscoder, CaloTPGRecord >	transcoderToken_

bool	useDebugTpgScales

bool	useEcal

bool	useHcal

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

src/L1RCTSaveInput/src/L1RCTSaveInput.cc

Description: Saves the input event from TPGs for loading simulated events in hardware

Implementation: Kind of kludgy – should think of a better way in future

Definition at line 55 of file L1RCTSaveInput.h.

Constructor & Destructor Documentation

◆ L1RCTSaveInput()

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

explicit

Definition at line 26 of file L1RCTSaveInput.cc.

References Exception, fileName, and ofs.

     : fileName(conf.getUntrackedParameter<std::string>("rctTestInputFile")),
       rctLookupTables(new L1RCTLookupTables),
       rct(new L1RCT(rctLookupTables)),
       useEcal(conf.getParameter<bool>("useEcal")),
       useHcal(conf.getParameter<bool>("useHcal")),
       ecalDigisToken_(consumes(conf.getParameter<edm::InputTag>("ecalDigisLabel"))),
       hcalDigisToken_(consumes(conf.getParameter<edm::InputTag>("hcalDigisLabel"))),
       rctParametersToken_(esConsumes<L1RCTParameters, L1RCTParametersRcd>()),
       channelMaskToken_(esConsumes<L1RCTChannelMask, L1RCTChannelMaskRcd>()),
       emScaleToken_(esConsumes<L1CaloEtScale, L1EmEtScaleRcd>()),
       transcoderToken_(esConsumes<CaloTPGTranscoder, CaloTPGRecord>()),
       hcalScaleToken_(esConsumes<L1CaloHcalScale, L1CaloHcalScaleRcd>()),
       ecalScaleToken_(esConsumes<L1CaloEcalScale, L1CaloEcalScaleRcd>()),
       useDebugTpgScales(conf.getParameter<bool>("useDebugTpgScales")),
       tokens_(consumesCollector()) {
   ofs.open(fileName.c_str(), std::ios::app);
   if (!ofs) {
     throw cms::Exception("FailedFileOpen") << "Could not create " << fileName << std::endl;
   }
 }

◆ ~L1RCTSaveInput()

L1RCTSaveInput::~L1RCTSaveInput ( )

override

Definition at line 48 of file L1RCTSaveInput.cc.

References rct, and rctLookupTables.

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

Member Function Documentation

◆ analyze()

void L1RCTSaveInput::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 55 of file L1RCTSaveInput.cc.

References HltBtagPostValidation_cff::c, channelMaskToken_, gather_cfg::cout, TauDecayModes::dec, L1RCT::digiInput(), MillePedeFileConverter_cfg::e, bsc_activity_cfg::ecal, EcalBarrel, L1RCT::ecalCompressedET(), ecalDigisToken_, EcalEndcap, L1RCT::ecalFineGrainBit(), ecalScaleToken_, emScaleToken_, options_cfi::eventSetup, EcalTPGScale::getTPGInGeV(), h, hltEgammaHLTExtra_cfi::hcal, L1RCT::hcalCompressedET(), hcalDigisToken_, CaloTPGTranscoder::hcaletValue(), L1RCT::hcalFineGrainBit(), hcalScaleToken_, hcalRecHitTable_cff::ieta, L1RCTLookupTables::lookup(), L1CaloEcalScale::nBinEta, L1CaloHcalScale::nBinEta, L1CaloEcalScale::nBinRank, L1CaloHcalScale::nBinRank, createIOVlist::nEvents, ofs, edm::ESHandle< T >::product(), alignCSCRings::r, rct, rctLookupTables, rctParametersToken_, alignCSCRings::s, L1CaloHcalScale::setBin(), L1CaloEcalScale::setBin(), L1RCTLookupTables::setChannelMask(), L1RCTLookupTables::setEcalScale(), L1RCTLookupTables::setHcalScale(), L1RCTLookupTables::setL1CaloEtScale(), L1RCTLookupTables::setRCTParameters(), tokens_, transcoderToken_, and useDebugTpgScales.

                                                                                    {
   edm::ESHandle<L1RCTParameters> rctParameters = eventSetup.getHandle(rctParametersToken_);
   const L1RCTParameters *r = rctParameters.product();
   edm::ESHandle<L1RCTChannelMask> channelMask = eventSetup.getHandle(channelMaskToken_);
   const L1RCTChannelMask *c = channelMask.product();
   edm::ESHandle<L1CaloEtScale> emScale = eventSetup.getHandle(emScaleToken_);
   const L1CaloEtScale *s = emScale.product();
 
   rctLookupTables->setRCTParameters(r);
   rctLookupTables->setChannelMask(c);
   rctLookupTables->setL1CaloEtScale(s);
 
   if (useDebugTpgScales) {
     // use old-style scales
     edm::ESHandle<CaloTPGTranscoder> transcoder = eventSetup.getHandle(transcoderToken_);
     const CaloTPGTranscoder *h_tpg = transcoder.product();
 
     EcalTPGScale e_tpg(tokens_, eventSetup);
 
     L1CaloEcalScale *ecalScale = new L1CaloEcalScale();
     L1CaloHcalScale *hcalScale = new L1CaloHcalScale();
 
     // create input scales, werner's code
     // ECAL
     std::cout << "ECAL Pos " << L1CaloEcalScale::nBinRank << std::endl;
     for (unsigned short ieta = 1; ieta <= L1CaloEcalScale::nBinEta; ++ieta) {
       for (unsigned short irank = 0; irank < L1CaloEcalScale::nBinRank; ++irank) {
         std::cout << ieta << " " << irank;
         EcalSubdetector subdet = (ieta <= 17) ? EcalBarrel : EcalEndcap;
         double etGeVPos = e_tpg.getTPGInGeV(irank,
                                             EcalTrigTowerDetId(1,  // +ve eta
                                                                subdet,
                                                                ieta,
                                                                1));  // dummy phi value
         ecalScale->setBin(irank, ieta, 1, etGeVPos);
         std::cout << etGeVPos << ", ";
       }
       std::cout << std::endl;
     }
     std::cout << std::endl;
 
     std::cout << "ECAL Neg" << std::endl;
     for (unsigned short ieta = 1; ieta <= L1CaloEcalScale::nBinEta; ++ieta) {
       for (unsigned short irank = 0; irank < L1CaloEcalScale::nBinRank; ++irank) {
         EcalSubdetector subdet = (ieta <= 17) ? EcalBarrel : EcalEndcap;
 
         std::cout << ieta << " " << irank;
         double etGeVNeg = e_tpg.getTPGInGeV(irank,
                                             EcalTrigTowerDetId(-1,  // -ve eta
                                                                subdet,
                                                                ieta,
                                                                2));  // dummy phi value
         ecalScale->setBin(irank, ieta, -1, etGeVNeg);
         std::cout << etGeVNeg << ", ";
       }
       std::cout << std::endl;
     }
     std::cout << std::endl;
 
     // HCAL
     std::cout << "HCAL" << std::endl;
     for (unsigned short ieta = 1; ieta <= L1CaloHcalScale::nBinEta; ++ieta) {
       for (unsigned short irank = 0; irank < L1CaloHcalScale::nBinRank; ++irank) {
         double etGeV = h_tpg->hcaletValue(ieta, irank);
 
         hcalScale->setBin(irank, ieta, 1, etGeV);
         hcalScale->setBin(irank, ieta, -1, etGeV);
         std::cout << etGeV << ", ";
         std::cout << std::endl;
       }
       std::cout << std::endl;
     }
 
     // set the input scales
     rctLookupTables->setEcalScale(ecalScale);
     rctLookupTables->setHcalScale(hcalScale);
 
   } else {
     edm::ESHandle<L1CaloHcalScale> hcalScale = eventSetup.getHandle(hcalScaleToken_);
     const L1CaloHcalScale *h = hcalScale.product();
     edm::ESHandle<L1CaloEcalScale> ecalScale = eventSetup.getHandle(ecalScaleToken_);
     const L1CaloEcalScale *e = ecalScale.product();
 
     rctLookupTables->setHcalScale(h);
     rctLookupTables->setEcalScale(e);
   }
 
   edm::Handle<EcalTrigPrimDigiCollection> ecal = event.getHandle(ecalDigisToken_);
   edm::Handle<HcalTrigPrimDigiCollection> hcal = event.getHandle(hcalDigisToken_);
   EcalTrigPrimDigiCollection ecalColl;
   HcalTrigPrimDigiCollection hcalColl;
   if (ecal.isValid()) {
     ecalColl = *ecal;
   }
   if (hcal.isValid()) {
     hcalColl = *hcal;
   }
   rct->digiInput(ecalColl, hcalColl);
   static int nEvents = 0;
   if (nEvents == 0) {
     ofs << "Crate = 0-17" << std::endl
         << "Card = 0-7 within the crate" << std::endl
         << "Tower = 0-31 covers 4 x 8 covered by the card" << std::endl
         << "EMAddr(0:8) = EMFGBit(0:0)+CompressedEMET(1:8)" << std::endl
         << "HDAddr(0:8) = HDFGBit(0:0)+CompressedHDET(1:8) - note: "
            "HDFGBit(0:0) is not part of the hardware LUT address"
         << std::endl
         << "LutOut(0:17)= "
            "LinearEMET(0:6)+HoEFGVetoBit(7:7)+LinearJetET(8:16)+ActivityBit(17:"
            "17)"
         << std::endl
         << "Event"
         << "\t"
         << "Crate"
         << "\t"
         << "Card"
         << "\t"
         << "Tower"
         << "\t"
         << "EMAddr"
         << "\t"
         << "HDAddr"
         << "\t"
         << "LUTOut" << std::endl;
   }
   if (nEvents < 64) {
     for (unsigned short iCrate = 0; iCrate < 18; iCrate++) {
       for (unsigned short iCard = 0; iCard < 7; iCard++) {
         // tower numbered from 0-31
         for (unsigned short iTower = 0; iTower < 32; iTower++) {
           unsigned short ecal = rct->ecalCompressedET(iCrate, iCard, iTower);
           unsigned short hcal = rct->hcalCompressedET(iCrate, iCard, iTower);
           unsigned short fgbit = rct->ecalFineGrainBit(iCrate, iCard, iTower);
           unsigned short mubit = rct->hcalFineGrainBit(iCrate, iCard, iTower);
           unsigned long lutOutput = rctLookupTables->lookup(ecal, hcal, fgbit, iCrate, iCard, iTower);
           ofs << std::hex << nEvents << "\t" << iCrate << "\t" << iCard << "\t" << iTower << "\t" << ecal * 2 + fgbit
               << "\t" << hcal * 2 + mubit << "\t" << lutOutput << std::dec << std::endl;
         }
       }
     }
   }
   nEvents++;
 }

Member Data Documentation

◆ channelMaskToken_

edm::ESGetToken<L1RCTChannelMask, L1RCTChannelMaskRcd> L1RCTSaveInput::channelMaskToken_

private

Definition at line 70 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ ecalDigisToken_

edm::EDGetTokenT<EcalTrigPrimDigiCollection> L1RCTSaveInput::ecalDigisToken_

private

Definition at line 67 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ ecalScaleToken_

edm::ESGetToken<L1CaloEcalScale, L1CaloEcalScaleRcd> L1RCTSaveInput::ecalScaleToken_

private

Definition at line 74 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ emScaleToken_

edm::ESGetToken<L1CaloEtScale, L1EmEtScaleRcd> L1RCTSaveInput::emScaleToken_

private

Definition at line 71 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ fileName

std::string L1RCTSaveInput::fileName

private

Definition at line 62 of file L1RCTSaveInput.h.

Referenced by L1RCTSaveInput().

◆ hcalDigisToken_

edm::EDGetTokenT<HcalTrigPrimDigiCollection> L1RCTSaveInput::hcalDigisToken_

private

Definition at line 68 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ hcalScaleToken_

edm::ESGetToken<L1CaloHcalScale, L1CaloHcalScaleRcd> L1RCTSaveInput::hcalScaleToken_

private

Definition at line 73 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ ofs

std::ofstream L1RCTSaveInput::ofs

private

Definition at line 76 of file L1RCTSaveInput.h.

Referenced by analyze(), and L1RCTSaveInput().

◆ rct

L1RCT* L1RCTSaveInput::rct

private

Definition at line 64 of file L1RCTSaveInput.h.

Referenced by analyze(), and ~L1RCTSaveInput().

◆ rctLookupTables

L1RCTLookupTables* L1RCTSaveInput::rctLookupTables

private

Definition at line 63 of file L1RCTSaveInput.h.

Referenced by analyze(), and ~L1RCTSaveInput().

◆ rctParametersToken_

edm::ESGetToken<L1RCTParameters, L1RCTParametersRcd> L1RCTSaveInput::rctParametersToken_

private

Definition at line 69 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ tokens_

EcalTPGScale::Tokens L1RCTSaveInput::tokens_

private

Definition at line 77 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ transcoderToken_

edm::ESGetToken<CaloTPGTranscoder, CaloTPGRecord> L1RCTSaveInput::transcoderToken_

private

Definition at line 72 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ useDebugTpgScales

bool L1RCTSaveInput::useDebugTpgScales

private

Definition at line 75 of file L1RCTSaveInput.h.

Referenced by analyze().

◆ useEcal

bool L1RCTSaveInput::useEcal

private

Definition at line 65 of file L1RCTSaveInput.h.

◆ useHcal

bool L1RCTSaveInput::useHcal

private

Definition at line 66 of file L1RCTSaveInput.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ L1RCTSaveInput()

◆ ~L1RCTSaveInput()

Member Function Documentation

◆ analyze()

Member Data Documentation

◆ channelMaskToken_

◆ ecalDigisToken_

◆ ecalScaleToken_

◆ emScaleToken_

◆ fileName

◆ hcalDigisToken_

◆ hcalScaleToken_

◆ ofs

◆ rct

◆ rctLookupTables

◆ rctParametersToken_

◆ tokens_

◆ transcoderToken_

◆ useDebugTpgScales

◆ useEcal

◆ useHcal