EcalSimRawData.cc

Go to the documentation of this file.

#include "SimCalorimetry/EcalElectronicsEmulation/interface/EcalSimRawData.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include <memory>

#include <fstream> //used for debugging
#include <iostream>
#include <iomanip>
#include <cmath>
#include "DataFormats/EcalDigi/interface/EcalMGPASample.h"
#include "FWCore/Utilities/interface/Exception.h"

using namespace std;
using namespace edm;

const int EcalSimRawData::ttType[nEbTtEta] = {
  0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,1, //EE-
  0,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1//EE+
};

const int EcalSimRawData::stripCh2Phi[nTtTypes][ttEdge][ttEdge] = {
  //TT type 0:
  /*ch-->*/
  {{4,3,2,1,0}, /*strip*/       
   {0,1,2,3,4}, /*|*/
   {4,3,2,1,0}, /*|*/
   {0,1,2,3,4}, /*|*/
   {4,3,2,1,0}},/*V*/
  //TT type 1:
  {{0,1,2,3,4},
   {4,3,2,1,0},
   {0,1,2,3,4},
   {4,3,2,1,0},
   {0,1,2,3,4}}
};

const int EcalSimRawData::strip2Eta[nTtTypes][ttEdge] = {
  {4,3,2,1,0}, //TT type 0
  {0,1,2,3,4}  //TT type 1
};

EcalSimRawData::EcalSimRawData(const edm::ParameterSet& params){
  //sets up parameters:
  digiProducer_ = params.getParameter<string>("unsuppressedDigiProducer");
  ebDigiCollection_ = params.getParameter<std::string>("EBdigiCollection");
  eeDigiCollection_ = params.getParameter<std::string>("EEdigiCollection");
  srDigiProducer_ = params.getParameter<string>("srProducer");
  ebSrFlagCollection_ = params.getParameter<std::string>("EBSrFlagCollection");
  eeSrFlagCollection_ = params.getParameter<std::string>("EESrFlagCollection");
  tpDigiCollection_
    = params.getParameter<std::string>("trigPrimDigiCollection");
  tcpDigiCollection_ = params.getParameter<std::string>("tcpDigiCollection");
  tpProducer_ = params.getParameter<string>("trigPrimProducer");
  xtalVerbose_ = params.getUntrackedParameter<bool>("xtalVerbose", false);
  tpVerbose_ = params.getUntrackedParameter<bool>("tpVerbose", false);
  tcc2dcc_ = params.getUntrackedParameter<bool>("tcc2dccData", true);
  srp2dcc_ = params.getUntrackedParameter<bool>("srp2dccData", true);
  fe2dcc_ = params.getUntrackedParameter<bool>("fe2dccData", true);
  fe2tcc_ = params.getUntrackedParameter<bool>("fe2tccData", true);
  dccNum_ = params.getUntrackedParameter<int>("dccNum", -1);
  tccNum_ = params.getUntrackedParameter<int>("tccNum", -1);
  tccInDefaultVal_ = params.getUntrackedParameter<int>("tccInDefaultVal", 0xffff);
  basename_ = params.getUntrackedParameter<std::string>("outputBaseName");


  string writeMode = params.getParameter<string>("writeMode");

  if(writeMode==string("littleEndian")){
    writeMode_ = littleEndian;
  } else if(writeMode==string("bigEndian")){
    writeMode_ = bigEndian;
  } else{
    writeMode_ = ascii;
  }
}

