#include <HcalTBWriter.h>

Inheritance diagram for HcalTBWriter:

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

void	endJob () override

	HcalTBWriter (const edm::ParameterSet &pset)

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

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 Member Functions
void	buildTree (const FEDRawDataCollection &raw)

void	extractEventInfo (const FEDRawDataCollection &raw, const edm::EventID &id)

Private Attributes
std::map< int, std::string >	blockToName_

CDFChunk *	chunkList_ [1024]

std::map< int, int >	chunkMap_

CDFEventInfo *	eventInfo_

TFile *	file_

std::string	namePattern_

CDFRunInfo	ri_

edm::EDGetTokenT< FEDRawDataCollection >	tok_raw_

TTree *	tree_

int	trigChunk_

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	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

Writes HCAL-style ROOT files from the RawData block

Author: J. Mans - Minnesota

Definition at line 25 of file HcalTBWriter.h.

Constructor & Destructor Documentation

HcalTBWriter::HcalTBWriter ( const edm::ParameterSet & pset )

Definition at line 11 of file HcalTBWriter.cc.

References blockToName_, eventInfo_, file_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), dqmiolumiharvest::j, Skims_PA_cff::name, names, EgammaValidation_cff::num, AlCaHLTBitMon_QueryRunRegistry::string, tok_raw_, and tree_.

     : namePattern_(pset.getUntrackedParameter<std::string>("FilenamePattern", "/tmp/HTB_%06d.root")) {
   tok_raw_ = consumes<FEDRawDataCollection>(pset.getParameter<edm::InputTag>("fedRawDataCollectionTag"));
 
   std::vector<edm::ParameterSet> names = pset.getUntrackedParameter<std::vector<edm::ParameterSet> >("ChunkNames");
   std::vector<edm::ParameterSet>::iterator j;
   for (j = names.begin(); j != names.end(); j++) {
     std::string name = j->getUntrackedParameter<std::string>("Name");
     int num = j->getUntrackedParameter<int>("Number");
     blockToName_[num] = name;
   }
 
   file_ = nullptr;
   tree_ = nullptr;
   eventInfo_ = nullptr;
 }

Member Function Documentation

void HcalTBWriter::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	es
	)

override

Definition at line 46 of file HcalTBWriter.cc.

References CDFChunk::adoptBuffer(), buildTree(), HltBtagPostValidation_cff::c, chunkList_, chunkMap_, FEDRawData::data(), extractEventInfo(), FEDRawDataCollection::FEDData(), file_, alignmentValidation::fname, edm::Event::getByToken(), mps_fire::i, edm::EventBase::id(), namePattern_, CDFChunk::releaseBuffer(), ri_, edm::EventID::run(), CDFRunInfo::setInfo(), FEDRawData::size(), tok_raw_, and tree_.

                                                                      {
   edm::Handle<FEDRawDataCollection> raw;
   e.getByToken(tok_raw_, raw);
 
   if (file_ == nullptr) {
     char fname[4096];
     snprintf(fname, 4096, namePattern_.c_str(), e.id().run());
     edm::LogInfo("HCAL") << "Opening " << fname << " for writing HCAL-format file.";
     file_ = new TFile(fname, "RECREATE");
     ri_.setInfo("OriginalFile", fname);
     buildTree(*raw);
   }
 
   // adopt the buffers for writing
   for (std::map<int, int>::const_iterator i = chunkMap_.begin(); i != chunkMap_.end(); i++) {
     CDFChunk* c = chunkList_[i->second];
     const FEDRawData& frd = raw->FEDData(i->first);
     c->adoptBuffer((ULong64_t*)frd.data(), frd.size() / 8);
   }
 
   // copy the event info bits
   extractEventInfo(*raw, e.id());
 
   // fill the tree
   tree_->Fill();
   // release all the buffers
   for (std::map<int, int>::const_iterator i = chunkMap_.begin(); i != chunkMap_.end(); i++) {
     CDFChunk* c = chunkList_[i->second];
     c->releaseBuffer();
   }
 }

