EcalSelectiveReadoutValidation Class Reference

#include <EcalSelectiveReadoutValidation.h>

Inheritance diagram for EcalSelectiveReadoutValidation:

Classes
struct	energiesEb_t
struct	energiesEe_t

Public Member Functions
void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
	EcalSelectiveReadoutValidation (const edm::ParameterSet &ps)
	Constructor. More...
	~EcalSelectiveReadoutValidation () override
	Destructor. More...
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &c) override
	Analyzes the event. More...
void	endRun (const edm::Run &r, const edm::EventSetup &c) override

Private Types
typedef EcalRecHit	RecHit
typedef EcalRecHitCollection	RecHitCollection
enum	subdet_t { EB, EE }
	distinguishes barral and endcap of ECAL. More...

Private Member Functions
template<class T , class U >
void	anaDigi (const T &frame, const U &srFlagColl)
void	anaDigiInit ()
void	analyzeDataVolume (const edm::Event &e, const edm::EventSetup &es)
void	analyzeEB (const edm::Event &event, const edm::EventSetup &es)
void	analyzeEE (const edm::Event &event, const edm::EventSetup &es)
void	analyzeTP (const edm::Event &event, const edm::EventSetup &es)
template<class T >
void	checkSrApplication (const edm::Event &event, T &srfs)
int	cIndex2iEta (int i) const
int	cIndex2iPhi (int i) const
int	cIndex2iTtEta (int i) const
int	cIndex2iTtPhi (int i) const
int	cIndex2iXY (int iX0) const
template<class T >
void	compareSrfColl (const edm::Event &event, T &srfFromData, T &computedSrf)
void	configFirWeights (const std::vector< double > &weightsForZsFIR)
std::pair< int, int >	dccCh (const DetId &xtalId) const
double	frame2Energy (const EcalDataFrame &frame) const
template<class T >
double	frame2EnergyForTp (const T &frame, int offset=0) const
double	getBytesPerCrystal () const
int	getCrystalCount (int iDcc, int iDccCh)
double	getDccEventSize (int iDcc0, double nReadXtals) const
double	getDccOverhead (subdet_t subdet) const
double	getDccSrDependentPayload (int iDcc0, double nReadRus, double nReadXtals) const
double	getEbEventSize (double nReadXtals) const
double	getEeEventSize (double nReadXtals) const
double	getL1aRate () const
int	getRuCount (int iDcc0) const
int	iEta2cIndex (int iEta) const
void	initAsciiFile ()
int	iPhi2cIndex (int iPhi) const
int	iTtEta2cIndex (int iEta) const
int	iTtPhi2cIndex (int iPhi) const
int	iXY2cIndex (int iX) const
void	printAvailableHists ()
void	readAllCollections (const edm::Event &e)
bool	registerHist (const std::string &name, const std::string &title)
void	selectFedsForLog ()
void	setTtEtSums (const edm::EventSetup &es, const EBDigiCollection &ebDigis, const EEDigiCollection &eeDigis)
void	updateL1aRate (const edm::Event &event)
EcalTrigTowerDetId	readOutUnitOf (const EBDetId &xtalId) const
EcalScDetId	readOutUnitOf (const EEDetId &xtalId) const
MonitorElement *	bookFloat (DQMStore::IBooker &, const std::string &name)
MonitorElement *	book1D (DQMStore::IBooker &, const std::string &name, const std::string &title, int nbins, double xmin, double xmax)
MonitorElement *	book2D (DQMStore::IBooker &, const std::string &name, const std::string &title, int nxbins, double xmin, double xmax, int nybins, double ymin, double ymax)
MonitorElement *	bookProfile (DQMStore::IBooker &, const std::string &name, const std::string &title, int nbins, double xmin, double xmax)
MonitorElement *	bookProfile2D (DQMStore::IBooker &, const std::string &name, const std::string &title, int nbinx, double xmin, double xmax, int nbiny, double ymin, double ymax, const char *option="")
void	fill (MonitorElement *me, float x)
void	fill (MonitorElement *me, float x, float yw)
void	fill (MonitorElement *me, float x, float y, float zw)
void	fill (MonitorElement *me, float x, float y, float z, float w)
int	ruGraphX (const EcalScDetId &id) const
int	ruGraphY (const EcalScDetId &id) const
int	ruGraphX (const EcalTrigTowerDetId &id) const
int	ruGraphY (const EcalTrigTowerDetId &id) const
int	xtalGraphX (const EEDetId &id) const
int	xtalGraphY (const EEDetId &id) const
int	xtalGraphX (const EBDetId &id) const
int	xtalGraphY (const EBDetId &id) const
int	dccId (const EcalScDetId &detId) const
int	dccId (const EcalTrigTowerDetId &detId) const

Static Private Member Functions
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)

Private Attributes
bool	allHists_
std::map< std::string, std::string >	availableHistList_
bool	collNotFoundWarn_
	Switch for collection-not-found warning. More...
energiesEb_t	ebEnergies [nEbEta][nEbPhi]
int	ebZsThr_
energiesEe_t	eeEnergies [nEndcaps][nEeX][nEeY]
int	eeZsThr_
const EcalElectronicsMapping *	elecMap_
int	firstFIRSample_
std::vector< int >	firWeights_
std::string	histDir_
std::set< std::string >	histList_
int	ievt_
bool	isRuComplete_ [nDccs_][nDccChs_]
bool	localReco_
std::vector< bool >	logErrForDccs_
bool	logSrApplicationErrors_
bool	logSrpAlgoErrors_
MonitorElement *	meL1aRate_
int	nCompleteZS_
	Counter of ZS-flagged RU fully read out. More...
int	nDroppedFRO_
	Counter of FRO-flagged RU dropped from data. More...
int	nEb_
int	nEbFROCnt_
	Counter of EB FRO-flagged RUs. More...
int	nEbHI_
int	nEbLI_
int	nEbZsErrors_
	Counter of EB ZS errors (LI channel below ZS threshold) More...
int	nEbZsErrorsType1_
int	nEe_
int	nEeFROCnt_
	Counter of EE FRO-flagged RUs. More...
int	nEeHI_
int	nEeLI_
int	nEeZsErrors_
	Counter of EE ZS errors (LI channel below ZS threshold) More...
int	nEeZsErrorsType1_
int	nHiPerDcc_ [nDccs_]
int	nHiRuPerDcc_ [nDccs_]
int	nIncompleteFRO_
	Counter of FRO-flagged RU only partial data. More...
int	nLiPerDcc_ [nDccs_]
int	nLiRuPerDcc_ [nDccs_]
int	nPerDcc_ [nDccs_]
int	nPerRu_ [nDccs_][nDccChs_]
int	nRuPerDcc_ [nDccs_]
std::string	outputFile_
	Output file for histograms. More...
bool	SkipInnerSC_
std::ofstream	srApplicationErrorLog_
	Output ascii file for unconsistency between Xtals and RU Flags. More...
std::string	srApplicationErrorLogFileName_
std::ofstream	srpAlgoErrorLog_
	Output ascii file for unconsistency on SR flags. More...
std::string	srpAlgoErrorLogFileName_
bool	tpInGeV_
const EcalTrigTowerConstituentsMap *	triggerTowerMap_
double	ttEtSums [nTtEta][nTtPhi]
bool	useEventRate_
bool	verbose_
	Verbosity switch. More...
std::vector< double >	weights_
bool	withEbSimHit_
bool	withEeSimHit_
std::ofstream	zsErrorLog_
	File to log ZS and other errors. More...
CollHandle< EBDigiCollection >	ebDigis_
CollHandle< EEDigiCollection >	eeDigis_
CollHandle< EBDigiCollection >	ebNoZsDigis_
CollHandle< EEDigiCollection >	eeNoZsDigis_
CollHandle< EBSrFlagCollection >	ebSrFlags_
CollHandle< EESrFlagCollection >	eeSrFlags_
CollHandle< EBSrFlagCollection >	ebComputedSrFlags_
CollHandle< EESrFlagCollection >	eeComputedSrFlags_
CollHandle< std::vector< PCaloHit > >	ebSimHits_
CollHandle< std::vector< PCaloHit > >	eeSimHits_
CollHandle< EcalTrigPrimDigiCollection >	tps_
CollHandle< RecHitCollection >	ebRecHits_
CollHandle< RecHitCollection >	eeRecHits_
CollHandle< FEDRawDataCollection >	fedRaw_
int64_t	tmax
int64_t	tmin
int64_t	l1aOfTmin
int64_t	l1aOfTmax
bool	l1aRateErr
MonitorElement *	meDccVol_
MonitorElement *	meDccLiVol_
MonitorElement *	meDccHiVol_
MonitorElement *	meDccVolFromData_
MonitorElement *	meVol_
MonitorElement *	meVolB_
MonitorElement *	meVolE_
MonitorElement *	meVolBLI_
MonitorElement *	meVolELI_
MonitorElement *	meVolLI_
MonitorElement *	meVolBHI_
MonitorElement *	meVolEHI_
MonitorElement *	meVolHI_
MonitorElement *	meChOcc_
MonitorElement *	meTp_
MonitorElement *	meTtf_
MonitorElement *	meTtfVsTp_
MonitorElement *	meTtfVsEtSum_
MonitorElement *	meTpVsEtSum_
MonitorElement *	meEbRecE_
MonitorElement *	meEbEMean_
MonitorElement *	meEbNoise_
MonitorElement *	meEbSimE_
MonitorElement *	meEbRecEHitXtal_
MonitorElement *	meEbRecVsSimE_
MonitorElement *	meEbNoZsRecVsSimE_
MonitorElement *	meEeRecE_
MonitorElement *	meEeEMean_
MonitorElement *	meEeNoise_
MonitorElement *	meEeSimE_
MonitorElement *	meEeRecEHitXtal_
MonitorElement *	meEeRecVsSimE_
MonitorElement *	meEeNoZsRecVsSimE_
MonitorElement *	meFullRoRu_
MonitorElement *	meZs1Ru_
MonitorElement *	meForcedRu_
MonitorElement *	meLiTtf_
MonitorElement *	meMiTtf_
MonitorElement *	meHiTtf_
MonitorElement *	meForcedTtf_
MonitorElement *	meTpMap_
MonitorElement *	meFullRoCnt_
MonitorElement *	meEbFullRoCnt_
MonitorElement *	meEeFullRoCnt_
MonitorElement *	meEbLiZsFir_
MonitorElement *	meEbHiZsFir_
MonitorElement *	meEbIncompleteRUZsFir_
MonitorElement *	meEeLiZsFir_
MonitorElement *	meEeHiZsFir_
MonitorElement *	meSRFlagsFromData_
MonitorElement *	meSRFlagsComputed_
MonitorElement *	meSRFlagsConsistency_
MonitorElement *	meIncompleteFRO_
MonitorElement *	meDroppedFRO_
MonitorElement *	meCompleteZS_
MonitorElement *	meIncompleteFROMap_
MonitorElement *	meDroppedFROMap_
MonitorElement *	meCompleteZSMap_
MonitorElement *	meIncompleteFRORateMap_
MonitorElement *	meDroppedFRORateMap_
MonitorElement *	meCompleteZSRateMap_
MonitorElement *	meIncompleteFROCnt_
MonitorElement *	meDroppedFROCnt_
MonitorElement *	meCompleteZSCnt_
MonitorElement *	meEbZsErrCnt_
MonitorElement *	meEeZsErrCnt_
MonitorElement *	meZsErrCnt_
MonitorElement *	meEbZsErrType1Cnt_
MonitorElement *	meEeZsErrType1Cnt_
MonitorElement *	meZsErrType1Cnt_
MonitorElement *	meEbFixedPayload_
MonitorElement *	meEeFixedPayload_
MonitorElement *	meFixedPayload_
bool	ebRuActive_ [nEbEta/ebTtEdge][nEbPhi/ebTtEdge]
bool	eeRuActive_ [nEndcaps][nEeX/scEdge][nEeY/scEdge]

Static Private Attributes
static const int	ebTtEdge = 5
	Number of crystals along an EB TT. More...
static const int	kByte_ = 1024
	number of bytes in 1 kByte: More...
static const int	maxDccId_ = minDccId_ + nDccs_ -1
static const int	minDccId_ = 1
static const unsigned	nDccChs_ = 68
	Number of input channels of a DCC. More...
static const int	nDccRus_ [nDccs_]
	number of RUs for each DCC More...
static const unsigned	nDccs_ = 54
	Total number of DCCs. More...
static const int	nEbDccs = 36
	number of DCCs for EB More...
static const int	nEbEta = 170
	number of crystals along Eta in EB More...
static const int	nEbPhi = 360
	number of crystals along Phi in EB More...
static const int	nEbRus = 36*68
	number of RUs for EB More...
static const int	nEbTtEta = 34
	Number of Trigger Towers in barrel along Eta. More...
static const int	nEeDccs = 18
	number of DCCs for EE More...
static const int	nEeRus = 2*(34+32+33+33+32+34+33+34+33)
	number of RUs for EE More...
static const int	nEeX = 100
	EE crystal grid size along X. More...
static const int	nEeY = 100
	EE crystal grid size along Y. More...
static const int	nEndcaps = 2
	number of endcaps More...
static const int	nMaxXtalPerRu = 25
	Number of crystals per Readout Unit excepted partial SCs. More...
static const int	nOneEeTtEta = 11
	Number of Trigger Towers in an endcap along Eta. More...
static const int	nTtEta = 2*nOneEeTtEta + nEbTtEta
	Number of Trigger Towers along Eta. More...
static const int	nTtPhi = 72
	Number of Trigger Towers along Phi. More...
static const double	rad2deg = 45./atan(1.)
	Conversion factor from radian to degree. More...
static const int	scEdge = 5
	Number of crystals along a supercrystal edge. More...

Detailed Description

Definition at line 34 of file EcalSelectiveReadoutValidation.h.

Member Typedef Documentation

typedef EcalRecHit EcalSelectiveReadoutValidation::RecHit

private

Definition at line 37 of file EcalSelectiveReadoutValidation.h.

typedef EcalRecHitCollection EcalSelectiveReadoutValidation::RecHitCollection

private

Definition at line 36 of file EcalSelectiveReadoutValidation.h.

Member Enumeration Documentation

enum EcalSelectiveReadoutValidation::subdet_t

private

distinguishes barral and endcap of ECAL.

Enumerator
EB
EE

Definition at line 58 of file EcalSelectiveReadoutValidation.h.

{EB, EE};

Constructor & Destructor Documentation

EcalSelectiveReadoutValidation::EcalSelectiveReadoutValidation

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 56 of file EcalSelectiveReadoutValidation.cc.

