#include <EcalDumpRaw.h>

Inheritance diagram for EcalDumpRaw:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
void	analyzeEB (const edm::Event &, const edm::EventSetup &) const
void	analyzeEE (const edm::Event &, const edm::EventSetup &) const
	EcalDumpRaw (const edm::ParameterSet &)
void	endJob ()
	~EcalDumpRaw ()
Private Types
enum	{ inDaqHeader, inDccHeader, inTccBlock, inSrBlock, inTowerBlock }
Private Member Functions
void	analyzeApd ()
void	analyzeFed (int fedId)
bool	decode (const uint32_t *data, int iWord32, std::ostream &out)
double	max (std::vector< double > a, unsigned &pos)
double	min (std::vector< double > a)
std::string	srRange (int offset) const
std::string	toNth (int n)
template<class T >
std::string	toString (T val)
std::string	tpgTag (int tccType, unsigned iSeq) const
std::string	ttfTag (int tccType, unsigned iSeq) const
Static Private Member Functions
static int	lme (int dccId1, int side)
static int	lmodOfRu (int ru1)
static int	modOfRu (int ru1)
static int	sideOfRu (int ru1)
Private Attributes
std::vector< double >	adc_
double	amplCut_
int	beg_fed_id_
int	bx_
int	dccCh_
unsigned	dccId_
enum EcalDumpRaw:: { ... }	decodeState_
int	detailedTrigType_
bool	dump_
bool	dumpAdc_
std::ofstream	dumpFile_
int	end_fed_id_
unsigned	eventId_
std::vector< unsigned >	eventList_
std::vector< int >	feBx_
unsigned	fedId_
std::vector< int >	feL1a_
std::vector< int >	feRuId_
std::string	filename_
int	first_event_
int	iEvent_
int	iRu_
unsigned	iSrWord64_
int	iTcc_
unsigned	iTccWord64_
int	iTow_
unsigned	iTowerWord64_
int	l1a_
int	l1amax_
int	l1aMaxX_
int	l1amin_
int	l1aMinX_
std::vector< std::vector < uint32_t > >	l1as_
int	last_event_
std::vector< uint32_t >	lastOrbit_
unsigned	maxEventId_
int	maxEvt_
unsigned	minEventId_
std::vector< int >	nTpgs_
int	nTts_
unsigned	orbit0_
bool	orbit0Set_
uint32_t	orbit_
std::vector< std::vector < uint32_t > >	orbits_
int	profileFedId_
int	profileRuId_
bool	pulsePerLme_
bool	pulsePerLmod_
bool	pulsePerRu_
unsigned	side_
int	simpleTrigType_
int	srpBx_
int	srpL1a_
int	tccBlockLen64_
int	tccBx_
int	tccId_
int	tccL1a_
int	tccType_
	type of TCC currently parsed
size_t	towerBlockLength_
std::vector< std::vector< int > >	tpg_
int	verbosity_
bool	writeDcc_
Static Private Attributes
static const unsigned	fedStart_ = 601
static const int	maxTccsPerDcc_ = 4
static const int	maxTpgsPerTcc_ = 68
static const unsigned	nDccs_ = 54
static const int	nRu_ = 70
static const int	nSamples = 10
static const int	ttId_ [nTccTypes_][maxTpgsPerTcc_]

static const int	ebmTcc_ = 0
static const int	ebpTcc_ = 1
static const int	eeInnerTcc_ = 2
static const int	eeOuterTcc_ = 3
static const int	nTccTypes_ = 4

Detailed Description

Utility to dump ECAL Raw data. Hexadecimal dump is accompagned with a side by data interpretention.

The script test/dumpRaw can be used to run this module. E. g.: dumpRaw /store/..../data_file.root Run dumpRaw -h to get help on this script.

Author: Ph. Gras CEA/IRFU Saclay

Definition at line 32 of file EcalDumpRaw.h.

Member Enumeration Documentation

anonymous enum [private]

Enumerator:

