EBSelectiveReadoutTask.h

Go to the documentation of this file.

#ifndef EBSelectiveReadoutTask_H
#define EBSelectiveReadoutTask_H

/*
 * \file EBSelectiveReadoutTask.h
 *
 * \author P. Gras
 * \author E. Di Marco
 *
*/

#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"

class MonitorElement;
class DQMStore;
class EcalSRSettings;

class EBSelectiveReadoutTask: public edm::EDAnalyzer{

public:

EBSelectiveReadoutTask(const edm::ParameterSet& ps);

virtual ~EBSelectiveReadoutTask();

protected:

void analyze(const edm::Event& e, const edm::EventSetup& c);

void beginJob(void);

void endJob(void);

void beginRun(const edm::Run& r, const edm::EventSetup& c);

void endRun(const edm::Run& r, const edm::EventSetup& c);

void reset(void);

void setup(void);

void cleanup(void);

static const int nTTEta = 34;
static const int nTTPhi = 72;
static const int nECALDcc = 54;
static const int nEBDcc = 36;
static const int kByte = 1024;

static const int nDccChs = 68;

static const int nEbRus = 36*68;

static const int nEbEta = 170;

static const int nEbPhi = 360;

static const int ebTtEdge = 5;
 
static const int scEdge = 5;

static const int nTtEta = 56;

static const int nTtPhi = 72;

static const int bytesPerCrystal = 24;

int nCryTower[72][34];

int nEvtFullReadout[72][34];

int nEvtRUForced[72][34];

int nEvtAnyReadout[72][34];

int nEvtZS1Readout[72][34];

int nEvtZSReadout[72][34];

int nEvtCompleteReadoutIfZS[72][34];

int nEvtDroppedReadoutIfFR[72][34];

int nEvtHighInterest[72][34];

int nEvtMediumInterest[72][34];

int nEvtLowInterest[72][34];

int nEvtAnyInterest[72][34];

private:

enum subdet_t {EB, EE};

void anaDigi(const EBDataFrame& frame, const EBSrFlagCollection& srFlagColl, uint16_t statusCode);

void anaDigiInit();

unsigned dccNum(const DetId& xtalId) const;

int iEta2cIndex(int iEta) const{
  return (iEta<0)?iEta+85:iEta+84;
}

int iPhi2cIndex(int iPhi) const{
  return iPhi-1;
}

EcalTrigTowerDetId readOutUnitOf(const EBDetId& xtalId) const;
  
int iXY2cIndex(int iX) const{
  return iX-1;
}

int cIndex2iXY(int iX0) const{
  return iX0+1;
}

int cIndex2iEta(int i) const{
  return (i<85)?i-85:i-84;
}

int cIndex2iPhi(int i) const {
  return i+1;
}

double getEbEventSize(double nReadXtals) const;

double getDccOverhead(subdet_t subdet) const{
  //  return (subdet==EB?34:25)*8;
  return (subdet==EB?34:52)*8;
}

double getDccEventSize(int iDcc0, double nReadXtals) const{
  subdet_t subdet;
  if(iDcc0<9 || iDcc0>=45){
    subdet = EE;
  } else{
    subdet = EB;
  }
  return getDccOverhead(subdet)+nReadXtals*bytesPerCrystal+nRuPerDcc_[iDcc0]*8;
}

int dccPhiIndexOfRU(int i, int j) const;

int dccPhiIndex(int i, int j) const {
  return dccPhiIndexOfRU(i/5, j/5);
}

int dccIndex(int i, int j) const;

void configFirWeights(const std::vector<double>& weightsForZsFIR);

static int dccZsFIR(const EcalDataFrame& frame,
                    const std::vector<int>& firWeights,
                    int firstFIRSample,
                    bool* saturated = 0);
  

static std::vector<int> getFIRWeights(const std::vector<double>&
                                      normalizedWeights);

int getCrystalCount() { return 25; }


int nEb_;

int nEbLI_;

int nEbHI_;

int nPerDcc_[nECALDcc];

int nPerRu_[nECALDcc][nDccChs];   

int nRuPerDcc_[nECALDcc];

bool ebRuActive_[nEbEta/ebTtEdge][nEbPhi/ebTtEdge];

std::vector<int> firWeights_;

int firstFIRSample_;

int ievt_;

DQMStore* dqmStore_;

std::string prefixME_;

bool enableCleanup_;

bool mergeRuns_;

edm::EDGetTokenT<EBDigiCollection> EBDigiCollection_;
edm::EDGetTokenT<EcalRecHitCollection> EcalRecHitCollection_;
edm::EDGetTokenT<EBSrFlagCollection> EBSRFlagCollection_;
edm::EDGetTokenT<EcalTrigPrimDigiCollection> EcalTrigPrimDigiCollection_;
edm::EDGetTokenT<FEDRawDataCollection> FEDRawDataCollection_;

float xbins[37];
float ybins[89];

MonitorElement* EBDccEventSize_;
MonitorElement* EBDccEventSizeMap_;
MonitorElement* EBTowerSize_;
MonitorElement* EBTTFMismatch_;
MonitorElement* EBReadoutUnitForcedBitMap_;
MonitorElement* EBFullReadoutSRFlagMap_;
MonitorElement* EBFullReadoutSRFlagCount_;
MonitorElement* EBZeroSuppression1SRFlagMap_;
MonitorElement* EBHighInterestTriggerTowerFlagMap_;
MonitorElement* EBMediumInterestTriggerTowerFlagMap_;
MonitorElement* EBLowInterestTriggerTowerFlagMap_;
MonitorElement* EBTTFlags_;
MonitorElement* EBCompleteZSMap_;
MonitorElement* EBCompleteZSCount_;
MonitorElement* EBDroppedFRMap_;
MonitorElement* EBDroppedFRCount_;
MonitorElement* EBEventSize_;
MonitorElement* EBHighInterestPayload_;
MonitorElement* EBLowInterestPayload_;
MonitorElement* EBHighInterestZsFIR_;
MonitorElement* EBLowInterestZsFIR_;

bool init_;

bool useCondDb_;
const EcalSRSettings* settings_;

};

#endif