void
EcalSimRawData::analyze(const edm::Event& event,
            const edm::EventSetup& es){
  //Event counter:
  static int iEvent = 0;
  ++iEvent; 
    
  if(xtalVerbose_ | tpVerbose_){
    cout << "======================================================================\n"
     << " Event " << iEvent << "\n"
     << "----------------------------------------------------------------------\n";
  }

  if(fe2dcc_){
    vector<uint16_t> adc[nEbEta][nEbPhi];
    getEbDigi(event, adc);
    genFeData(basename_, iEvent, adc);
  }

  if(fe2tcc_){
    int tcp[nTtEta][nTtPhi]={{0}};
    getTp(event, tcpDigiCollection_, tcp);
    genTccIn(basename_, iEvent, tcp);
  }
  
  if(tcc2dcc_){
    int tp[nTtEta][nTtPhi]={{0}};
    getTp(event, tpDigiCollection_, tp); 
    genTccOut(basename_, iEvent, tp);
  }
  
  //SR flags:
  int ebSrf[nTtEta][nTtPhi];
  int eeSrf[nEndcaps][nScX][nScY];

  if(srp2dcc_){
    getSrfs(event, ebSrf, eeSrf);
    genSrData(basename_, iEvent, ebSrf);
  }
  
}

void EcalSimRawData::elec2GeomNum(int ittEta0, int ittPhi0, int strip1,
                  int ch1, int& iEta0, int& iPhi0) const{
  assert(0<=ittEta0 && ittEta0<nEbTtEta);
  assert(0<=ittPhi0 && ittPhi0<nTtPhi);
  assert(1<=strip1&& strip1<=ttEdge);
  assert(1<=ch1 && ch1<=ttEdge);
  const int type = ttType[ittEta0];
  iEta0 = ittEta0*ttEdge + strip2Eta[type][strip1-1];
  iPhi0 = ittPhi0*ttEdge + stripCh2Phi[type][strip1-1][ch1-1];
  assert(0<=iEta0 && iEta0<nEbEta);
  assert(0<=iPhi0 && iPhi0<nEbPhi);
}

void EcalSimRawData::fwrite(ofstream& f, uint16_t data,
                int& iWord, bool hpar) const{

  if(hpar){
    //set horizontal odd parity bit:
    setHParity(data);
  }
  
  switch(writeMode_){
  case littleEndian:
    {
      char c = data & 0x00FF;
      f.write(&c, sizeof(c));
      c = (data >>8) & 0x00FF;
      f.write(&c, sizeof(c));
    }
    break;
  case bigEndian:
    {
      char c = (data >>8) & 0x00FF; 
      f.write(&c, sizeof(c));
      c = data & 0x00FF;
      f.write(&c, sizeof(c));
    }
    break;
  case ascii:
    f << ((iWord%8==0&&iWord!=0)?"\n":"")
      << "0x" << setfill('0') << setw(4) << hex << data << "\t"
      << dec << setfill(' ');
    break;
  }
  ++iWord;
}

string EcalSimRawData::getExt() const{
  switch(writeMode_){
  case littleEndian:
    return ".le";
  case bigEndian:
    return ".be";
  case ascii:
    return ".txt";
  default:
    return".?";
  }  
}