inDaqHeader
inDccHeader
inTccBlock
inSrBlock
inTowerBlock

Definition at line 101 of file EcalDumpRaw.h.

{inDaqHeader, inDccHeader, inTccBlock, inSrBlock, inTowerBlock}

Constructor & Destructor Documentation

EcalDumpRaw::EcalDumpRaw ( const edm::ParameterSet & ps ) [explicit]

Definition at line 130 of file EcalDumpRaw.cc.

References beg_fed_id_, gather_cfg::cout, dumpFile_, end_fed_id_, filename_, first_event_, edm::ParameterSet::getUntrackedParameter(), last_event_, max(), verbosity_, and writeDcc_.

                                                 :
  iEvent_(0),
  adc_(nSamples, 0.),
  amplCut_(ps.getUntrackedParameter<double>("amplCut", 5.)),
  dump_(ps.getUntrackedParameter<bool>("dump", true)),
  dumpAdc_(ps.getUntrackedParameter<bool>("dumpAdc", true)),
  //  doHisto_(ps.getUntrackedParameter<bool>("doHisto", true)),
  maxEvt_(ps.getUntrackedParameter<int>("maxEvt", 10000)),
  profileFedId_(ps.getUntrackedParameter<int>("profileFedId", 0)),
  profileRuId_(ps.getUntrackedParameter<int>("profileRuId", 1)),
  l1aMinX_(ps.getUntrackedParameter<int>("l1aMinX", 1)),
  l1aMaxX_(ps.getUntrackedParameter<int>("l1aMaxX", 601)),
  lastOrbit_(nDccs_, numeric_limits<uint32_t>::max()),
  eventId_(numeric_limits<unsigned>::max()),
  eventList_(ps.getUntrackedParameter<vector<unsigned> >("eventList", vector<unsigned>())),
  minEventId_(999999),
  maxEventId_(0),
  orbit0_(0),
  orbit0Set_(false),
  bx_(-1),
  l1a_(-1),
  l1amin_(numeric_limits<int>::max()),
  l1amax_(-numeric_limits<int>::min()),
  simpleTrigType_(-1),
  detailedTrigType_(-1),
  //  histo_("hist.root", "RECREATE"),
  l1as_(36+2),
  orbits_(36+2),
  tpg_(maxTccsPerDcc_, std::vector<int>(maxTpgsPerTcc_)),
  nTpgs_(maxTccsPerDcc_, 0),
  srpL1a_(-1),
  tccL1a_(-1),
  nTts_(-1),
  tccBlockLen64_(19),
  feL1a_(nRu_,-1),
  srpBx_(-1),
  tccBx_(-1),
  tccType_(0),
  feBx_(nRu_,-1),
  feRuId_(nRu_,-1),
  iTow_(0),
  pulsePerRu_(ps.getUntrackedParameter<bool>("pulsePerRu", true)),
  pulsePerLmod_(ps.getUntrackedParameter<bool>("pulsePerLmod", true)),
  pulsePerLme_(ps.getUntrackedParameter<bool>("pulsePerLme", true)),
  tccId_(0)
{
  verbosity_= ps.getUntrackedParameter<int>("verbosity",1);

  beg_fed_id_= ps.getUntrackedParameter<int>("beg_fed_id",601);
  end_fed_id_= ps.getUntrackedParameter<int>("end_fed_id",654);


  first_event_ = ps.getUntrackedParameter<int>("first_event",1);
  last_event_  = ps.getUntrackedParameter<int>("last_event",
                                               numeric_limits<int>::max());

  writeDcc_ = ps.getUntrackedParameter<bool>("writeDCC",false);
  filename_  = ps.getUntrackedParameter<string>("filename","dump.bin");
  if(writeDcc_){
    dumpFile_.open(filename_.c_str());
    if(dumpFile_.bad()){
      /*edm::LogError("EcalDumpRaw")*/ std::cout << "Failed to open file '"
                               << filename_.c_str() << "' specified by "
                               << "parameter filename for writing. DCC data "
        " dump will be disabled.";
      writeDcc_ = false;
    }
  }
}

