#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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

edm::InputTag	ecalDigisLabel

edm::ESGetToken < L1CaloEcalScale, L1CaloEcalScaleRcd >	ecalScaleToken_

edm::ESGetToken< L1CaloEtScale, L1EmEtScaleRcd >	emScaleToken_

std::string	fileName

edm::InputTag	hcalDigisLabel

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

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)

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 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<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 53 of file L1RCTSaveInput.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 27 of file L1RCTSaveInput.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, beamvalidation::exit(), 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")),
       ecalDigisLabel(conf.getParameter<edm::InputTag>("ecalDigisLabel")),
       hcalDigisLabel(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) {
     std::cerr << "Could not create " << fileName << std::endl;
     exit(1);
   }
 }

L1RCTSaveInput::~L1RCTSaveInput ( )

override

Definition at line 50 of file L1RCTSaveInput.cc.

References rct, and rctLookupTables.

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

Member Function Documentation

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 57 of file L1RCTSaveInput.cc.

References c, channelMaskToken_, gather_cfg::cout, TauDecayModes::dec, L1RCT::digiInput(), alignCSCRings::e, digitizers_cfi::ecal, EcalBarrel, L1RCT::ecalCompressedET(), ecalDigisLabel, EcalEndcap, L1RCT::ecalFineGrainBit(), ecalScaleToken_, emScaleToken_, edm::EventSetup::getHandle(), EcalTPGScale::getTPGInGeV(), h, digitizers_cfi::hcal, L1RCT::hcalCompressedET(), hcalDigisLabel, CaloTPGTranscoder::hcaletValue(), L1RCT::hcalFineGrainBit(), hcalScaleToken_, edm::HandleBase::isValid(), L1RCTLookupTables::lookup(), L1CaloHcalScale::nBinEta, L1CaloEcalScale::nBinEta, L1CaloEcalScale::nBinRank, L1CaloHcalScale::nBinRank, submitPVValidationJobs::nEvents, ofs, edm::ESHandle< class >::product(), alignCSCRings::r, rct, rctLookupTables, rctParametersToken_, alignCSCRings::s, L1CaloEcalScale::setBin(), L1CaloHcalScale::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;
   edm::Handle<HcalTrigPrimDigiCollection> hcal;
   event.getByLabel(ecalDigisLabel, ecal);
   event.getByLabel(hcalDigisLabel, hcal);
   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++;
 }