CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Types | Private Member Functions | Static Private Member Functions | Private Attributes | Static Private Attributes
EcalDumpRaw Class Reference

#include <EcalDumpRaw.h>

Inheritance diagram for EcalDumpRaw:
edm::EDAnalyzer edm::EDConsumerBase

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 ()
 
- Public Member Functions inherited from edm::EDAnalyzer
 EDAnalyzer ()
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndex indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

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 (const std::vector< double > &a, unsigned &pos)
 
double min (const 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 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_
 
std::vector< int > dccChStatus_
 
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_
 
edm::InputTag fedRawDataCollectionTag_
 
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_
 
edm::InputTag l1AcceptBunchCrossingCollectionTag_
 
bool l1aHistory_
 
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 More...
 
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
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
typedef WorkerT< EDAnalyzerWorkerType
 
- 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
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
CurrentProcessingContext const * currentContext () const
 
- 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

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 33 of file EcalDumpRaw.h.

Member Enumeration Documentation

anonymous enum
private
Enumerator
inDaqHeader 
inDccHeader 
inTccBlock 
inSrBlock 
inTowerBlock 

Definition at line 102 of file EcalDumpRaw.h.

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_.

130  :
131  iEvent_(0),
132  adc_(nSamples, 0.),
133  amplCut_(ps.getUntrackedParameter<double>("amplCut", 5.)),
134  dump_(ps.getUntrackedParameter<bool>("dump", true)),
135  dumpAdc_(ps.getUntrackedParameter<bool>("dumpAdc", true)),
136  l1aHistory_(ps.getUntrackedParameter<bool>("l1aHistory", true)),
137  // doHisto_(ps.getUntrackedParameter<bool>("doHisto", true)),
138  maxEvt_(ps.getUntrackedParameter<int>("maxEvt", 10000)),
139  profileFedId_(ps.getUntrackedParameter<int>("profileFedId", 0)),
140  profileRuId_(ps.getUntrackedParameter<int>("profileRuId", 1)),
141  l1aMinX_(ps.getUntrackedParameter<int>("l1aMinX", 1)),
142  l1aMaxX_(ps.getUntrackedParameter<int>("l1aMaxX", 601)),
145  eventList_(ps.getUntrackedParameter<vector<unsigned> >("eventList", vector<unsigned>())),
146  minEventId_(999999),
147  maxEventId_(0),
148  orbit0_(0),
149  orbit0Set_(false),
150  bx_(-1),
151  l1a_(-1),
154  simpleTrigType_(-1),
155  detailedTrigType_(-1),
156  // histo_("hist.root", "RECREATE"),
157  l1as_(36+2),
158  orbits_(36+2),
159  tpg_(maxTccsPerDcc_, std::vector<int>(maxTpgsPerTcc_)),
161  dccChStatus_(70, 0),
162  srpL1a_(-1),
163  tccL1a_(-1),
164  nTts_(-1),
165  tccBlockLen64_(19),
166  feL1a_(nRu_,-1),
167  srpBx_(-1),
168  tccBx_(-1),
169  tccType_(0),
170  feBx_(nRu_,-1),
171  feRuId_(nRu_,-1),
172  iTow_(0),
173  pulsePerRu_(ps.getUntrackedParameter<bool>("pulsePerRu", true)),
174  pulsePerLmod_(ps.getUntrackedParameter<bool>("pulsePerLmod", true)),
175  pulsePerLme_(ps.getUntrackedParameter<bool>("pulsePerLme", true)),
176  tccId_(0),
177  fedRawDataCollectionTag_(ps.getParameter<edm::InputTag>("fedRawDataCollectionTag")),
178  l1AcceptBunchCrossingCollectionTag_(ps.getParameter<edm::InputTag>("l1AcceptBunchCrossingCollectionTag"))
179 {
180  verbosity_= ps.getUntrackedParameter<int>("verbosity",1);
181 
182  beg_fed_id_= ps.getUntrackedParameter<int>("beg_fed_id",601);
183  end_fed_id_= ps.getUntrackedParameter<int>("end_fed_id",654);
184 
185 
186  first_event_ = ps.getUntrackedParameter<int>("first_event",1);
187  last_event_ = ps.getUntrackedParameter<int>("last_event",
189 
190  writeDcc_ = ps.getUntrackedParameter<bool>("writeDCC",false);
191  filename_ = ps.getUntrackedParameter<string>("filename","dump.bin");
192  if(writeDcc_){
193  dumpFile_.open(filename_.c_str());
194  if(dumpFile_.bad()){
195  /*edm::LogError("EcalDumpRaw")*/ std::cout << "Failed to open file '"
196  << filename_.c_str() << "' specified by "
197  << "parameter filename for writing. DCC data "
198  " dump will be disabled.";
199  writeDcc_ = false;
200  }
201  }
202 }
T getParameter(std::string const &) const
int tccBlockLen64_
Definition: EcalDumpRaw.h:169
T getUntrackedParameter(std::string const &, T const &) const
int end_fed_id_
Definition: EcalDumpRaw.h:93
bool pulsePerLme_
Definition: EcalDumpRaw.h:182
static const int nRu_
Definition: EcalDumpRaw.h:170
bool orbit0Set_
Definition: EcalDumpRaw.h:150
int detailedTrigType_
Definition: EcalDumpRaw.h:156
std::vector< double > adc_
Definition: EcalDumpRaw.h:106
#define min(a, b)
Definition: mlp_lapack.h:161
bool l1aHistory_
Definition: EcalDumpRaw.h:112
int first_event_
Definition: EcalDumpRaw.h:94
double amplCut_
Definition: EcalDumpRaw.h:109
unsigned minEventId_
Definition: EcalDumpRaw.h:146
bool writeDcc_
Definition: EcalDumpRaw.h:91
std::vector< uint32_t > lastOrbit_
Definition: EcalDumpRaw.h:120
bool pulsePerLmod_
Definition: EcalDumpRaw.h:181
std::vector< unsigned > eventList_
Definition: EcalDumpRaw.h:145
std::vector< int > feL1a_
Definition: EcalDumpRaw.h:171
std::vector< std::vector< uint32_t > > orbits_
Definition: EcalDumpRaw.h:159
std::vector< int > feBx_
Definition: EcalDumpRaw.h:176
const T & max(const T &a, const T &b)
unsigned maxEventId_
Definition: EcalDumpRaw.h:147
int profileRuId_
Definition: EcalDumpRaw.h:116
std::ofstream dumpFile_
Definition: EcalDumpRaw.h:179
edm::InputTag fedRawDataCollectionTag_
Definition: EcalDumpRaw.h:186
std::string filename_
Definition: EcalDumpRaw.h:96
std::vector< int > nTpgs_
Definition: EcalDumpRaw.h:161
static const int nSamples
Definition: EcalDumpRaw.h:108
std::vector< std::vector< int > > tpg_
Definition: EcalDumpRaw.h:160
bool pulsePerRu_
Definition: EcalDumpRaw.h:180
int tccType_
type of TCC currently parsed
Definition: EcalDumpRaw.h:175
int simpleTrigType_
Definition: EcalDumpRaw.h:155
std::vector< int > dccChStatus_
Definition: EcalDumpRaw.h:162
std::vector< int > feRuId_
Definition: EcalDumpRaw.h:177
int profileFedId_
Definition: EcalDumpRaw.h:115
std::vector< std::vector< uint32_t > > l1as_
Definition: EcalDumpRaw.h:158
static const int maxTpgsPerTcc_
Definition: EcalDumpRaw.h:123
edm::InputTag l1AcceptBunchCrossingCollectionTag_
Definition: EcalDumpRaw.h:187
unsigned eventId_
Definition: EcalDumpRaw.h:144
unsigned orbit0_
Definition: EcalDumpRaw.h:148
tuple cout
Definition: gather_cfg.py:121
int beg_fed_id_
Definition: EcalDumpRaw.h:92
static const unsigned nDccs_
Definition: EcalDumpRaw.h:121
static const int maxTccsPerDcc_
Definition: EcalDumpRaw.h:124
int verbosity_
Definition: EcalDumpRaw.h:90
int last_event_
Definition: EcalDumpRaw.h:95
EcalDumpRaw::~EcalDumpRaw ( )

Definition at line 207 of file EcalDumpRaw.cc.

207  {
208 }

Member Function Documentation

void EcalDumpRaw::analyze ( const edm::Event event,
const edm::EventSetup es 
)
virtual

Implements edm::EDAnalyzer.

Definition at line 212 of file EcalDumpRaw.cc.

References beg_fed_id_, bx_, dtNoiseDBValidation_cfg::cerr, gather_cfg::cout, FEDRawData::data(), data, dccChStatus_, dccId_, decode(), dump_, dumpFile_, end_fed_id_, eventId_, eventList_, feBx_, feBxOffset, fedRawDataCollectionTag_, feL1a_, feRuId_, lumiContext::fill, spr::find(), first_event_, i, iEvent_, iRu_, edm::HandleBase::isValid(), iTcc_, iTow_, l1a_, l1AcceptBunchCrossingCollectionTag_, l1aHistory_, l1amax_, l1amin_, last_event_, FEDNumbering::lastFEDId(), maxEventId_, minEventId_, nRu_, nTpgs_, nTts_, lumiPlot::rawdata, alignCSCRings::s, FEDRawData::size(), srpBx_, srpL1a_, dqm_diff::start, tccBx_, tccL1a_, tccType_, toNth(), trigNames, and writeDcc_.

212  {
213  ++iEvent_;
214  eventId_ = event.id().event();
215 
216  if(eventList_.size()!=0 && find(eventList_.begin(), eventList_.end(),
217  eventId_) == eventList_.end()){
218  cout << "Skipping event " << eventId_ << ".\n";
219  return;
220  }
221 
222  if ((first_event_ > 0 && iEvent_ < first_event_) ||
223  (last_event_ > 0 && last_event_ < iEvent_)) return;
224  timeval start;
225  timeval stop;
226  gettimeofday(&start, 0);
227 
229  event.getByLabel(fedRawDataCollectionTag_, rawdata);
230 
231  if(dump_ || l1aHistory_) cout << "\n======================================================================\n"
232  << toNth(iEvent_)
233  << " read event. "
234  << "Event id: "
235  << " " << eventId_
236  << "\n----------------------------------------------------------------------\n";
237 
238  if(l1aHistory_){
240  event.getByLabel(l1AcceptBunchCrossingCollectionTag_, l1aHist);
241  if(!l1aHist.isValid()) {
242  cout << "L1A history not found.\n";
243  } else if (l1aHist->size() == 0) {
244  cout << "L1A history is empty.\n";
245  } else{
246  cout << "L1A history: \n";
247  for(L1AcceptBunchCrossingCollection::const_iterator it = l1aHist->begin();
248  it != l1aHist->end();
249  ++it){
250  cout << "L1A offset: " << it->l1AcceptOffset() << "\t"
251  << "BX: " << it->bunchCrossing() << "\t"
252  << "Orbit ID: " << it->orbitNumber() << "\t"
253  << "Trigger type: " << it->eventType() << " ("
254  << trigNames[it->eventType()&0xF] << ")\n";
255  }
256  }
257  cout << "----------------------------------------------------------------------\n";
258  }
259 
260  if(eventId_ < minEventId_) minEventId_ = eventId_;
261  if(eventId_ > maxEventId_) maxEventId_ = eventId_;
262 
263 #if 1
264 
265  bool dccIdErr = false;
266  unsigned iFed = 0;
267  unsigned refDccId = 0;
268  // static bool recordNextPhys = false;
269  //static int bxCalib = -1;
270  //x static int nCalib = 0;
271 
272  for (int id = 0; id<=FEDNumbering::lastFEDId(); ++id){
273 
274  if (id < beg_fed_id_ || end_fed_id_ < id) continue;
275 
276  const FEDRawData& data = rawdata->FEDData(id);
277 
278  if (data.size()>4){
279  ++iFed;
280  if ((data.size() % 8) !=0){
281  cout << "***********************************************\n";
282  cout << " Fed size in bits not multiple of 64, strange.\n";
283  cout << "***********************************************\n";
284  }
285 
286 
287  size_t nWord32 = data.size()/4;
288  const uint32_t * pData = ( reinterpret_cast<uint32_t*>(const_cast<unsigned char*> ( data.data())));
289  stringstream s;
290  srpL1a_ = -1;
291  tccL1a_ = -1;
292  srpBx_ = -1;
293  tccBx_ = -1;
294  iTow_ = 0;
295  iRu_ = 0;
296  nTts_ = -1;
297  iTcc_ = 0;
298  tccType_ = 0;
299 
300  for(int i = 0; i < nRu_; ++i){
301  feL1a_[i] = -1;
302  feBx_[i] = -1;
303  feRuId_[i] = -1;
304  }
305 
306  fill(nTpgs_.begin(), nTpgs_.end(), 0);
307 
308  fill(dccChStatus_.begin(), dccChStatus_.end(), 0);
309 
310  bool rc;
311  for(size_t iWord32=0; iWord32 < nWord32; iWord32+=2){
312  s.str("");
313  if(id>=601 && id<=654){// ECAL DCC data
314  rc = decode(pData+iWord32, iWord32/2, s);
315  } else{
316  rc = true;
317  }
318  if(rc && dump_){
319  cout << setfill('0') << hex
320  << "[" << setw(8) << iWord32*4 << "] "
321  << setw(4) << (pData[iWord32+1]>>16 & 0xFFFF) << " "
322  << setw(4) << (pData[iWord32+1]>>0 & 0xFFFF) << " "
323  << setw(4) << (pData[iWord32]>>16 & 0xFFFF) << " "
324  << setw(4) << (pData[iWord32]>>0 & 0xFFFF) << " "
325  << setfill(' ') << dec
326  << s.str() << "\n";
327  }
328  }
329 
330  if(iFed==1){
331  refDccId = dccId_;
332  } else{
333  if(dccId_!=refDccId){
334  dccIdErr = true;
335  }
336  }
337 
338  if(dump_) cout << flush; //flushing cout before writing to cerr
339 
340  if(srpBx_!=-1 && srpBx_!=bx_){
341  cerr << "Bx discrepancy between SRP and DCC, Bx(SRP) = "
342  << srpBx_ << ", Bx(DCC) = " << bx_
343  << " in " << toNth(iEvent_) << " event, FED "
344  << id << endl;
345  }
346 
347  if(tccBx_!=-1 && tccBx_!=bx_){
348  cerr << "Bx discrepancy between TCC and DCC, Bx(TCC) = "
349  << srpBx_ << ", Bx(DCC) = " << bx_
350  << " in " << toNth(iEvent_) << " event, FED "
351  << id << endl;
352  }
353 
354  bool feBxErr = false;
355  for(int i=0; i < nRu_; ++i){
356  int expectedFeBx;
357  if(feBxOffset==0){
358  expectedFeBx = bx_ - 1;
359  } else{
360  expectedFeBx = (bx_==3564) ? 0 : bx_;
361  }
362  if(feBx_[i]!=-1 && feBx_[i]!=expectedFeBx){
363  cerr << "BX error for " << toNth(i+1) << " RU, RU ID "
364  << feRuId_[i];
365  if((unsigned) feRuId_[i] <= dccChStatus_.size()){
366  bool detected = (dccChStatus_[feRuId_[i]-1] == 10 || dccChStatus_[feRuId_[i]-1] == 11);
367  cerr << (detected?" ":" not ") << "detected by DCC (ch status: "
368  << dccChStatus_[feRuId_[i]-1] << ")";
369  }
370  cerr << " in " << toNth(iEvent_) << " event, FED "
371  << id << "." << endl;
372 
373  feBxErr = true;
374  }
375  }
376  if(feBxErr) cerr << "Bx discrepancy between DCC and at least one FE"
377  << " in " << toNth(iEvent_) << " event, FED "
378  << id << "\n";
379 
380 
381  int localL1a = l1a_ & 0xFFF;
382  if(srpL1a_!=-1 && srpL1a_!=localL1a){
383  cerr << "Discrepancy between SRP and DCC L1a counter, L1a(SRP) = "
384  << srpL1a_ << ", L1a(DCC) & 0xFFF = " << localL1a
385  << " in " << toNth(iEvent_) << " event, FED "
386  << id << endl;
387 
388  }
389 
390  if(tccL1a_!=-1 && tccL1a_!=localL1a){
391  cerr << "Discrepancy between TCC and DCC L1a counter, L1a(TCC) = "
392  << srpL1a_ << ", L1a(DCC) & 0xFFF = " << localL1a
393  << " in " << toNth(iEvent_) << " event, FED "
394  << id << endl;
395 
396  }
397 
398  bool feL1aErr = false;
399  for(int i=0; i < nRu_; ++i){
400  if(feL1a_[i] != -1 && feL1a_[i] != ((localL1a-1) & 0xFFF)){
401  cerr << "FE L1A error for " << toNth(i+1) << " RU, RU ID "
402  << feRuId_[i];
403  if((unsigned) feRuId_[i] <= dccChStatus_.size()){
404  bool detected = (dccChStatus_[feRuId_[i]-1] == 9 || dccChStatus_[feRuId_[i]-1] == 11);
405  cerr << (detected?" ":" not ") << "detected by DCC (ch status: "
406  << dccChStatus_[feRuId_[i]-1] << ")";
407  }
408  cerr << " in " << toNth(iEvent_) << " event, FED "
409  << id << "." << endl;
410  feL1aErr = true;
411  }
412  }
413  if(feL1aErr) cerr << "Discrepancy in L1a counter between DCC "
414  "and at least one FE (L1A(DCC) & 0xFFF = " << localL1a << ")"
415  << " in " << toNth(iEvent_) << " event, FED "
416  << id << "\n";
417 
418 
419  if(iTow_>0 && iTow_< nRu_ && feRuId_[iTow_] < feRuId_[iTow_-1]){
420  cerr << "Error in RU ID (TT/SC ID)"
421  << " in " << toNth(iEvent_) << " event, FED "
422  << id << endl;
423  }
424 
425  if (beg_fed_id_ <= id && id <= end_fed_id_ && writeDcc_){
426  dumpFile_.write( reinterpret_cast <const char *> (pData), nWord32*4);
427  }
428 
429  if(dump_) cout << "\n";
430  } else{
431  // cout << "No data for FED " << id << ". Size = "
432  // << data.size() << " byte(s).\n";
433  }
434  } //next fed
435 
436  if(dump_) cout << "Number of selected FEDs with a data block: "
437  << iFed << "\n";
438 
439  if(dccIdErr){
440  cout << flush;
441  cerr << "DCC ID discrepancy in detailed trigger type "
442  << " of " << toNth(iEvent_) << " event." << endl;
443  }
444 
445  if(l1a_>0 && l1a_< l1amin_) l1amin_ = l1a_;
446  if(l1a_>l1amax_) l1amax_ = l1a_;
447 
448 
449 #endif
450 
451  gettimeofday(&stop, 0);
452  // double dt = (stop.tv_sec-start.tv_sec)*1.e3
453  // + (stop.tv_usec-start.tv_usec)*1.e-3;
454  // histo_.fillD("hCodeTime", "Code execution time;Duration (ms);Event count",
455  // PGXAxis(100, 0, 100),
456  // dt);
457 }
int i
Definition: DBlmapReader.cc:9
std::string toNth(int n)
Definition: EcalDumpRaw.cc:459
tuple start
Check for commandline option errors.
Definition: dqm_diff.py:58
int end_fed_id_
Definition: EcalDumpRaw.h:93
string fill
Definition: lumiContext.py:319
static const int nRu_
Definition: EcalDumpRaw.h:170
bool l1aHistory_
Definition: EcalDumpRaw.h:112
int first_event_
Definition: EcalDumpRaw.h:94
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:49
unsigned minEventId_
Definition: EcalDumpRaw.h:146
bool writeDcc_
Definition: EcalDumpRaw.h:91
bool decode(const uint32_t *data, int iWord32, std::ostream &out)
Definition: EcalDumpRaw.cc:483
std::vector< unsigned > eventList_
Definition: EcalDumpRaw.h:145
std::vector< int > feL1a_
Definition: EcalDumpRaw.h:171
std::vector< int > feBx_
Definition: EcalDumpRaw.h:176
unsigned maxEventId_
Definition: EcalDumpRaw.h:147
std::ofstream dumpFile_
Definition: EcalDumpRaw.h:179
edm::InputTag fedRawDataCollectionTag_
Definition: EcalDumpRaw.h:186
std::vector< int > nTpgs_
Definition: EcalDumpRaw.h:161
bool isValid() const
Definition: HandleBase.h:76
static int lastFEDId()
Definition: FEDNumbering.cc:19
static const char *const trigNames[]
Definition: EcalDumpRaw.cc:75
int tccType_
type of TCC currently parsed
Definition: EcalDumpRaw.h:175
std::vector< int > dccChStatus_
Definition: EcalDumpRaw.h:162
std::vector< int > feRuId_
Definition: EcalDumpRaw.h:177
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
const int feBxOffset
Definition: EcalDumpRaw.cc:34
edm::InputTag l1AcceptBunchCrossingCollectionTag_
Definition: EcalDumpRaw.h:187
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:29
unsigned eventId_
Definition: EcalDumpRaw.h:144
tuple cout
Definition: gather_cfg.py:121
dictionary rawdata
Definition: lumiPlot.py:393
int beg_fed_id_
Definition: EcalDumpRaw.h:92
unsigned dccId_
Definition: EcalDumpRaw.h:142
int last_event_
Definition: EcalDumpRaw.h:95
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

Definition at line 483 of file EcalDumpRaw.cc.

References adc_, amplCut_, bx_, dtNoiseDBValidation_cfg::cerr, colorNames, gather_cfg::cout, pyrootRender::da, dtTPAnalyzer_cfg::dataType, dccCh_, dccChStatus_, dccId_, detailedTrigNames, detailedTrigType_, dump_, dumpAdc_, ebmTcc_, ebpTcc_, eeInnerTcc_, eeOuterTcc_, feBx_, fedId_, fedStart_, feL1a_, feRuId_, lumiContext::fill, g, i, iRu_, iSrWord64_, iTcc_, iTccWord64_, iTowerWord64_, l1a_, lastOrbit_, max(), max(), maxTccsPerDcc_, maxTpgsPerTcc_, mgpaGainFactors, min(), nDccs_, nRu_, nTpgs_, nTts_, orbit0_, orbit0Set_, orbit_, alignCSCRings::s, side_, sideOfRu(), simpleTrigType_, srpBx_, srpL1a_, srRange(), tccBlockLen64_, tccId_, tccL1a_, tccType_, towerBlockLength_, tpg_, tpgTag(), trigNames, ttfTag(), and ttsNames.

Referenced by BeautifulSoup.Tag::_invert(), analyze(), and BeautifulSoup.NavigableString::encode().

483  {
484  bool rc = true;
485  const bool d = dump_;
486  if(iWord64==0){//start of event
487  iSrWord64_ = 0;
488  iTccWord64_ = 0;
489  iTowerWord64_ = 0;
490  }
491  int dataType = (data[1] >>28) & 0xF;
492  const int boe = 5;
493  const int eoe = 10;
494  if(dataType==boe){//Begin of Event header
495  /**********************************************************************
496  * DAQ header
497  *
498  **********************************************************************/
499  simpleTrigType_ = (data[1] >>24) & 0xF;
500  l1a_ = (data[1]>>0 ) & 0xFFffFF;
501  bx_ = (data[0] >>20) & 0xFFF;
502  fedId_ = (data[0] >>8 ) & 0xFFF;
503  if(d) out << "Trigger type: " << simpleTrigType_
504  << "(" << trigNames[(data[1]>>24) & 0xF] << ")"
505  << " L1A: " << l1a_
506  << " BX: " << bx_
507  << " FED ID: " << fedId_
508  << " FOV: " << ((data[0] >>4 ) & 0xF)
509  << " H: " << ((data[0] >>3 ) & 0x1);
510  } else if((dataType>>2)==0){//DCC header
511  /**********************************************************************
512  * ECAL DCC header
513  *
514  **********************************************************************/
515  int dccHeaderId = (data[1] >>24) & 0x3F;
516  switch(dccHeaderId){
517  case 1:
518  if(d) out << "Run #: " << ((data[1] >>0 ) & 0xFFFFFF)
519  << " DCC Err: " << ((data[0] >>24) & 0xFF)
520  << " Evt Len: " << ((data[0] >>0 ) & 0xFFFFFF);
521  break;
522  case 2:
523  side_ = (data[1] >>11) & 0x1;
524  detailedTrigType_ = (data[1] >>8 ) & 0x7;
525  dccId_ = (data[1] >>0 ) & 0x3F;
526  if(d) out << "DCC FOV: " << ((data[1] >>16) & 0xF)
527  << " Side: " << side_
528  << " Trig.: " << detailedTrigType_
529  << " (" << detailedTrigNames[(data[1]>>8)&0x7] << ")"
530  << " Color: " << ((data[1] >>6 ) & 0x3)
531  << " (" << colorNames[(data[1]>>6)&0x3] << ")"
532  << " DCC ID: " << dccId_;
533  break;
534  case 3:
535  {
536  if(d) out << "TCC Status ch<4..1>: 0x"
537  << hex << ((data[1]>>8) & 0xFFFF) << dec
538  << " SR status: " << ((data[1] >>4 ) & 0xF)
539  << " TZS: " << ((data[1] >>2 ) & 0x1)
540  << " ZS: " << ((data[1] >>1 ) & 0x1)
541  << " SR: " << ((data[1] >>0 ) & 0x1);
542  orbit_ = data[0];
543  if(d) out << " Orbit: " << orbit_;
544  if(!orbit0Set_){
545  orbit0_ = orbit_;
546  orbit0Set_ = true;
547  }
548  int iDcc0 = fedId_-fedStart_;
549  if((unsigned)iDcc0<nDccs_){
551  if(d) out << " (+" << (int)orbit_-(int)lastOrbit_[iDcc0] <<")";
552  }
553  lastOrbit_[iDcc0] = orbit_;
554  }
555  }
556  break;
557  case 4:
558  case 5:
559  case 6:
560  case 7:
561  case 8:
562  {
563  int chOffset = (dccHeaderId-4)*14;
564  dccChStatus_[13+chOffset] = ((data[1] >>20) & 0xF);
565  dccChStatus_[12+chOffset] = ((data[1] >>16) & 0xF);
566  dccChStatus_[11+chOffset] = ((data[1] >>12) & 0xF);
567  dccChStatus_[10+chOffset] = ((data[1] >>8 ) & 0xF);
568  dccChStatus_[ 9+chOffset] = ((data[1] >>4 ) & 0xF);
569  dccChStatus_[ 8+chOffset] = ((data[1] >>0) & 0xF);
570  dccChStatus_[ 7+chOffset] = ((data[0] >>28) & 0xF);
571  dccChStatus_[ 6+chOffset] = ((data[0] >>24) & 0xF);
572  dccChStatus_[ 5+chOffset] = ((data[0] >>20) & 0xF);
573  dccChStatus_[ 4+chOffset] = ((data[0] >>16) & 0xF);
574  dccChStatus_[ 3+chOffset] = ((data[0] >>12) & 0xF);
575  dccChStatus_[ 2+chOffset] = ((data[0] >>8 ) & 0xF);
576  dccChStatus_[ 1+chOffset] = ((data[0] >>4 ) & 0xF);
577  dccChStatus_[ 0+chOffset] = ((data[0] >>0 ) & 0xF);
578 
579  if(d){
580  out << "FE CH status:";
581  for(int i = chOffset; i < chOffset + 14; ++i){
582  out << " #" << (i+1) << ":" << dccChStatus_[i];
583  }
584  }
585  }
586  break;
587  default:
588  if(d) out << " bits<63..62>=0 (DCC header) bits<61..56>=" << dccHeaderId
589  << "(unknown=>ERROR?)";
590  }
591  } else if((dataType>>1)==3){//TCC block
592  /**********************************************************************
593  * TCC block
594  *
595  **********************************************************************/
596  if(iTccWord64_==0){
597  //header
598  tccL1a_ = (data[1] >>0 ) & 0xFFF;
599  tccId_ = ((data[0] >>0 ) & 0xFF);
600  nTts_ = ((data[1] >>16) & 0x7F);
602  ++iTcc_;
603  if(d) out << "LE1: " << ((data[1] >>28) & 0x1)
604  << " LE0: " << ((data[1] >>27) & 0x1)
605  << " N_samples: " << ((data[1] >>23) & 0x1F)
606  << " N_TTs: " << nTts_
607  << " E1: " << ((data[1] >>12) & 0x1)
608  << " L1A: " << tccL1a_
609  << " '3': " << ((data[0] >>29) & 0x7)
610  << " E0: " << ((data[0] >>28) & 0x1)
611  << " Bx: " << ((data[0] >>16) & 0xFFF)
612  << " TTC ID: " << tccId_;
613  if(nTts_==68){ //EB TCC (TCC68)
614  if(fedId_ < 628) tccType_ = ebmTcc_;
615  else tccType_ = ebpTcc_;
616  } else if(nTts_ == 16){//Inner EE TCC (TCC48)
618  } else if(nTts_ == 28){//Outer EE TCC (TCC48)
620  } else {
621  cout << flush;
622  cerr << "Error in #TT field of TCC block."
623  "This field is normally used to determine type of TCC "
624  "(TCC48 or TCC68). Type of TCC will be deduced from the TCC ID.\n";
625  if(tccId_ < 19) tccType_ = eeInnerTcc_;
626  else if(tccId_ < 37) tccType_ = eeOuterTcc_;
627  else if(tccId_ < 55) tccType_ = ebmTcc_;
628  else if(tccId_ < 73) tccType_ = ebpTcc_;
629  else if(tccId_ < 91) tccType_ = eeOuterTcc_;
630  else if(tccId_ < 109) tccType_ = eeInnerTcc_;
631  else{
632  cerr << "TCC ID is also invalid. EB- TCC type will be assumed.\n";
633  tccType_ = ebmTcc_;
634  }
635  cerr << flush;
636  }
637  tccBlockLen64_ = (tccType_==ebmTcc_ || tccType_==ebpTcc_) ? 18 : 9;
638  } else{// if(iTccWord64_<18){
639  int tpgOffset = (iTccWord64_-1)*4;
640  if(iTcc_ > maxTccsPerDcc_){
641  out << "Too many TCC blocks";
642  } else if(tpgOffset > (maxTpgsPerTcc_ - 4)){
643  out << "Too many TPG in one TCC block";
644  } else{
645  tpg_[iTcc_-1][3+tpgOffset] = (data[1] >>16) & 0x1FF;
646  tpg_[iTcc_-1][2+tpgOffset] = (data[1] >>0 ) & 0x1FF;
647  tpg_[iTcc_-1][1+tpgOffset] = (data[0] >>16) & 0x1FF;
648  tpg_[iTcc_-1][0+tpgOffset] = (data[0] >>0 ) & 0x1FF;
649  //int n[2][4] = {{1,2,3,4},
650  // {4,3,2,1}};
651  //int iorder = (628<=fedId_ && fedId_<=645)?1:0;
652  if(d) out << ttfTag(tccType_, 3+tpgOffset) << ":" //"TTF# " << setw(2) << ttId_[3 + tpgOffset] << ":"
653  << ((data[1] >>25) & 0x7) << " "
654  << tpgTag(tccType_, 3+tpgOffset) << ":" //" TPG# "<< setw(2) << ttId_[3 + tpgOffset] << ":"
655  << setw(3) << tpg_[iTcc_-1][3+tpgOffset] << " "
656  << ttfTag(tccType_, 2+tpgOffset) << ":" //" TTF# "<< setw(2) << ttId_[2 + tpgOffset] << ":"
657  << ((data[1] >>9 ) & 0x7) << " "
658  << tpgTag(tccType_, 2+tpgOffset) << ":" //" TPG# "<< setw(2) << ttId_[2 + tpgOffset] << ":"
659  << setw(3) << tpg_[iTcc_-1][2+tpgOffset] << " "
660  << " '3': " << ((data[0] >>29) & 0x7) << " "
661  << ttfTag(tccType_, 1+tpgOffset) << ":" //" TTF# "<< setw(2) << ttId_[1 + tpgOffset] << ":"
662  << ((data[0] >>25) & 0x7) << " "
663  << setw(3) << tpgTag(tccType_, 1+tpgOffset) << ": "//" TPG# "<< setw(2) << ttId_[1 + tpgOffset] << ":"
664  << tpg_[iTcc_-1][1+tpgOffset] << " "
665  << ttfTag(tccType_, 0+tpgOffset) << ":" //" TTF# "<< setw(2) << ttId_[0 + tpgOffset] << ":"
666  << ((data[0] >>9 ) & 0x7) << " "
667  << setw(3) << tpgTag(tccType_, 0+tpgOffset) << ":" //" TPG# "<< setw(2) << ttId_[0 + tpgOffset] << ":"
668  << tpg_[iTcc_-1][0+tpgOffset];
669  }
670  }// else{
671  // if(d) out << "ERROR";
672  //}
673  ++iTccWord64_;
674  if(iTccWord64_ >= (unsigned)tccBlockLen64_) iTccWord64_ = 0;
675  } else if((dataType>>1)==4){//SRP block
676  /**********************************************************************
677  * SRP block
678  *
679  **********************************************************************/
680  if(iSrWord64_==0){//header
681  srpL1a_ = (data[1] >>0 ) & 0xFFF;
682  srpBx_ = (data[0] >>16) & 0xFFF;
683  if(d) out << "LE1: " << ((data[1] >>28) & 0x1)
684  << " LE0: " << ((data[1] >>27) & 0x1)
685  << " N_SRFs: " << ((data[1] >>16) & 0x7F)
686  << " E1: " << ((data[1] >>12) & 0x1)
687  << " L1A: " << srpL1a_
688  << " '4': " << ((data[0] >>29) & 0x7)
689  << " E0: " << ((data[0] >>28) & 0x1)
690  << " Bx: " << srpBx_
691  << " SRP ID: " << ((data[0] >>0 ) & 0xFF);
692  } else if(iSrWord64_<6){
693  int ttfOffset = (iSrWord64_-1)*16;
694  if(d){
695  if(iSrWord64_<5){
696  out <<"SRF# " << setw(6) << right << srRange(12+ttfOffset)/*16+ttfOffset << "..#" << 13+ttfOffset*/ << ": "
697  << oct << ((data[1] >>16) & 0xFFF) << dec
698  << " SRF# " << srRange(8+ttfOffset) /*12+ttfOffset << "..#" << 9+ttfOffset*/ << ": "
699  << oct << ((data[1] >>0 ) & 0xFFF) << dec
700  << " '4':" << ((data[0] >>29) & 0x7)
701  << " SRF# " << srRange(4+ttfOffset) /*8+ttfOffset << "..#" << 5+ttfOffset*/ << ": "
702  << oct << ((data[0] >>16) & 0xFFF) << dec;
703  } else{//last 64-bit word has only 4 SRFs.
704  out << " ";
705  }
706  out << " SRF# " << srRange(ttfOffset) /*4+ttfOffset << "..#" << 1+ttfOffset*/ << ": "
707  << oct << ((data[0] >>0 ) & 0xFFF) << dec;
708  }
709  } else{
710  if(d) out << "ERROR";
711  }
712  ++iSrWord64_;
713  } else if((dataType>>2)==3){//Tower block
714  /**********************************************************************
715  * "Tower" block (crystal channel data from a RU (=1 FE cards))
716  *
717  **********************************************************************/
718  if(iTowerWord64_==0){//header
719  towerBlockLength_ = (data[1]>>16) & 0x1FF;
720  int l1a;
721  int bx;
722  l1a = (data[1] >>0 ) & 0xFFF;
723  bx = (data[0] >>16) & 0xFFF;
724  dccCh_=(data[0] >>0 ) & 0xFF;
725  if(d) out << "Block Len: " << towerBlockLength_
726  << " E1: " << ((data[1] >>12) & 0x1)
727  << " L1A: " << l1a
728  << " '3': " << ((data[0] >>30) & 0x3)
729  << " E0: " << ((data[0] >>28) & 0x1)
730  << " Bx: " << bx
731  << " N_samples: " << ((data[0] >>8 ) & 0x7F)
732  << " RU ID: " << dccCh_;
733  if(iRu_ < nRu_){
734  feL1a_[iRu_] = l1a;
735  feBx_[iRu_] = bx;
736  feRuId_[iRu_] = dccCh_;
737  ++iRu_;
738  }
739  } else if((unsigned)iTowerWord64_<towerBlockLength_){
740  if(!dumpAdc_){
741  //no output.
742  rc = false;
743  }
744  const bool da = dumpAdc_ && dump_;
745  switch((iTowerWord64_-1)%3){
746  int s[4];
747  int g[4];
748  case 0:
749  s[0]=(data[0] >>16) & 0xFFF;
750  g[0]=(data[0] >>28) & 0x3;
751  s[1]=(data[1] >>0 ) & 0xFFF;
752  g[1]=(data[1] >>12) & 0x3;
753  s[2]=(data[1] >>16) & 0xFFF;
754  g[2]=(data[1] >>28) & 0x3;
755  fill(adc_.begin(), adc_.end(), 0.);
756  if(da) out << "GMF: " << ((data[0] >>11) & 0x1)
757  << " SMF: " << ((data[0] >>9 ) & 0x1)
758  << " M: " << ((data[0] >>8 ) & 0x1)
759  << " XTAL: " << ((data[0] >>4 ) & 0x7)
760  << " STRIP: " << ((data[0] >>0 ) & 0x7)
761  << " " << setw(4) << s[0]
762  << "G" << g[0]
763  << " " << setw(4) << s[1]
764  << "G" << g[1]
765  << " " << setw(4) << s[2]
766  << "G" << g[2];
767  for(int i=0; i<3; ++i) adc_[i] = s[i]*mgpaGainFactors[g[i]];
768  break;
769  case 1:
770  s[0]=(data[0] >>0 ) & 0xFFF;
771  g[0]=(data[0] >>12) & 0x3;
772  s[1]=(data[0] >>16) & 0xFFF;
773  g[1]=(data[0] >>28) & 0x3;
774  s[2]=(data[1] >>0 ) & 0xFFF;
775  g[2]=(data[1] >>12) & 0x3;
776  s[3]=(data[1] >>16) & 0xFFF;
777  g[3]=(data[1] >>28) & 0x3;
778  if(da) out << " "
779  << " " << setw(4) << s[0]
780  << "G" << g[0]
781  << " " << setw(4) << s[1]
782  << "G" << g[1]
783  << " " << setw(4) << s[2]
784  << "G" << g[2]
785  << " " << setw(4) << s[3]
786  << "G" << g[3];
787  for(int i=0; i<4; ++i) adc_[i+3] = s[i]*mgpaGainFactors[g[i]];
788  break;
789  case 2:
790  if(da) out << "TZS: " << ((data[1] >>14) & 0x1);
791 
792  s[0]=(data[0] >>0 ) & 0xFFF;
793  g[0]=(data[0] >>12) & 0x3;
794  s[1]=(data[0] >>16) & 0xFFF;
795  g[1]=(data[0] >>28) & 0x3;
796  s[2]=(data[1] >>0 ) & 0xFFF;
797  g[2]=(data[1] >>12) & 0x3 ;
798 
799  for(int i=0; i<3; ++i) adc_[i+7] = s[i]*mgpaGainFactors[g[i]];
800  if(dccCh_<=68){
801  unsigned bom0; //Bin of Maximum, starting counting from 0
802  double ampl = max(adc_, bom0)-min(adc_);
803  if(da) out << " Ampl: " << setw(4) << ampl
804  << (ampl>amplCut_?"*":" ")
805  << " BoM:" << setw(2) << (bom0+1)
806  << " ";
807  if(fedId_ == dccId_ + 600 //block of the read-out SM
808  //if laser, only one side:
809  && (detailedTrigType_!=4 || sideOfRu(dccCh_)==(int)side_)
810  ){
811  }
812  } else{
813  if(da) out << setw(29) << "";
814  }
815  if(da) out << " " << setw(4) << s[0]
816  << "G" << g[0]
817  << " " << setw(4) << s[1]
818  << "G" << g[1]
819  << " " << setw(4) << s[2]
820  << "G" << g[2];
821  break;
822  default:
823  assert(false);
824  }
825  } else {
826  if(d) out << "ERROR";
827  }
828  ++iTowerWord64_;
831  ++dccCh_;
832  }
833  } else if(dataType==eoe){//End of event trailer
834  /**********************************************************************
835  * Event DAQ trailer
836  *
837  **********************************************************************/
838  int tts = (data[0] >>4) & 0xF;
839  if(d) out << "Evt Len.: " << ((data[1] >>0 ) & 0xFFFFFF)
840  << " CRC16: " << ((data[0] >>16) & 0xFFFF)
841  << " Evt Status: " << ((data[0] >>8 ) & 0xF)
842  << " TTS: " << tts
843  << " (" << ttsNames[tts] << ")"
844  << " T:" << ((data[0] >>3) & 0x1);
845  } else{
846  if(d) out << " incorrect 64-bit word type marker (see MSBs)";
847  }
848  return rc;
849 }
int tccBlockLen64_
Definition: EcalDumpRaw.h:169
static int sideOfRu(int ru1)
Definition: EcalDumpRaw.cc:886
int i
Definition: DBlmapReader.cc:9
std::string tpgTag(int tccType, unsigned iSeq) const
Definition: EcalDumpRaw.cc:951
string fill
Definition: lumiContext.py:319
unsigned iTowerWord64_
Definition: EcalDumpRaw.h:99
static const int nRu_
Definition: EcalDumpRaw.h:170
static const unsigned fedStart_
Definition: EcalDumpRaw.h:122
bool orbit0Set_
Definition: EcalDumpRaw.h:150
size_t towerBlockLength_
Definition: EcalDumpRaw.h:104
std::string ttfTag(int tccType, unsigned iSeq) const
Definition: EcalDumpRaw.cc:934
int detailedTrigType_
Definition: EcalDumpRaw.h:156
double min(const std::vector< double > &a)
Definition: EcalDumpRaw.h:60
std::vector< double > adc_
Definition: EcalDumpRaw.h:106
double amplCut_
Definition: EcalDumpRaw.h:109
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e g
Definition: Activities.doc:4
std::vector< uint32_t > lastOrbit_
Definition: EcalDumpRaw.h:120
unsigned side_
Definition: EcalDumpRaw.h:143
std::vector< int > feL1a_
Definition: EcalDumpRaw.h:171
static const char *const ttsNames[]
Definition: EcalDumpRaw.cc:104
static const int eeOuterTcc_
Definition: EcalDumpRaw.h:133
uint32_t orbit_
Definition: EcalDumpRaw.h:149
std::vector< int > feBx_
Definition: EcalDumpRaw.h:176
const T & max(const T &a, const T &b)
static const int ebmTcc_
Definition: EcalDumpRaw.h:130
unsigned iTccWord64_
Definition: EcalDumpRaw.h:101
std::vector< int > nTpgs_
Definition: EcalDumpRaw.h:161
static const char *const detailedTrigNames[]
Definition: EcalDumpRaw.cc:86
unsigned iSrWord64_
Definition: EcalDumpRaw.h:100
std::vector< std::vector< int > > tpg_
Definition: EcalDumpRaw.h:160
tuple out
Definition: dbtoconf.py:99
unsigned fedId_
Definition: EcalDumpRaw.h:141
static const char *const trigNames[]
Definition: EcalDumpRaw.cc:75
std::string srRange(int offset) const
Definition: EcalDumpRaw.cc:917
int tccType_
type of TCC currently parsed
Definition: EcalDumpRaw.h:175
int simpleTrigType_
Definition: EcalDumpRaw.h:155
std::vector< int > dccChStatus_
Definition: EcalDumpRaw.h:162
static const char *const colorNames[]
Definition: EcalDumpRaw.cc:97
static const int eeInnerTcc_
Definition: EcalDumpRaw.h:132
std::vector< int > feRuId_
Definition: EcalDumpRaw.h:177
static const int maxTpgsPerTcc_
Definition: EcalDumpRaw.h:123
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
unsigned orbit0_
Definition: EcalDumpRaw.h:148
tuple cout
Definition: gather_cfg.py:121
static const unsigned nDccs_
Definition: EcalDumpRaw.h:121
static const int maxTccsPerDcc_
Definition: EcalDumpRaw.h:124
double mgpaGainFactors[]
Definition: EcalDumpRaw.cc:127
double max(const std::vector< double > &a, unsigned &pos)
Definition: EcalDumpRaw.h:52
unsigned dccId_
Definition: EcalDumpRaw.h:142
static const int ebpTcc_
Definition: EcalDumpRaw.h:131
void EcalDumpRaw::endJob ( void  )
virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 204 of file EcalDumpRaw.cc.