EcalDumpRaw::~EcalDumpRaw ( )

Definition at line 203 of file EcalDumpRaw.cc.

{
}

Member Function Documentation

void EcalDumpRaw::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	es
	)		`[virtual]`

Implements edm::EDAnalyzer.

Definition at line 208 of file EcalDumpRaw.cc.

References beg_fed_id_, bx_, ExpressReco_HICollisions_FallBack::cerr, gather_cfg::cout, FEDRawData::data(), runTheMatrix::data, dccId_, decode(), dump_, dumpFile_, end_fed_id_, eventId_, eventList_, feBx_, feBxOffset, feL1a_, feRuId_, spr::find(), first_event_, i, ExpressReco_HICollisions_FallBack::id, iEvent_, iRu_, iTcc_, iTow_, l1a_, l1amax_, l1amin_, last_event_, FEDNumbering::lastFEDId(), maxEventId_, minEventId_, nRu_, nTpgs_, nTts_, asciidump::s, simpleTrigType_, FEDRawData::size(), srpBx_, srpL1a_, tccBx_, tccL1a_, tccType_, toNth(), and writeDcc_.

                                                                  {
  ++iEvent_;
  eventId_ = event.id().event();

  if(eventList_.size()!=0 && find(eventList_.begin(), eventList_.end(),
                                  eventId_) == eventList_.end()){
    cout << "Skipping event " << eventId_ << ".\n";
    return;
  }
  
  if ((first_event_ > 0 && iEvent_ < first_event_) ||
      (last_event_ > 0 && last_event_ < iEvent_)) return;
  timeval start;
  timeval stop;
  gettimeofday(&start, 0);

  edm::Handle<FEDRawDataCollection> rawdata;
  event.getByType(rawdata);


  if(dump_) cout << "\n----------------------------------------------------------------------\n"
                 << toNth(iEvent_)
                 << " read event. "
                 << "Event id: "
                 << " " << eventId_
                 << "\n----------------------------------------------------------------------\n";

  if(eventId_ < minEventId_) minEventId_ = eventId_;
  if(eventId_ > maxEventId_) maxEventId_ = eventId_;

#if 1

  bool dccIdErr = false;
  bool simpleTrigTypeErr = false;
  bool outOfSync = false;
  unsigned iFed = 0;
  unsigned refDccId = 0;
  int refSimpleTrigType = -1;
  int refBx = -1;
  //  static bool recordNextPhys = false;
  //static int bxCalib = -1;
  //x static int nCalib = 0;

  for (int id = 0; id<=FEDNumbering::lastFEDId(); ++id){

    if (id < beg_fed_id_ || end_fed_id_ < id) continue;

    const FEDRawData& data = rawdata->FEDData(id);

    if (data.size()>4){
      ++iFed;
      if ((data.size() %16) !=0){
        cout << "***********************************************\n";
        cout << " Fed size in bits not multiple of 64, strange.\n";
        cout << "***********************************************\n";
      }


      size_t nWord32 = data.size()/4;
      const uint32_t * pData = ( reinterpret_cast<uint32_t*>(const_cast<unsigned char*> ( data.data())));
      stringstream s;
      srpL1a_ = -1;
      tccL1a_ = -1;
      srpBx_ = -1;
      tccBx_ = -1;
      iTow_ = 0;
      iRu_ = 0;
      nTts_ = -1;
      iTcc_ = 0;
      tccType_ = 0;

      for(int i = 0; i < nRu_; ++i){
        feL1a_[i]  = -1;
        feBx_[i]   = -1;
        feRuId_[i] = -1;
      }

      fill(nTpgs_.begin(), nTpgs_.end(), 0);

      bool rc;
      for(size_t iWord32=0; iWord32 < nWord32; iWord32+=2){
        s.str("");
        if(id>=601 && id<=654){// ECAL DCC data
          rc = decode(pData+iWord32, iWord32/2, s);
        } else{
          rc = true;
        }
        if(rc && dump_){
          cout << setfill('0') << hex
               << "[" << setw(8) << iWord32*4 << "] "
               <<        setw(4) << (pData[iWord32+1]>>16 & 0xFFFF) << " "
               <<        setw(4) << (pData[iWord32+1]>>0  & 0xFFFF) << " "
               <<        setw(4) << (pData[iWord32]>>16 & 0xFFFF) << " "
               <<        setw(4) << (pData[iWord32]>>0  & 0xFFFF) << " "
               << setfill(' ') << dec
               << s.str() << "\n";
        }
      }
      if(dump_) cout << "\n";

      if(iFed==1){
        refDccId = dccId_;
        refSimpleTrigType = simpleTrigType_;
        refBx = bx_;
      } else{
        if(dccId_!=refDccId){
          dccIdErr = true;
        }
        if(simpleTrigType_!=refSimpleTrigType){
          simpleTrigTypeErr = true;
        }
        if(refBx!=bx_){
          outOfSync = true;
        }
      }

      if(dump_) cout << flush; //flushing cout before writing to cerr

      if(srpBx_!=-1 && srpBx_!=bx_){
        cerr << "Bx discrepancy between SRP and DCC, Bx(SRP) = "
             << srpBx_ << ", Bx(DCC) = " << bx_
             << " in " << toNth(iEvent_) << " event, FED "
             << id << "\n";
      }

      if(tccBx_!=-1 && tccBx_!=bx_){
        cerr << "Bx discrepancy between TCC and DCC, Bx(TCC) = "
             << srpBx_ << ", Bx(DCC) = " << bx_
             << " in " << toNth(iEvent_) << " event, FED "
             << id << "\n";
      }

      bool feBxErr = false;
      for(int i=0; i < nRu_; ++i){
        int expectedFeBx;
        if(feBxOffset==0){
          expectedFeBx = bx_ - 1;
        } else{
          expectedFeBx = (bx_==3564) ? 0 : bx_;
        }
        if(feBx_[i]!=-1 && feBx_[i]!=bx_-1+feBxOffset) feBxErr = true;
      }
      if(feBxErr) cerr << "Bx discrepancy between DCC and at least one FE"
                       << " in " << toNth(iEvent_) << " event, FED "
                       << id << "\n";


      int localL1a = l1a_ & 0xFFF;
      if(srpL1a_!=-1 && srpL1a_!=localL1a){
        cerr << "Discrepancy between SRP and DCC L1a counter, L1a(SRP) = "
             << srpL1a_ << ", L1a(DCC) & 0xFFF = " << localL1a
             << " in " << toNth(iEvent_) << " event, FED "
             << id << "\n";

      }

      if(tccL1a_!=-1 && tccL1a_!=localL1a){
        cerr << "Discrepancy between TCC and DCC L1a counter, L1a(TCC) = "
             << srpL1a_ << ", L1a(DCC) & 0xFFF = " << localL1a
             << " in " << toNth(iEvent_) << " event, FED "
             << id << "\n";

      }

      bool feL1aErr = false;
      for(int i=0; i < nRu_; ++i){
        if(feL1a_[i]!=-1 && feL1a_[i]!=localL1a-1){
          cout << "FE L1A error for RU " << (i+1) << endl;
          feL1aErr = true;
        }
      }
      if(feL1aErr) cerr << "Discrepancy in L1a counter between DCC "
                     "and at least one FE (L1A(DCC) & 0xFFF = " << localL1a << ")"
                        << " in " << toNth(iEvent_) << " event, FED "
                        << id << "\n";
      

      if(iTow_>0 && iTow_< nRu_ && feRuId_[iTow_] < feRuId_[iTow_-1]){
        cerr << "Error in RU ID (TT/SC ID)"
             << " in " << toNth(iEvent_) << " event, FED "
             << id << "\n";
      }

      if (beg_fed_id_ <= id && id <= end_fed_id_ && writeDcc_){
        dumpFile_.write( reinterpret_cast <const char *> (pData), nWord32*4);
      }
    } else{
      //      cout << "No data for FED " <<  id << ". Size = "
      //     << data.size() << " byte(s).\n";
    }
  } //next fed

  if(dump_) cout << "Number of selected FEDs with a data block: "
                 << iFed << "\n";

  if(dccIdErr){
    cerr << "DCC ID discrepancy in detailed trigger type "
         << " of " << toNth(iEvent_) << " event.\n";
  }
  int bx = -1;
  if(!outOfSync){
    bx = bx_;
  }

  if(l1a_>0 && l1a_< l1amin_) l1amin_ = l1a_;
  if(l1a_>l1amax_) l1amax_ = l1a_;


#endif

  gettimeofday(&stop, 0);
  //  double dt  = (stop.tv_sec-start.tv_sec)*1.e3
  //  + (stop.tv_usec-start.tv_usec)*1.e-3;
  //  histo_.fillD("hCodeTime", "Code execution time;Duration (ms);Event count",
  //             PGXAxis(100, 0, 100),
  //             dt);
}