References allHists_, configFirWeights(), ebComputedSrFlags_, ebDigis_, ebNoZsDigis_, ebRecHits_, ebSimHits_, ebSrFlags_, ebZsThr_, eeComputedSrFlags_, eeDigis_, eeNoZsDigis_, eeRecHits_, eeSimHits_, eeSrFlags_, eeZsThr_, fedRaw_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), histList_, compare::hists, logSrApplicationErrors_, logSrpAlgoErrors_, outputFile_, CollHandle< T >::setToken(), srApplicationErrorLogFileName_, srpAlgoErrorLogFileName_, tps_, and verbose_.

                                                                                    :
  collNotFoundWarn_(ps.getUntrackedParameter<bool>("warnIfCollectionNotFound", true)),
  ebDigis_(ps.getParameter<edm::InputTag>("EbDigiCollection"), false,
       collNotFoundWarn_),
  eeDigis_(ps.getParameter<edm::InputTag>("EeDigiCollection"), false,
       collNotFoundWarn_),
  ebNoZsDigis_(ps.getParameter<edm::InputTag>("EbUnsuppressedDigiCollection"),
           false, false/*collNotFoundWarn_*/),
  eeNoZsDigis_(ps.getParameter<edm::InputTag>("EeUnsuppressedDigiCollection"),
           false, false/*collNotFoundWarn_*/),
  ebSrFlags_(ps.getParameter<edm::InputTag>("EbSrFlagCollection"), false,
         collNotFoundWarn_),
  eeSrFlags_(ps.getParameter<edm::InputTag>("EeSrFlagCollection"), false,
         collNotFoundWarn_),
  ebComputedSrFlags_(ps.getParameter<edm::InputTag>("EbSrFlagFromTTCollection"), false,
             false/*collNotFoundWarn_*/),
  eeComputedSrFlags_(ps.getParameter<edm::InputTag>("EeSrFlagFromTTCollection"), false,
             false/*collNotFoundWarn_*/),
  ebSimHits_(ps.getParameter<edm::InputTag>("EbSimHitCollection"), false,
         false/*collNotFoundWarn_*/),
  eeSimHits_(ps.getParameter<edm::InputTag>("EeSimHitCollection"), false,
         false/*collNotFoundWarn_*/),
  tps_(ps.getParameter<edm::InputTag>("TrigPrimCollection"), false,
       collNotFoundWarn_),
  ebRecHits_(ps.getParameter<edm::InputTag>("EbRecHitCollection"), false,
         false/*collNotFoundWarn_*/),
  eeRecHits_(ps.getParameter<edm::InputTag>("EeRecHitCollection"), false,
         false/*collNotFoundWarn_*/),
  fedRaw_(ps.getParameter<edm::InputTag>("FEDRawCollection"), false,
      false/*collNotFoundWarn_*/),
  tmax(0),
  tmin(numeric_limits<int64_t>::max()),
  l1aOfTmin(0),
  l1aOfTmax(0),
  triggerTowerMap_(nullptr),
  localReco_(ps.getParameter<bool>("LocalReco")),
  weights_(ps.getParameter<vector<double> >("weights")),
  tpInGeV_(ps.getParameter<bool>("tpInGeV")),
  firstFIRSample_(ps.getParameter<int>("ecalDccZs1stSample")),
  useEventRate_(ps.getParameter<bool>("useEventRate")),
  logErrForDccs_(nDccs_, false),
  ievt_(0),
  allHists_(false),
  histDir_(ps.getParameter<string>("histDir")),
  withEeSimHit_(false),
  withEbSimHit_(false){

  edm::ConsumesCollector collector(consumesCollector());
  ebDigis_.setToken(collector);
  eeDigis_.setToken(collector);
  ebNoZsDigis_.setToken(collector);
  eeNoZsDigis_.setToken(collector);
  ebSrFlags_.setToken(collector);
  eeSrFlags_.setToken(collector);
  ebComputedSrFlags_.setToken(collector);
  eeComputedSrFlags_.setToken(collector);
  ebSimHits_.setToken(collector);
  eeSimHits_.setToken(collector);
  tps_.setToken(collector);
  ebRecHits_.setToken(collector);
  eeRecHits_.setToken(collector);
  fedRaw_.setToken(collector);

  double ebZsThr = ps.getParameter<double>("ebZsThrADCCount");
  double eeZsThr = ps.getParameter<double>("eeZsThrADCCount");

  ebZsThr_ = lround(ebZsThr*4);
  eeZsThr_ = lround(eeZsThr*4);

  //File to log SRP algorithem inconsistency
  srpAlgoErrorLogFileName_
    = ps.getUntrackedParameter<string>("srpAlgoErrorLogFile","");
  logSrpAlgoErrors_ = (!srpAlgoErrorLogFileName_.empty());

  //File to log SRP decision application inconsistency
  srApplicationErrorLogFileName_
    = ps.getUntrackedParameter<string>("srApplicationErrorLogFile","");
  logSrApplicationErrors_ = (!srApplicationErrorLogFileName_.empty());

  //FIR ZS weights
  configFirWeights(ps.getParameter<vector<double> >("dccWeights"));

  // DQM ROOT output
  outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");

  if(!outputFile_.empty()){
    LogInfo("OutputInfo") << " Ecal Digi Task histograms will be saved to '"
              << outputFile_.c_str() << "'";
  } else{
    LogInfo("OutputInfo") << " Ecal Digi Task histograms will NOT be saved";
  }

  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);

  // get hold of back-end interface

  vector<string>
    hists(ps.getUntrackedParameter<vector<string> >("histograms",
                            vector<string>(1, "all")));

  for(vector<string>::iterator it = hists.begin();
      it!=hists.end(); ++it) histList_.insert(*it);
  if(histList_.find("all") != histList_.end()) allHists_ = true;

}

EcalSelectiveReadoutValidation::~EcalSelectiveReadoutValidation ( )

override

Destructor.

Definition at line 830 of file EcalSelectiveReadoutValidation.cc.

                                                               {
}

Member Function Documentation

template<class T , class U >

void EcalSelectiveReadoutValidation::anaDigi ( const T &  frame, const U &  srFlagColl  )

private

Accumulates statitics for data volume analysis. To be called for each ECAL digi. See anaDigiInit().

Definition at line 1498 of file EcalSelectiveReadoutValidation.cc.

References Reference_intrackfit_cff::barrel, dccCh(), ebRuActive_, ebTtEdge, EcalBarrel, eeRuActive_, Exception, RemoveAddSevLevel::flag, iEta2cIndex(), iPhi2cIndex(), iXY2cIndex(), minDccId_, nEb_, nEbHI_, nEbLI_, nEe_, nEeHI_, nEeLI_, nHiPerDcc_, nHiRuPerDcc_, nLiPerDcc_, nLiRuPerDcc_, nPerDcc_, nPerRu_, nRuPerDcc_, DetId::rawId(), readOutUnitOf(), scEdge, EcalSrFlag::SRF_FULL, DetId::subdetId(), and ecaldqm::zside().

Referenced by analyzeDataVolume().

                                             {
  const DetId& xtalId = frame.id();
  typedef typename U::key_type RuDetId;
  const RuDetId& ruId = readOutUnitOf(frame.id());
  typename U::const_iterator srf = srFlagColl.find(ruId);

  bool highInterest = false;
  int flag = 0;
  
  if(srf != srFlagColl.end()){
    flag = srf->value() & ~EcalSrFlag::SRF_FORCED_MASK;
    
    highInterest = (flag == EcalSrFlag::SRF_FULL);
    
  }

  bool barrel = (xtalId.subdetId()==EcalBarrel);

  pair<int,int> ch = dccCh(xtalId);
  
  if(barrel){
    ++nEb_;
    if(highInterest){
      ++nEbHI_;
    } else{//low interest
      ++nEbLI_;
    }
    int iEta0 = iEta2cIndex(static_cast<const EBDetId&>(xtalId).ieta());
    int iPhi0 = iPhi2cIndex(static_cast<const EBDetId&>(xtalId).iphi());
    if(!ebRuActive_[iEta0/ebTtEdge][iPhi0/ebTtEdge]){
      ++nRuPerDcc_[ch.first-minDccId_];
      if(highInterest){
        ++nHiRuPerDcc_[ch.first-minDccId_];
      } else{
        ++nLiRuPerDcc_[ch.first-minDccId_];
      }

      ebRuActive_[iEta0/ebTtEdge][iPhi0/ebTtEdge] = true;
    }
  } else{//endcap
    ++nEe_;
    if(highInterest){
      ++nEeHI_;
    } else{//low interest
      ++nEeLI_;
    }
    int iX0 = iXY2cIndex(static_cast<const EEDetId&>(frame.id()).ix());
    int iY0 = iXY2cIndex(static_cast<const EEDetId&>(frame.id()).iy());
    int iZ0 = static_cast<const EEDetId&>(frame.id()).zside()>0?1:0;

    if(!eeRuActive_[iZ0][iX0/scEdge][iY0/scEdge]){
      ++nRuPerDcc_[ch.first-minDccId_];
      if(highInterest){
        ++nHiRuPerDcc_[ch.first-minDccId_];
      } else{
        ++nLiRuPerDcc_[ch.first-minDccId_];
      }

      eeRuActive_[iZ0][iX0/scEdge][iY0/scEdge] = true;
    }
  }

  if(ch.second < 1 || ch.second > 68){
    throw cms::Exception("EcalSelectiveReadoutValidation")
      << "Error in DCC channel retrieval for crystal with detId "
      << xtalId.rawId() << "DCC channel out of allowed range [1..68]\n";
  }
  ++nPerDcc_[ch.first-minDccId_];
  ++nPerRu_[ch.first-minDccId_][ch.second-1];
  if(highInterest){
    ++nHiPerDcc_[ch.first-minDccId_];
  } else{//low interest channel
    ++nLiPerDcc_[ch.first-minDccId_];
  }
}

void EcalSelectiveReadoutValidation::anaDigiInit ( )

private

Initializes statistics accumalator for data volume analysis. To be call at start of each event analysis.

Definition at line 1575 of file EcalSelectiveReadoutValidation.cc.

References ebRuActive_, eeRuActive_, nEb_, nEbHI_, nEbLI_, nEe_, nEeHI_, nEeLI_, nHiPerDcc_, nHiRuPerDcc_, nLiPerDcc_, nLiRuPerDcc_, nPerDcc_, nPerRu_, and nRuPerDcc_.

Referenced by analyzeDataVolume().

                                                {
  nEb_ = 0;
  nEe_ = 0;
  nEeLI_ = 0;
  nEeHI_ = 0;
  nEbLI_ = 0;
  nEbHI_ = 0;
  bzero(nPerDcc_, sizeof(nPerDcc_));
  bzero(nLiPerDcc_, sizeof(nLiPerDcc_));
  bzero(nHiPerDcc_, sizeof(nHiPerDcc_));
  bzero(nRuPerDcc_, sizeof(nRuPerDcc_));
  bzero(ebRuActive_, sizeof(ebRuActive_));
  bzero(eeRuActive_, sizeof(eeRuActive_));
  bzero(nPerRu_, sizeof(nPerRu_));
  bzero(nLiRuPerDcc_, sizeof(nLiRuPerDcc_));
  bzero(nHiRuPerDcc_, sizeof(nHiRuPerDcc_));
}

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

overrideprotected

Analyzes the event.

Definition at line 189 of file EcalSelectiveReadoutValidation.cc.

References analyzeDataVolume(), analyzeEB(), analyzeEE(), analyzeTP(), checkSrApplication(), compareSrfColl(), ebComputedSrFlags_, ebNoZsDigis_, ebSimHits_, ebSrFlags_, eeComputedSrFlags_, eeNoZsDigis_, eeSimHits_, eeSrFlags_, fill(), ievt_, meCompleteZSCnt_, meDroppedFROCnt_, meEbFullRoCnt_, meEbZsErrCnt_, meEbZsErrType1Cnt_, meEeFullRoCnt_, meEeZsErrCnt_, meEeZsErrType1Cnt_, meFullRoCnt_, meIncompleteFROCnt_, meZsErrCnt_, meZsErrType1Cnt_, nCompleteZS_, nDroppedFRO_, nEbFROCnt_, nEbZsErrors_, nEbZsErrorsType1_, nEeFROCnt_, nEeZsErrors_, nEeZsErrorsType1_, nIncompleteFRO_, readAllCollections(), selectFedsForLog(), setTtEtSums(), tps_, updateL1aRate(), withEbSimHit_, and withEeSimHit_.

                                               {

  updateL1aRate(event);
  
  //retrieves event products:
  readAllCollections(event);

  withEeSimHit_ = (!eeSimHits_->empty());
  withEbSimHit_ = (!ebSimHits_->empty());

  if(ievt_<10){
    edm::LogInfo("EcalSrValid") << "Size of TP collection: " << tps_->size() << std::endl
                << "Size of EB SRF collection read from data: "
                << ebSrFlags_->size() << std::endl
                << "Size of EB SRF collection computed from data TTFs: "
                << ebComputedSrFlags_->size() << std::endl
                << "Size of EE SRF collection read from data: "
                << eeSrFlags_->size() << std::endl
                << "Size of EE SRF collection computed from data TTFs: "
                << eeComputedSrFlags_->size() << std::endl;
  }
  
  if(ievt_==0){
    selectFedsForLog(); //note: must be called after readAllCollection
  }

  //computes Et sum trigger tower crystals:
  setTtEtSums(es, *ebNoZsDigis_, *eeNoZsDigis_);

  //Data Volume
  analyzeDataVolume(event, es);

  //EB digis
  //must be called after analyzeDataVolume because it uses
  //isRuComplete_ array that this method fills
  analyzeEB(event, es);

  //EE digis
  //must be called after analyzeDataVolume because it uses
  //isRuComplete_ array that this method fills
  analyzeEE(event, es);

  fill(meFullRoCnt_, nEeFROCnt_+nEbFROCnt_);
  fill(meEbFullRoCnt_, nEbFROCnt_);
  fill(meEeFullRoCnt_, nEeFROCnt_);

  fill(meEbZsErrCnt_, nEbZsErrors_);
  fill(meEeZsErrCnt_, nEeZsErrors_);
  fill(meZsErrCnt_, nEbZsErrors_ + nEeZsErrors_);

  fill(meEbZsErrType1Cnt_, nEbZsErrorsType1_);
  fill(meEeZsErrType1Cnt_, nEeZsErrorsType1_);
  fill(meZsErrType1Cnt_, nEbZsErrorsType1_ + nEeZsErrorsType1_);

  //TP
  analyzeTP(event, es);

  if(!ebComputedSrFlags_->empty()){
    compareSrfColl(event, *ebSrFlags_, *ebComputedSrFlags_);
  }
  if(!eeComputedSrFlags_->empty()){
    compareSrfColl(event, *eeSrFlags_, *eeComputedSrFlags_);
  }
  nDroppedFRO_ = 0;
  nIncompleteFRO_ = 0;
  nCompleteZS_ = 0;
  checkSrApplication(event, *ebSrFlags_);
  checkSrApplication(event, *eeSrFlags_);
  fill(meDroppedFROCnt_, nDroppedFRO_);
  fill(meIncompleteFROCnt_, nIncompleteFRO_);
  fill(meCompleteZSCnt_, nCompleteZS_);
  ++ievt_;
}

void EcalSelectiveReadoutValidation::analyzeDataVolume ( const edm::Event &  e, const edm::EventSetup &  es  )

private

Data volume analysis. To be called for each event.

Parameters

event	EDM event
es	event setup

Definition at line 1432 of file EcalSelectiveReadoutValidation.cc.

References a, anaDigi(), anaDigiInit(), b, ebDigis_, ebSrFlags_, eeDigis_, eeSrFlags_, FEDRawDataCollection::FEDData(), fedRaw_, fill(), getBytesPerCrystal(), getCrystalCount(), getDccEventSize(), getDccSrDependentPayload(), getEbEventSize(), getEeEventSize(), isRuComplete_, kByte_, maxDccId_, meDccHiVol_, meDccLiVol_, meDccVol_, meDccVolFromData_, meVol_, meVolB_, meVolBHI_, meVolBLI_, meVolE_, meVolEHI_, meVolELI_, meVolHI_, meVolLI_, minDccId_, nDccRus_, nDccs_, nEb_, nEbHI_, nEbLI_, nEe_, nEeHI_, nEeLI_, nHiPerDcc_, nHiRuPerDcc_, nLiPerDcc_, nLiRuPerDcc_, nPerDcc_, nPerRu_, FEDRawData::size(), and edm::DataFrameContainer::size().