void EcalSimRawData::genFeData(string basename, int iEvent,
                   const vector<uint16_t> adcCount[nEbEta][nEbPhi]
                   ) const{
  int smf = 0;
  int gmf = 0;
  int nPendingEvt = 0;
  int monitorFlag = 0;
  int chFrameLen = adcCount[0][0].size() + 1;
  
  int iWord = 0;
  
  for(int iZ0 = 0; iZ0<2; ++iZ0){
    for(int iDccPhi0 = 0; iDccPhi0<nDccInPhi; ++iDccPhi0){
      int iDcc1 = iDccPhi0 + iZ0*nDccInPhi + nDccEndcap + 1;

      if(dccNum_!=-1  && dccNum_!=iDcc1) continue;
      
      stringstream s;
      s.str("");
      const string& ext = getExt();
      s << basename << "_fe2dcc" << setfill('0') << setw(2) << iDcc1
    << setfill(' ') << ext;
      ofstream f(s.str().c_str(), (iEvent==1?ios::ate:ios::app));

      if(!f) return;


      if(writeMode_==ascii){
    f << (iEvent==1?"":"\n") << "[Event:" << iEvent << "]\n";
      }
      
      for(int iTtEtaInSm0 = 0; iTtEtaInSm0 < nTtSmEta; ++iTtEtaInSm0){
    int iTtEta0 = iZ0*nTtSmEta + iTtEtaInSm0;
    for(int iTtPhiInSm0 = 0; iTtPhiInSm0 < nTtSmPhi; ++iTtPhiInSm0){
      //phi=0deg at middle of 1st barrel DCC:
      int iTtPhi0 = -nTtPhisPerEbDcc/2 + iDccPhi0*nTtPhisPerEbDcc
        + iTtPhiInSm0;
      if(iTtPhi0<0) iTtPhi0 += nTtPhi;
      for(int stripId1 = 1; stripId1 <= ttEdge; ++stripId1){
        uint16_t stripHeader =
          0xF << 11
          | (nPendingEvt & 0x3F) << 5
          | (gmf & 0x1) << 4
          | (smf & 0x1) << 3
          | (stripId1 & 0x7);
        fwrite(f,stripHeader, iWord);

        for(int xtalId1 = 1; xtalId1 <= ttEdge; ++xtalId1){
          
          uint16_t crystalHeader =
        1 <<14
        | (chFrameLen & 0xFF) <<4
        | (monitorFlag & 0x1) <<3
        | (xtalId1 & 0x7);
          //          crystalHeader |=parity(crystalHeader) << 15;
          fwrite(f, crystalHeader, iWord);
          
          int iEta0;
          int iPhi0;
          elec2GeomNum(iTtEta0, iTtPhi0, stripId1, xtalId1,
               iEta0, iPhi0);
          if(xtalVerbose_){
        cout << dec
             << "iDcc1 = " << iDcc1 << "\t"
             << "iEbTtEta0 = " << iTtEta0 << "\t"
             << "iEbTtPhi0 = " << iTtPhi0 << "\t"
             << "stripId1 = " << stripId1 << "\t"
             << "xtalId1 = " << xtalId1 << "\t"
             << "iEta0 = " << iEta0 << "\t"
             << "iPhi0 = " << iPhi0 << "\t"
             << "adc[5] = 0x" << hex << adcCount[iEta0][iPhi0][5]
             << dec << "\n";
          }
          

          const vector<uint16_t>& adc = adcCount[iEta0][iPhi0];
          for(unsigned iSample=0; iSample  < adc.size(); ++iSample){
        uint16_t data = adc[iSample] & 0x3FFF;
        //      data |= parity(data);
        fwrite(f, data, iWord);
          } //next time sample
        } //next crystal in strip
      } //next strip in TT
    } //next TT along phi
      } //next TT along eta
    } //next DCC
  } //next half-barrel
}