204  {
205 }
int EcalDumpRaw::lmodOfRu ( int  ru1)
staticprivate

Definition at line 904 of file EcalDumpRaw.cc.

904  {
905  int iEta0 = (ru1-1)/4;
906  int iPhi0 = (ru1-1)%4;
907  int rs;
908  if(iEta0==0){
909  rs = 1;
910  } else{
911  rs = 2 + ((iEta0-1)/4)*2 + (iPhi0%4)/2;
912  }
913  // cout << "ru1 = " << ru1 << " -> lmod = " << rs << "\n";
914  return rs;
915 }
double EcalDumpRaw::max ( const std::vector< double > &  a,
unsigned &  pos 
)
inlineprivate

Definition at line 52 of file EcalDumpRaw.h.

References i, m, and pos.

Referenced by decode(), and srRange().

52  {
53  pos = 0;
54  double m = a[pos];
55  for(unsigned i = 1; i < a.size(); ++i){
56  if(a[i]>m){ m = a[i]; pos = i;}
57  }
58  return m;
59  }
int i
Definition: DBlmapReader.cc:9
double a
Definition: hdecay.h:121
double EcalDumpRaw::min ( const std::vector< double > &  a)
inlineprivate

Definition at line 60 of file EcalDumpRaw.h.

References i, and m.

