#include <RctDigiToRctText.h>

Inheritance diagram for RctDigiToRctText:

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

	~RctDigiToRctText () 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
void	analyze (const edm::Event &, const edm::EventSetup &) override

Private Attributes
std::ofstream	fdebug
	handle for debug file More...

std::ofstream	m_file [NUM_RCT_CRATES]
	handles for output files More...

bool	m_hexUpperCase
	write upper case hex words More...

edm::InputTag	m_rctInputLabel
	label for RCT digis More...

std::string	m_textFileName
	basename for output files More...

int	nevt = -1

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

Definition at line 30 of file RctDigiToRctText.h.

Constructor & Destructor Documentation

◆ RctDigiToRctText()

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

explicit

open output text files

open info|debug file

Definition at line 11 of file RctDigiToRctText.cc.

References Exception, fdebug, MillePedeFileConverter_cfg::fileName, caHitNtupletGeneratorKernels::good, mps_fire::i, m_file, m_textFileName, NUM_RCT_CRATES, MillePedeFileConverter_cfg::out, and AlCaHLTBitMon_QueryRunRegistry::string.

     : m_rctInputLabel(iConfig.getParameter<edm::InputTag>("RctInputLabel")),
       m_textFileName(iConfig.getParameter<std::string>("TextFileName")),
       m_hexUpperCase(iConfig.getParameter<bool>("HexUpperCase")) {
   for (unsigned i = 0; i < NUM_RCT_CRATES; i++) {
     std::stringstream fileStream;
     fileStream << m_textFileName << std::setw(2) << std::setfill('0') << i << ".txt";
     std::string fileName(fileStream.str());
     m_file[i].open(fileName.c_str(), std::ios::out);
 
     if (!m_file[i].good()) {
       throw cms::Exception("RctDigiToRctTextTextFileOpenError")
           << "RctDigiToRctText::RctDigiToRctText : "
           << " couldn't create the file " << fileName << std::endl;
     }
   }
 
   fdebug.open("rctdigitorcttext_debug.txt", std::ios::out);
 }

◆ ~RctDigiToRctText()

RctDigiToRctText::~RctDigiToRctText ( )

override

close files

Definition at line 33 of file RctDigiToRctText.cc.