void EcalSimRawData::genSrData(string basename, int iEvent,
                   int srf[nEbTtEta][nTtPhi]) const{
  for(int iZ0 = 0; iZ0<2; ++iZ0){
    for(int iDccPhi0 = 0; iDccPhi0<nDccInPhi; ++iDccPhi0){
      int iDcc1 = iDccPhi0 + iZ0*nDccInPhi + nDccEndcap + 1;
      if(dccNum_!=-1  && dccNum_!=iDcc1) continue;
      stringstream s;
      s.str("");
      s << basename << "_ab2dcc" << setfill('0') << setw(2) << iDcc1
    << setfill(' ') << getExt();
      ofstream f(s.str().c_str(), (iEvent==1?ios::ate:ios::app));
      
      if(!f) throw cms::Exception(string("Cannot create/open file ")
                  + s.str() + ".");
      
      int iWord = 0;
  
      if(writeMode_==ascii){
    f << (iEvent==1?"":"\n")<< "[Event:" << iEvent << "]\n";
      }

      const uint16_t le1 = 0;
      const uint16_t le0 = 0;
      const uint16_t h1 = 1;
      const uint16_t nFlags = 68;
      uint16_t data =  (h1 & 0x1)<< 14
    | (le1 & 0x1) << 12
    | (le0 & 0x1) << 11
    | (nFlags & 0x7F);
      
      fwrite(f, data, iWord, true);
      
      int iFlag = 0;
      data = 0;

      for(int iTtEtaInSm0 = 0; iTtEtaInSm0 < nTtSmEta; ++iTtEtaInSm0){
    //  int iTtEbEta0 = iZ0*nTtSmEta + iTtEtaInSm0;
    int iTtEta0 = nEeTtEta + iZ0*nTtSmEta + iTtEtaInSm0;
    for(int iTtPhiInSm0 = 0; iTtPhiInSm0 < nTtSmPhi; ++iTtPhiInSm0){
      //phi=0deg at middle of 1st barrel DCC:
      int iTtPhi0 = -nTtPhisPerEbDcc/2 + iDccPhi0*nTtPhisPerEbDcc
        + iTtPhiInSm0;
      if(iTtPhi0<0) iTtPhi0 += nTtPhi;
      //flags are packed by four:
      //|15 |14 |13-12 |11      9|8      6|5      3|2      0| 
      //| P | 0 | X  X |  srf i+3| srf i+2| srf i+1| srf i  |
      //|   |   |      | field 3 |field 2 | field 1| field 0|
      const int field = iFlag%4;
      //cout << "TtEta0: " << iTtEta0 << "\tTtPhi0: " << iTtPhi0 << "\n";
      //cout << "#" << oct << (int)srf[iTtEta0][iTtPhi0] << "o ****> #" << oct << (srf[iTtEta0][iTtPhi0] << (field*3)) << "o\n" << dec;
      
      data |= srf[iTtEta0][iTtPhi0] << (field*3);

      if(field==3){
        //cout <<  srf[iTtEta0][iTtPhi0] << "----> 0x" << hex << data << "\n";
        fwrite(f, data, iWord, true);
        data = 0;
      }
      ++iFlag;
    } //next TT along phi
      } //next TT along eta
    } //next DCC
  } //next half-barrel
}


void EcalSimRawData::genTccIn(string basename, int iEvent,
                  const int tcp[nTtEta][nTtPhi]) const{
  for(int iZ0 = 0; iZ0<2; ++iZ0){
    for(int iTccPhi0 = 0; iTccPhi0<nTccInPhi; ++iTccPhi0){
      int iTcc1 = iTccPhi0 + iZ0*nTccInPhi + nTccEndcap + 1;
      
      if(tccNum_!=-1  && tccNum_!=iTcc1) continue;
      
      stringstream s;
      s.str("");
      const char* ext = ".txt"; //only ascii mode supported for TCP

      s << basename << "_tcc" << setfill('0') << setw(2) << iTcc1
    << setfill(' ') << ext;
      ofstream fe2tcc(s.str().c_str(), (iEvent==1?ios::ate:ios::app));

      if(!fe2tcc) throw cms::Exception(string("Failed to create file ")
                       + s.str() + ".");
      
      int memPos = iEvent-1;
      int iCh1 = 1;
      for(int iTtEtaInSm0 = 0; iTtEtaInSm0 < nTtSmEta; ++iTtEtaInSm0){
    int iTtEta0 = (iZ0==0) ? 27 - iTtEtaInSm0 : 28 + iTtEtaInSm0; 
    for(int iTtPhiInSm0 = 0; iTtPhiInSm0 < nTtSmPhi; ++iTtPhiInSm0){
      //phi=0deg at middle of 1st barrel DCC:
      int iTtPhi0 = -nTtPhisPerEbTcc/2 + iTccPhi0*nTtPhisPerEbTcc
        + iTtPhiInSm0;
      iTtPhi0 += nTtPhisPerEbTcc*iTccPhi0;
      if(iTtPhi0<0) iTtPhi0 += nTtPhi;
      uint16_t tp_fe2tcc = (tcp[iTtEta0][iTtPhi0] & 0x7ff) ; //keep only Et (9:0) and FineGrain (10)
      
      if(tpVerbose_ && tp_fe2tcc!=0){
        cout << dec
         << "iTcc1 = " << iTcc1 << "\t"
         << "iTtEta0 = " << iTtEta0 << "\t"
         << "iTtPhi0 = " << iTtPhi0 << "\t"
         << "iCh1 = " << iCh1 << "\t"
         << "memPos = " << memPos << "\t" 
         << "tp = 0x" << setfill('0') << hex << setw(3)
         << tp_fe2tcc
         << dec << setfill(' ') << "\n";
      }
      fe2tcc << iCh1 << "\t"
         << memPos << "\t"
         << setfill('0') << hex
         << "0x" << setw(4) << tp_fe2tcc << "\t"
         << "0"
         << dec << setfill(' ') << "\n";
      ++iCh1;
    } //next TT along phi
      } //next TT along eta
      fe2tcc << std::flush;
      fe2tcc.close();
    } //next TCC
  } //next half-barrel
}