Referenced by analyze().

                                                    {

  anaDigiInit();


  //Complete RU, i.e. RU actually fully readout
  for(int iDcc = minDccId_; iDcc <= maxDccId_; ++iDcc){
    for(int iCh = 1; iCh < nDccRus_[iDcc-minDccId_]; ++iCh){
      isRuComplete_[iDcc-minDccId_][iCh-1]
        = (nPerRu_[iDcc-minDccId_][iCh-1]==getCrystalCount(iDcc, iCh));
    }
  }


  //Barrel
  for (unsigned int digis=0; digis<ebDigis_->size(); ++digis){
    EBDataFrame ebdf = (*ebDigis_)[digis];
    anaDigi(ebdf, *ebSrFlags_);
  }

  // Endcap
  for (unsigned int digis=0; digis<eeDigis_->size(); ++digis){
    EEDataFrame eedf = (*eeDigis_)[digis];
    anaDigi(eedf, *eeSrFlags_);
  }

  //histos
  for(unsigned iDcc0 = 0; iDcc0 <  nDccs_; ++iDcc0){
    fill(meDccVol_, iDcc0+1, getDccEventSize(iDcc0, nPerDcc_[iDcc0])/kByte_);
    fill(meDccLiVol_, iDcc0+1,
         getDccSrDependentPayload(iDcc0, nLiRuPerDcc_[iDcc0],
                                  nLiPerDcc_[iDcc0])/kByte_);
    fill(meDccHiVol_, iDcc0+1,
         getDccSrDependentPayload(iDcc0, nHiRuPerDcc_[iDcc0],
                                  nHiPerDcc_[iDcc0])/kByte_);
    const FEDRawDataCollection& raw = *fedRaw_;
    fill(meDccVolFromData_, iDcc0+1,
     ((double)raw.FEDData(601+iDcc0).size())/kByte_);
  }


  //low interesest channels:
  double a = nEbLI_*getBytesPerCrystal()/kByte_; //getEbEventSize(nEbLI_)/kByte_;
  fill(meVolBLI_, a);
  double b = nEeLI_*getBytesPerCrystal()/kByte_; //getEeEventSize(nEeLI_)/kByte_;
  fill(meVolELI_, b);
  fill(meVolLI_, a+b);

  //high interest chanels:
  a = nEbHI_*getBytesPerCrystal()/kByte_; //getEbEventSize(nEbHI_)/kByte_;
  fill(meVolBHI_, a);
  b = nEeHI_*getBytesPerCrystal()/kByte_; //getEeEventSize(nEeHI_)/kByte_;
  fill(meVolEHI_, b);
  fill(meVolHI_, a+b);

  //any-interest channels:
  a = getEbEventSize(nEb_)/kByte_;
  fill(meVolB_, a);
  b = getEeEventSize(nEe_)/kByte_;
  fill(meVolE_, b);
  fill(meVol_, a+b);
}

void EcalSelectiveReadoutValidation::analyzeEB ( const edm::Event &  event, const edm::EventSetup &  es  )

private

ECAL barrel data analysis. To be called for each event.

Parameters

event	EDM event
es	event setup

Definition at line 543 of file EcalSelectiveReadoutValidation.cc.

References edm::DataFrameContainer::begin(), gather_cfg::cout, dccCh(), dccZsFIR(), ebDigis_, ebEnergies, ebNoZsDigis_, ebRecHits_, ebSimHits_, ebSrFlags_, ebZsThr_, DetId::Ecal, EcalBarrel, edm::DataFrameContainer::end(), PCaloHit::energy(), EcalRecHit::energy(), PV3DBase< T, PVType, FrameType >::eta(), EcalSelectiveReadoutValidation::energiesEb_t::eta, Exception, fill(), plotBeamSpotDB::first, firstFIRSample_, firWeights_, RemoveAddSevLevel::flag, frame2Energy(), EcalSelectiveReadoutValidation::energiesEb_t::gain12, ecalMGPA::gainId(), edm::EventSetup::get(), CaloSubdetectorGeometry::getGeometry(), mps_fire::i, EBDataFrame::id(), EBSrFlag::id(), PCaloHit::id(), EcalRecHit::id(), EcalTrigTowerDetId::ieta(), EBDetId::ieta(), iEta2cIndex(), EcalTrigTowerDetId::ietaAbs(), ievt_, EcalTrigTowerDetId::iphi(), iPhi2cIndex(), isRuComplete_, localReco_, SiStripPI::max, meChOcc_, meEbEMean_, meEbHiZsFir_, meEbLiZsFir_, meEbNoise_, meEbNoZsRecVsSimE_, meEbRecE_, meEbRecEHitXtal_, meEbRecVsSimE_, meEbSimE_, meForcedRu_, meFullRoRu_, meZs1Ru_, nEbEta, nEbFROCnt_, nEbPhi, nEbZsErrors_, nEbZsErrorsType1_, EcalSelectiveReadoutValidation::energiesEb_t::noZsRecE, PV3DBase< T, PVType, FrameType >::phi(), EcalSelectiveReadoutValidation::energiesEb_t::phi, rad2deg, readOutUnitOf(), EcalSelectiveReadoutValidation::energiesEb_t::recE, alignCSCRings::s, edm::second(), EcalSelectiveReadoutValidation::energiesEb_t::simE, rpcPointValidation_cfi::simHit, EcalSelectiveReadoutValidation::energiesEb_t::simHit, EcalDataFrame::size(), srApplicationErrorLog_, EcalSrFlag::SRF_FORCED_MASK, EcalSrFlag::SRF_FULL, EcalSrFlag::SRF_ZS1, findQualityFiles::v, EcalSrFlag::value(), withEbSimHit_, x, xtalGraphX(), xtalGraphY(), y, and EcalTrigTowerDetId::zside().

Referenced by analyze().

                                              {

  bool eventError = false;
  nEbZsErrors_ = 0;
  nEbZsErrorsType1_ = 0;
  vector<pair<int,int> > xtalEtaPhi;

  xtalEtaPhi.reserve(nEbPhi*nEbEta);
  for(int iEta0=0; iEta0<nEbEta; ++iEta0){
    for(int iPhi0=0; iPhi0<nEbPhi; ++iPhi0){
      ebEnergies[iEta0][iPhi0].noZsRecE = -numeric_limits<double>::max();
      ebEnergies[iEta0][iPhi0].recE = -numeric_limits<double>::max();
      ebEnergies[iEta0][iPhi0].simE = 0; //must be zero.
      ebEnergies[iEta0][iPhi0].simHit = 0;
      ebEnergies[iEta0][iPhi0].gain12 = false;
      xtalEtaPhi.push_back(pair<int,int>(iEta0, iPhi0));
    }
  }

    // get the barrel geometry:
  edm::ESHandle<CaloGeometry> geoHandle;

  es.get<CaloGeometryRecord>().get(geoHandle);
  const CaloSubdetectorGeometry *geometry_p
    = (*geoHandle).getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
//CaloSubdetectorGeometry const& geometry = *geometry_p;

  //EB unsuppressed digis:
  for(EBDigiCollection::const_iterator it = ebNoZsDigis_->begin();
      it != ebNoZsDigis_->end(); ++it){
    const EBDataFrame& frame = *it;
    int iEta0 = iEta2cIndex(static_cast<const EBDetId&>(frame.id()).ieta());
    int iPhi0 = iPhi2cIndex(static_cast<const EBDetId&>(frame.id()).iphi());
    if(iEta0<0 || iEta0>=nEbEta){
      stringstream s;
      s << "EcalSelectiveReadoutValidation: "
        << "iEta0 (= " << iEta0 << ") is out of range ("
        << "[0," << nEbEta -1 << "]\n";
      throw cms::Exception(s.str());
    }
    if(iPhi0<0 || iPhi0>=nEbPhi){
      stringstream s;
      s << "EcalSelectiveReadoutValidation: "
        << "iPhi0 (= " << iPhi0 << ") is out of range ("
        << "[0," << nEbPhi -1 << "]\n";
      throw cms::Exception(s.str());
    }

    ebEnergies[iEta0][iPhi0].noZsRecE = frame2Energy(frame);
    ebEnergies[iEta0][iPhi0].gain12 = true;
    for(int i = 0; i< frame.size(); ++i){
      const int gain12Code = 0x1;
      if(frame[i].gainId()!=gain12Code) ebEnergies[iEta0][iPhi0].gain12 =  false;
    }

    const GlobalPoint xtalPos
      = geometry_p->getGeometry(frame.id())->getPosition();

    ebEnergies[iEta0][iPhi0].phi = rad2deg*((double)xtalPos.phi());
    ebEnergies[iEta0][iPhi0].eta = xtalPos.eta();
  } //next non-zs digi

    //EB sim hits
  for(vector<PCaloHit>::const_iterator it = ebSimHits_->begin();
      it != ebSimHits_->end(); ++it){
    const PCaloHit& simHit = *it;
    EBDetId detId(simHit.id());
    int iEta = detId.ieta();
    int iEta0 =iEta2cIndex(iEta);
    int iPhi = detId.iphi();
    int iPhi0 = iPhi2cIndex(iPhi);
    ebEnergies[iEta0][iPhi0].simE += simHit.energy();
    ++ebEnergies[iEta0][iPhi0].simHit;
  }

  bool crystalShot[nEbEta][nEbPhi];
  pair<int,int> EBxtalCoor[nEbEta][nEbPhi]; 

  for(int iEta0=0; iEta0<nEbEta; ++iEta0){
    for(int iPhi0=0; iPhi0<nEbPhi; ++iPhi0){
      crystalShot[iEta0][iPhi0] = false;
      EBxtalCoor[iEta0][iPhi0] = make_pair(0,0);
    }
  }

  int nEbDigi = 0;

  for(EBDigiCollection::const_iterator it = ebDigis_->begin();
      it != ebDigis_->end(); ++it){
    ++nEbDigi;
    const EBDataFrame& frame = *it;
    int iEta = static_cast<const EBDetId&>(frame.id()).ieta();
    int iPhi = static_cast<const EBDetId&>(frame.id()).iphi();
    int iEta0 = iEta2cIndex(iEta);
    int iPhi0 = iPhi2cIndex(iPhi);
    if(iEta0<0 || iEta0>=nEbEta){
      throw (cms::Exception("EcalSelectiveReadoutValidation")
             << "iEta0 (= " << iEta0 << ") is out of range ("
             << "[0," << nEbEta -1 << "]");
    }
    if(iPhi0<0 || iPhi0>=nEbPhi){
      throw (cms::Exception("EcalSelectiveReadoutValidation")
             << "iPhi0 (= " << iPhi0 << ") is out of range ("
             << "[0," << nEbPhi -1 << "]");
    }
    assert(iEta0>=0 && iEta0<nEbEta);
    assert(iPhi0>=0 && iPhi0<nEbPhi);
    if(!crystalShot[iEta0][iPhi0]){
      crystalShot[iEta0][iPhi0] = true;
      EBxtalCoor[iEta0][iPhi0] = make_pair(xtalGraphX(frame.id()), xtalGraphY(frame.id()));
    } else{
      cout << "Error: several digi for same crystal!";
      abort();
    }
    if(localReco_){
      ebEnergies[iEta0][iPhi0].recE = frame2Energy(frame);
    }

    ebEnergies[iEta0][iPhi0].gain12 = true;
    for(int i = 0; i< frame.size(); ++i){
      const int gain12Code = 0x1;
      if(frame[i].gainId()!=gain12Code){
        ebEnergies[iEta0][iPhi0].gain12 =  false;
      }
    }

    EBSrFlagCollection::const_iterator srf
      = ebSrFlags_->find(readOutUnitOf(frame.id()));

    bool highInterest = false;

    // if(srf == ebSrFlags_->end()){
    //       throw cms::Exception("EcalSelectiveReadoutValidation")
    //  << __FILE__ << ":" << __LINE__ << ": SR flag not found";
    //}

    if(srf != ebSrFlags_->end()){
      highInterest = ((srf->value() & ~EcalSrFlag::SRF_FORCED_MASK)
                      == EcalSrFlag::SRF_FULL);
    }


    if(highInterest){
      fill(meEbHiZsFir_, dccZsFIR(frame, firWeights_, firstFIRSample_, nullptr));
    } else{
      int v = dccZsFIR(frame, firWeights_, firstFIRSample_, nullptr);
      fill(meEbLiZsFir_, v);
      if(v < ebZsThr_){
        eventError = true;
        ++nEbZsErrors_;
        pair<int,int> ru = dccCh(frame.id());
        if(isRuComplete_[ru.first][ru.second-1]) ++nEbZsErrorsType1_;
        if(nEbZsErrors_ < 3){
          srApplicationErrorLog_ << event.id() << ", "
                                 << "RU " << frame.id() << ", "
                                 << "DCC " << ru.first
                                 << " Ch : " << ru.second << ": "
                                 << "LI channel under ZS threshold.\n";
        }
        if(nEbZsErrors_==3){
          srApplicationErrorLog_ << event.id() << ": "
                                 << "more ZS errors for this event...\n";
        }
      }
    }
  } //next EB digi

  for(int iEta0=0; iEta0<nEbEta; ++iEta0){
    for(int iPhi0=0; iPhi0<nEbPhi; ++iPhi0) 
      fill(meChOcc_, EBxtalCoor[iEta0][iPhi0].first, EBxtalCoor[iEta0][iPhi0].second,crystalShot[iEta0][iPhi0] ? 1. : 0.);
  }

  if(!localReco_){
    for(RecHitCollection::const_iterator it
          = ebRecHits_->begin();
        it != ebRecHits_->end(); ++it){
      ++nEbDigi;
      const RecHit& hit = *it;
      int iEta = static_cast<const EBDetId&>(hit.id()).ieta();
      int iPhi = static_cast<const EBDetId&>(hit.id()).iphi();
      int iEta0 = iEta2cIndex(iEta);
      int iPhi0 = iPhi2cIndex(iPhi);
      if(iEta0<0 || iEta0>=nEbEta){
        LogError("EcalSrValid") << "iEta0 (= " << iEta0 << ") is out of range ("
                                << "[0," << nEbEta -1 << "]\n";
      }
      if(iPhi0<0 || iPhi0>=nEbPhi){
        LogError("EcalSrValid") << "iPhi0 (= " << iPhi0 << ") is out of range ("
                                << "[0," << nEbPhi -1 << "]\n";
      }
      ebEnergies[iEta0][iPhi0].recE = hit.energy();
    }
  }

  for(unsigned int i=0; i<xtalEtaPhi.size(); ++i){
    int iEta0 = xtalEtaPhi[i].first;
    int iPhi0=  xtalEtaPhi[i].second;
    energiesEb_t& energies = ebEnergies[iEta0][iPhi0];

    double recE = energies.recE;
    if(recE!=-numeric_limits<double>::max()){//not zero suppressed
      fill(meEbRecE_, ebEnergies[iEta0][iPhi0].recE);
      fill(meEbEMean_, ievt_+1, recE);
    } //not zero suppressed

    if(withEbSimHit_){
      if(!energies.simHit){//noise only crystal channel
        fill(meEbNoise_, energies.noZsRecE);
      } else{
        fill(meEbSimE_, energies.simE);
        fill(meEbRecEHitXtal_, energies.recE);
      }
      fill(meEbRecVsSimE_, energies.simE, energies.recE);
      fill(meEbNoZsRecVsSimE_, energies.simE, energies.noZsRecE);
    }
  }

  int EBZs1RuCount[2][17][72];
  int EBFullRuCount[2][17][72];
  int EBForcedRuCount[2][17][72];
  std::pair<int, int> EBtowerCoor[2][17][72];
  for(int iZ(0); iZ < 2; iZ++){
    for(int iEta(0);iEta < 17; iEta++){
      for(int iPhi(0);iPhi < 72; iPhi++){
        EBZs1RuCount[iZ][iEta][iPhi] = 0;
        EBFullRuCount[iZ][iEta][iPhi] = 0;
        EBForcedRuCount[iZ][iEta][iPhi] = 0;
      }
    }
  }

  //SRF
  nEbFROCnt_ = 0;
  char ebSrfMark[2][17][72];
  bzero(ebSrfMark, sizeof(ebSrfMark));
  //      int idbg = 0;
  for(EBSrFlagCollection::const_iterator it = ebSrFlags_->begin();
      it != ebSrFlags_->end(); ++it){
    const EBSrFlag& srf = *it;
    int iEtaAbs = srf.id().ietaAbs();
    int iPhi = srf.id().iphi();
    int iZ = srf.id().zside();

    //  cout << "--> " << ++idbg << iEtaAbs << " " << iPhi << " "  << iZ
    //       << " " << srf.id() << "\n";
    
    if(iEtaAbs < 1 || iEtaAbs > 17
       || iPhi < 1 || iPhi > 72) throw cms::Exception("EcalSelectiveReadoutValidation")
                                   << "Found a barrel SRF with an invalid det ID: " << srf.id() << ".\n";
    ++ebSrfMark[iZ>0?1:0][iEtaAbs-1][iPhi-1];
    if(ebSrfMark[iZ>0?1:0][iEtaAbs-1][iPhi-1] > 1) throw cms::Exception("EcalSelectiveReadoutValidation")
                                                     << "Duplicate SRF for RU " << srf.id() << ".\n";

    EBtowerCoor[iZ > 0 ? 1 : 0][iEtaAbs - 1][iPhi - 1] = std::pair<int, int>(srf.id().ieta(), srf.id().iphi());

    int flag = srf.value() & ~EcalSrFlag::SRF_FORCED_MASK;
    if(flag == EcalSrFlag::SRF_ZS1){
      EBZs1RuCount[iZ>0?1:0][iEtaAbs-1][iPhi-1] +=1;
    }
    if(flag == EcalSrFlag::SRF_FULL){
      EBFullRuCount[iZ>0?1:0][iEtaAbs-1][iPhi-1] +=1;
      ++nEbFROCnt_;
    }
    if(srf.value() & EcalSrFlag::SRF_FORCED_MASK){
      EBForcedRuCount[iZ>0?1:0][iEtaAbs-1][iPhi-1] +=1;
    }
  }
  for(int iZ(0); iZ < 2; iZ++){
    for(int iEta(0);iEta < 17; iEta++){
      for(int iPhi(0);iPhi < 72; iPhi++){
        float x(EBtowerCoor[iZ][iEta][iPhi].first);
        float y(EBtowerCoor[iZ][iEta][iPhi].second);
        fill(meZs1Ru_, x, y, EBZs1RuCount[iZ][iEta][iPhi]);
        fill(meFullRoRu_, x, y, EBFullRuCount[iZ][iEta][iPhi]);
        fill(meForcedRu_, x, y, EBForcedRuCount[iZ][iEta][iPhi]);
      }
    }
  }

  if(eventError) srApplicationErrorLog_ << event.id()
                                        << ": " << nEbZsErrors_
                                        << " ZS-flagged EB channels under "
                   "the ZS threshold, whose " << nEbZsErrorsType1_
                                        << " in a complete RU.\n";
}