void HcalTBWriter::buildTree ( const FEDRawDataCollection & raw )

private

Definition at line 78 of file HcalTBWriter.cc.

References blockToName_, HltBtagPostValidation_cff::c, chunkList_, chunkMap_, eventInfo_, FEDRawDataCollection::FEDData(), mps_fire::i, dqmiolumiharvest::j, Skims_PA_cff::name, FEDRawData::size(), AlCaHLTBitMon_QueryRunRegistry::string, tree_, and trigChunk_.

Referenced by analyze().

                                                             {
   tree_ = new TTree("CMSRAW", "CMS Common Data Format Tree");
   chunkMap_.clear();
   trigChunk_ = -1;
   int j = 0;
   for (int i = 0; i < 2048; i++) {
     const FEDRawData& frd = raw.FEDData(i);
     if (frd.size() < 16)
       continue;  // it's empty... like
 
     std::string name;
     if (blockToName_.find(i) != blockToName_.end())
       name = blockToName_[i];
     else {
       char sname[64];
       snprintf(sname, 64, "Chunk%03d", i);
       name = sname;
     }
 
     CDFChunk* c = new CDFChunk(name.c_str());
     chunkList_[j] = c;
     tree_->Branch(name.c_str(), "CDFChunk", &(chunkList_[j]));
     chunkMap_[i] = j;
 
     if (name == "HCAL_Trigger" || name == "SliceTest_Trigger")
       trigChunk_ = j;
 
     j++;
   }
   eventInfo_ = new CDFEventInfo();
   tree_->Branch("CDFEventInfo", "CDFEventInfo", &eventInfo_, 16000, 2);
 }

void HcalTBWriter::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 28 of file HcalTBWriter.cc.

References edmScanValgrind::buffer, chunkMap_, eventInfo_, file_, ri_, CDFRunInfo::setInfo(), CDFRunInfo::store(), and tree_.

                           {
   char buffer[1024];
   if (file_ != nullptr) {
     file_->Write();
 
     ri_.setInfo("DAQSofwareRelease", "UNKNOWN -- HcalTBWriter");
     gethostname(buffer, 1024);
     ri_.setInfo("WriterHostname", buffer);
     ri_.store(file_);
 
     file_->Close();
     file_ = nullptr;
     tree_ = nullptr;
     chunkMap_.clear();
     eventInfo_ = nullptr;
   }
 }

void HcalTBWriter::extractEventInfo	(	const FEDRawDataCollection &	raw,
		const edm::EventID &	id
	)

private

Definition at line 135 of file HcalTBWriter.cc.

References FEDHeader::bxID(), chunkList_, DEFINE_FWK_MODULE, eventInfo_, printContent_cfi::getData, FEDHeader::lvl1ID(), StandardTrgMsgBlkStruct::orbitNumber, visualization-live-secondInstance_cfg::runno, newExtendedTrgMsgBlkStruct::runNumberSequenceId, CDFEventInfo::Set(), newExtendedTrgMsgBlkStruct::stdBlock, and trigChunk_.

Referenced by analyze().

                                                                                          {
   int runno = id.run();
   const char* seqid = "";
   int eventNo = id.event();
   int l1aNo = eventNo;
   int orbitNo = 0;
   int bunchNo = 0;
 
   if (trigChunk_ >= 0) {
     const newExtendedTrgMsgBlk* tinfo = (const newExtendedTrgMsgBlk*)(chunkList_[trigChunk_]->getData() +
                                                                       2);  // assume 2 64-bit words for the CDF header
     orbitNo = tinfo->stdBlock.orbitNumber;
     seqid = tinfo->runNumberSequenceId;
     FEDHeader head((const unsigned char*)chunkList_[trigChunk_]->getData());
     bunchNo = head.bxID();
     l1aNo = head.lvl1ID();
   }
 
   eventInfo_->Set(runno, seqid, eventNo, l1aNo, orbitNo, bunchNo);
 }