void EcalSimRawData::genTccOut(string basename, int iEvent,
                   const int tps[nTtEta][nTtPhi]) const{
  int iDccWord = 0;

  for(int iZ0 = 0; iZ0<2; ++iZ0){
    for(int iTccPhi0 = 0; iTccPhi0<nTccInPhi; ++iTccPhi0){
      int iTcc1 = iTccPhi0 + iZ0*nTccInPhi + nTccEndcap + 1;

      if(tccNum_!=-1  && tccNum_!=iTcc1) continue;
      
      stringstream s;
      s.str("");
      const char* ext = ".txt"; //only ascii mode supported for TCP

      s << basename << "_tcc" << setfill('0') << setw(2) << iTcc1
    << setfill(' ') << ext;

      s.str("");
      s << basename << "_tcc2dcc" << setfill('0') << setw(2) << iTcc1
    << setfill(' ') << getExt();
      ofstream dccF(s.str().c_str(), (iEvent==1?ios::ate:ios::app));
      
      if(!dccF){
    cout << "Warning: failed to create or open file " << s << ".\n";
    return;
      }
      
      const uint16_t h1 = 1;
      const uint16_t le1 = 0;
      const uint16_t le0 = 0;
      const uint16_t nSamples = 1;
      const uint16_t nTts = 68;
      const uint16_t data = (h1 & 0x1) << 14
    | (le1 & 0x1) << 12
    | (le0 & 0x1) << 11
    | (nSamples & 0xF) << 7
    | (nTts & 0x7F);    
      dccF << (iEvent==1?"":"\n") << "[Event:" << iEvent << "]\n";
      fwrite(dccF, data, iDccWord, false);
      
      int memPos = iEvent-1;
      int iCh1 = 1;
      for(int iTtEtaInSm0 = 0; iTtEtaInSm0 < nTtSmEta; ++iTtEtaInSm0){
    int iTtEta0 = nEeTtEta + iZ0*nTtSmEta + iTtEtaInSm0;
    for(int iTtPhiInSm0 = 0; iTtPhiInSm0 < nTtSmPhi; ++iTtPhiInSm0){
      //phi=0deg at middle of 1st barrel DCC:
      int iTtPhi0 = -nTtPhisPerEbTcc/2 + iTccPhi0*nTtPhisPerEbTcc
        + iTtPhiInSm0;
      if(iTtPhi0<0) iTtPhi0 += nTtPhi;
      
      if(tpVerbose_){
        cout << dec
         << "iTcc1 = " << iTcc1 << "\t"
         << "iTtEta0 = " << iTtEta0 << "\t"
         << "iTtPhi0 = " << iTtPhi0 << "\t"
         << "iCh1 = " << iCh1 << "\t"
         << "memPos = " << memPos << "\t" 
         << "tp = 0x" << hex << tps[iTtEta0][iTtPhi0]
         << dec << "\n";
      }
      fwrite(dccF, tps[iTtEta0][iTtPhi0], iDccWord, false);
      ++iCh1;
    } //next TT along phi
      } //next TT along eta
    } //next TCC
  } //next half-barrel
}


void EcalSimRawData::setHParity(uint16_t& a) const{
  const int odd = 1 <<15;
  const int even = 0;
  //parity bit of numbers from 0x0 to 0xF:
  //                    0   1   2    3   4    5    6   7   8    9    A   B    C   D   E    F           
  const int p[16] = {even,odd,odd,even,odd,even,even,odd,odd,even,even,odd,even,odd,odd,even};
  //inverts parity bit (LSB) of 'a' in case of even parity:
  a ^= p[a&0xF] ^ p[(a>>4)&0xF] ^ p[(a>>8)&0xF] ^ p[a>>12&0xF] ^ odd;   
}