void EcalSelectiveReadoutValidation::analyzeEE ( const edm::Event &  event, const edm::EventSetup &  es  )

private

ECAL endcap data analysis. To be called for each event.

Parameters

event	EDM event
es	event setup

Definition at line 265 of file EcalSelectiveReadoutValidation.cc.

References edm::DataFrameContainer::begin(), gather_cfg::cout, dccCh(), dccZsFIR(), DetId::Ecal, EcalEndcap, eeDigis_, eeEnergies, eeNoZsDigis_, eeRecHits_, eeSimHits_, eeSrFlags_, eeZsThr_, edm::DataFrameContainer::end(), PCaloHit::energy(), EcalRecHit::energy(), PV3DBase< T, PVType, FrameType >::eta(), EcalSelectiveReadoutValidation::energiesEe_t::eta, Exception, fill(), plotBeamSpotDB::first, firstFIRSample_, firWeights_, RemoveAddSevLevel::flag, frame2Energy(), EcalSelectiveReadoutValidation::energiesEe_t::gain12, ecalMGPA::gainId(), edm::EventSetup::get(), CaloSubdetectorGeometry::getGeometry(), mps_fire::i, EEDataFrame::id(), EESrFlag::id(), PCaloHit::id(), EcalRecHit::id(), ievt_, isRuComplete_, EcalScDetId::ix(), EEDetId::ix(), iXY2cIndex(), EcalScDetId::iy(), EEDetId::iy(), localReco_, SiStripPI::max, meChOcc_, meEeEMean_, meEeHiZsFir_, meEeLiZsFir_, meEeNoise_, meEeNoZsRecVsSimE_, meEeRecE_, meEeRecEHitXtal_, meEeRecVsSimE_, meEeSimE_, meForcedRu_, meFullRoRu_, meZs1Ru_, nEeFROCnt_, nEeX, nEeY, nEeZsErrors_, nEeZsErrorsType1_, nEndcaps, EcalSelectiveReadoutValidation::energiesEe_t::noZsRecE, PV3DBase< T, PVType, FrameType >::phi(), EcalSelectiveReadoutValidation::energiesEe_t::phi, rad2deg, readOutUnitOf(), EcalSelectiveReadoutValidation::energiesEe_t::recE, edm::second(), EcalSelectiveReadoutValidation::energiesEe_t::simE, rpcPointValidation_cfi::simHit, EcalSelectiveReadoutValidation::energiesEe_t::simHit, EcalDataFrame::size(), edm::DataFrameContainer::size(), srApplicationErrorLog_, EcalSrFlag::SRF_FORCED_MASK, EcalSrFlag::SRF_FULL, EcalSrFlag::SRF_ZS1, findQualityFiles::v, EcalSrFlag::value(), withEeSimHit_, phase1PixelTopology::xOffset, xtalGraphX(), xtalGraphY(), EcalScDetId::zside(), EEDetId::zside(), and ecaldqm::zside().

Referenced by analyze().

                                                   {
  bool eventError = false;
  nEeZsErrors_ = 0;
  nEeZsErrorsType1_ = 0;

  for(int iZ0=0; iZ0<nEndcaps; ++iZ0){
    for(int iX0=0; iX0<nEeX; ++iX0){
      for(int iY0=0; iY0<nEeY; ++iY0){
        eeEnergies[iZ0][iX0][iY0].noZsRecE = -numeric_limits<double>::max();
        eeEnergies[iZ0][iX0][iY0].recE = -numeric_limits<double>::max();
        eeEnergies[iZ0][iX0][iY0].simE = 0; //must be set to zero.
        eeEnergies[iZ0][iX0][iY0].simHit = 0;
        eeEnergies[iZ0][iX0][iY0].gain12   = false;
      }
    }
  }

  // gets the endcap geometry:
  edm::ESHandle<CaloGeometry> geoHandle;
  es.get<CaloGeometryRecord>().get(geoHandle);
  const CaloSubdetectorGeometry *geometry_p
    = (*geoHandle).getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
//CaloSubdetectorGeometry const& geometry = *geometry_p;

  //EE unsupressed digis:
  for (unsigned int digis=0; digis<eeNoZsDigis_->size(); ++digis){

    EEDataFrame frame = (*eeNoZsDigis_)[digis];
    int iX0 = iXY2cIndex(frame.id().ix());
    int iY0 = iXY2cIndex(frame.id().iy());
    int iZ0 = frame.id().zside()>0?1:0;

    if(iX0<0 || iX0>=nEeX){
      edm::LogError("EcalSrValid") << "iX0 (= " << iX0 << ") is out of range ("
                                   << "[0," << nEeX -1 << "]\n";
    }
    if(iY0<0 || iY0>=nEeY){
      edm::LogError("EcalSrValid") << "iY0 (= " << iY0 << ") is out of range ("
                                   << "[0," << nEeY -1 << "]\n";
    }
    //    cout << "EE no ZS energy computation..." ;
    eeEnergies[iZ0][iX0][iY0].noZsRecE = frame2Energy(frame);

    eeEnergies[iZ0][iX0][iY0].gain12 = true;
    for(int i = 0; i< frame.size(); ++i){
      const int gain12Code = 0x1;
      if(frame[i].gainId()!=gain12Code) eeEnergies[iZ0][iX0][iY0].gain12 =  false;
    }

    const GlobalPoint xtalPos
      = geometry_p->getGeometry(frame.id())->getPosition();

    eeEnergies[iZ0][iX0][iY0].phi = rad2deg*((double)xtalPos.phi());
    eeEnergies[iZ0][iX0][iY0].eta = xtalPos.eta();
  }

  //EE rec hits:
  if(!localReco_){
    for(RecHitCollection::const_iterator it
          = eeRecHits_->begin();
        it != eeRecHits_->end(); ++it){
      const RecHit& hit = *it;
      int iX0 = iXY2cIndex(static_cast<const EEDetId&>(hit.id()).ix());
      int iY0 = iXY2cIndex(static_cast<const EEDetId&>(hit.id()).iy());
      int iZ0 = static_cast<const EEDetId&>(hit.id()).zside()>0?1:0;

      if(iX0<0 || iX0>=nEeX){
        LogError("EcalSrValid") << "iX0 (= " << iX0 << ") is out of range ("
                                << "[0," << nEeX -1 << "]\n";
      }
      if(iY0<0 || iY0>=nEeY){
        LogError("EcalSrValid") << "iY0 (= " << iY0 << ") is out of range ("
                                << "[0," << nEeY -1 << "]\n";
      }
      //    cout << "EE no ZS energy computation..." ;
      eeEnergies[iZ0][iX0][iY0].recE = hit.energy();
    }
  }

  //EE sim hits:
  for(vector<PCaloHit>::const_iterator it = eeSimHits_->begin();
      it != eeSimHits_->end(); ++it){
    const PCaloHit& simHit = *it;
    EEDetId detId(simHit.id());
    int iX = detId.ix();
    int iX0 =iXY2cIndex(iX);
    int iY = detId.iy();
    int iY0 = iXY2cIndex(iY);
    int iZ0 = detId.zside()>0?1:0;
    eeEnergies[iZ0][iX0][iY0].simE += simHit.energy();
    ++eeEnergies[iZ0][iX0][iY0].simHit;
  }

  bool EEcrystalShot[nEeX][nEeY][2];
  pair<int,int> EExtalCoor[nEeX][nEeY][2];  

  for(int iEeZ=0; iEeZ<2; ++iEeZ){
    for(int iEeX=0; iEeX<nEeX; ++iEeX){
      for(int iEeY=0; iEeY<nEeY; ++iEeY){
        EEcrystalShot[iEeX][iEeY][iEeZ] = false;
        EExtalCoor[iEeX][iEeY][iEeZ] = make_pair(0,0); 
      }
    }
  }

  //EE suppressed digis
  for(EEDigiCollection::const_iterator it = eeDigis_->begin();
      it != eeDigis_->end(); ++it){
    const EEDataFrame& frame = *it;
    int iX0 = iXY2cIndex(static_cast<const EEDetId&>(frame.id()).ix());
    int iY0 = iXY2cIndex(static_cast<const EEDetId&>(frame.id()).iy());
    int iZ0 = static_cast<const EEDetId&>(frame.id()).zside()>0?1:0;
    if(iX0<0 || iX0>=nEeX){
      LogError("EcalSrValid") << "iX0 (= " << iX0 << ") is out of range ("
                              << "[0," << nEeX -1 << "]\n";
    }
    if(iY0<0 || iY0>=nEeY){
      LogError("EcalSrValid") << "iY0 (= " << iY0 << ") is out of range ("
                              << "[0," << nEeY -1 << "]\n";
    }

    if(!EEcrystalShot[iX0][iY0][iZ0]){
      EEcrystalShot[iX0][iY0][iZ0] = true;
      EExtalCoor[iX0][iY0][iZ0] = make_pair(xtalGraphX(frame.id()), xtalGraphY(frame.id())); 
    } else{
      cout << "Error: several digi for same crystal!";
      abort();
    }

    if(localReco_){
      eeEnergies[iZ0][iX0][iY0].recE = frame2Energy(frame);
    }

    eeEnergies[iZ0][iX0][iY0].gain12 = true;
    for(int i = 0; i< frame.size(); ++i){
      const int gain12Code = 0x1;
      if(frame[i].gainId()!=gain12Code){
        eeEnergies[iZ0][iX0][iY0].gain12 =  false;
      }
    }

    EESrFlagCollection::const_iterator srf
      = eeSrFlags_->find(readOutUnitOf(frame.id()));

    bool highInterest = false;


    if(srf==eeSrFlags_->end()) continue;

    if(srf!=eeSrFlags_->end()){
      highInterest = ((srf->value() & ~EcalSrFlag::SRF_FORCED_MASK)
                      == EcalSrFlag::SRF_FULL);
    }

    if(highInterest){
      fill(meEeHiZsFir_, dccZsFIR(frame, firWeights_, firstFIRSample_, nullptr));
    } else{
      int v = dccZsFIR(frame, firWeights_, firstFIRSample_, nullptr);
      fill(meEeLiZsFir_, v);
      if(v < eeZsThr_){
        eventError = true;
        ++nEeZsErrors_;
        pair<int,int> ru = dccCh(frame.id());
        if(isRuComplete_[ru.first][ru.second-1]) ++nEeZsErrorsType1_;
        if(nEeZsErrors_ < 3){
          srApplicationErrorLog_ << event.id() << ", "
                                 << "RU " << frame.id() << ", "
                                 << "DCC " << ru.first
                                 << " Ch : " << ru.second << ": "
                                 << "LI channel under ZS threshold.\n";
        }
        if(nEeZsErrors_==3){
          srApplicationErrorLog_ << event.id() << ": "
                                 << "more ZS errors for this event...\n";
        }
      }
    }
  } //next ZS digi.

  for(int iEeZ=0; iEeZ<2; ++iEeZ){
    for(int iEeX=0; iEeX<nEeX; ++iEeX){
      for(int iEeY=0; iEeY<nEeY; ++iEeY){
        fill(meChOcc_, EExtalCoor[iEeX][iEeY][iEeZ].first, EExtalCoor[iEeX][iEeY][iEeZ].second,EEcrystalShot[iEeX][iEeY][iEeZ] ? 1 : 0);
      }
    }
  }

  for(int iZ0=0; iZ0<nEndcaps; ++iZ0){
    for(int iX0=0; iX0<nEeX; ++iX0){
      for(int iY0=0; iY0<nEeY; ++iY0){
        double recE = eeEnergies[iZ0][iX0][iY0].recE;
        if(recE==-numeric_limits<double>::max()) continue; //not a crystal or ZS
        fill(meEeRecE_, eeEnergies[iZ0][iX0][iY0].recE);

        fill(meEeEMean_, ievt_+1,
             eeEnergies[iZ0][iX0][iY0].recE);

        if(withEeSimHit_){
          if(!eeEnergies[iZ0][iX0][iY0].simHit){//noise only crystal channel
            fill(meEeNoise_, eeEnergies[iZ0][iX0][iY0].noZsRecE);
          } else{
            fill(meEeSimE_, eeEnergies[iZ0][iX0][iY0].simE);
            fill(meEeRecEHitXtal_, eeEnergies[iZ0][iX0][iY0].recE);
          }
          fill(meEeRecVsSimE_, eeEnergies[iZ0][iX0][iY0].simE,
               eeEnergies[iZ0][iX0][iY0].recE);
          fill(meEeNoZsRecVsSimE_, eeEnergies[iZ0][iX0][iY0].simE,
               eeEnergies[iZ0][iX0][iY0].noZsRecE);
        }
      }
    }
  }

  int EEZs1RuCount[2][20][20];
  int EEFullRuCount[2][20][20];
  int EEForcedRuCount[2][20][20];
  for(int iZ(0); iZ < 2; iZ++){
    for(int iX(0);iX < 20; iX++){
      for(int iY(0);iY < 20; iY++){
        EEZs1RuCount[iZ][iX][iY] = 0;
        EEFullRuCount[iZ][iX][iY] = 0;
        EEForcedRuCount[iZ][iX][iY] = 0;
      }
    }
  }

  nEeFROCnt_ = 0;
  char eeSrfMark[2][20][20];
  bzero(eeSrfMark, sizeof(eeSrfMark));
  //Filling RU histo
  for(EESrFlagCollection::const_iterator it = eeSrFlags_->begin();
      it != eeSrFlags_->end(); ++it){
    const EESrFlag& srf = *it;
    // srf.id() is EcalScDetId; 1 <= ix <= 20 1 <= iy <= 20
    int iX = srf.id().ix();
    int iY = srf.id().iy();
    int zside = srf.id().zside(); //-1 for EE-, +1 for EE+
    if(iX<1 || iY > 100) throw cms::Exception("EcalSelectiveReadoutValidation")
                           << "Found an endcap SRF with an invalid det ID: " << srf.id() << ".\n";
    ++eeSrfMark[zside>0?1:0][iX-1][iY-1];
    if(eeSrfMark[zside>0?1:0][iX-1][iY-1] > 1) throw cms::Exception("EcalSelectiveReadoutValidation")
                                                 << "Duplicate SRF for supercrystal " << srf.id() << ".\n";
    int flag = srf.value() & ~EcalSrFlag::SRF_FORCED_MASK;
    if(flag == EcalSrFlag::SRF_ZS1){
      EEZs1RuCount[zside>0?1:0][iX-1][iY-1] += 1;
    }

    if(flag == EcalSrFlag::SRF_FULL){
      EEFullRuCount[zside>0?1:0][iX-1][iY-1] += 1;
      ++nEeFROCnt_;
    }
      
    if(srf.value() & EcalSrFlag::SRF_FORCED_MASK){
      EEForcedRuCount[zside>0?1:0][iX-1][iY-1] += 1;
    }
  }
  for(int iZ(0); iZ < 2; iZ++){
    int xOffset(iZ == 0 ? -40 : 20);
    for(int iX(0);iX < 20; iX++){
      for(int iY(0);iY < 20; iY++){
        int GraphX = (iX + 1)+ xOffset;
        int GraphY = (iY+1);
        fill(meZs1Ru_, GraphX, GraphY, EEZs1RuCount[iZ][iX][iY]);
        fill(meFullRoRu_, GraphX, GraphY, EEFullRuCount[iZ][iX][iY]);
        fill(meForcedRu_, GraphX, GraphY, EEForcedRuCount[iZ][iX][iY]);
      }
    }
  }
  
  if(eventError) srApplicationErrorLog_ << event.id()
                                        << ": " << nEeZsErrors_
                                        << " ZS-flagged EE channels under "
                   "the ZS threshold, whose " << nEeZsErrorsType1_
                                        << " in a complete RU.\n";
} //end of analyzeEE

