#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

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)

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 ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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 13 of file HcalTBWriter.cc.

References blockToName_, eventInfo_, file_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), dataset::name, names, pileupDistInMC::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 50 of file HcalTBWriter.cc.

References CDFChunk::adoptBuffer(), buildTree(), EnergyCorrector::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 83 of file HcalTBWriter.cc.

References blockToName_, EnergyCorrector::c, chunkList_, chunkMap_, eventInfo_, FEDRawDataCollection::FEDData(), mps_fire::i, dataset::name, FEDRawData::size(), lumiContext::sname, 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 32 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 137 of file HcalTBWriter.cc.

References FEDHeader::bxID(), chunkList_, DEFINE_FWK_MODULE, eventInfo_, AlCaHLTBitMon_QueryRunRegistry::getData(), FEDHeader::lvl1ID(), StandardTrgMsgBlkStruct::orbitNumber, lumiregperbunch::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); 
 }