RctDigiToRctText Class Reference

#include <RctDigiToRctText.h>

Inheritance diagram for RctDigiToRctText:

Public Member Functions
	RctDigiToRctText (const edm::ParameterSet &)
	~RctDigiToRctText ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

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
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 32 of file RctDigiToRctText.h.

Constructor & Destructor Documentation

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

explicit

open output text files

open info|debug file

Definition at line 12 of file RctDigiToRctText.cc.

References edm::hlt::Exception, fdebug, convertXMLtoSQLite_cfg::fileName, i, m_file, m_textFileName, NUM_RCT_CRATES, and dbtoconf::out.

                                                                 : 
  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

(

)

close files

Definition at line 37 of file RctDigiToRctText.cc.

References fdebug, 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

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

privatevirtual

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 46 of file RctDigiToRctText.cc.

References data, edm::hlt::Exception, fdebug, edm::Event::getByLabel(), i, errorMatrix2Lands_multiChannel::id, j, LogDebug, m_file, m_hexUpperCase, m_rctInputLabel, nevt, NUM_RCT_CRATES, metsig::tau, and tmp.

                                                                             {

  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]={0};
  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

std::ofstream RctDigiToRctText::fdebug

private

handle for debug file

Definition at line 55 of file RctDigiToRctText.h.

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

std::ofstream RctDigiToRctText::m_file[NUM_RCT_CRATES]

private

handles for output files

Definition at line 52 of file RctDigiToRctText.h.

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

bool RctDigiToRctText::m_hexUpperCase

private

write upper case hex words

Definition at line 49 of file RctDigiToRctText.h.

Referenced by analyze().

edm::InputTag RctDigiToRctText::m_rctInputLabel

private

label for RCT digis

Definition at line 43 of file RctDigiToRctText.h.

Referenced by analyze().

std::string RctDigiToRctText::m_textFileName

private

basename for output files

Definition at line 46 of file RctDigiToRctText.h.

Referenced by RctDigiToRctText().