void EcalSelectiveReadoutValidation::analyzeTP ( const edm::Event &  event, const edm::EventSetup &  es  )

private

Trigger primitive analysis. To be called for each event.

Parameters

event	EDM event
es	event setup

Definition at line 1307 of file EcalSelectiveReadoutValidation.cc.

References cIndex2iTtEta(), cIndex2iTtPhi(), DetId::Ecal, EcalBarrel, EcalEndcap, fill(), edm::EventSetup::get(), EcalTPGLutIdMap::getMap(), EcalTPGGroups::getMap(), EcalTPGPhysicsConst::getMap(), MonitorElement::getTProfile(), mps_fire::i, iTtEta2cIndex(), iTtPhi2cIndex(), lumiPlot::lut, meForcedTtf_, meHiTtf_, meLiTtf_, meMiTtf_, meTp_, meTpMap_, meTpVsEtSum_, meTtf_, meTtfVsEtSum_, meTtfVsTp_, nTtEta, nTtPhi, edm::ESHandle< T >::product(), l1t::HGCalTowerID::rawId(), DataFormats::L1THGCal::towerId, tpInGeV_, tps_, and ttEtSums.

Referenced by analyze().

                                               {
  int TTFlagCount[8];
  int LiTTFlagCount[nTtEta][nTtPhi];
  int MiTTFlagCount[nTtEta][nTtPhi];
  int HiTTFlagCount[nTtEta][nTtPhi];
  for(int iTTFlag(0);iTTFlag <8; iTTFlag++){
    TTFlagCount[iTTFlag]=0;
  }
  for(int iTtEta(0); iTtEta < nTtEta; iTtEta++){
    for(int iTtPhi(0); iTtPhi < nTtPhi; iTtPhi++){
      LiTTFlagCount[iTtEta][iTtPhi] = 0;
      MiTTFlagCount[iTtEta][iTtPhi] = 0;
      HiTTFlagCount[iTtEta][iTtPhi] = 0;
    }
  }
  int tpEtCount[100];
  for(int iEt(0);iEt < 100; iEt++){
     tpEtCount[iEt] = 0;
  }

  edm::ESHandle<EcalTPGPhysicsConst> physHandle;
  es.get<EcalTPGPhysicsConstRcd>().get( physHandle );
  const EcalTPGPhysicsConstMap & physMap = physHandle.product()->getMap() ;

  edm::ESHandle<EcalTPGLutGroup> lutGrpHandle;
  es.get<EcalTPGLutGroupRcd>().get( lutGrpHandle );
  const EcalTPGGroups::EcalTPGGroupsMap & lutGrpMap = lutGrpHandle.product()->getMap() ;  

  edm::ESHandle<EcalTPGLutIdMap> lutMapHandle;
  es.get<EcalTPGLutIdMapRcd>().get( lutMapHandle );
  const EcalTPGLutIdMap::EcalTPGLutMap & lutMap = lutMapHandle.product()->getMap() ;  

  EcalTPGPhysicsConstMapIterator ebItr(physMap.find(DetId(DetId::Ecal,EcalBarrel).rawId()));
  double lsb10bitsEB(ebItr == physMap.end() ? 0. : ebItr->second.EtSat / 1024.);
  EcalTPGPhysicsConstMapIterator eeItr(physMap.find(DetId(DetId::Ecal,EcalEndcap).rawId()));
  double lsb10bitsEE(eeItr == physMap.end() ? 0. : eeItr->second.EtSat / 1024.);

  for(EcalTrigPrimDigiCollection::const_iterator it = tps_->begin(); it != tps_->end(); ++it){

    double tpEt;
    if(tpInGeV_){
      EcalTrigTowerDetId const& towerId(it->id());
      unsigned int ADC = it->compressedEt();

      double lsb10bits(0.);
      if(towerId.subDet() == EcalBarrel) lsb10bits = lsb10bitsEB;
      else if(towerId.subDet() == EcalEndcap) lsb10bits = lsb10bitsEE;

      int tpg10bits = 0 ;
      EcalTPGGroups::EcalTPGGroupsMapItr itgrp = lutGrpMap.find(towerId.rawId()) ;
      uint32_t lutGrp = 999 ;
      if (itgrp != lutGrpMap.end()) lutGrp = itgrp->second ;

      EcalTPGLutIdMap::EcalTPGLutMapItr itLut = lutMap.find(lutGrp) ;
      if (itLut != lutMap.end()) {
        const unsigned int * lut = (itLut->second).getLut() ;
        for (unsigned int i=0 ; i<1024 ; i++)
          if (ADC == (0xff & lut[i])) {
            tpg10bits = i ;
            break ;
          }
      }

      tpEt = lsb10bits * tpg10bits;
    }
    else{
      tpEt = it->compressedEt();
    }
    int iEta = it->id().ieta();
    int iEta0 = iTtEta2cIndex(iEta);
    int iPhi = it->id().iphi();
    int iPhi0 = iTtPhi2cIndex(iPhi);
    double etSum = ttEtSums[iEta0][iPhi0];

    int iE = meTp_->getTProfile()->FindFixBin(tpEt);
    if ((iE >= 0) && (iE < 100)) {
      ++tpEtCount[iE];
    } else {
      // FindFixBin might return an overflow bin (outside tpEtCount range).
      // To prevent a memory overflow / segfault, these values are ignored.
      //std::cout << "EcalSelectiveReadoutValidation: Invalid iE value: " << iE << std::endl;
    }

    fill(meTpVsEtSum_, etSum, tpEt);
    ++TTFlagCount[it->ttFlag()];
    if((it->ttFlag() & 0x3) == 0){
      LiTTFlagCount[iEta0][iPhi0] += 1;
    }
    else if((it->ttFlag() & 0x3) == 1){
      MiTTFlagCount[iEta0][iPhi0] += 1;
    }
    else if((it->ttFlag() & 0x3) == 3){
      HiTTFlagCount[iEta0][iPhi0] += 1;
    }
    if((it->ttFlag() & 0x4)){
      fill(meForcedTtf_, iEta, iPhi);
    }

    fill(meTtfVsTp_, tpEt, it->ttFlag());
    fill(meTtfVsEtSum_, etSum, it->ttFlag());
    fill(meTpMap_, iEta, iPhi, tpEt, 1.);
  }

  for(int ittflag(0); ittflag < 8; ittflag++){
    fill(meTtf_, ittflag, TTFlagCount[ittflag]);
  }
  for(int iTtEta(0); iTtEta < nTtEta; iTtEta++){
    for(int iTtPhi(0); iTtPhi < nTtPhi; iTtPhi++){ 
      fill(meLiTtf_, cIndex2iTtEta(iTtEta), cIndex2iTtPhi(iTtPhi), LiTTFlagCount[iTtEta][iTtPhi]);
      fill(meMiTtf_, cIndex2iTtEta(iTtEta), cIndex2iTtPhi(iTtPhi), MiTTFlagCount[iTtEta][iTtPhi]);
      fill(meHiTtf_, cIndex2iTtEta(iTtEta), cIndex2iTtPhi(iTtPhi), HiTTFlagCount[iTtEta][iTtPhi]);
    }
  }
  if(tpInGeV_){
    for(int iE(0);iE<100; iE++){
      fill(meTp_, iE, tpEtCount[iE]);
    }
  } else{
    for(int iE(0); iE<40; iE++){
      fill(meTp_, iE, tpEtCount[iE]);
    }
  }
}

MonitorElement * EcalSelectiveReadoutValidation::book1D ( DQMStore::IBooker &  ibook, const std::string &  name, const std::string &  title, int  nbins, double  xmin, double  xmax  )

private

Definition at line 1841 of file EcalSelectiveReadoutValidation.cc.

References DQMStore::IBooker::book1D(), Exception, dataset::name, registerHist(), and mps_fire::result.

Referenced by bookHistograms(), and cIndex2iTtPhi().

                                                                                                                                                                {
  if(!registerHist(name, title)) return nullptr; //this histo is disabled
  MonitorElement* result = ibook.book1D(name, title, nbins, xmin, xmax);
  if(result==nullptr){
    throw cms::Exception("Histo")
      << "Failed to book histogram " << name;
  }
  return result;
}

MonitorElement * EcalSelectiveReadoutValidation::book2D ( DQMStore::IBooker &  ibook, const std::string &  name, const std::string &  title, int  nxbins, double  xmin, double  xmax, int  nybins, double  ymin, double  ymax  )

private

Definition at line 1851 of file EcalSelectiveReadoutValidation.cc.

References DQMStore::IBooker::book2D(), Exception, dataset::name, registerHist(), and mps_fire::result.

Referenced by bookHistograms(), and cIndex2iTtPhi().

                                                                                                                                                                                                       {
  if(!registerHist(name, title)) return nullptr; //this histo is disabled
  MonitorElement* result = ibook.book2D(name, title, nxbins, xmin, xmax,
                    nybins, ymin, ymax);
  if(result==nullptr){
    throw cms::Exception("Histo")
      << "Failed to book histogram " << name;
  }
  return result;
}

MonitorElement * EcalSelectiveReadoutValidation::bookFloat ( DQMStore::IBooker &  ibook, const std::string &  name  )

private

Wrappers to the book methods of the DQMStore DQM histogramming interface.

Definition at line 1830 of file EcalSelectiveReadoutValidation.cc.

References DQMStore::IBooker::bookFloat(), Exception, dataset::name, registerHist(), and mps_fire::result.

Referenced by bookHistograms(), and cIndex2iTtPhi().

                                                                                                      {
  if(!registerHist(name, "")) return nullptr; //this histo is disabled
  MonitorElement* result = ibook.bookFloat(name);
  if(result==nullptr){
    throw cms::Exception("DQM")
      << "Failed to book integer DQM monitor element" << name;
  }
  return result;
}

void EcalSelectiveReadoutValidation::bookHistograms ( DQMStore::IBooker &  i, edm::Run const &  , edm::EventSetup const &    )

override

Definition at line 851 of file EcalSelectiveReadoutValidation.cc.

