#include <CalibCalorimetry/HcalLuttoDB/src/HcalLuttoDB.cc>

Inheritance diagram for HcalLuttoDB:

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

void	beginJob () override

void	endJob () override

	HcalLuttoDB (const edm::ParameterSet &)

	~HcalLuttoDB () 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 Member Functions
std::vector< unsigned char >	extractOutputLut (const CaloTPGTranscoder &coder, HcalTrigTowerDetId chan)

std::ostream *	openChecksums ()

std::ostream *	openPerCrate (int crate)

std::ostream *	openPerLut1 (HcalElectronicsId eid)

std::ostream *	openPerLut2 (HcalElectronicsId eid)

void	writeoutlut1 (HcalDetId id, HcalElectronicsId eid, const std::vector< unsigned short > &lut, std::ostream &os)

void	writeoutlut2 (HcalTrigTowerDetId id, HcalElectronicsId eid, const std::vector< unsigned char > &lut, std::ostream &os)

Private Attributes
std::string	creationstamp_

const std::string	creationtag_

const std::string	fileformat_

const bool	filePerCrate_

std::ostream *	oc_

const std::string	targetfirmware_

edm::ESGetToken< CaloTPGTranscoder, CaloTPGRecord >	tokcaloCode_

edm::ESGetToken< HcalDbService, HcalDbRecord >	tokdb_

edm::ESGetToken< HcalTPGCoder, HcalTPGRecord >	tokhcalCode_

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

Description: <one line="" class="" summary>=""> R Implementation: <Notes on="" implementation>="">

Definition at line 54 of file HcalLuttoDB.cc.

Constructor & Destructor Documentation

◆ HcalLuttoDB()

HcalLuttoDB::HcalLuttoDB ( const edm::ParameterSet & iConfig )

explicit

Definition at line 90 of file HcalLuttoDB.cc.

References tokcaloCode_, tokdb_, and tokhcalCode_.

     : creationtag_(iConfig.getParameter<std::string>("creationtag")),
       targetfirmware_(iConfig.getParameter<std::string>("targetfirmware")),
       filePerCrate_(iConfig.getUntrackedParameter<bool>("filePerCrate", true)),
       fileformat_(iConfig.getParameter<std::string>("filePrefix")) {
   tokdb_ = esConsumes<HcalDbService, HcalDbRecord>();
   tokhcalCode_ = esConsumes<HcalTPGCoder, HcalTPGRecord>();
   tokcaloCode_ = esConsumes<CaloTPGTranscoder, CaloTPGRecord>();
 }

◆ ~HcalLuttoDB()

HcalLuttoDB::~HcalLuttoDB ( )

override

Definition at line 100 of file HcalLuttoDB.cc.

100 {}

Member Function Documentation

◆ analyze()

void HcalLuttoDB::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 270 of file HcalLuttoDB.cc.

