GtPsbTextToDigi Class Reference

#include <GtPsbTextToDigi.h>

Inheritance diagram for GtPsbTextToDigi:

Public Member Functions
	GtPsbTextToDigi (const edm::ParameterSet &)
	~GtPsbTextToDigi () override
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () 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	endJob () override
void	produce (edm::Event &, const edm::EventSetup &) override
void	putEmptyDigi (edm::Event &)
	Create empty digi collection. More...

Private Attributes
int	m_bc0 [4]
	Hold detected BC0 signal position per file. More...
std::ifstream	m_file [4]
	File handle. More...
int	m_fileEventOffset
	Number of events to be offset wrt input. More...
int	m_nevt
	Event counter. More...
std::string	m_textFileName
	Name out input file. More...

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 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

Definition at line 34 of file GtPsbTextToDigi.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 15 of file GtPsbTextToDigi.cc.

References Exception, MillePedeFileConverter_cfg::fileName, mps_fire::i, compare_using_db::ifile, cuy::ii, recoMuon::in, m_bc0, m_file, m_textFileName, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                              :
  m_fileEventOffset(iConfig.getUntrackedParameter<int>("FileEventOffset",0)),
  m_textFileName(iConfig.getParameter<std::string>("TextFileName")),
  m_nevt(0) {
  // Produces collections
  produces<L1GctEmCandCollection>("isoEm");
  produces<L1GctEmCandCollection>("nonIsoEm");
  produces<L1GctJetCandCollection>("cenJets");
  produces<L1GctJetCandCollection>("forJets");
  produces<L1GctJetCandCollection>("tauJets");
  //produces<L1GctEtTotal>();

  // Open the input files
  for (int ifile=0; ifile<4; ifile++){
    // gct em cand coll: (noiso) 0<->0, 1<->1, (iso) 6<->2, 7<->3
    int ii = (ifile<2)?ifile:ifile+4;
    std::stringstream fileStream;
    fileStream << m_textFileName << ii << ".txt";
    std::string fileName(fileStream.str());
    m_file[ifile].open(fileName.c_str(),std::ios::in);
    if(!m_file[ifile].good()) {
      throw cms::Exception("GtPsbTextToDigiTextFileOpenError")
    //LogDebug("GtPsbTextToDigi")
    << "GtPsbTextToDigi::GtPsbTextToDigi : "
    << " couldn't open the file " << fileName 
    //<< "...skipping!" 
    << std::endl;
    }
    std::hex(m_file[ifile]); 
  }

  // Initialize bc0 position holder
  for(int i=0; i<4; i++)
    m_bc0[i]=-1;  
}

GtPsbTextToDigi::~GtPsbTextToDigi ( )

override

Definition at line 51 of file GtPsbTextToDigi.cc.

References mps_fire::i, and m_file.

                                  {
  // Close the input files
  for (unsigned i=0; i<4; i++){  
    m_file[i].close();
  }
}

Member Function Documentation

void GtPsbTextToDigi::endJob ( void  )

overrideprivatevirtual

Check BC0 signals consistency

Reimplemented from edm::EDProducer.

Definition at line 223 of file GtPsbTextToDigi.cc.

References mps_fire::i, LogDebug, m_bc0, and match().

                             {
  int nmem = 4;
  bool match = true;
  for(int i=0; i<nmem-1; i++)
    match &= (m_bc0[i]==m_bc0[i+1]);  
  LogDebug("GtPsbTextToDigi") << "[GtPsbTextToDigi::endJob()] ";
  if(!match) 
    LogDebug("GtPsbTextToDigi") << "did not find matching BC0 in all input files: ";
  else
    LogDebug("GtPsbTextToDigi") << "detected common BC0 in all input files: ";
  for(int i=0; i<nmem; i++)
    LogDebug("GtPsbTextToDigi") << " " << m_bc0[i];  
  LogDebug("GtPsbTextToDigi") << std::flush << std::endl;

}

void GtPsbTextToDigi::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

buffer

Read in file

cycle debuging (temporary)

Fill in digi collections

Put collections

Definition at line 84 of file GtPsbTextToDigi.cc.