References availableHistList_, book1D(), book2D(), bookFloat(), bookProfile(), bookProfile2D(), ebZsThr_, eeZsThr_, histList_, meChOcc_, meCompleteZSCnt_, meCompleteZSMap_, meCompleteZSRateMap_, meDccHiVol_, meDccLiVol_, meDccVol_, meDccVolFromData_, meDroppedFROCnt_, meDroppedFROMap_, meDroppedFRORateMap_, meEbEMean_, meEbFullRoCnt_, meEbHiZsFir_, meEbLiZsFir_, meEbNoise_, meEbNoZsRecVsSimE_, meEbRecE_, meEbRecEHitXtal_, meEbRecVsSimE_, meEbSimE_, meEbZsErrCnt_, meEbZsErrType1Cnt_, meEeEMean_, meEeFullRoCnt_, meEeHiZsFir_, meEeLiZsFir_, meEeNoise_, meEeNoZsRecVsSimE_, meEeRecE_, meEeRecEHitXtal_, meEeRecVsSimE_, meEeSimE_, meEeZsErrCnt_, meEeZsErrType1Cnt_, meForcedRu_, meForcedTtf_, meFullRoCnt_, meFullRoRu_, meHiTtf_, meIncompleteFROCnt_, meIncompleteFROMap_, meIncompleteFRORateMap_, meL1aRate_, meLiTtf_, meMiTtf_, meSRFlagsConsistency_, meTp_, meTpMap_, meTpVsEtSum_, meTtf_, meTtfVsEtSum_, meTtfVsTp_, meVol_, meVolB_, meVolBHI_, meVolBLI_, meVolE_, meVolEHI_, meVolELI_, meVolHI_, meVolLI_, meZs1Ru_, meZsErrCnt_, meZsErrType1Cnt_, nDccs_, printAvailableHists(), alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, and tpInGeV_.

                                                                                                                {

  meL1aRate_ = bookFloat(ibooker, "l1aRate_"); 

  meDccVol_ = bookProfile(ibooker, "hDccVol", //"EcalDccEventSizeComputed",
              "ECAL DCC event fragment size;Dcc id; "
              "<Event size> (kB)", nDccs_, .5, .5+nDccs_);

  meDccLiVol_ = bookProfile(ibooker, "hDccLiVol",
                            "LI channel payload per DCC;Dcc id; "
                            "<Event size> (kB)", nDccs_, .5, .5+nDccs_);

  meDccHiVol_ = bookProfile(ibooker, "hDccHiVol",
                            "HI channel payload per DCC;Dcc id; "
                            "<Event size> (kB)", nDccs_, .5, .5+nDccs_);

  meDccVolFromData_ = bookProfile(ibooker, "hDccVolFromData", //"EcalDccEventSize",
                  "ECAL DCC event fragment size;Dcc id; "
                  "<Event size> (kB)", nDccs_, .5, .5+nDccs_);

  meVolBLI_ = book1D(ibooker, "hVolBLI",// "EBLowInterestPayload",
             "ECAL Barrel low interest crystal data payload;"
             "Event size (kB);Nevts",
             100, 0., 200.);

  meVolELI_ = book1D(ibooker, "hVolELI", //"EELowInterestPayload",
             "Endcap low interest crystal data payload;"
             "Event size (kB);Nevts",
             100, 0., 200.);

  meVolLI_ = book1D(ibooker, "hVolLI", //"EcalLowInterestPayload",
            "ECAL low interest crystal data payload;"
            "Event size (kB);Nevts",
            100, 0., 200.);

  meVolBHI_ = book1D(ibooker, "hVolBHI", //"EBHighInterestPayload",
             "Barrel high interest crystal data payload;"
             "Event size (kB);Nevts",
             100, 0., 200.);

  meVolEHI_ = book1D(ibooker, "hVolEHI", //"EEHighInterestPayload",
             "Endcap high interest crystal data payload;"
             "Event size (kB);Nevts",
             100, 0., 200.);

  meVolHI_ = book1D(ibooker, "hVolHI", //"EcalHighInterestPayload",
            "ECAL high interest crystal data payload;"
            "Event size (kB);Nevts",
            100, 0., 200.);

  meVolB_ = book1D(ibooker, "hVolB", //"EBEventSize",
           "Barrel data volume;Event size (kB);Nevts",
           100, 0., 200.);

  meVolE_ = book1D(ibooker, "hVolE", //"EEEventSize",
           "Endcap data volume;Event size (kB);Nevts",
           100, 0., 200.);

  meVol_ = book1D(ibooker, "hVol", //"EcalEventSize",
          "ECAL data volume;Event size (kB);Nevts",
          100, 0., 200.);

  meChOcc_ = bookProfile2D(ibooker, "h2ChOcc", //"EcalChannelOccupancy",
            "ECAL crystal channel occupancy after zero suppression;"
            "iX -200 / iEta / iX + 100;"
            "iY / iPhi (starting from -10^{o}!);"
            "Event count rate",
            401, -200.5, 200.5,
            360, .5, 360.5);

  //TP
  string tpUnit;
  if(tpInGeV_) tpUnit = string("GeV"); else tpUnit = string("TP hw unit");
  string title;
  title = string("Trigger primitive TT E_{T};E_{T} ")
    + tpUnit + string(";Event Count");
  meTp_ = bookProfile(ibooker, "hTp", //"EcalTriggerPrimitiveEt",
         title,
         (tpInGeV_?100:40), 0., (tpInGeV_?10.:40.));

  meTtf_ = bookProfile(ibooker, "hTtf", //"EcalTriggerTowerFlag",
          "Trigger primitive TT flag;Flag number;Event count",
          8, -.5, 7.5);

  title = string("Trigger tower flag vs TP;E_{T}(TT) (")
    + tpUnit + string(");Flag number");
  meTtfVsTp_ = book2D(ibooker, "h2TtfVsTp",
              title,
              100, 0., (tpInGeV_?10.:40.),
              8, -.5, 7.5);

  meTtfVsEtSum_ = book2D(ibooker, "h2TtfVsEtSum",
             "Trigger tower flag vs #sumE_{T};"
             "E_{T}(TT) (GeV);"
             "TTF",
             100, 0., 10.,
             8, -.5, 7.5);
  title = string("Trigger primitive Et (TP) vs #sumE_{T};"
         "E_{T} (sum) (GeV);"
         "E_{T} (TP) (") + tpUnit + string (")");

  meTpVsEtSum_ = book2D(ibooker, "h2TpVsEtSum",
            title,
            100, 0., 10.,
            100, 0., (tpInGeV_?10.:40.));

  title = string("Trigger primitive E_{T};"
                 "iEta;"
                 "iPhi;"
                 "E_{T} (TP) (") + tpUnit + string (")");
  meTpMap_ = bookProfile2D(ibooker, "h2Tp",
                           title,
                           57, -28.5, 28.5,
               72, .5, 72.5);

  //SRF
  meFullRoRu_ = book2D(ibooker, "h2FRORu", //"EcalFullReadoutSRFlagMap",
               "Full Read-out readout unit;"
               "iX - 40 / iEta / iX + 20;"
               "iY / iPhi (iPhi = 1 at phi = 0 rad);"
               "Event count",
                       80, -39.5, 40.5,
               72, .5, 72.5);

  meFullRoCnt_ = book1D(ibooker, "hFROCnt",
                        "Number of Full-readout-flagged readout units;"
                        "FRO RU count;Event count",
                        300, -.5, 299.5);

  meEbFullRoCnt_ = book1D(ibooker, "hEbFROCnt",
                          "Number of EB Full-readout-flagged readout units;"
                          "FRO RU count;Event count",
                          200, -.5, 199.5);

  meEeFullRoCnt_ = book1D(ibooker, "hEeFROCnt",
                          "Number of EE Full-readout-flagged readout units;"
                          "FRO RU count;Event count",
                          200, -.5, 199.5);

  meZs1Ru_ = book2D(ibooker, "h2Zs1Ru", //"EbZeroSupp1SRFlagMap",
            "Readout unit with ZS-thr-1 flag;"
            "iX - 40 / iEta / iX + 20;"
            "iY0 / iPhi0 (iPhi = 1 at phi = 0 rad);"
            "Event count",
                    80, -39.5, 40.5,
            72, .5, 72.5);

  meForcedRu_ = book2D(ibooker, "h2ForcedRu", //"EcalReadoutUnitForcedBitMap",
               "ECAL readout unit with forced bit of SR flag on;"
               "iX - 40 / iEta / iX + 20;"
               "iY / iPhi (iPhi = 1 at phi = 0 rad);"
               "Event count",
                       80, -39.5, 40.5,
               72, .5, 72.5);

  meLiTtf_ = bookProfile2D(ibooker, "h2LiTtf", //"EcalLowInterestTriggerTowerFlagMap",
            "Low interest trigger tower flags;"
            "iEta;"
            "iPhi;"
            "Event count",
                    57, -28.5, 28.5,
            72, .5, 72.5);

  meMiTtf_ = bookProfile2D(ibooker, "h2MiTtf", //"EcalMidInterestTriggerTowerFlagMap",
            "Mid interest trigger tower flags;"
            "iEta;"
            "iPhi;"
            "Event count",
            57, -28.5, 28.5,
                    72, .5, 72.5);

  meHiTtf_ = bookProfile2D(ibooker, "h2HiTtf", //"EcalHighInterestTriggerTowerFlagMap",
            "High interest trigger tower flags;"
            "iEta;"
            "iPhi;"
            "Event count",
                    57, -28.5, 28.5,
            72, .5, 72.5);

  meForcedTtf_ = book2D(ibooker, "h2ForcedTtf", //"EcalTtfForcedBitMap",
                        "Trigger tower flags with forced bit set;"
                        "iEta;"
                        "iPhi;"
                        "Event count",
                        57, -28.5, 28.5,
                        72, .5, 72.5);


  const float ebMinNoise = -1.;
  const float ebMaxNoise = 1.;

  const float eeMinNoise = -1.;
  const float eeMaxNoise = 1.;

  const float ebMinE = ebMinNoise;
  const float ebMaxE = ebMaxNoise;

  const float eeMinE = eeMinNoise;
  const float eeMaxE = eeMaxNoise;

  const int evtMax = 500;

  meEbRecE_ = book1D(ibooker, "hEbRecE",
             "Crystal reconstructed energy;E (GeV);Event count",
             100, ebMinE, ebMaxE);

  meEbEMean_ = bookProfile(ibooker, "hEbEMean",
               "EE <E_hit>;event #;<E_hit> (GeV)",
               evtMax, .5, evtMax + .5);

  meEbNoise_ = book1D(ibooker, "hEbNoise",
              "Crystal noise "
              "(rec E of crystal without deposited energy)"
                      ";Rec E (GeV);Event count",
              100, ebMinNoise, ebMaxNoise);

  meEbLiZsFir_   = book1D(ibooker, "zsEbLiFIRemu",
              "Emulated ouput of ZS FIR filter for EB "
              "low interest crystals;"
              "ADC count*4;"
              "Event count",
              60, -30, 30);

  meEbHiZsFir_ = book1D(ibooker, "zsEbHiFIRemu",
                        "Emulated ouput of ZS FIR filter for EB "
                        "high interest crystals;"
                        "ADC count*4;"
                        "Event count",
                        60, -30, 30);

  //TODO: Fill this histogram...
//   meEbIncompleteRUZsFir_ = book1D(ibooker, "zsEbIncompleteRUFIRemu",
//                                   "Emulated ouput of ZS FIR filter for EB "
//                                   "incomplete FRO-flagged RU;"
//                                   "ADC count*4;"
//                                   "Event count",
//                                   60, -30, 30);

  meEbSimE_ = book1D(ibooker, "hEbSimE", "EB hit crystal simulated energy",
             100, ebMinE, ebMaxE);

  meEbRecEHitXtal_ = book1D(ibooker, "hEbRecEHitXtal",
                "EB rec energy of hit crystals",
                100, ebMinE, ebMaxE);

  meEbRecVsSimE_ = book2D(ibooker, "hEbRecVsSimE",
              "Crystal simulated vs reconstructed energy;"
              "Esim (GeV);Erec GeV);Event count",
              100, ebMinE, ebMaxE,
              100, ebMinE, ebMaxE);

  meEbNoZsRecVsSimE_ = book2D(ibooker, "hEbNoZsRecVsSimE",
                  "Crystal no-zs simulated vs reconstructed "
                  "energy;"
                  "Esim (GeV);Erec GeV);Event count",
                  100, ebMinE, ebMaxE,
                  100, ebMinE, ebMaxE);

  meEeRecE_ = book1D(ibooker, "hEeRecE",
             "EE crystal reconstructed energy;E (GeV);"
             "Event count",
             100, eeMinE, eeMaxE);

  meEeEMean_ = bookProfile(ibooker, "hEeEMean",
               "<E_{EE hit}>;event;<E_{hit}> (GeV)",
               evtMax, .5, evtMax + .5);


  meEeNoise_ = book1D(ibooker, "hEeNoise",
              "EE crystal noise "
              "(rec E of crystal without deposited energy);"
              "E (GeV);Event count",
              200, eeMinNoise, eeMaxNoise);

  meEeLiZsFir_   = book1D(ibooker, "zsEeLiFIRemu",
              "Emulated ouput of ZS FIR filter for EE "
              "low interest crystals;"
              "ADC count*4;"
              "Event count",
              60, -30, 30);

  meEeHiZsFir_   = book1D(ibooker, "zsEeHiFIRemu",
              "Emulated ouput of ZS FIR filter for EE "
              "high interest crystals;"
              "ADC count*4;"
              "Event count",
              60, -30, 30);

  //TODO: Fill this histogram...
//   meEeIncompleteRUZsFir_ = book1D(ibooker, "zsEeIncompleteRUFIRemu",
//                                     "Emulated ouput of ZS FIR filter for EE "
//                                     "incomplete FRO-flagged RU;"
//                                     "ADC count*4;"
//                                     "Event count",
//                                     60, -30, 30);


  meEeSimE_ = book1D(ibooker, "hEeSimE", "EE hit crystal simulated energy",
             100, eeMinE, eeMaxE);

  meEeRecEHitXtal_ = book1D(ibooker, "hEeRecEHitXtal",
                "EE rec energy of hit crystals",
                100, eeMinE, eeMaxE);

  meEeRecVsSimE_ = book2D(ibooker, "hEeRecVsSimE",
              "EE crystal simulated vs reconstructed energy;"
              "Esim (GeV);Erec GeV);Event count",
              100, eeMinE, eeMaxE,
              100, eeMinE, eeMaxE);

  meEeNoZsRecVsSimE_ = book2D(ibooker, "hEeNoZsRecVsSimE",
                  "EE crystal no-zs simulated vs "
                  "reconstructed "
                  "energy;Esim (GeV);Erec GeV);Event count",
                  100, eeMinE, eeMaxE,
                  100, eeMinE, eeMaxE);

  meSRFlagsConsistency_ = book2D(ibooker, "hSRAlgoErrorMap",
                 "TTFlags and SR Flags mismatch;"
                 "iX - 40 / iEta / iX + 20;"
                 "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                 "Event count",
                 80, -39.5, 40.5,
                 72, .5, 72.5);

  //Readout Units histos (interest/Ncrystals)
  meIncompleteFROMap_ = book2D(ibooker, "hIncompleteFROMap",
                               "Incomplete full-readout-flagged readout units;"
                               "iX - 40 / iEta / iX + 20;"
                               "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                               "Event count",
                               80, -39.5, 40.5,
                               72, .5, 72.5);

  meIncompleteFROCnt_ = book1D(ibooker, "hIncompleteFROCnt",
                               "Number of incomplete full-readout-flagged "
                               "readout units;"
                               "Number of RUs;Event count;",
                               200, -.5, 199.5);

  meIncompleteFRORateMap_
    = bookProfile2D(ibooker, "hIncompleteFRORateMap",
            "Incomplete full-readout-flagged readout units;"
            "iX - 40 / iEta / iX + 20;"
            "iY / iPhi (iPhi = 1 at phi = 0 rad);"
            "Incomplete error rate",
                    80, -39.5, 40.5,
            72, .5, 72.5);


  meDroppedFROMap_ = book2D(ibooker, "hDroppedFROMap",
                            "Dropped full-readout-flagged readout units;"
                            "iX - 40 / iEta / iX + 20;"
                            "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                            "Event count",
                            80, -39.5, 40.5,
                            72, .5, 72.5);

  meDroppedFROCnt_ = book1D(ibooker, "hDroppedFROCnt",
                            "Number of dropped full-readout-flagged "
                            "RU count;RU count;Event count",
                            200, -.5, 199.5);

  meCompleteZSCnt_ = book1D(ibooker, "hCompleteZsCnt",
                            "Number of zero-suppressed-flagged RU fully "
                            "readout;"
                            "RU count;Event count",
                            200, -.5, 199.5);

  stringstream buf;
  buf << "Number of LI EB channels below the " << ebZsThr_/4. << " ADC count ZS threshold;"
    "Channel count;Event count",
  meEbZsErrCnt_ = book1D(ibooker, "hEbZsErrCnt",
                         buf.str(),
                         200, -.5, 199.5);

  buf.str("");
  buf << "Number of LI EE channels below the " << eeZsThr_/4. << " ADC count ZS theshold;"
    "Channel count;Event count",
  meEeZsErrCnt_ = book1D(ibooker, "hEeZsErrCnt",
                         buf.str(),
                         200, -.5, 199.5);

  meZsErrCnt_ = book1D(ibooker, "hZsErrCnt",
                       "Number of LI channels below the ZS threshold;"
                       "Channel count;Event count",
                       200, -.5, 199.5);

  meEbZsErrType1Cnt_ = book1D(ibooker, "hEbZsErrType1Cnt",
                              "Number of EB channels below the ZS "
                              "threshold in a LI but fully readout RU;"
                              "Channel count;Event count;",
                              200, -.5, 199.5);

  meEeZsErrType1Cnt_ = book1D(ibooker, "hEeZsErrType1Cnt",
                              "Number EE channels below the ZS threshold"
                              " in a LI but fully readout RU;"
                              "Channel count;Event count",
                              200, -.5, 199.5);

  meZsErrType1Cnt_ = book1D(ibooker, "hZsErrType1Cnt",
                            "Number of LI channels below the ZS threshold "
                            "in a LI but fully readout RU;"
                            "Channel count;Event count",
                            200, -.5, 199.5);


  meDroppedFRORateMap_
    = bookProfile2D(ibooker, "hDroppedFRORateMap",
            "Dropped full-readout-flagged readout units"
                    "iX - 40 / iEta / iX + 20;"
                    "iY / iPhi (iPhi = 1 at phi = 0 rad);"
                    "Dropping rate",
                    80, -39.5, 40.5,
                    72, .5, 72.5);

  meCompleteZSMap_ = book2D(ibooker, "hCompleteZSMap",
             "Complete zero-suppressed-flagged readout units;"
             "iX - 40 / iEta / iX + 20;"
             "iY / iPhi (iPhi = 1 at phi = 0 rad);"
             "Event count",
             80, -39.5, 40.5,
             72, .5, 72.5);

  meCompleteZSRateMap_
    = bookProfile2D(ibooker, "hCompleteZSRate",
            "Complete zero-suppressed-flagged readout units;"
            "iX - 40 / iEta / iX + 20;"
            "iY / iPhi (iPhi = 1 at phi = 0 rad);"
            "Completeness rate",
            80, -39.5, 40.5,
            72, .5, 72.5);

  //print list of available histograms (must be called after
  //the bookXX methods):
  printAvailableHists();

  //check the histList parameter:
  stringstream s;
  for(set<string>::iterator it = histList_.begin();
      it!=histList_.end();
      ++it){
    if(*it!=string("all")
       && availableHistList_.find(*it)==availableHistList_.end()){
      s << (s.str().empty()?"":", ") << *it;
    }
  }
  if(!s.str().empty()){
    LogWarning("Configuration")
      << "Parameter 'histList' contains some unknown histogram(s). "
      "Check spelling. Following name were not found: "
      << s.str();
  }
}

