#include <L1GtPatternGenerator.h>

Inheritance diagram for L1GtPatternGenerator:

Public Member Functions
	L1GtPatternGenerator (const edm::ParameterSet &)

	~L1GtPatternGenerator ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
void	extractGlobalTriggerData (const edm::Event &iEvent, L1GtPatternMap &patterns)

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

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
	analyze each event More...

virtual void	beginJob ()
	analyze More...

virtual void	endJob ()
	end of job More...

void	packHfRecords (const std::string &resultName, L1GtPatternMap &allPatterns)

Static Private Member Functions
static uint32_t	packEtMiss (uint32_t rawValue)

static uint32_t	packRegionalMuons (uint32_t rawValue)

Private Attributes
std::vector< int >	m_bx

std::vector< uint32_t >	m_columnDefaults

std::vector< uint32_t >	m_columnLengths

std::vector< std::string >	m_columnNames

edm::InputTag	m_cscTag

bool	m_debug

std::string	m_destPath
	an algorithm and a condition in that algorithm to test the object maps More...

edm::InputTag	m_dtTag
	input tags for regional muon data More...

std::string	m_fileName
	output file name More...

std::ofstream	m_fileStream

std::string	m_footer

edm::InputTag	m_gctTag
	input tag for GCT data More...

edm::InputTag	m_gmtTag
	input tag for GMT data More...

edm::InputTag	m_gtTag
	input tag for GT data More...

std::string	m_header
	formatting instructions More...

edm::InputTag	m_rpcbTag

edm::InputTag	m_rpcfTag

std::auto_ptr< L1GtPatternWriter >	m_writer

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

Detailed Description

Description: A generator of pattern files for L1 GT hardware testing.

Implementation: <TODO: enter implementation details>

Author: : Thomas Themel - HEPHY Vienna

Description: see header file.

Implementation: <TODO: enter implementation details>

Author: : Thomas Themel - HEPHY Vienna

Definition at line 37 of file L1GtPatternGenerator.h.

Constructor & Destructor Documentation

L1GtPatternGenerator::L1GtPatternGenerator ( const edm::ParameterSet & parSet )

explicit

Definition at line 37 of file L1GtPatternGenerator.cc.

References edm::ParameterSet::getParameter(), LogDebug, m_bx, m_columnDefaults, m_columnLengths, m_columnNames, m_cscTag, m_debug, m_dtTag, m_fileName, m_footer, m_gctTag, m_gmtTag, m_gtTag, m_header, m_rpcbTag, m_rpcfTag, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
     // input tags for trigger records
     m_gctTag = parSet.getParameter<edm::InputTag>("GctInputTag");
     m_gmtTag = parSet.getParameter<edm::InputTag>("GmtInputTag");
     m_gtTag = parSet.getParameter<edm::InputTag>("GtInputTag");
     m_dtTag = parSet.getParameter<edm::InputTag>("DtInputTag");
     m_cscTag = parSet.getParameter<edm::InputTag>("CscInputTag");
     m_rpcbTag = parSet.getParameter<edm::InputTag>("RpcbInputTag");
     m_rpcfTag = parSet.getParameter<edm::InputTag>("RpcfInputTag");
 
     // output formatting stuff
     m_header = parSet.getParameter<std::string>("PatternFileHeader");
     m_footer = parSet.getParameter<std::string>("PatternFileFooter");
     m_columnNames = parSet.getParameter<std::vector<std::string> >("PatternFileColumns");
     m_columnLengths = parSet.getParameter<std::vector<uint32_t> >("PatternFileLengths");
     m_columnDefaults = parSet.getParameter<std::vector<uint32_t> >("PatternFileDefaultValues");
     m_fileName = parSet.getParameter<std::string>("PatternFileName");
     m_bx = parSet.getParameter<std::vector<int> >("bx");
     m_debug = parSet.getParameter<bool>("DebugOutput");
 
 
     if(m_columnLengths.size() != m_columnNames.size()) { 
       edm::LogWarning("L1GtPatternGenerator") 
     << "Length of PatternFileColumns does not match length of PatternFileLenghts, " <<
     m_columnNames.size() << " vs " << m_columnLengths.size() << std::endl;
     }
 
     LogDebug("L1GtPatternGenerator")
       << "\nL1 GCT  record:            "
       << m_gctTag 
       << "\nL1 GMT record:             "
       << m_gmtTag
       << "\nL1 GT record:              "
       << m_gtTag << std::endl;
 }

L1GtPatternGenerator::~L1GtPatternGenerator ( )

Definition at line 75 of file L1GtPatternGenerator.cc.

76 {}

Member Function Documentation

void L1GtPatternGenerator::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	evSetup
	)