References HcalElectronicsMap::allElectronicsId(), CaloTPGTranscoder_cfi::CaloTPGTranscoder, extractOutputLut(), filePerCrate_, HcalGenericDetId::genericSubdet(), edm::EventSetup::getData(), HcalDbService::getHcalMapping(), HcalTPGCoder::getLinearizationLUT(), HcalGenericDetId::HcalGenBarrel, HcalGenericDetId::HcalGenEndcap, HcalGenericDetId::HcalGenForward, HcalElectronicsMap::lookup(), HcalElectronicsMap::lookupTrigger(), DetId::null(), oc_, openChecksums(), openPerCrate(), openPerLut1(), openPerLut2(), tokcaloCode_, tokdb_, tokhcalCode_, writeoutlut1(), and writeoutlut2().

                                                                              {
   //using namespace edm;
   //using namespace std;
 
   edm::LogInfo("Hcal") << "Beginning dump of Hcal TPG LUTS (this may take a minute or two)";
 
   const HcalElectronicsMap* Map_;
   const HcalDbService* pSetup = &iSetup.getData(tokdb_);
   Map_ = pSetup->getHcalMapping();
   // get the conditions, for the decoding
   const HcalTPGCoder* inputCoder = &iSetup.getData(tokhcalCode_);
   const CaloTPGTranscoder* outTranscoder = &iSetup.getData(tokcaloCode_);
 
   std::vector<HcalElectronicsId> allEID = Map_->allElectronicsId();
   std::vector<HcalElectronicsId>::iterator itreid;
 
   std::ostream* pfile = nullptr;
   oc_ = openChecksums();
 
   for (int crate = 0; crate < 20; crate++) {
     edm::LogInfo("Hcal") << "Beginning crate " << crate;
     for (itreid = allEID.begin(); itreid != allEID.end(); ++itreid) {
       if (itreid->readoutVMECrateId() != crate)
         continue;
       if (itreid->isTriggerChainId()) {  // lut2
         HcalTrigTowerDetId tid = Map_->lookupTrigger(*itreid);
         if (tid.null())
           continue;
 
         if (filePerCrate_ && pfile == nullptr)
           pfile = openPerCrate(crate);
         else if (pfile == nullptr)
           pfile = openPerLut2(*itreid);
 
         std::vector<unsigned char> lut = extractOutputLut(*outTranscoder, tid);
         writeoutlut2(tid, *itreid, lut, *pfile);
         if (!filePerCrate_) {
           delete pfile;
           pfile = nullptr;
         }
       } else {  // lut1
         HcalGenericDetId gid = Map_->lookup(*itreid);
         if (gid.null() || !(gid.genericSubdet() == HcalGenericDetId::HcalGenBarrel ||
                             gid.genericSubdet() == HcalGenericDetId::HcalGenEndcap ||
                             gid.genericSubdet() == HcalGenericDetId::HcalGenForward))
           continue;
 
         if (filePerCrate_ && pfile == nullptr)
           pfile = openPerCrate(crate);
         else if (pfile == nullptr)
           pfile = openPerLut1(*itreid);
 
         std::vector<unsigned short> lut = inputCoder->getLinearizationLUT(HcalDetId(gid));
         writeoutlut1(HcalDetId(gid), *itreid, lut, *pfile);
         if (!filePerCrate_) {
           delete pfile;
           pfile = nullptr;
         }
       }
     }
     if (pfile != nullptr) {
       if (filePerCrate_)
         *pfile << "</CFGBrickSet>\n";
       delete pfile;
       pfile = nullptr;
     }
   }
   *oc_ << "</CFGBrick>\n";
   delete oc_;
 }

◆ beginJob()

void HcalLuttoDB::beginJob ( )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 342 of file HcalLuttoDB.cc.

References edmScanValgrind::buffer, creationstamp_, submitPVValidationJobs::now, and hcalRecHitTable_cff::time.

                            {
   char buffer[120];
   time_t now = time(nullptr);
   struct tm* tm = localtime(&now);
   strftime(buffer, 120, "%F %T", tm);
   creationstamp_ = buffer;
 }

◆ endJob()

void HcalLuttoDB::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 351 of file HcalLuttoDB.cc.

Referenced by o2olib.O2ORunMgr::executeJob().

351 {}

◆ extractOutputLut()

std::vector< unsigned char > HcalLuttoDB::extractOutputLut	(	const CaloTPGTranscoder &	coder,
		HcalTrigTowerDetId	chan
	)

private

Definition at line 260 of file HcalLuttoDB.cc.

References officialStyle::chan, mps_fire::i, and alignCSCRings::s.

Referenced by analyze().

                                                                                                               {
   std::vector<unsigned char> lut;
   for (int i = 0; i < 1024; i++) {
     HcalTriggerPrimitiveSample s = coder.hcalCompress(chan, i, false);
     lut.push_back(s.compressedEt());
   }
   return lut;
 }

◆ openChecksums()

std::ostream * HcalLuttoDB::openChecksums ( )

private

Definition at line 106 of file HcalLuttoDB.cc.

References fileformat_, and alignmentValidation::fname.

Referenced by analyze().

                                        {
   char fname[1024];
   snprintf(fname, 1024, "%s_checksums.xml", fileformat_.c_str());
   std::ostream* os = new std::ofstream(fname);
   (*os) << "<?xml version=\"1.0\"?>\n<CFGBrick>\n";
   return os;
 }

◆ openPerCrate()

std::ostream * HcalLuttoDB::openPerCrate ( int crate )

private

Definition at line 114 of file HcalLuttoDB.cc.

References fileformat_, and alignmentValidation::fname.

Referenced by analyze().

                                                {
   char fname[1024];
   snprintf(fname, 1024, "%s_%d.xml", fileformat_.c_str(), crate);
   std::ostream* os = new std::ofstream(fname);
   (*os) << "<?xml version=\"1.0\"?>\n<CFGBrickSet>\n";
   return os;
 }

◆ openPerLut1()

std::ostream * HcalLuttoDB::openPerLut1 ( HcalElectronicsId eid )

private

Definition at line 122 of file HcalLuttoDB.cc.

References runTauDisplay::eid, fileformat_, and alignmentValidation::fname.