References funct::abs(), gather_cfg::cout, data, TauDecayModes::dec, Exception, mps_fire::i, compare_using_db::ifile, cuy::ii, LogDebug, m_bc0, m_file, m_fileEventOffset, m_nevt, m_textFileName, eostools::move(), edm::Event::put(), putEmptyDigi(), tmp, and globals_cff::x1.

                                                                           {

  // specify clock cycle bit sequence 1 0 1 0... or 0 1 0 1...
  unsigned short cbs[2] = {1,0};

  // Skip event if required
  if (m_nevt < m_fileEventOffset){ 
    putEmptyDigi(iEvent);
    LogDebug("GtPsbTextToDigi")
      << "[GtPsbTextToDigi::produce()] skipping event " << m_nevt 
      << std::endl;
    m_nevt++;
    return;
  } else if (m_nevt==0 && m_fileEventOffset<0) {
    // skip first fileEventOffset input events
    unsigned long int buff;
    for (int ievt=0; ievt<abs(m_fileEventOffset); ievt++){  
      for (int ifile=0; ifile<4; ifile++){  
    for (int cycle=0; cycle<2; cycle++){  
      std::hex(m_file[ifile]);
      m_file[ifile] >> buff;
      unsigned tmp = (buff>>15)&0x1;
      if(tmp!=cbs[cycle]) {
        if(m_bc0[ifile]==-1 && cycle==1 && tmp==1)
          m_bc0[ifile] = ievt;
        else 
          throw cms::Exception("GtPsbTextToDigiTextFileFormatError") 
        //std::cout << "GtPsbTextToDigiTextFileFormatError" 
        << "GtPsbTextToDigi::produce : "
        << " found format inconsistency in file #" << ifile 
        << "\n in skipped line:" << ievt*2+1
        << " cycle:" << tmp << " is different from " << cbs[cycle]
        << std::endl;
      }
    } 
      }
      LogDebug("GtPsbTextToDigi")
    << "[GtPsbTextToDigi::produce()] skipping input " << ievt 
    << std::endl;
    }
  }
  m_nevt++;
  
  // New collections
  std::unique_ptr<L1GctEmCandCollection>  gctIsolaEm( new L1GctEmCandCollection () );
  std::unique_ptr<L1GctEmCandCollection>  gctNoIsoEm( new L1GctEmCandCollection () );
  std::unique_ptr<L1GctJetCandCollection> gctCenJets( new L1GctJetCandCollection() ); 
  std::unique_ptr<L1GctJetCandCollection> gctForJets( new L1GctJetCandCollection() ); 
  std::unique_ptr<L1GctJetCandCollection> gctTauJets( new L1GctJetCandCollection() );
  //std::unique_ptr<L1GctEtTotal>           gctEtTotal( new L1GctEtTotal          () );

  uint16_t data[4][2]= {{0}};
  for (int i=0; i<4; i++)
    for (int j=0; j<2; j++)
      data[i][j]=0;

  // Loop over files
  for (int ifile=0; ifile<4; ifile++){  
    int ii = (ifile<2)?ifile:ifile+4;
    
    // Check we're not at the end of the file
    if(m_file[ifile].eof()) {
      LogDebug("GtPsbTextToDigi")
    << "GtPsbTextToDigi::produce : "
    << " unexpected end of file " << m_textFileName << ii << ".txt"
    << std::endl;
      putEmptyDigi(iEvent);
      continue;
    }      
    
    if(!m_file[ifile].good()) {
      LogDebug("GtPsbTextToDigi")
    << "GtPsbTextToDigi::produce : "
    << " problem reading file " << m_textFileName << ii << ".txt"
    << std::endl;
      putEmptyDigi(iEvent);
      continue;
    }
    
    unsigned long int uLongBuffer;
    
    for(unsigned cycle=0; cycle<2; cycle++) {
      m_file[ifile] >> uLongBuffer;
      unsigned tmp = (uLongBuffer>>15)&0x1;
      
      if(false && tmp!=cbs[cycle]) 
        std::cout << "[GtPsbTextToDigi::produce()] asserting " 
              << " evt:"   << m_nevt
              << " ifile:" << ifile
              << " cycle:" << cbs[cycle]
              << std::hex
              << " buffer:"<< uLongBuffer
              << " tmp: "  << tmp
              << std::dec 
              << "\n\n" << std::flush;

      if(tmp!=cbs[cycle]){
    if(m_bc0[ifile]==-1 && cycle==1 && tmp==1){
      m_bc0[ifile] = (m_nevt-m_fileEventOffset);
    }else{
      throw cms::Exception("GtPsbTextToDigiTextFileFormatError") 
        //std::cout << "GtPsbTextToDigiTextFileFormatError " 
        << "GtPsbTextToDigi::produce : "
        << " found format inconsistency in file #" << ifile 
        << "\n in line:" << (m_nevt-m_fileEventOffset)*2-1  
        << " cycle:" << tmp << " is different from " << cbs[cycle]
        << std::endl;
        }
      }
      data[ifile][cycle] = (uLongBuffer&0x7fff);
    } //cycle loop
  } //ifile loop
  
  unsigned iIsola, iNoIso;
  for (unsigned cycle=0; cycle<2; cycle++){  
    for (unsigned i=0; i<2; i++){  
      iIsola = i+2;
      iNoIso = i;
      gctIsolaEm->push_back(L1GctEmCand(data[iIsola][cycle]&0x7fff,true));
      gctNoIsoEm->push_back(L1GctEmCand(data[iNoIso][cycle]&0x7fff,false));
      L1GctEmCand candI(data[iIsola][cycle],true);
      L1GctEmCand candN(data[iNoIso][cycle],false);
    }
  }

  iEvent.put(std::move(gctIsolaEm), "isoEm");
  iEvent.put(std::move(gctNoIsoEm), "nonIsoEm");
  iEvent.put(std::move(gctCenJets), "cenJets");
  iEvent.put(std::move(gctForJets), "forJets");
  iEvent.put(std::move(gctTauJets), "tauJets");
  //iEvent.put(std::move(gctEtTotal));

}

