#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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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 FCDTask_cfi::crate, data, TauDecayModes::dec, l1tnanotables_cff::et, Exception, fdebug, tpScales_cff::HF, mps_fire::i, EcalPhiSymFlatTableProducers_cfi::id, l1tnanotables_cff::iem, iEvent, genparticles_cff::iso, 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