RctDigiToRctText Class Reference

#include <RctDigiToRctText.h>

Inheritance diagram for RctDigiToRctText:

Public Member Functions
	RctDigiToRctText (const edm::ParameterSet &)
	~RctDigiToRctText () override
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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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	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...

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	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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.

    : 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);
}

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

~RctDigiToRctText()

RctDigiToRctText::~RctDigiToRctText ( )

override

close files

Definition at line 33 of file RctDigiToRctText.cc.

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

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

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::EDAnalyzer.

Definition at line 40 of file RctDigiToRctText.cc.

                                                                                  {
  static int nevt = -1;
  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;
}

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

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