Referenced by decode(), and srRange().

60  {
61  double m = a[0];
62  for(unsigned i = 1; i < a.size(); ++i){
63  if(a[i]<m) m = a[i];
64  }
65  return m;
66  }
int i
Definition: DBlmapReader.cc:9
double a
Definition: hdecay.h:121
int EcalDumpRaw::modOfRu ( int  ru1)
staticprivate

Definition at line 895 of file EcalDumpRaw.cc.

895  {
896  int iEta0 = (ru1-1)/4;
897  if(iEta0<5){
898  return 1;
899  } else{
900  return 2 + (iEta0-5)/4;
901  }
902 }
int EcalDumpRaw::sideOfRu ( int  ru1)
staticprivate

Definition at line 886 of file EcalDumpRaw.cc.

Referenced by decode().

886  {
887  if(ru1 < 5 || (ru1-5)%4 >= 2){
888  return 0;
889  } else{
890  return 1;
891  }
892 }
std::string EcalDumpRaw::srRange ( int  offset) const
private

Definition at line 917 of file EcalDumpRaw.cc.

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

Referenced by decode().

917  {
918  int min = offset+1;
919  int max = offset+4;
920  stringstream buf;
921  if(628 <= fedId_ && fedId_ <= 646){//EB+
922  buf << right << min << ".."
923  << left << max;
924  } else{
925  buf << right << max << ".."
926  << left << min;
927  }
928  string s = buf.str();
929  buf.str("");
930  buf << setw(6) << right << s;
931  return buf.str();
932 }
double min(const std::vector< double > &a)
Definition: EcalDumpRaw.h:60
unsigned int offset(bool)
unsigned fedId_
Definition: EcalDumpRaw.h:141
double max(const std::vector< double > &a, unsigned &pos)
Definition: EcalDumpRaw.h:52
string EcalDumpRaw::toNth ( int  n)
private