void EcalDumpRaw::analyzeApd ( ) [private]

void EcalDumpRaw::analyzeEB	(	const edm::Event &	,
		const edm::EventSetup &
	)		const

void EcalDumpRaw::analyzeEE	(	const edm::Event &	,
		const edm::EventSetup &
	)		const

void EcalDumpRaw::analyzeFed ( int fedId ) [private]

bool EcalDumpRaw::decode	(	const uint32_t *	data,
		int	iWord32,
		std::ostream &	out
	)		`[private]`

Referenced by analyze().

void EcalDumpRaw::endJob ( void ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 200 of file EcalDumpRaw.cc.

{
}

int EcalDumpRaw::lme	(	int	dccId1,
		int	side
	)		`[static, private]`

Definition at line 816 of file EcalDumpRaw.cc.

References min().

Referenced by ESElectronicsSimFast::decode().

                                      {
  int fedid = ((dcc1-1)%600) + 600; //to handle both FED and DCC id.
   vector<int> lmes;
   // EE -
   if( fedid <= 609 ) {
     if ( fedid <= 607 ) {
       lmes.push_back(fedid-601+83);
     } else if ( fedid == 608 ) {
       lmes.push_back(90);
       lmes.push_back(91);
     } else if ( fedid == 609 ) {
       lmes.push_back(92);
     }
   } //EB
   else if ( fedid >= 610  && fedid <= 645 ) {
     lmes.push_back(2*(fedid-610)+1);
     lmes.push_back(lmes[0]+1);
   } // EE+
   else if ( fedid >= 646 ) {
     if ( fedid <= 652 ) {
       lmes.push_back(fedid-646+73);
     } else if ( fedid == 653 ) {
       lmes.push_back(80);
       lmes.push_back(81);
     } else if ( fedid == 654 ) {
       lmes.push_back(82);
     }
   }
   return lmes.size()==0?-1:lmes[std::min(lmes.size(), (size_t)side)];
}

int EcalDumpRaw::lmodOfRu ( int ru1 ) [static, private]

Definition at line 866 of file EcalDumpRaw.cc.

                                {
  int iEta0 = (ru1-1)/4;
  int iPhi0 = (ru1-1)%4;
  int rs;
  if(iEta0==0){
    rs =  1;
  } else{
    rs = 2 + ((iEta0-1)/4)*2 + (iPhi0%4)/2;
  }
  //  cout << "ru1 = " << ru1 << " -> lmod = " << rs << "\n";
  return rs;
}

double EcalDumpRaw::max	(	std::vector< double >	a,
		unsigned &	pos
	)		`[inline, private]`

Definition at line 51 of file EcalDumpRaw.h.

References i, m, and pos.

Referenced by ESElectronicsSimFast::decode(), EcalDumpRaw(), and srRange().

                                                {
    pos = 0;
    double m = a[pos];
    for(unsigned i = 1; i < a.size(); ++i){
      if(a[i]>m){ m = a[i]; pos = i;}
    }
    return m;
  }

double EcalDumpRaw::min ( std::vector< double > a ) [inline, private]

Definition at line 59 of file EcalDumpRaw.h.

References i, and m.

Referenced by ESElectronicsSimFast::decode(), lme(), and srRange().

                                 {
    double m = a[0];
    for(unsigned i = 1; i < a.size(); ++i){
      if(a[i]<m) m = a[i];
    }
    return m;
  }

int EcalDumpRaw::modOfRu ( int ru1 ) [static, private]

Definition at line 857 of file EcalDumpRaw.cc.

                               {
  int iEta0 = (ru1-1)/4;
  if(iEta0<5){
    return 1;
  } else{
    return 2 + (iEta0-5)/4;
  }
}

int EcalDumpRaw::sideOfRu ( int ru1 ) [static, private]

Definition at line 848 of file EcalDumpRaw.cc.

Referenced by ESElectronicsSimFast::decode().

                                {
  if(ru1 < 5 || (ru1-5)%4 >= 2){
    return 0;
  } else{
    return 1;
  }
}

std::string EcalDumpRaw::srRange ( int offset ) const [private]

Definition at line 879 of file EcalDumpRaw.cc.

References fedId_, max(), min(), and asciidump::s.

Referenced by ESElectronicsSimFast::decode().

                                              {
  int min = offset+1;
  int max = offset+4;
  stringstream buf;
  if(628 <= fedId_ && fedId_ <= 646){//EB+
    buf << right << min << ".."
        << left  << max;
  } else{
    buf << right << max << ".."
        << left  << min;
  }
  string s = buf.str();
  buf.str("");
  buf << setw(6) << right << s;
  return buf.str();
}

string EcalDumpRaw::toNth ( int n ) [private]

Definition at line 426 of file EcalDumpRaw.cc.

References n, and asciidump::s.

Referenced by analyze().

                              {
  stringstream s;
  s << n;
  if(n%100<10 || n%100>20){
    switch(n%10){
    case 1:
      s << "st";
      break;
    case 2:
      s << "nd";
      break;
    case 3:
      s << "rd";
      break;
    default:
      s << "th";
    }
  } else{
    s << "th";
  }
  return s.str();
}

template<class T >

std::string EcalDumpRaw::toString ( T val ) [inline, private]

Definition at line 69 of file EcalDumpRaw.h.

References asciidump::s.

                           {
    std::stringstream s;
    s << val;
    return s.str();
  }

std::string EcalDumpRaw::tpgTag	(	int	tccType,
		unsigned	iSeq
	)		const `[private]`

Definition at line 913 of file EcalDumpRaw.cc.

References Exception, and ttId_.

Referenced by ESElectronicsSimFast::decode().

                                                             {
  if((unsigned)iSeq > sizeof(ttId_))
    throw cms::Exception("OutOfRange")
      << __FILE__ << ":"  << __LINE__ << ": "
      << "parameter out of range\n";
                             
  const int ttId = ttId_[tccType][iSeq];
  stringstream buf;
  buf.str("");
  if(ttId==0){
    buf << "    '0'";
  } else{
    buf << "TPG# " << setw(2) << ttId;
  }
  return buf.str();
}

std::string EcalDumpRaw::ttfTag	(	int	tccType,
		unsigned	iSeq
	)		const `[private]`

Definition at line 896 of file EcalDumpRaw.cc.

References Exception, and ttId_.

Referenced by ESElectronicsSimFast::decode().

                                                             {
  if((unsigned)iSeq > sizeof(ttId_))
    throw cms::Exception("OutOfRange")
      << __FILE__ << ":"  << __LINE__ << ": "
      << "parameter out of range\n";
                             
  const int ttId = ttId_[tccType][iSeq];
  stringstream buf;
  buf.str("");
  if(ttId==0){
    buf << "    '0'";
  } else{
    buf << "TTF# " << setw(2) << ttId;
  }
  return buf.str();
}