Referenced by analyze().

                                                           {
   char fname[1024];
   snprintf(fname,
            1024,
            "%s_%d_%d%c_%d_%d_1.xml",
            fileformat_.c_str(),
            eid.readoutVMECrateId(),
            eid.htrSlot(),
            ((eid.htrTopBottom()) ? ('t') : ('b')),
            eid.fiberIndex(),
            eid.fiberChanId());
   std::ostream* os = new std::ofstream(fname);
   (*os) << "<?xml version=\"1.0\"?>\n";
   return os;
 }

◆ openPerLut2()

std::ostream * HcalLuttoDB::openPerLut2 ( HcalElectronicsId eid )

private

Definition at line 138 of file HcalLuttoDB.cc.

References runTauDisplay::eid, fileformat_, and alignmentValidation::fname.

Referenced by analyze().

                                                           {
   char fname[1024];
   snprintf(fname,
            1024,
            "%s_%d_%d%c_%d_%d_2.xml",
            fileformat_.c_str(),
            eid.readoutVMECrateId(),
            eid.htrSlot(),
            ((eid.htrTopBottom()) ? ('t') : ('b')),
            eid.slbSiteNumber(),
            eid.slbChannelIndex());
   std::ostream* os = new std::ofstream(fname);
   (*os) << "<?xml version=\"1.0\"?>\n";
   return os;
 }

◆ writeoutlut1()

void HcalLuttoDB::writeoutlut1	(	HcalDetId	id,
		HcalElectronicsId	eid,
		const std::vector< unsigned short > &	lut,
		std::ostream &	os
	)

private

Definition at line 154 of file HcalLuttoDB.cc.

References creationstamp_, creationtag_, TauDecayModes::dec, runTauDisplay::eid, formatRevision_, HcalForward, mps_fire::i, oc_, and targetfirmware_.

Referenced by analyze().

                                                {
   os << "<CFGBrick> " << std::endl;
   os << " <Parameter name='IETA' type='int'>" << id.ieta() << "</Parameter>" << std::endl;
   os << " <Parameter name='IPHI' type='int'>" << id.iphi() << "</Parameter>" << std::endl;
   os << " <Parameter name='DEPTH' type='int'>" << id.depth() << "</Parameter>" << std::endl;
   os << " <Parameter name='CRATE' type='int'>" << eid.readoutVMECrateId() << "</Parameter>" << std::endl;
   os << " <Parameter name='SLOT' type='int'>" << eid.htrSlot() << "</Parameter>" << std::endl;
   os << " <Parameter name='TOPBOTTOM' type='int'>" << eid.htrTopBottom() << "</Parameter>" << std::endl;
   os << " <Parameter name='FIBER' type='int'>" << eid.fiberIndex() << "</Parameter>" << std::endl;
   os << " <Parameter name='FIBERCHAN' type='int'>" << eid.fiberChanId() << "</Parameter>" << std::endl;
   os << " <Parameter name='LUT_TYPE' type='int'>1</Parameter>" << std::endl;
   os << " <Parameter name='CREATIONTAG' type='string'>" << creationtag_ << "</Parameter>" << std::endl;
   os << " <Parameter name='CREATIONSTAMP' type='string'>" << creationstamp_ << "</Parameter>" << std::endl;
   os << " <Parameter name='FORMATREVISION' type='string'>" << formatRevision_ << "</Parameter>" << std::endl;
   os << " <Parameter name='TARGETFIRMWARE' type='string'>" << targetfirmware_ << "</Parameter>" << std::endl;
   int generalizedIndex = id.ietaAbs() + 1000 * id.depth() + 10000 * id.iphi() + ((id.ieta() < 0) ? (0) : (100)) +
                          ((id.subdet() == HcalForward && id.ietaAbs() == 29) ? (4 * 10000) : (0));
 
   os << " <Parameter name='GENERALIZEDINDEX' type='int'>" << generalizedIndex << "</Parameter>" << std::endl;
   // do checksum
   md5_state_t md5er;
   md5_byte_t digest[16];
   unsigned char tool[2];
   md5_init(&md5er);
   for (int i = 0; i < 128; i++) {
     tool[0] = lut[i] & 0xFF;
     tool[1] = (lut[i] >> 8) & 0xFF;
     md5_append(&md5er, tool, 2);
   }
   md5_finish(&md5er, digest);
   os << " <Parameter name='CHECKSUM' type='string'>";
   for (int i = 0; i < 16; i++)
     os << std::hex << (((int)(digest[i])) & 0xFF);
   os << "</Parameter>\n";
 
   *oc_ << "  <Data crate='" << eid.readoutVMECrateId() << "' slot='" << eid.htrSlot() << "' fpga='"
        << eid.htrTopBottom() << "' fiber='" << eid.fiberIndex() << "' fiberchan='" << eid.fiberChanId()
        << "' luttype='1' elements='1' encoding='hex'>";
   for (int i = 0; i < 16; i++)
     *oc_ << std::hex << (((int)(digest[i])) & 0xFF);
   *oc_ << "</Data>\n";
 
   os << " <Data elements='128' encoding='hex'> " << std::endl;
   os << std::hex;
   for (int initr2 = 0; initr2 < 128; initr2++) {
     os << lut[initr2] << " ";
   }
   os << std::dec;
   os << std::endl;
   os << " </Data> " << std::endl;
   os << "</CFGBrick> " << std::endl;
 }