Definition at line 459 of file EcalDumpRaw.cc.

References n, and alignCSCRings::s.

Referenced by analyze().

459  {
460  stringstream s;
461  s << n;
462  if(n%100<10 || n%100>20){
463  switch(n%10){
464  case 1:
465  s << "st";
466  break;
467  case 2:
468  s << "nd";
469  break;
470  case 3:
471  s << "rd";
472  break;
473  default:
474  s << "th";
475  }
476  } else{
477  s << "th";
478  }
479  return s.str();
480 }
template<class T >
std::string EcalDumpRaw::toString ( T  val)
inlineprivate

Definition at line 70 of file EcalDumpRaw.h.

References alignCSCRings::s.

70  {
71  std::stringstream s;
72  s << val;
73  return s.str();
74  }
std::string EcalDumpRaw::tpgTag ( int  tccType,
unsigned  iSeq 
) const
private

Definition at line 951 of file EcalDumpRaw.cc.

References edm::hlt::Exception, ecaldqm::ttId(), and ttId_.

Referenced by decode().

951  {
952  if((unsigned)iSeq > sizeof(ttId_))
953  throw cms::Exception("OutOfRange")
954  << __FILE__ << ":" << __LINE__ << ": "
955  << "parameter out of range\n";
956 
957  const int ttId = ttId_[tccType][iSeq];
958  stringstream buf;
959  buf.str("");
960  if(ttId==0){
961  buf << " '0'";
962  } else{
963  buf << "TPG# " << setw(2) << ttId;
964  }
965  return buf.str();
966 }
static const int ttId_[nTccTypes_][maxTpgsPerTcc_]
Definition: EcalDumpRaw.h:139
unsigned ttId(const DetId &)
std::string EcalDumpRaw::ttfTag ( int  tccType,
unsigned  iSeq 
) const
private