privatevirtual

analyze each event

Analyze each event:

Extract the input records that interest us from the event
Split them into pattern file lines according to their bx number
Format the lines and write them to the file.

Implements edm::EDAnalyzer.

Definition at line 230 of file L1GtPatternGenerator.cc.

References L1GctEtMiss::et(), edm::EventID::event(), extractGlobalTriggerData(), extractRecordData(), L1MuGMTCand::getDataWord(), L1MuRegionalCand::getDataWord(), edm::EventBase::id(), edm::InputTag::instance(), edm::InputTag::label(), LogTrace, edm::EventBase::luminosityBlock(), m_cscTag, m_dtTag, m_gctTag, m_gmtTag, m_rpcbTag, m_rpcfTag, m_writer, packHfRecords(), packRegionalMuons(), L1GctEtMiss::phi(), L1GctEtHad::raw(), L1GctEtTotal::raw(), L1GctJetCand::raw(), L1GctHtMiss::raw(), L1GctEmCand::raw(), L1GctHFRingEtSums::raw(), L1GctHFBitCounts::raw(), edm::Event::run(), and convertSQLiteXML::runNumber.

 {
     // debug information
     const unsigned int runNumber = iEvent.run();
     const unsigned int lsNumber = iEvent.luminosityBlock();
     const unsigned int eventNumber = iEvent.id().event();
 
     LogTrace("L1GtPatternGenerator") << "\n\nL1GtPatternGenerator::analyze: Run: " << runNumber << " LS: " << lsNumber << " Event: "
             << eventNumber << "\n\n" << std::endl;
     
 
   L1GtPatternMap allPatterns;
 
   // GMT muon candidates
   extractRecordData(iEvent, allPatterns, m_gmtTag.label(), m_gmtTag.instance(), &L1MuGMTCand::getDataWord, "gmtMuon");
 
   // regional muon candidates
   extractRecordData(iEvent, allPatterns, m_cscTag.label(), m_cscTag.instance(), &L1MuRegionalCand::getDataWord, "cscMuon", packRegionalMuons);  
   extractRecordData(iEvent, allPatterns, m_dtTag.label(),  m_dtTag.instance(), &L1MuRegionalCand::getDataWord, "dtMuon", packRegionalMuons);
   extractRecordData(iEvent, allPatterns, m_rpcfTag.label(), m_rpcfTag.instance(), &L1MuRegionalCand::getDataWord, "fwdMuon", packRegionalMuons);
   extractRecordData(iEvent, allPatterns, m_rpcbTag.label(), m_rpcbTag.instance(), &L1MuRegionalCand::getDataWord, "brlMuon", packRegionalMuons);
 
   // GCT objects
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "nonIsoEm", &L1GctEmCand::raw, "gctEm");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "isoEm", &L1GctEmCand::raw, "gctIsoEm");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "", &L1GctEtMiss::et, "etMiss");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "", &L1GctEtMiss::phi, "etMissPhi");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "", &L1GctHtMiss::raw, "htMiss");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "", &L1GctEtHad::raw, "etHad");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "", &L1GctEtTotal::raw, "etTotal");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "cenJets", &L1GctJetCand::raw, "cenJet");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "forJets", &L1GctJetCand::raw, "forJet");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "tauJets", &L1GctJetCand::raw, "tauJet");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "", &L1GctHFBitCounts::raw, "hfBitCounts");
   extractRecordData(iEvent, allPatterns, m_gctTag.label(), "", &L1GctHFRingEtSums::raw, "hfRingEtSums");
 
   // Post processing: 
   // HFBitCounts/HFRingEtSums need to be mangled to PSB values
   packHfRecords("hfPsbValue", allPatterns);
 
   // GT objects
   extractGlobalTriggerData(iEvent, allPatterns);
 
   // Output
   m_writer->writePatterns(allPatterns);
 }

void L1GtPatternGenerator::beginJob ( void )

privatevirtual

analyze

Method called once each job just before starting event loop.

Initialize the output file and the writer object.

Reimplemented from edm::EDAnalyzer.

Definition at line 280 of file L1GtPatternGenerator.cc.

References m_bx, m_columnDefaults, m_columnLengths, m_columnNames, m_debug, m_fileName, m_fileStream, m_footer, m_header, and m_writer.

 {
   m_fileStream.open(m_fileName.c_str());
 
   if(!m_fileStream) { 
     edm::LogError("L1GtPatternGenerator") <<  "Failed to open output file " << m_fileName;
   }
 
   m_writer = std::auto_ptr<L1GtPatternWriter>(new L1GtPatternWriter(m_fileStream, m_header, m_footer, m_columnNames, m_columnLengths, m_columnDefaults, m_bx, m_debug));
 }