MonitorElement * EcalSelectiveReadoutValidation::bookProfile ( DQMStore::IBooker &  ibook, const std::string &  name, const std::string &  title, int  nbins, double  xmin, double  xmax  )

private

Definition at line 1862 of file EcalSelectiveReadoutValidation.cc.

References DQMStore::IBooker::bookProfile(), Exception, dataset::name, registerHist(), and mps_fire::result.

Referenced by bookHistograms(), and cIndex2iTtPhi().

                                                                                                                                                                     {
  if(!registerHist(name, title)) return nullptr; //this histo is disabled
  MonitorElement* result = ibook.bookProfile(name, title, nbins, xmin, xmax,
                         0, 0, 0);
  if(result==nullptr){
    throw cms::Exception("Histo")
      << "Failed to book histogram " << name;
  }
  return result;
}

MonitorElement * EcalSelectiveReadoutValidation::bookProfile2D ( DQMStore::IBooker &  ibook, const std::string &  name, const std::string &  title, int  nbinx, double  xmin, double  xmax, int  nbiny, double  ymin, double  ymax, const char *  option = ""  )

private

Definition at line 1873 of file EcalSelectiveReadoutValidation.cc.

References DQMStore::IBooker::bookProfile2D(), Exception, dataset::name, registerHist(), and mps_fire::result.

Referenced by bookHistograms(), and cIndex2iTtPhi().

                                                                                                                                                                                                                                {
  if(!registerHist(name, title)) return nullptr; //this histo is disabled
  MonitorElement* result
    = ibook.bookProfile2D(name,
              title,
              nbinx, xmin, xmax,
              nbiny, ymin, ymax,
              0, 0, 0,
              option);
  if(result==nullptr){
    throw cms::Exception("Histo")
      << "Failed to book histogram " << name;
  }
  return result;
}

template<class T >

void EcalSelectiveReadoutValidation::checkSrApplication ( const edm::Event &  event, T &  srfs  )

private

Checks application of SR decision by the DCC.

Parameters

event	event currently analyzed.
srfs	Selective readou flags

Definition at line 2197 of file EcalSelectiveReadoutValidation.cc.

References dccCh(), fill(), RemoveAddSevLevel::flag, getCrystalCount(), meCompleteZSMap_, meCompleteZSRateMap_, meDroppedFROMap_, meDroppedFRORateMap_, meIncompleteFROMap_, meIncompleteFRORateMap_, minDccId_, nCompleteZS_, nDroppedFRO_, nIncompleteFRO_, nPerRu_, ruGraphX(), ruGraphY(), srApplicationErrorLog_, EcalSrFlag::SRF_FULL, EcalSrFlag::SRF_ZS1, and EcalSrFlag::SRF_ZS2.