Definition at line 934 of file EcalDumpRaw.cc.

References edm::hlt::Exception, ecaldqm::ttId(), and ttId_.

Referenced by decode().

934  {
935  if((unsigned)iSeq > sizeof(ttId_))
936  throw cms::Exception("OutOfRange")
937  << __FILE__ << ":" << __LINE__ << ": "
938  << "parameter out of range\n";
939 
940  const int ttId = ttId_[tccType][iSeq];
941  stringstream buf;
942  buf.str("");
943  if(ttId==0){
944  buf << " '0'";
945  } else{
946  buf << "TTF# " << setw(2) << ttId;
947  }
948  return buf.str();
949 }
static const int ttId_[nTccTypes_][maxTpgsPerTcc_]
Definition: EcalDumpRaw.h:139
unsigned ttId(const DetId &)

Member Data Documentation

std::vector<double> EcalDumpRaw::adc_
private

Definition at line 106 of file EcalDumpRaw.h.

Referenced by decode().

double EcalDumpRaw::amplCut_
private

Definition at line 109 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::beg_fed_id_
private

Definition at line 92 of file EcalDumpRaw.h.

Referenced by analyze(), and EcalDumpRaw().

int EcalDumpRaw::bx_
private

Definition at line 151 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

int EcalDumpRaw::dccCh_
private

