BeSlinkAna.hh

#ifndef BeSlinkAna_HH
#define BeSlinkAna_HH

#include <iostream>
#include <iomanip>
#include <cassert>

#include "TH1D.h"
#include "TH2D.h"
#include "TH1I.h"

#include "Buffering/BeSlinkSim.hh"


class BeSlinkAna {
public:
  BeSlinkAna() {
  }

  void boj(const BeSlinkSim &b, unsigned nHistory=100000) {
    fBeSlinkSim=&b;
    fPreviousThrottleActive=false;
    fRaisedThrottleBx=-1;
    fEndedThrottleBx=-1;
    
    std::ostringstream sout;
    sout << "BeSlinkAna" << std::setw(2) << std::setfill('0')
	 << fBeSlinkSim->slinkNumber();
    
    std::string sHgcroc[18]={"00","01","02","03","04","05","06","07","08",
			     "09","10","11","12","13","14","15","16","17"};

    hTotalEventWords=new TH1D((sout.str()+"TotalEventWords").c_str(),
			      ";Total number of bytes per event",250,0,4000);
    
    const SlinkArch &slinkArch(*(fBeSlinkSim->fSlinkArch));
    for(unsigned e(0);e<slinkArch.numberOfEcons();e++) {
      const EconArch &econArch(slinkArch.econArch(e));
      for(unsigned b(0);b<econArch.numberOfBlocks();b++) {
	hEventBufferUsed[e][b]=new TH1D((sout.str()+"Econ"+sHgcroc[e]+"Block"+sHgcroc[b]+"EventBufferUsed").c_str(),
					";Number of event buffer words used",1000,0,10000);
	hEventBufferUsedHistory[e][b]=new TH1I((sout.str()+"Econ"+sHgcroc[e]+"Block"+sHgcroc[b]+"EventBufferUsedHistory").c_str(),
					       ";BX;Number of event buffer words used",nHistory,0,nHistory);
	hEventBufferExceeded[e][b]=new TH1D((sout.str()+"Econ"+sHgcroc[e]+"Block"+sHgcroc[b]+"EventBufferExceeded").c_str(),
					";Number of event buffer words used",1000,0,10000);
      }
    }

    hThrottleActiveHistory=new TH1I((sout.str()+"ThrottleActiveHistory").c_str(),
				    ";BX;Throttle active",nHistory,0,nHistory);
    hRaisedThrottleToL1Accept=new TH1D((sout.str()+"RaisedThrottleToL1Accept").c_str(),
				    ";L1Accept - ended throttle (BX)",1000,0,1000);
    hEndedThrottleToL1Accept=new TH1D((sout.str()+"EndedThrottleToL1Accept").c_str(),
				    ";L1Accept - ended throttle (BX)",1000,0,1000);
  }
  
  void analyseL1Accept(uint64_t bx) {
    if(fRaisedThrottleBx!=-1) hRaisedThrottleToL1Accept->Fill(bx-fRaisedThrottleBx);
    if(fEndedThrottleBx!=-1) hEndedThrottleToL1Accept->Fill(bx-fEndedThrottleBx);
  }
  
  void analyseBx(uint64_t bx) {
    const SlinkArch &slinkArch(*(fBeSlinkSim->fSlinkArch));
    for(unsigned e(0);e<slinkArch.numberOfEcons();e++) {
      const EconArch &econArch(slinkArch.econArch(e));
      for(unsigned b(0);b<econArch.numberOfBlocks();b++) {
	hEventBufferUsed[e][b]->Fill(fBeSlinkSim->fEventBuffer[e][b].used());
	hEventBufferUsedHistory[e][b]->Fill(bx,fBeSlinkSim->fEventBuffer[e][b].used());
	if(fBeSlinkSim->exceededBuffer()) hEventBufferExceeded[e][b]->Fill(fBeSlinkSim->fEventBuffer[e][b].used());
      }
    }

    if(fBeSlinkSim->fL1AcceptEoe) {
      hTotalEventWords->Fill(16*fBeSlinkSim->fTotalNumberOfEventWords);
    }

    if(fBeSlinkSim->fThrottleActive) {
      if(!fPreviousThrottleActive) fRaisedThrottleBx=bx;
      hThrottleActiveHistory->Fill(bx);
    } else {
      if(fPreviousThrottleActive) fEndedThrottleBx=bx;
    }

    fPreviousThrottleActive=fBeSlinkSim->fThrottleActive;
  }
  
  void eoj() {
  }
  
private:
  const BeSlinkSim *fBeSlinkSim;

  bool fPreviousThrottleActive;
  uint64_t fRaisedThrottleBx;
  uint64_t fEndedThrottleBx;
  
  TH1D *hTotalEventWords;
  TH1I *hThrottleActiveHistory;
  TH1D *hRaisedThrottleToL1Accept;
  TH1D *hEndedThrottleToL1Accept;
  
  TH1D *hEventBufferUsed[Dimensions::MaxNumberOfEconsPerSlink][Dimensions::MaxNumberOfBlocksPerEcon];
  TH1I *hEventBufferUsedHistory[Dimensions::MaxNumberOfEconsPerSlink][Dimensions::MaxNumberOfBlocksPerEcon];
  TH1D *hEventBufferExceeded[Dimensions::MaxNumberOfEconsPerSlink][Dimensions::MaxNumberOfBlocksPerEcon];
};

#endif