◆ writeoutlut2()

void HcalLuttoDB::writeoutlut2	(	HcalTrigTowerDetId	id,
		HcalElectronicsId	eid,
		const std::vector< unsigned char > &	lut,
		std::ostream &	os
	)

private

Definition at line 210 of file HcalLuttoDB.cc.

References creationstamp_, creationtag_, TauDecayModes::dec, runTauDisplay::eid, formatRevision_, mps_fire::i, createfilelist::int, oc_, and targetfirmware_.

Referenced by analyze().

                                                {
   os << "<CFGBrick> " << std::endl;
   os << " <Parameter name='IETA' type='int'>" << id.ieta() << "</Parameter>" << std::endl;
   os << " <Parameter name='IPHI' type='int'>" << id.iphi() << "</Parameter>" << std::endl;
   os << " <Parameter name='CRATE' type='int'>" << eid.readoutVMECrateId() << "</Parameter>" << std::endl;
   os << " <Parameter name='SLOT' type='int'>" << eid.htrSlot() << "</Parameter>" << std::endl;
   os << " <Parameter name='TOPBOTTOM' type='int'>" << eid.htrTopBottom() << "</Parameter>" << std::endl;
   os << " <Parameter name='SLB' type='int'>" << eid.slbSiteNumber() << "</Parameter>" << std::endl;
   os << " <Parameter name='SLBCHAN' type='int'>" << eid.slbChannelIndex() << "</Parameter>" << std::endl;
   os << " <Parameter name='LUT_TYPE' type='int'>2</Parameter>" << std::endl;
   os << " <Parameter name='CREATIONTAG' type='string'>" << creationtag_ << "</Parameter>" << std::endl;
   os << " <Parameter name='CREATIONSTAMP' type='string'>" << creationstamp_ << "</Parameter>" << std::endl;
   os << " <Parameter name='FORMATREVISION' type='string'>" << formatRevision_ << "</Parameter>" << std::endl;
   os << " <Parameter name='TARGETFIRMWARE' type='string'>" << targetfirmware_ << "</Parameter>" << std::endl;
   int generalizedIndex = id.ietaAbs() + 10000 * id.iphi() + ((id.ieta() < 0) ? (0) : (100));
 
   os << " <Parameter name='GENERALIZEDINDEX' type='int'>" << generalizedIndex << "</Parameter>" << std::endl;
 
   // do checksum
   md5_state_t md5er;
   md5_byte_t digest[16];
   md5_init(&md5er);
   md5_append(&md5er, &(lut[0]), 1024);
   md5_finish(&md5er, digest);
   os << " <Parameter name='CHECKSUM' type='string'>";
   for (int i = 0; i < 16; i++)
     os << std::hex << (((int)(digest[i])) & 0xFF);
   os << "</Parameter>\n";
 
   *oc_ << "  <Data crate='" << eid.readoutVMECrateId() << "' slot='" << eid.htrSlot() << "' fpga='"
        << eid.htrTopBottom() << "' slb='" << eid.slbSiteNumber() << "' slbchan='" << eid.slbChannelIndex()
        << "' luttype='2' elements='1' encoding='hex'>";
   for (int i = 0; i < 16; i++)
     *oc_ << std::hex << (((int)(digest[i])) & 0xFF);
   *oc_ << "</Data>\n";
 
   os << " <Data elements='1024' encoding='hex'> " << std::endl;
   os << std::hex;
   for (int initr2 = 0; initr2 < 1024; initr2++) {
     os << (int(lut[initr2]) & 0xFF) << " ";
   }
   os << std::dec;
   os << std::endl;
   os << " </Data> " << std::endl;
   os << "</CFGBrick> " << std::endl;
 }