Definition at line 119 of file EcalDumpRaw.h.

Referenced by decode().

std::vector<int> EcalDumpRaw::dccChStatus_
private

Definition at line 162 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

unsigned EcalDumpRaw::dccId_
private

Definition at line 142 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

enum { ... } EcalDumpRaw::decodeState_
int EcalDumpRaw::detailedTrigType_
private

Definition at line 156 of file EcalDumpRaw.h.

Referenced by decode().

bool EcalDumpRaw::dump_
private

Definition at line 110 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

bool EcalDumpRaw::dumpAdc_
private

Definition at line 111 of file EcalDumpRaw.h.

Referenced by decode().

std::ofstream EcalDumpRaw::dumpFile_
private

Definition at line 179 of file EcalDumpRaw.h.

Referenced by analyze(), and EcalDumpRaw().

const int EcalDumpRaw::ebmTcc_ = 0
staticprivate

TCC types

Definition at line 130 of file EcalDumpRaw.h.

Referenced by decode().

const int EcalDumpRaw::ebpTcc_ = 1
staticprivate

Definition at line 131 of file EcalDumpRaw.h.

Referenced by decode().

const int EcalDumpRaw::eeInnerTcc_ = 2
staticprivate

Definition at line 132 of file EcalDumpRaw.h.