void EcalSimRawData::getSrfs(const edm::Event& event,
                 int ebSrf[nTtEta][nTtPhi],
                 int eeSrf[nEndcaps][nScX][nScY]) const{

  //EE
  edm::Handle<EESrFlagCollection> hEeSrFlags;
  event.getByLabel(srDigiProducer_, eeSrFlagCollection_, hEeSrFlags);
  for(size_t i=0; i < sizeof(eeSrf)/sizeof(int); ((int*)eeSrf)[i++] = -1){};
  if(hEeSrFlags.isValid()){
    for(EESrFlagCollection::const_iterator it = hEeSrFlags->begin();
    it != hEeSrFlags->end(); ++it){
      const EESrFlag& flag = *it;
      int iZ0 = flag.id().zside()>0?1:0;
      int iX0 = flag.id().ix()-1;
      int iY0 = flag.id().iy()-1;
      assert(iZ0>=0 && iZ0<nEndcaps);
      assert(iX0>=0 && iX0<nScX);
      assert(iY0>=0 && iY0<nScY);
      eeSrf[iZ0][iX0][iY0] = flag.value();
    }
  } else{
    LogWarning("EcalSimRawData") << "EE SR flag not found ("
                 << "Product label: " << srDigiProducer_
                 << "Producet instance: "
                 << eeSrFlagCollection_ << ")";
  }

  //EB
  edm::Handle<EBSrFlagCollection> hEbSrFlags;
  event.getByLabel(srDigiProducer_, ebSrFlagCollection_, hEbSrFlags);
  for(size_t i=0; i<sizeof(ebSrf)/sizeof(int); ((int*)ebSrf)[i++] = -1){};
  if(hEbSrFlags.isValid()){
    for(EBSrFlagCollection::const_iterator it = hEbSrFlags->begin();
    it != hEbSrFlags->end(); ++it){
      
      const EBSrFlag& flag = *it;
      int iEta = flag.id().ieta();
      int iEta0 = iEta + nTtEta/2 - (iEta>=0?1:0); //0->55 from eta=-3 to eta=3
      int iEbEta0 = iEta0 - nEeTtEta;//0->33 from eta=-1.48 to eta=1.48
      int iPhi0 = flag.id().iphi() - 1;

      assert(iEbEta0>=0 && iEbEta0<nEbTtEta);
      assert(iPhi0>=0 && iPhi0<nTtPhi);
      
      ebSrf[iEbEta0][iPhi0] = flag.value();
    }
  } else{
    LogWarning("EcalSimRawData") << "EB SR flag not found ("
                 << "Product label: " << srDigiProducer_
                 << "Producet instance: "
                 << ebSrFlagCollection_ << ")";
  }
}