Referenced by analyze().

                                                                {
  typedef typename T::const_iterator SrFlagCollectionConstIt;
  typedef typename T::key_type MyRuDetIdType;

  for(SrFlagCollectionConstIt itSrf = srfs.begin();
      itSrf != srfs.end(); ++itSrf){
    int flag = itSrf->value() & ~EcalSrFlag::SRF_FORCED_MASK;
    pair<int,int> ru = dccCh(itSrf->id());

    if(flag == EcalSrFlag::SRF_FULL){
      if(nPerRu_[ru.first-minDccId_][ru.second-1]==getCrystalCount(ru.first, ru.second)){ //no error
    fill(meIncompleteFRORateMap_, ruGraphX(itSrf->id()),
             ruGraphY(itSrf->id()), 0);
    fill(meDroppedFRORateMap_,
             ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
      } else if(nPerRu_[ru.first-minDccId_][ru.second-1]==0) {//tower dropped!
    fill(meIncompleteFRORateMap_,
             ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
    fill(meDroppedFRORateMap_,
             ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
    fill(meDroppedFROMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
        ++nDroppedFRO_;
    srApplicationErrorLog_ << event.id() << ": Flag of RU "
                   << itSrf->id() << " (DCC " << ru.first
                               << " ch " << ru.second << ") is 'Full readout' "
                   << "while none of its channel was read out\n";
      } else{ //tower partially read out
    fill(meIncompleteFRORateMap_,
             ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
    fill(meDroppedFRORateMap_,
             ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
    fill(meIncompleteFROMap_,
             ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);
        ++nIncompleteFRO_;
    srApplicationErrorLog_ << event.id() << ": Flag of RU"
                   << itSrf->id() << " (DCC " << ru.first
                               << " ch " << ru.second << ") is 'Full readout' "
                   << "while only "
                               << nPerRu_[ru.first-minDccId_][ru.second-1]
                   << " / " << getCrystalCount(ru.first, ru.second)
                               << " channels were read out.\n";
      }
    }

    if(flag == EcalSrFlag::SRF_ZS1 || flag == EcalSrFlag::SRF_ZS2){
      if(nPerRu_[ru.first-minDccId_][ru.second-1]
         ==getCrystalCount(ru.first, ru.second)){
    //ZS readout unit whose every channel was read

        fill(meCompleteZSMap_, ruGraphX(itSrf->id()), ruGraphY(itSrf->id()));
    fill(meCompleteZSRateMap_,
             ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 1);

        ++nCompleteZS_;
      } else{
    fill(meCompleteZSRateMap_,
             ruGraphX(itSrf->id()), ruGraphY(itSrf->id()), 0);
      }
    }
  }
}

int EcalSelectiveReadoutValidation::cIndex2iEta ( int  i ) const

inlineprivate

converse of iEta2cIndex() method.

Definition at line 251 of file EcalSelectiveReadoutValidation.h.

Referenced by dccCh().

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

int EcalSelectiveReadoutValidation::cIndex2iPhi ( int  i ) const

inlineprivate

converse of iPhi2cIndex() method.

Definition at line 258 of file EcalSelectiveReadoutValidation.h.

Referenced by dccCh().

                               {
    return (i+11) % 360;
  }

int EcalSelectiveReadoutValidation::cIndex2iTtEta ( int  i ) const

inlineprivate

converse of iTtEta2cIndex() method.

Definition at line 285 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeTP().

                                {
    return (i<27)?i-28:i-27;
  }

int EcalSelectiveReadoutValidation::cIndex2iTtPhi ( int  i ) const

inlineprivate

converse of iTtPhi2cIndex() method.

Definition at line 291 of file EcalSelectiveReadoutValidation.h.

References book1D(), book2D(), bookFloat(), bookProfile(), bookProfile2D(), dccZsFIR(), getFIRWeights(), dataset::name, pileupCalc::nbins, TSGForRoadSearch_cfi::option, readOutUnitOf(), AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, Phase2TrackerMonitorDigi_cff::ymax, and Phase2TrackerMonitorDigi_cff::ymin.

Referenced by analyzeTP().

                                {
    return i + 1;
  }

int EcalSelectiveReadoutValidation::cIndex2iXY ( int  iX0 ) const

inlineprivate

converse of iXY2cIndex() method.

Definition at line 245 of file EcalSelectiveReadoutValidation.h.

                               {
    return iX0+1;
  }

template<class T >

void EcalSelectiveReadoutValidation::compareSrfColl ( const edm::Event &  event, T &  srfFromData, T &  computedSrf  )

private

Compares two SR flag collection, flags read from data and computed flags. Descripencies are recorded in relevant histogram and log file.

Template Parameters

T	collection type. Must be either an EESrFlagCollection or an EBSrFlagCollection.

Parameters

event	event currently analyzed. Used in logs.
srfFromData	SR flag collection read from data
compareSrf	SR flag collection computed from TTF by SRP emulation

Definition at line 2085 of file EcalSelectiveReadoutValidation.cc.

References dccCh(), fill(), logErrForDccs_, meSRFlagsConsistency_, minDccId_, ruGraphX(), ruGraphY(), and srpAlgoErrorLog_.

Referenced by analyze().

                                                                                                        {
  typedef typename T::const_iterator SrFlagCollectionConstIt;
  typedef typename T::key_type MyRuDetIdType;
  SrFlagCollectionConstIt itSrfFromData = srfFromData.begin();
  SrFlagCollectionConstIt itComputedSr = computedSrf.begin();

  while(itSrfFromData != srfFromData.end()
        || itComputedSr != computedSrf.end()){
    MyRuDetIdType inconsistentRu = 0;
    bool inconsistent = false;
    if(itComputedSr == computedSrf.end() ||
       (itSrfFromData != srfFromData.end()
        && itSrfFromData->id() < itComputedSr->id())){
      //computedSrf is missig a detid found in srfFromData
      pair<int, int> ch = dccCh(itSrfFromData->id());
      srpAlgoErrorLog_ << event.id() << ": " << itSrfFromData->id()
                       << ", DCC " << ch.first << " ch " << ch.second
                       << " found in data (SRF:" << itSrfFromData->flagName()
                       << ") but not in the set of SRFs computed from the data TTF.\n";
      inconsistentRu = itSrfFromData->id();
      inconsistent = true;
      ++itSrfFromData;
    } else if(itSrfFromData==srfFromData.end() ||
              (itComputedSr != computedSrf.end()
               && itComputedSr->id() < itSrfFromData->id())){
      //ebSrFlags is missing a detid found in computedSrf
      pair<int, int> ch = dccCh(itComputedSr->id());
      if(logErrForDccs_[ch.first-minDccId_]){
        srpAlgoErrorLog_ << event.id() << ": " << itComputedSr->id()
                         << ", DCC " << ch.first << " ch " << ch.second
                         << " not found in data. Computed SRF: "
                         << itComputedSr->flagName() << ".\n";
        inconsistentRu = itComputedSr->id();
        inconsistent = true;
      }
      ++itComputedSr;
    } else{
      //*itSrfFromData and *itComputedSr has same detid
      if(itComputedSr->value()!=itSrfFromData->value()){
        pair<int, int> ch = dccCh(itSrfFromData->id());
        srpAlgoErrorLog_ << event.id() << ", "
                         << itSrfFromData->id()
                         << ", DCC " << ch.first << " ch " << ch.second
                         << ", SRF inconsistency: "
                         << "from data: " << itSrfFromData->flagName()
                         << ", computed from TTF: "
                         << itComputedSr->flagName()
                         << "\n";
        inconsistentRu = itComputedSr->id();
        inconsistent = true;
      }
      if(itComputedSr != computedSrf.end()) ++itComputedSr;
      if(itSrfFromData != srfFromData.end()) ++itSrfFromData;
    }

    if(inconsistent) fill(meSRFlagsConsistency_, ruGraphX(inconsistentRu),
                          ruGraphY(inconsistentRu));
  }
}

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

private

Configure DCC ZS FIR weights. Heuristic is used to determine if input weights are normalized weights or integer weights in the hardware representation.

Parameters

weightsForZsFIR

weights from configuration file

Definition at line 2014 of file EcalSelectiveReadoutValidation.cc.

References Exception, firstFIRSample_, firWeights_, getFIRWeights(), mps_fire::i, createfilelist::int, cmsBatch::log, indexGen::s2, and mathSSE::sqrt().

Referenced by EcalSelectiveReadoutValidation().

                                                                                     {
  bool notNormalized  = false;
  bool notInt = false;
  for(unsigned i=0; i < weightsForZsFIR.size(); ++i){
    if(weightsForZsFIR[i] > 1.) notNormalized = true;
    if((int)weightsForZsFIR[i]!=weightsForZsFIR[i]) notInt = true;
  }
  if(notInt && notNormalized){
    throw cms::Exception("InvalidParameter")
      << "weigtsForZsFIR paramater values are not valid: they "
      << "must either be integer and uses the hardware representation "
      << "of the weights or less or equal than 1 and used the normalized "
      << "representation.";
  }
  LogInfo log("DccFir");
  if(notNormalized){
    firWeights_ = vector<int>(weightsForZsFIR.size());
    for(unsigned i = 0; i< weightsForZsFIR.size(); ++i){
      firWeights_[i] = (int)weightsForZsFIR[i];
    }
  } else{
    firWeights_ = getFIRWeights(weightsForZsFIR);
  }

  log << "Input weights for FIR: ";
  for(unsigned i = 0; i < weightsForZsFIR.size(); ++i){
    log << weightsForZsFIR[i] << "\t";
  }

  double s2 = 0.;
  log << "\nActual FIR weights: ";
  for(unsigned i = 0; i < firWeights_.size(); ++i){
    log << firWeights_[i] << "\t";
    s2 += firWeights_[i]*firWeights_[i];
  }

  s2 = sqrt(s2);
  log << "\nNormalized FIR weights after hw representation rounding: ";
  for(unsigned i = 0; i < firWeights_.size(); ++i){
    log << firWeights_[i] / (double)(1<<10) << "\t";
  }

  log <<"\nFirst FIR sample: " << firstFIRSample_;
}

pair< int, int > EcalSelectiveReadoutValidation::dccCh ( const DetId &  xtalId ) const

private

Retrieves the logical number of the DCC reading a given crystal channel.

Parameters

xtarId

crystal channel identifier

Returns

the DCC logical number starting from 1. Retrieves the DCC channel reading out a crystal, the crystals of a barrel trigger tower or the crystals, of an endcap supercrystal.

Parameters

xtarId

crystal channel, barrel trigger tower or endcap supercrystal identifier

Returns

pair of (DCC ID, DCC channel)

Definition at line 1631 of file EcalSelectiveReadoutValidation.cc.

References cIndex2iEta(), cIndex2iPhi(), EcalElectronicsId::dccId(), TauDecayModes::dec, DetId::det(), DetId::Ecal, EcalBarrel, EcalEndcap, EcalTriggerTower, elecMap_, Exception, EcalElectronicsMapping::getDCCandSC(), EcalElectronicsMapping::getElectronicsId(), EcalTrigTowerDetId::ieta(), EcalTrigTowerDetId::iphi(), iTtEta2cIndex(), iTtPhi2cIndex(), maxDccId_, minDccId_, nOneEeTtEta, DetId::rawId(), mps_fire::result, DetId::subdetId(), EcalElectronicsId::towerId(), and groupFilesInBlocks::tt.

Referenced by anaDigi(), analyzeEB(), analyzeEE(), checkSrApplication(), compareSrfColl(), dccId(), and getDccSrDependentPayload().

                                                                           {
  if(detId.det()!=DetId::Ecal){
    throw cms::Exception("InvalidParameter")
      << "Wrong type of DetId passed to the "
      "EcalSelectiveReadoutValidation::dccCh(const DetId&). "
      "An ECAL DetId was expected.\n";
  }
  
  DetId xtalId;
  switch(detId.subdetId()){
  case EcalTriggerTower:       //Trigger tower
    {
      const EcalTrigTowerDetId tt = detId;
      //pick up one crystal of the trigger tower: they are however all readout by
      //the same DCC channel in the barrel.
      //Arithmetic is easier on the "c" indices:
      const int iTtPhi0 = iTtPhi2cIndex(tt.iphi());
      const int iTtEta0 = iTtEta2cIndex(tt.ieta());
      const int oneXtalPhi0 = iTtPhi0 * 5;
      const int oneXtalEta0 = (iTtEta0 - nOneEeTtEta) * 5;

      xtalId = EBDetId(cIndex2iEta(oneXtalEta0),
               cIndex2iPhi(oneXtalPhi0));
    }
    break;
  case EcalEndcap:
    if(detId.rawId() & 0x8000){ //Supercrystal
      return elecMap_->getDCCandSC(EcalScDetId(detId));
    } else {                    //EE crystal
      xtalId = detId;
    }
    break;
  case EcalBarrel:              //EB crystal
    xtalId = detId;
    break;
  default:
    throw cms::Exception("InvalidParameter")
      << "Wrong type of DetId passed to the method "
      "EcalSelectiveReadoutValidation::dccCh(const DetId&). "
      "A valid EcalTriggerTower, EcalBarrel or EcalEndcap DetId was expected. "
      "detid = " << xtalId.rawId() << ".\n";
  }

  const EcalElectronicsId& EcalElecId = elecMap_->getElectronicsId(xtalId);

  pair<int,int> result;
  result.first = EcalElecId.dccId();

  if(result.first < minDccId_ || result.second > maxDccId_){
    throw cms::Exception("OutOfRange")
      << "Got an invalid DCC ID, DCCID = " << result.first
      << " for DetId 0x" << hex << detId.rawId()
      << " and 0x" << xtalId.rawId() << dec << "\n";
  }

  result.second = EcalElecId.towerId();

  if(result.second < 1 || result.second > 68){
    throw cms::Exception("OutOfRange")
      << "Got an invalid DCC channel ID, DCC_CH = " << result.second
      << " for DetId 0x" << hex << detId.rawId()
      << " and 0x"  << xtalId.rawId() << dec << "\n";
  }

  return result;
}

int EcalSelectiveReadoutValidation::dccId ( const EcalScDetId &  detId ) const

private

Retrieves the ID of the DCC reading a readout unit

Parameters

detId

detid of the readout unit

Definition at line 2146 of file EcalSelectiveReadoutValidation.cc.

References elecMap_, and EcalElectronicsMapping::getDCCandSC().

Referenced by selectFedsForLog(), and xtalGraphY().

                                                                       {
  return elecMap_->getDCCandSC(detId).first;
}

int EcalSelectiveReadoutValidation::dccId ( const EcalTrigTowerDetId &  detId ) const

private

Definition at line 2150 of file EcalSelectiveReadoutValidation.cc.

References dccCh(), Exception, and EcalTrigTowerDetId::ietaAbs().

                                                                              {
  if(detId.ietaAbs()>17){
    throw cms::Exception("InvalidArgument")
      << "Argument of EcalSelectiveReadoutValidation::dccId(const EcalTrigTowerDetId&) "
       << "must be a barrel trigger tower Id\n";
  }
  return dccCh(detId).first;
}

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

staticprivate

Emulates the DCC zero suppression FIR filter. If one of the time sample is not in gain 12, numeric_limits<int>::max() is returned.

Parameters

frame	data frame
firWeights	TAP weights
firstFIRSample	index (starting from 1) of the first time sample to be used in the filter
saturated	if not null, *saturated is set to true if all the time sample are not in gain 12 and set to false otherwise.

Returns

FIR output or numeric_limits<int>::max().

Definition at line 1950 of file EcalSelectiveReadoutValidation.cc.

References EcalMGPASample::adc(), EcalMGPASample::gainId(), LogTrace, SiStripPI::max, simplePhotonAnalyzer_cfi::sample, EcalDataFrame::size(), and w.

Referenced by analyzeEB(), analyzeEE(), and cIndex2iTtPhi().

                                     {
  const int nFIRTaps = 6;
  //FIR filter weights:
  const vector<int>& w = firWeights;

  //accumulator used to compute weighted sum of samples
  int acc = 0;
  bool gain12saturated = false;
  const int gain12 = 0x01;
  const int lastFIRSample = firstFIRSample + nFIRTaps - 1;
  //LogDebug("DccFir") << "DCC FIR operation: ";
  int iWeight = 0;
  for(int iSample=firstFIRSample-1;
      iSample<lastFIRSample; ++iSample, ++iWeight){
    if(iSample>=0 && iSample < frame.size()){
      EcalMGPASample sample(frame[iSample]);
      if(sample.gainId()!=gain12) gain12saturated = true;
      LogTrace("DccFir") << (iSample>=firstFIRSample?"+":"") << sample.adc()
             << "*(" << w[iWeight] << ")";
      acc+=sample.adc()*w[iWeight];
    } else{
      edm::LogWarning("DccFir") << __FILE__ << ":" << __LINE__ <<
    ": Not enough samples in data frame or 'ecalDccZs1stSample' module "
    "parameter is not valid...";
    }
  }
  LogTrace("DccFir") << "\n";
  //discards the 8 LSBs
  //(shift operator cannot be used on negative numbers because
  // the result depends on compilator implementation)
  acc = (acc>=0)?(acc >> 8):-(-acc >> 8);
  //ZS passed if weighted sum acc above ZS threshold or if
  //one sample has a lower gain than gain 12 (that is gain 12 output
  //is saturated)

  LogTrace("DccFir") << "acc: " << acc << "\n"
             << "saturated: " << (gain12saturated?"yes":"no") << "\n";

  if(saturated){
    *saturated = gain12saturated;
  }

  return gain12saturated?numeric_limits<int>::max():acc;
}

void EcalSelectiveReadoutValidation::dqmBeginRun ( edm::Run const &  r, edm::EventSetup const &  es  )

override

Definition at line 833 of file EcalSelectiveReadoutValidation.cc.

References elecMap_, edm::EventSetup::get(), initAsciiFile(), edm::ESHandle< T >::product(), and triggerTowerMap_.

                                                                                        {
  // endcap mapping
  edm::ESHandle<EcalTrigTowerConstituentsMap> hTriggerTowerMap;
  es.get<IdealGeometryRecord>().get(hTriggerTowerMap);
  triggerTowerMap_ = hTriggerTowerMap.product();

  //electronics map
  ESHandle< EcalElectronicsMapping > ecalmapping;
  es.get< EcalMappingRcd >().get(ecalmapping);
  elecMap_ = ecalmapping.product();

  initAsciiFile();
}

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

overrideprotected

Definition at line 847 of file EcalSelectiveReadoutValidation.cc.

References MonitorElement::Fill(), getL1aRate(), and meL1aRate_.

                                                                                   {
  meL1aRate_->Fill(getL1aRate());
}

void EcalSelectiveReadoutValidation::fill ( MonitorElement *  me, float  x  )

inlineprivate

Wrapper to fill methods of DQM monitor elements.

Definition at line 358 of file EcalSelectiveReadoutValidation.h.

References MonitorElement::Fill().

Referenced by analyze(), analyzeDataVolume(), analyzeEB(), analyzeEE(), analyzeTP(), checkSrApplication(), and compareSrfColl().

                                        {
    if(me) me->Fill(x);
  }

void EcalSelectiveReadoutValidation::fill ( MonitorElement *  me, float  x, float  yw  )

inlineprivate

Definition at line 361 of file EcalSelectiveReadoutValidation.h.

References MonitorElement::Fill().

                                                 {
    if(me) me->Fill(x, yw);
  }

void EcalSelectiveReadoutValidation::fill ( MonitorElement *  me, float  x, float  y, float  zw  )

inlineprivate

Definition at line 364 of file EcalSelectiveReadoutValidation.h.

References MonitorElement::Fill().

                                                          {
    if(me) me->Fill(x, y, zw);
  }

void EcalSelectiveReadoutValidation::fill ( MonitorElement *  me, float  x, float  y, float  z, float  w  )

inlineprivate

Definition at line 367 of file EcalSelectiveReadoutValidation.h.

References MonitorElement::Fill(), getL1aRate(), initAsciiFile(), and updateL1aRate().

                                                                  {
    if(me) me->Fill(x, y, z, w);
  }

double EcalSelectiveReadoutValidation::frame2Energy ( const EcalDataFrame &  frame ) const

private

Selective Readout decisions Validation

Parameters

event	EDM event
es	event setup Energy reconstruction from ADC samples.
frame	the ADC sample of an ECA channel

Definition at line 1593 of file EcalSelectiveReadoutValidation.cc.

References mps_fire::i, createfilelist::int, min(), gen::n, EcalDataFrame::size(), and weights_.

Referenced by analyzeEB(), and analyzeEE().

                                                                                   {
  static std::atomic<bool> firstCall {true};
  bool expected = true;
  if( firstCall.compare_exchange_strong(expected, false) ) {
    stringstream buf;
    buf << "Weights:";
    for(unsigned i=0; i<weights_.size();++i){
      buf << "\t" << weights_[i];
    }
    edm::LogInfo("EcalSrValid") << buf.str() << "\n";
    firstCall = false;
  }
  double adc2GeV = 0.;

  if(typeid(EBDataFrame)==typeid(frame)){//barrel APD
    adc2GeV = .035;
  } else if(typeid(EEDataFrame)==typeid(frame)){//endcap VPT
    adc2GeV = 0.06;
  } else{
    assert(false);
  }

  double acc = 0;

  const int n = min(frame.size(), (int)weights_.size());

  double gainInv[] = {12., 1., 6., 12.};

  for(int i=0; i < n; ++i){
    acc += weights_[i]*frame[i].adc()*gainInv[frame[i].gainId()]*adc2GeV;
  }
  return acc;
}

template<class T >

double EcalSelectiveReadoutValidation::frame2EnergyForTp ( const T &  frame, int  offset = 0  ) const

private

Energy reconstruction from ADC samples to be used for trigger primitive estimate.

Parameters

frame	the ADC sample of an ECA channel
offset	time offset. To be used to evaluate energy of the event previous (offset=-1) and next (offset=+1) to the triggered one.

Definition at line 1795 of file EcalSelectiveReadoutValidation.cc.

References Exception, mps_fire::i, gen::n, PFRecoTauDiscriminationByIsolation_cfi::offset, and HGCalRecHit_cfi::weights.

Referenced by setTtEtSums().

                                              {
  //we have to start by 0 in order to handle offset=-1
  //(however Fenix FIR has AFAK only 5 taps)
  double weights[] = {0., -1/3., -1/3., -1/3., 0., 1.};

  double adc2GeV = 0.;
  if(typeid(frame) == typeid(EBDataFrame)){
    adc2GeV = 0.035;
  } else if(typeid(frame) == typeid(EEDataFrame)){
    adc2GeV = 0.060;
  } else{ //T is an invalid type!
    //TODO: replace message by a cms exception
    throw cms::Exception("Severe Error")
      << __FILE__ << ":" << __LINE__ << ": "
      << "this is a bug. Please report it.\n";
  }

  double acc = 0;

  const int n = min<int>(frame.size(), sizeof(weights)/sizeof(weights[0]));

  double gainInv[] = {12., 1., 6., 12};

  for(int i=offset; i < n; ++i){
    int iframe = i + offset;
    if(iframe>=0 && iframe<frame.size()){
      acc += weights[i]*frame[iframe].adc()
    *gainInv[frame[iframe].gainId()]*adc2GeV;
    }
  }
  //cout << "\n";
  return acc;
}

double EcalSelectiveReadoutValidation::getBytesPerCrystal ( ) const

inlineprivate

Gets the number of bytes per crystal channel of the event part depending on the number of read crystal channels.

Returns

the number of bytes.

Definition at line 145 of file EcalSelectiveReadoutValidation.h.

Referenced by analyzeDataVolume(), getDccSrDependentPayload(), getEbEventSize(), and getEeEventSize().

                                   {
    return 3*8;
  }

int EcalSelectiveReadoutValidation::getCrystalCount ( int  iDcc, int  iDccCh  )

private

Retrieves number of crystal channel read out by a DCC channel

Parameters

iDcc	DCC ID starting from 1
iDccCh	DCC channel starting from 1

Returns

crystal count

Definition at line 2260 of file EcalSelectiveReadoutValidation.cc.

References maxDccId_.

Referenced by analyzeDataVolume(), checkSrApplication(), and xtalGraphY().

                                                                       {
  if(iDcc < minDccId_ || iDcc > maxDccId_){ //invalid DCC
    return 0;
  } else if (10 <= iDcc && iDcc <= 45) {//EB
    return 25;
  } else { //EE
    int iDccPhi;
    if(iDcc < 10) iDccPhi = iDcc;
    else iDccPhi = iDcc - 45;
    switch(iDccPhi*100+iDccCh){
    case 110:
    case 232:
    case 312:
    case 412:
    case 532:
    case 610:
    case 830:
    case 806:
      //inner partials at 12, 3, and 9 o'clock
      return 20;
    case 134:
    case 634:
    case 827:
    case 803:
      return 10;
    case 330:
    case 430:
      return 20;
    case 203:
    case 503:
    case 721:
    case 921:
      return 21;
    default:
      return 25;
    }
  }
}

double EcalSelectiveReadoutValidation::getDccEventSize ( int  iDcc0, double  nReadXtals  ) const

inlineprivate

Gets the size of an DCC event fragment.

Parameters

iDcc0	the DCC logical number starting from 0.
nReadXtals	number of read crystal channels.

Returns

the DCC event fragment size in bytes.

Definition at line 154 of file EcalSelectiveReadoutValidation.h.

References EB, EE, getDccOverhead(), getDccSrDependentPayload(), and getRuCount().

Referenced by analyzeDataVolume().

                                                            {
    subdet_t subdet;
     if(iDcc0<9 || iDcc0>=45){
       subdet = EE;
     } else{
       subdet = EB;
     }
     //     return getDccOverhead(subdet)+nReadXtals*getBytesPerCrystal()
     //       + getRuCount(iDcc0)*8;
    return getDccOverhead(subdet)
      + getDccSrDependentPayload(iDcc0, getRuCount(iDcc0), nReadXtals);
  }

double EcalSelectiveReadoutValidation::getDccOverhead ( subdet_t  subdet ) const

inlineprivate

Gets the size in bytes fixed-size part of a DCC event fragment.

Returns

the fixed size in bytes.

Definition at line 136 of file EcalSelectiveReadoutValidation.h.

References EB.

Referenced by getDccEventSize(), getEbEventSize(), and getEeEventSize().

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

double EcalSelectiveReadoutValidation::getDccSrDependentPayload ( int  iDcc0, double  nReadRus, double  nReadXtals  ) const

inlineprivate

Gets DCC event fragment payload depending on the channel selection made by the selective readout.

Parameters

iDcc0	the DCC logical number starting from 0.
nReadRus	number of read-out RUs
nReadXtals	number of read-out crystal channels.

Returns

the DCC event fragment payload in bytes.

Definition at line 174 of file EcalSelectiveReadoutValidation.h.

References dccCh(), getBytesPerCrystal(), getRuCount(), readAllCollections(), and setTtEtSums().

Referenced by analyzeDataVolume(), and getDccEventSize().

                                                          {
    return nReadXtals*getBytesPerCrystal() + nReadRus*8;
  }

double EcalSelectiveReadoutValidation::getEbEventSize ( double  nReadXtals ) const

private

Computes the size of an ECAL barrel event fragment.

Parameters

nReadXtals

number of read crystal channels

Returns

the event fragment size in bytes

Definition at line 1924 of file EcalSelectiveReadoutValidation.cc.

References EB, getBytesPerCrystal(), getDccOverhead(), getRuCount(), nEbDccs, and nEeDccs.

Referenced by analyzeDataVolume().

                                                                            {
  double ruHeaderPayload = 0.;
  const int firstEbDcc0 = nEeDccs/2;
  for(int iDcc0 = firstEbDcc0; iDcc0 < firstEbDcc0 + nEbDccs; ++iDcc0){
    ruHeaderPayload += getRuCount(iDcc0)*8.;
  }

  return getDccOverhead(EB)*nEbDccs + nReadXtals*getBytesPerCrystal()
    + ruHeaderPayload;
}