void L1GtPatternGenerator::endJob ( void )

privatevirtual

end of job

Method called once each job just after ending the event loop.

Close the output file stream.

Reimplemented from edm::EDAnalyzer.

Definition at line 294 of file L1GtPatternGenerator.cc.

References m_fileStream, and m_writer.

 {
   m_writer->close();
   m_fileStream.close();
 }

void L1GtPatternGenerator::extractGlobalTriggerData	(	const edm::Event &	iEvent,
		L1GtPatternMap &	patterns
	)

protected

Definition at line 167 of file L1GtPatternGenerator.cc.

References edm::EventID::event(), Exception, extractGlobalTriggerWord(), edm::Event::getByLabel(), L1GtPatternMap::getLine(), patZpeak::handle, edm::EventBase::id(), edm::InputTag::instance(), edm::HandleBase::isValid(), edm::InputTag::label(), geometryCSVtoXML::line, m_gtTag, and L1GtPatternLine::push().

Referenced by analyze().

                                                                                                     {
 
   // extract global trigger readout record
   edm::Handle<L1GlobalTriggerReadoutRecord> handle;
   iEvent.getByLabel(m_gtTag, handle);
 
   // continue if it's present
   if(!handle.isValid()) { 
     throw cms::Exception(__func__) << "Failed to extract GT readout record labeled " 
                    << m_gtTag.label() << ", instance " << m_gtTag.instance();
   }
 
   edm::EventNumber_t eventNr = iEvent.id().event();
 
   // for each FDL word...
   const std::vector<L1GtFdlWord>& fdlWords = handle->gtFdlVector();
   for(std::vector<L1GtFdlWord>::const_iterator it = fdlWords.begin();
       it != fdlWords.end() ; ++it) { 
     // extract relevant data
     int bx = it->bxInEvent();
 
     // find matching pattern file line
     L1GtPatternLine& line = patterns.getLine(eventNr, bx);
 
     extractGlobalTriggerWord(it->gtDecisionWord(), line, "gtDecision");
     extractGlobalTriggerWord(it->gtDecisionWordExtended(), line, "gtDecisionExt");
     extractGlobalTriggerWord(it->gtTechnicalTriggerWord(), line, "gtTechTrigger");
 
     line.push("gtFinalOr", it->finalOR());
   }
 }

static uint32_t L1GtPatternGenerator::packEtMiss ( uint32_t rawValue )

staticprivate

Post-processing for etMissing

void L1GtPatternGenerator::packHfRecords	(	const std::string &	resultName,
		L1GtPatternMap &	allPatterns
	)

private

Post-processing for complex mapping of HF records to PSB values

The mapping from hfBitCounts/hfRingEtSums raw data to the PSBs is non-trivial, see http://wwwhephy.oeaw.ac.at/p3w/electronic1/GlobalTrigger/doc/InterfaceDesc/update_CMS_NOTE_2002_069.pdf

Definition at line 202 of file L1GtPatternGenerator.cc.

References L1GtPatternMap::begin(), and L1GtPatternMap::end().

Referenced by analyze().

 {
   // iterate over each pattern line
   for(L1GtPatternMap::LineMap::iterator it = allPatterns.begin(); 
       it != allPatterns.end(); ++it) { 
     // Get the HF bit counts and ring sums 
     uint32_t counts = it->second.get("hfBitCounts1");
     uint32_t sums = it->second.get("hfRingEtSums1");
     
     
     // Bits 0..11 -> 4 bit counts
     uint32_t hfPsbValue = (counts & 0xFFF) |
         // Bit 12..14 ring 1 pos. rap. HF Et sum
         (sums & 0x7) << 12 | 
         // Bits 16.. rest of the ET sums
         (sums >> 3) << 16;
     // TODO: Spec states non-data values for Bits 15, 31, 47 and 63.
 
     // Export computed value to pattern writer. */
     it->second.push(resultName, hfPsbValue);
   }
 }

uint32_t L1GtPatternGenerator::packRegionalMuons ( uint32_t rawData )

staticprivate

Post-processing for regional muon trigger inputs

Bits 8..15 (5 bits Pt, 3 bits quality) need to be inverted on the GMT inputs. See http://wwwhephy.oeaw.ac.at/p3w/cms/trigger/globalMuonTrigger/notes/in04_022.pdf

Definition at line 159 of file L1GtPatternGenerator.cc.

Referenced by analyze().

                                                                  {
     uint32_t invertMask = 0x0000FF00;
     uint32_t toKeep =   rawData & (~invertMask);
     return toKeep | (~rawData & invertMask);
 }