void EcalSimRawData::getEbDigi(const edm::Event& event,
                   vector<uint16_t> adc[nEbEta][nEbPhi]) const{

  edm::Handle<EBDigiCollection> hEbDigis;
  event.getByLabel(digiProducer_, ebDigiCollection_, hEbDigis);

  int nSamples = 0;
  if(hEbDigis.isValid() && hEbDigis->size()>0){//there is at least one digi
    nSamples  = hEbDigis->begin()->size();//gets the sample count from 1st digi
  }


  const uint16_t suppressed = 0xFFFF;

  adc[0][0] = vector<uint16_t>(nSamples, suppressed);
  
  for(int iEbEta=0; iEbEta<nEbEta; ++iEbEta){
    for(int iEbPhi=0; iEbPhi<nEbPhi; ++iEbPhi){
      adc[iEbEta][iEbPhi] = adc[0][0];
    }
  }
  if(hEbDigis.isValid()){
    if(xtalVerbose_) cout << setfill('0');
    for(EBDigiCollection::const_iterator it = hEbDigis->begin();
    it != hEbDigis->end(); ++it){
      const EBDataFrame& frame = *it;
      
      int iEta0 = iEta2cIndex((frame.id()).ieta());
      int iPhi0 = iPhi2cIndex((frame.id()).iphi());
      
      //     cout << "xtl indices conv: (" << frame.id().ieta() << ","
      //     << frame.id().iphi() << ") -> ("
      //     << iEta0 << "," << iPhi0 << ")\n";
    
      if(iEta0<0 || iEta0>=nEbEta){
    cout << "iEta0 (= " << iEta0 << ") is out of range ("
         << "[0," << nEbEta -1 << "])\n";
      }
      if(iPhi0<0 || iPhi0>=nEbPhi){
    cout << "iPhi0 (= " << iPhi0 << ") is out of range ("
         << "[0," << nEbPhi -1 << "])\n";
      }
    
      if(xtalVerbose_){
    cout << iEta0 << "\t" << iPhi0 << ":\t";
    cout << hex;
      }
      
      if(nSamples!=frame.size()){
    throw cms::Exception("EcalSimRawData", "Found EB digis with different sample count! This is not supported by EcalSimRawData.");
      }
    
      for(int iSample=0; iSample<nSamples; ++iSample){
    const EcalMGPASample& sample = frame.sample(iSample);
    uint16_t encodedAdc = sample.raw();
    adc[iEta0][iPhi0][iSample] = encodedAdc;  
    if(xtalVerbose_){
      cout << (iSample>0?" ":"") << "0x" << setw(4) 
           << encodedAdc;
    }
      }
      if(xtalVerbose_) cout << "\n" << dec;
    }
    if(xtalVerbose_) cout << setfill(' ');
  }
}
  

void EcalSimRawData::getTp(const edm::Event& event,
               const std::string& collName,
               int tcp[nTtEta][nTtPhi]) const{
  edm::Handle<EcalTrigPrimDigiCollection> hTpDigis;
  event.getByLabel(tpProducer_, collName, hTpDigis);
  if(hTpDigis.isValid() && hTpDigis->size()>0){
    const EcalTrigPrimDigiCollection& tpDigis = *hTpDigis.product();

    //    EcalSelectiveReadout::ttFlag_t ttf[nTtEta][nTtPhi];
    for(int iTtEta0=0; iTtEta0 < nTtEta; ++iTtEta0){
      for(int iTtPhi0=0; iTtPhi0 < nTtPhi; ++iTtPhi0){
    tcp[iTtEta0][iTtPhi0] = tccInDefaultVal_ ;
      }
    }
    if(tpVerbose_){
      cout << setfill('0');
    }
    for(EcalTrigPrimDigiCollection::const_iterator it = tpDigis.begin();
    it != tpDigis.end(); ++it){
      const EcalTriggerPrimitiveDigi& tp = *it;
      int iTtEta0 = iTtEta2cIndex(tp.id().ieta());
      int iTtPhi0 = iTtPhi2cIndex(tp.id().iphi());
      if(iTtEta0<0 || iTtEta0>=nTtEta){
    cout << "iTtEta0 (= " << iTtEta0 << ") is out of range ("
         << "[0," << nEbTtEta -1 << "])\n";
      }
      if(iTtPhi0<0 || iTtPhi0>=nTtPhi){
    cout << "iTtPhi0 (= " << iTtPhi0 << ") is out of range ("
         << "[0," << nTtPhi -1 << "])\n";
      }

      tcp[iTtEta0][iTtPhi0] = tp[tp.sampleOfInterest()].raw();

      if(tpVerbose_){
    if(tcp[iTtEta0][iTtPhi0]!=0) //print non-zero values only
    cout << collName << (collName.size()==0?"":" ")
         << "TP(" << setw(2) << iTtEta0 << "," << iTtPhi0 << ") = "
         << "0x" << setw(4) 
         << tcp[iTtEta0][iTtPhi0]
         << "\tcmssw indices: "
         << tp.id().ieta() << " " << tp.id().iphi() << "\n";
      }
    }//next TP
    if(tpVerbose_) cout << setfill(' ');
  }
}