Referenced by decode().

const int EcalDumpRaw::eeOuterTcc_ = 3
staticprivate

Definition at line 133 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::end_fed_id_
private

Definition at line 93 of file EcalDumpRaw.h.

Referenced by analyze(), and EcalDumpRaw().

unsigned EcalDumpRaw::eventId_
private

Definition at line 144 of file EcalDumpRaw.h.

Referenced by analyze().

std::vector<unsigned> EcalDumpRaw::eventList_
private

Definition at line 145 of file EcalDumpRaw.h.

Referenced by analyze().

std::vector<int> EcalDumpRaw::feBx_
private

Definition at line 176 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

unsigned EcalDumpRaw::fedId_
private

Definition at line 141 of file EcalDumpRaw.h.

Referenced by decode(), and srRange().

edm::InputTag EcalDumpRaw::fedRawDataCollectionTag_
private

Definition at line 186 of file EcalDumpRaw.h.

Referenced by analyze().

const unsigned EcalDumpRaw::fedStart_ = 601
staticprivate

Definition at line 122 of file EcalDumpRaw.h.

Referenced by decode().

std::vector<int> EcalDumpRaw::feL1a_
private

Definition at line 171 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

std::vector<int> EcalDumpRaw::feRuId_
private

Definition at line 177 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