References fdebug, mps_fire::i, m_file, and NUM_RCT_CRATES.

                                     {
   for (unsigned i = 0; i < NUM_RCT_CRATES; i++)
     m_file[i].close();
   fdebug.close();
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

count bunch crossing

get the RCT data

debug flags and stream

— ELECTRONS —

debug

— REGIONS —

print electrons

print regions

mip bits

quiet bits

region info

HF

flush data to files

Implements edm::one::EDAnalyzerBase.

Definition at line 40 of file RctDigiToRctText.cc.

References data, TauDecayModes::dec, EgHLTOffHistBins_cfi::et, Exception, fdebug, tpScales_cff::HF, mps_fire::i, l1ctLayer2EG_cff::id, iEvent, dqmiolumiharvest::j, LogDebug, m_file, m_hexUpperCase, m_rctInputLabel, nevt, NUM_RCT_CRATES, and createJobs::tmp.

                                                                                   {
   nevt++;
 
   edm::Handle<L1CaloEmCollection> em;
   edm::Handle<L1CaloRegionCollection> rgn;
   iEvent.getByLabel(m_rctInputLabel, em);
   iEvent.getByLabel(m_rctInputLabel, rgn);
 
   bool ldebug = false;
   bool debug_NOTEMPTY[18] = {false};
   for (int i = 0; i < 18; i++)
     debug_NOTEMPTY[i] = false;
   std::stringstream dstrm;
 
 
   unsigned long int CAND[18][8];
   int n_iso[18] = {0};
   int n_niso[18] = {0};
   bool iso;
   int id;
 
   for (L1CaloEmCollection::const_iterator iem = em->begin(); iem != em->end(); iem++) {
     int crate = iem->rctCrate();
     iso = iem->isolated();
     unsigned data = iem->raw();
 
     id = iso ? n_iso[crate]++ : 4 + n_niso[crate]++;
 
     CAND[crate][id] = data;
 
     if (crate > 17 || id > 7)
       throw cms::Exception("RctDigiToRctTextElectronIndexOutBounds")
           << "out of bounds indices  crate:" << crate << "id:" << id << std::endl;
     if (ldebug && data != 0)
       debug_NOTEMPTY[crate] = true;
     dstrm.str("");
     dstrm << "electron "
           << " bx:" << nevt << " crate:" << crate << " iso:" << iso << " raw:" << data << " \t cand:" << *iem;
     if (debug_NOTEMPTY[crate])
       fdebug << dstrm.str() << std::endl;
   }
 
 
   unsigned short MIPbits[18][7][2] = {{{0}}};
   unsigned short QIEbits[18][7][2] = {{{0}}};
   unsigned short RC[18][7][2] = {{{0}}};
   unsigned short RCof[18][7][2] = {{{0}}};
   unsigned short RCtau[18][7][2] = {{{0}}};
   unsigned short HF[18][4][2] = {{{0}}};
 
   for (L1CaloRegionCollection::const_iterator irgn = rgn->begin(); irgn != rgn->end(); irgn++) {
     int crate = irgn->rctCrate();
     int card = irgn->rctCard();
     int rgnidx = irgn->rctRegionIndex();
 
     dstrm.str("");
     if (!irgn->id().isHf()) {
       RC[crate][card][rgnidx] = irgn->et();
       RCof[crate][card][rgnidx] = irgn->overFlow();
       RCtau[crate][card][rgnidx] = irgn->tauVeto();
       MIPbits[crate][card][rgnidx] = irgn->mip();
       QIEbits[crate][card][rgnidx] = irgn->quiet();
       // debug info
       dstrm << hex << "Et=" << irgn->et() << " OverFlow=" << irgn->overFlow() << " tauVeto=" << irgn->tauVeto()
             << " mip=" << irgn->mip() << " quiet=" << irgn->quiet() << " Card=" << irgn->rctCard()
             << " Region=" << irgn->rctRegionIndex() << " Crate=" << irgn->rctCrate() << dec;
       if (ldebug)
         LogDebug("Regions") << dstrm.str();
     } else {
       HF[crate][irgn->id().rctEta() - 7][irgn->id().rctPhi()] = irgn->et();
       // debug info
       dstrm << hex << "Et=" << irgn->et() << " FGrain=" << irgn->fineGrain() << " Eta=" << irgn->id().rctEta()
             << " Phi=" << irgn->id().rctPhi() << " Crate=" << irgn->rctCrate() << dec;
       if (ldebug)
         LogDebug("HFRegions") << dstrm.str();
     }
 
     if (ldebug && irgn->et() != 0)
       debug_NOTEMPTY[crate] = true;  // debug
     if (debug_NOTEMPTY[crate]) {
       fdebug << "region"
              << " bx:" << nevt << " crate:" << crate << "\t";
       fdebug << dstrm.str() << std::endl;
     }
   }
 
   std::stringstream sstrm;
   if (m_hexUpperCase)
     sstrm << std::uppercase;
   else
     sstrm.unsetf(std::ios::uppercase);
 
 
   for (unsigned crate = 0; crate < NUM_RCT_CRATES; crate++) {
     sstrm.str("");
     sstrm << "Crossing " << nevt << std::endl;
 
     for (int j = 0; j < 8; j++) {
       sstrm << setw(3) << setfill('0') << hex << (CAND[crate][j] & 0x3ff);
       if (j < 7)
         sstrm << " ";
     }
     sstrm << setfill(' ') << dec;
     m_file[crate] << sstrm.str() << std::endl;
 
     // debug
     if (debug_NOTEMPTY[crate])
       fdebug << sstrm.str() << std::endl;
     if (ldebug)
       LogDebug("Electrons") << sstrm.str() << std::endl;
   }
 
 
   for (unsigned crate = 0; crate < NUM_RCT_CRATES; crate++) {
     sstrm.str("");
     for (int card = 0; card < 7; card++) {
       for (int j = 0; j < 2; j++) {
         sstrm << " " << MIPbits[crate][card][j];
       }
     }
     m_file[crate] << sstrm.str() << std::endl;
     if (debug_NOTEMPTY[crate])
       fdebug << sstrm.str() << std::endl;  // debug
 
     sstrm.str("");
     for (int card = 0; card < 7; card++) {
       for (int j = 0; j < 2; j++) {
         sstrm << " " << QIEbits[crate][card][j];
       }
     }
     m_file[crate] << sstrm.str() << std::endl;
     if (debug_NOTEMPTY[crate])
       fdebug << sstrm.str() << std::endl;  // debug
 
     sstrm.str("");
     for (int card = 0; card < 7; card++) {
       for (int j = 0; j < 2; j++) {
         unsigned long int tmp;
         unsigned et = RC[crate][card][j];
         unsigned ovf = RCof[crate][card][j];
         unsigned tau = RCtau[crate][card][j];
         // ovf = ovf || (et>=0x400);
         tmp = ((tau & 0x1) << 11) | ((ovf & 0x1) << 10) | ((et & 0x3ff));
         sstrm << " " << setw(3) << setfill('0') << hex << tmp;
       }
     }
     m_file[crate] << sstrm.str() << std::endl;
     if (debug_NOTEMPTY[crate])
       fdebug << sstrm.str() << std::endl << std::endl;  // debug
 
     sstrm.str("");
     for (int ip = 0; ip < 2; ip++) {
       for (int ie = 0; ie < 4; ie++) {
         unsigned et = HF[crate][ie][ip] & 0xff;
         sstrm << " " << setw(2) << setfill('0') << hex << et;
       }
     }
     m_file[crate] << sstrm.str() << std::endl;
     if (debug_NOTEMPTY[crate])
       fdebug << sstrm.str() << std::endl;  // debug
     sstrm << setfill(' ') << dec;
 
   }  // end crate loop
 
   for (unsigned i = 0; i < NUM_RCT_CRATES; i++)
     m_file[i] << std::flush;
 
   fdebug << std::flush;
 }

Member Data Documentation

◆ fdebug

std::ofstream RctDigiToRctText::fdebug

private

handle for debug file

Definition at line 51 of file RctDigiToRctText.h.

Referenced by analyze(), RctDigiToRctText(), and ~RctDigiToRctText().

◆ m_file

std::ofstream RctDigiToRctText::m_file[NUM_RCT_CRATES]

private

handles for output files

Definition at line 48 of file RctDigiToRctText.h.

Referenced by analyze(), RctDigiToRctText(), and ~RctDigiToRctText().

◆ m_hexUpperCase

bool RctDigiToRctText::m_hexUpperCase

private

write upper case hex words

Definition at line 45 of file RctDigiToRctText.h.

Referenced by analyze().

◆ m_rctInputLabel

edm::InputTag RctDigiToRctText::m_rctInputLabel

private

label for RCT digis

Definition at line 39 of file RctDigiToRctText.h.

Referenced by analyze().

◆ m_textFileName

std::string RctDigiToRctText::m_textFileName

private

basename for output files

Definition at line 42 of file RctDigiToRctText.h.

Referenced by RctDigiToRctText().

◆ nevt

int RctDigiToRctText::nevt = -1

private

Definition at line 53 of file RctDigiToRctText.h.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RctDigiToRctText()

◆ ~RctDigiToRctText()

Member Function Documentation

◆ analyze()

Member Data Documentation

◆ fdebug

◆ m_file

◆ m_hexUpperCase

◆ m_rctInputLabel

◆ m_textFileName

◆ nevt