void GtPsbTextToDigi::putEmptyDigi ( edm::Event &  iEvent )

private

Create empty digi collection.

Append empty digi collection.

Definition at line 59 of file GtPsbTextToDigi.cc.

References mps_fire::i, LogDebug, m_nevt, eostools::move(), and edm::Event::put().

Referenced by produce().

                                                   {
  std::unique_ptr<L1GctEmCandCollection>  gctIsolaEm( new L1GctEmCandCollection () );
  std::unique_ptr<L1GctEmCandCollection>  gctNoIsoEm( new L1GctEmCandCollection () );
  std::unique_ptr<L1GctJetCandCollection> gctCenJets( new L1GctJetCandCollection() );
  std::unique_ptr<L1GctJetCandCollection> gctForJets( new L1GctJetCandCollection() );
  std::unique_ptr<L1GctJetCandCollection> gctTauJets( new L1GctJetCandCollection() );
  //std::unique_ptr<L1GctEtTotal>           gctEtTotal( new L1GctEtTotal          () );
  for (int i=0; i<4; i++){  
    gctIsolaEm->push_back(L1GctEmCand (0,true));
    gctNoIsoEm->push_back(L1GctEmCand (0,false));
    gctCenJets->push_back(L1GctJetCand(0,false,false));
    gctForJets->push_back(L1GctJetCand(0,false,true));
    gctTauJets->push_back(L1GctJetCand(0,true,false));
    //gctEtTotal->push_back(());
  }
  iEvent.put(std::move(gctIsolaEm), "isoEm");
  iEvent.put(std::move(gctNoIsoEm), "nonIsoEm");
  iEvent.put(std::move(gctCenJets), "cenJets");
  iEvent.put(std::move(gctForJets), "forJets");
  iEvent.put(std::move(gctTauJets), "tauJets");
  //iEvent.put(std::move(gctEtTotal));

  LogDebug("GtPsbTextToDigi") << "putting empty digi (evt:" << m_nevt << ")\n"; 
}

Member Data Documentation

int GtPsbTextToDigi::m_bc0[4]

private

Hold detected BC0 signal position per file.

Definition at line 60 of file GtPsbTextToDigi.h.

Referenced by endJob(), GtPsbTextToDigi(), and produce().

std::ifstream GtPsbTextToDigi::m_file[4]

private

File handle.

Definition at line 57 of file GtPsbTextToDigi.h.

Referenced by GtPsbTextToDigi(), produce(), and ~GtPsbTextToDigi().

int GtPsbTextToDigi::m_fileEventOffset

private

Number of events to be offset wrt input.

Definition at line 48 of file GtPsbTextToDigi.h.

Referenced by produce().

int GtPsbTextToDigi::m_nevt

private

Event counter.

Definition at line 54 of file GtPsbTextToDigi.h.

Referenced by produce(), and putEmptyDigi().

std::string GtPsbTextToDigi::m_textFileName

private

Name out input file.

Definition at line 51 of file GtPsbTextToDigi.h.

Referenced by GtPsbTextToDigi(), and produce().