std::string EcalDumpRaw::filename_
private

Definition at line 96 of file EcalDumpRaw.h.

Referenced by EcalDumpRaw().

int EcalDumpRaw::first_event_
private

Definition at line 94 of file EcalDumpRaw.h.

Referenced by analyze(), and EcalDumpRaw().

int EcalDumpRaw::iEvent_
private

Definition at line 97 of file EcalDumpRaw.h.

Referenced by analyze().

int EcalDumpRaw::iRu_
private

Definition at line 163 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

unsigned EcalDumpRaw::iSrWord64_
private

Definition at line 100 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::iTcc_
private

Definition at line 185 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

unsigned EcalDumpRaw::iTccWord64_
private

Definition at line 101 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::iTow_
private

Definition at line 178 of file EcalDumpRaw.h.

Referenced by analyze().

unsigned EcalDumpRaw::iTowerWord64_
private

Definition at line 99 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::l1a_
private

Definition at line 152 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

edm::InputTag EcalDumpRaw::l1AcceptBunchCrossingCollectionTag_
private

Definition at line 187 of file EcalDumpRaw.h.

Referenced by analyze().

bool EcalDumpRaw::l1aHistory_
private

Definition at line 112 of file EcalDumpRaw.h.

Referenced by analyze().

int EcalDumpRaw::l1amax_
private

Definition at line 154 of file EcalDumpRaw.h.

Referenced by analyze().

int EcalDumpRaw::l1aMaxX_
private

Definition at line 118 of file EcalDumpRaw.h.

int EcalDumpRaw::l1amin_
private

Definition at line 153 of file EcalDumpRaw.h.

Referenced by analyze().

int EcalDumpRaw::l1aMinX_
private

Definition at line 117 of file EcalDumpRaw.h.

std::vector<std::vector<uint32_t> > EcalDumpRaw::l1as_
private

Definition at line 158 of file EcalDumpRaw.h.

int EcalDumpRaw::last_event_
private

Definition at line 95 of file EcalDumpRaw.h.

Referenced by analyze(), and EcalDumpRaw().

std::vector<uint32_t> EcalDumpRaw::lastOrbit_
private

Definition at line 120 of file EcalDumpRaw.h.

Referenced by decode().

unsigned EcalDumpRaw::maxEventId_
private

Definition at line 147 of file EcalDumpRaw.h.

Referenced by analyze().

int EcalDumpRaw::maxEvt_
private

Definition at line 114 of file EcalDumpRaw.h.

const int EcalDumpRaw::maxTccsPerDcc_ = 4
staticprivate

Definition at line 124 of file EcalDumpRaw.h.

Referenced by decode().

const int EcalDumpRaw::maxTpgsPerTcc_ = 68
staticprivate

Definition at line 123 of file EcalDumpRaw.h.

Referenced by decode().

unsigned EcalDumpRaw::minEventId_
private

Definition at line 146 of file EcalDumpRaw.h.

Referenced by analyze().

const unsigned EcalDumpRaw::nDccs_ = 54
staticprivate

Definition at line 121 of file EcalDumpRaw.h.

Referenced by decode().

const int EcalDumpRaw::nRu_ = 70
staticprivate

Definition at line 170 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

const int EcalDumpRaw::nSamples = 10
staticprivate

Definition at line 108 of file EcalDumpRaw.h.

const int EcalDumpRaw::nTccTypes_ = 4
staticprivate

Definition at line 134 of file EcalDumpRaw.h.

std::vector<int> EcalDumpRaw::nTpgs_
private

Definition at line 161 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

int EcalDumpRaw::nTts_
private

Definition at line 167 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

unsigned EcalDumpRaw::orbit0_
private

Definition at line 148 of file EcalDumpRaw.h.

Referenced by decode().

bool EcalDumpRaw::orbit0Set_
private

Definition at line 150 of file EcalDumpRaw.h.

Referenced by decode().

uint32_t EcalDumpRaw::orbit_
private

Definition at line 149 of file EcalDumpRaw.h.

Referenced by decode().

std::vector<std::vector<uint32_t> > EcalDumpRaw::orbits_
private

Definition at line 159 of file EcalDumpRaw.h.

int EcalDumpRaw::profileFedId_
private

Definition at line 115 of file EcalDumpRaw.h.

int EcalDumpRaw::profileRuId_
private

Definition at line 116 of file EcalDumpRaw.h.

bool EcalDumpRaw::pulsePerLme_
private

Definition at line 182 of file EcalDumpRaw.h.

bool EcalDumpRaw::pulsePerLmod_
private

Definition at line 181 of file EcalDumpRaw.h.

bool EcalDumpRaw::pulsePerRu_
private

Definition at line 180 of file EcalDumpRaw.h.

unsigned EcalDumpRaw::side_
private

Definition at line 143 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::simpleTrigType_
private

Definition at line 155 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::srpBx_
private

Definition at line 172 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

int EcalDumpRaw::srpL1a_
private

Definition at line 164 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

int EcalDumpRaw::tccBlockLen64_
private

Definition at line 169 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::tccBx_
private

Definition at line 173 of file EcalDumpRaw.h.

Referenced by analyze().

int EcalDumpRaw::tccId_
private

Definition at line 183 of file EcalDumpRaw.h.

Referenced by decode().

int EcalDumpRaw::tccL1a_
private

Definition at line 165 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

int EcalDumpRaw::tccType_
private

type of TCC currently parsed

Definition at line 175 of file EcalDumpRaw.h.

Referenced by analyze(), and decode().

size_t EcalDumpRaw::towerBlockLength_
private

Definition at line 104 of file EcalDumpRaw.h.

Referenced by decode().

std::vector<std::vector<int> > EcalDumpRaw::tpg_
private

Definition at line 160 of file EcalDumpRaw.h.

Referenced by decode().

const int EcalDumpRaw::ttId_
staticprivate

TT ID in the order the TPG appears in the data

Definition at line 139 of file EcalDumpRaw.h.

Referenced by tpgTag(), and ttfTag().

int EcalDumpRaw::verbosity_
private

Definition at line 90 of file EcalDumpRaw.h.

Referenced by EcalDumpRaw().

bool EcalDumpRaw::writeDcc_
private

Definition at line 91 of file EcalDumpRaw.h.

Referenced by analyze(), and EcalDumpRaw().