RctTextToRctDigi.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    RctTextToRctDigi
// Class:      RctTextToRctDigi
// 
//
// Original Author:  Alex Tapper
//         Created:  Fri Mar  9 19:11:51 CET 2007

// Rct Input File Format 
// Line 1: Crossing no as "Crossing x" (2)     
// Line 2: isoe0 isoe1 isoe2 isoe3 nonIsoe0 nonIsoe1 nonIso2 nonIso3 (8) 
// Line 3: RC0mip0 RC0mip1 RC1mip0 RC1mip1 RC2mip0 RC2mip1 RC3mip0 RC3mip1 RC4mip0 RC4mip1 
//         RC5mip0 RC5mip1 RC6mip0 RC6mip1 (14)
// Line 4: RC0qt0 RCqt1 RC1qt0 RC1qt1 RC2qt0 RC2qt1 RC3qt0 RC3qt1 RC4qt0 RC4qt1 
//         RC5qt0 RC5qt1 RC6qt0 RC6qt1 (14)
// Line 5: RC0reg0 RC0reg1 RC1reg0 RC1reg1 RC2reg0 RC2reg1 RC3reg0 RC3reg1 RC4reg0 RC4reg1
//         RC5reg0 RC5reg1 RC6reg0 RC6reg1 (14)
// Line 6: HF0eta0 HF0eta1 HF0eta2 HF0eta3 HF1eta0 HF1eta1 HF1eta2 HF1eta3 (8)
//
// NOTE:  CMS IN 2004/009 specifies that cable four provides 8 Quiet (fineGrain) bits for the HF.  These are not
//        detailed in the fileformat above, and are not currently dealt with in any way. Set to true.
// 

#include "RctTextToRctDigi.h"
#include "FWCore/ServiceRegistry/interface/Service.h" // Framework services
#include "FWCore/MessageLogger/interface/MessageLogger.h" // Logger
#include <iomanip>

using namespace edm;
using namespace std;

// Set constant
const static unsigned NUM_RCT_CRATES = 18;

RctTextToRctDigi::RctTextToRctDigi(const edm::ParameterSet& iConfig):
  m_textFileName(iConfig.getParameter<std::string>("TextFileName")),
  m_fileEventOffset(iConfig.getParameter<int>("FileEventOffset")),
  m_nevt(0)
{
  // Produces collections
  produces<L1CaloEmCollection>();
  produces<L1CaloRegionCollection>();

  // Open the input files
  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::in);

    if(!m_file[i].good())
      {
        //throw cms::Exception("RctTextToRctDigiTextFileOpenError")
    LogDebug("RctTextToRctDigi")
          << "RctTextToRctDigi::RctTextToRctDigi : "
          << " couldn't open the file " << fileName << "...skipping!" << std::endl;
      }
  }
}

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

void RctTextToRctDigi::putEmptyDigi(edm::Event& iEvent) {
  std::auto_ptr<L1CaloEmCollection> em (new L1CaloEmCollection);
  std::auto_ptr<L1CaloRegionCollection> rgn (new L1CaloRegionCollection);
    for (unsigned i=0; i<NUM_RCT_CRATES; i++){  
      for (unsigned j=0; j<4; j++) {
    em->push_back(L1CaloEmCand(0, i, true));
    em->push_back(L1CaloEmCand(0, i, false));
      }
      for (unsigned j=0; j<14; j++)
    rgn->push_back(L1CaloRegion(0,false,false,false,false,i,j/2,j%2));
      for (unsigned j=0; j<8; j++)
    rgn->push_back(L1CaloRegion(0,true,i,j));
    }
    iEvent.put(em);
    iEvent.put(rgn);
}

void RctTextToRctDigi::bxSynchro(int &bx,int crate) {
  std::string tmp;
  // bypass bx input until correct bx is reached
  while(bx<m_nevt+m_fileEventOffset) {
    for (int j=0; j<6; j++){
      getline(m_file[crate],tmp);
    }
    m_file[crate] >> tmp >> bx;
    if(tmp!="Crossing")
      throw cms::Exception("RctTextToRctDigiTextFileReadError")
    << "RctTextToRctDigi::bxSynchro : "
    << " something screwy happened Crossing!=" << tmp << std::endl;
  }
}

// ------------ method called to produce the data  ------------
void RctTextToRctDigi::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
{

  // Skip event if required
  if (m_nevt < m_fileEventOffset){ 
    //string tmp;
    //for (int i=0; i<6; i++){
    //  getline(m_file[i],tmp);
    //} 
    putEmptyDigi(iEvent);
    m_nevt++;
    return;
  }

  // New collections
  std::auto_ptr<L1CaloEmCollection> em (new L1CaloEmCollection);
  std::auto_ptr<L1CaloRegionCollection> rgn (new L1CaloRegionCollection);

  // Loop over RCT crates
  for (unsigned i=0; i<NUM_RCT_CRATES; i++){  

    if(!m_file[i].good()) {
      continue;
    }

    // Check we're not at the end of the file
    if(m_file[i].eof())
      {
        //throw cms::Exception("RctTextToRctDigiTextFileReadError")
    LogDebug("RctTextToRctDigi")
          << "RctTextToRctDigi::produce : "
          << " unexpected end of file " << m_textFileName << i 
      << " adding empty collection for event !"
      << std::endl;
    putEmptyDigi(iEvent);
    continue;
      }      
    
    // Check we're at the start of an event
    std::string tmp;
    m_file[i]>> tmp;
    if(tmp!="Crossing")
      {
        throw cms::Exception("RctTextToRctDigiTextFileReadError")
          << "RctTextToRctDigi::produce : "
          << " something screwy happened Crossing!=" << tmp << std::endl;
      }      

    // Read BX number
    dec(m_file[i]);
    int BXNum;
    m_file[i]>>BXNum;
    
    bxSynchro(BXNum,i);
 
    if(BXNum!=m_nevt+m_fileEventOffset)
      throw cms::Exception("RctTextToRctDigiTextSyncError")
    << "RctTextToRctDigi::produce : "
    << " something screwy happened "
    << "evt:" << m_nevt << " != bx:" << BXNum << " + " << m_fileEventOffset 
    << std::endl;
    
    // Buffers
    unsigned long int uLongBuffer;
    bool mipBitBuffer[14],qBitBuffer[14];

    // All in hex from now on
    hex(m_file[i]); 

    // Isolated electrons
    for (unsigned j=0; j<4; j++){
      m_file[i] >> uLongBuffer;
      em->push_back(L1CaloEmCand(uLongBuffer, i, true, j,BXNum,0));
    }

    // Non-isolated electrons
    for (unsigned j=0; j<4; j++){
      m_file[i] >> uLongBuffer;
      em->push_back(L1CaloEmCand(uLongBuffer, i, false, j,BXNum,0));
    }      
    
    // MIP bits 
    for (unsigned j=0; j<14; j++){
      m_file[i] >> mipBitBuffer[j];
    }   

    // Quiet bits 
    for (unsigned j=0; j<14; j++){
      m_file[i] >> qBitBuffer[j];
    }     

    // Barrel and endcap regions
    for (unsigned j=0; j<14; j++){
      m_file[i] >> uLongBuffer;

      unsigned et = uLongBuffer & 0x3ff;  // put the first 10 bits of rawData into the Et
      uLongBuffer >>= 10;  // shift the remaining bits down to remove the 10 bits of Et
      
      bool overFlow = ((uLongBuffer & 0x1)        != 0); //LSB is now overflow bit
      bool tauVeto  = (((uLongBuffer & 0x2) >> 1) != 0); //2nd bit is tauveto      

      rgn->push_back(L1CaloRegion(et,overFlow,tauVeto,mipBitBuffer[j],qBitBuffer[j],i,j/2,j%2));
    }      
    
    // HF
    for (unsigned j=0; j<8; j++){
      m_file[i] >> uLongBuffer;

      unsigned et = uLongBuffer & 0xff;  // put the first 8 bits into the Et

      rgn->push_back(L1CaloRegion(et,true,i,j));
    }        

    dec(m_file[i]); 
  }
  
  iEvent.put(em);
  iEvent.put(rgn);

  m_nevt++;
}