CSCDQM_Summary.h

Go to the documentation of this file.

/*
 * =====================================================================================
 *
 *       Filename:  Summary.h
 *
 *    Description:  CSC summary map and appropriate functions.
 *
 *        Version:  1.0
 *        Created:  05/19/2008 10:52:21 AM
 *       Revision:  none
 *       Compiler:  gcc
 *
 *         Author:  Valdas Rapsevicius (VR), Valdas.Rapsevicius@cern.ch
 *        Company:  CERN, CH
 *
 * =====================================================================================
 */

#ifndef CSCDQM_Summary_H
#define CSCDQM_Summary_H

#include <TH2.h>
#include <cmath>
#include <vector>
#include <bitset>
#include <iostream>

#include "CSCDQM_Detector.h"
#include "CSCDQM_Logger.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"

#define HWSTATUSBITSETSIZE    14
#define HWSTATUSERRORBITS     0xffc
#define HWSTATUSEQUALS(s, m)  (((std::bitset<HWSTATUSBITSETSIZE>) m & s) == m)
#define HWSTATUSANY(s, m)     (((std::bitset<HWSTATUSBITSETSIZE>) m & s).any())
#define HWSTATUSANYERROR(s)   (HWSTATUSANY(s, HWSTATUSERRORBITS))

#define NTICS                 100

namespace cscdqm {

enum HWStatusBit {

  DATA,         
  MASKED,       

  HOT,          
  COLD,         

  FORMAT_ERR,   
  L1SYNC_ERR,   
  FIFOFULL_ERR, 
  INPUTTO_ERR,  

  NODATA_ALCT,  
  NODATA_CLCT,  
  NODATA_CFEB,  
  CFEB_BWORDS,  

  STANDBY,      
  WAS_ON        

};

typedef std::bitset<HWSTATUSBITSETSIZE> HWStatusBitSet;

class Summary {

  public:

    Summary();
    ~Summary();

    void Reset();

    const Detector getDetector() const { return detector; }

    void ReadReportingChambers(const TH2*& h2, const double threshold = 1.0);
    void ReadReportingChambersRef(const TH2*& h2, const TH2*& refh2, const double cold_coef = 0.1, const double cold_Sfail = 5.0, const double hot_coef = 2.0, const double hot_Sfail = 5.0);
    void ReadErrorChambers(const TH2*& evs, const TH2*& err, const HWStatusBit bit, const double eps_max = 0.1, const double Sfail = 5.0);

    const unsigned int setMaskedHWElements(std::vector<std::string>& tokens);

    void Write(TH2*& h2, const unsigned int station) const;
    void WriteMap(TH2*& h2);
    void WriteChamberState(TH2*& h2, const int mask, const int value = 1, const bool reset = true, const bool op_any = false) const;

    void ReSetValue(const HWStatusBit bit);
    void ReSetValue(const Address& adr, const HWStatusBit bit);
    void SetValue(const HWStatusBit bit, const int value = 1);
    void SetValue(Address adr, const HWStatusBit bit, const int value = 1);

    const HWStatusBitSet GetValue(Address adr) const;
    const int IsPhysicsReady(const unsigned int px, const unsigned int py);
    //const int IsPhysicsReady(const float xmin, const float xmax, const float ymin, const float ymax) const;

    const double GetEfficiencyHW() const;
    const double GetEfficiencyHW(const unsigned int station) const;
    const double GetEfficiencyHW(Address adr) const; 
    const double GetEfficiencyArea(const unsigned int station) const; 
    const double GetEfficiencyArea(const Address& adr) const; 

    bool isChamberStandby(unsigned int side, unsigned int station, unsigned int ring, unsigned int chamber) const;
    bool isChamberStandby(CSCDetId cid) const;

  private:

    const bool ChamberCoordsToAddress(const unsigned int x, const unsigned int y, Address& adr) const;
    const bool ChamberAddressToCoords(const Address& adr, unsigned int& x, unsigned int& y) const;
    const double GetReportingArea(Address adr) const; 

    HWStatusBitSet map[N_SIDES][N_STATIONS][N_RINGS][N_CHAMBERS][N_LAYERS][N_CFEBS][N_HVS];

    Detector detector;

};

}

#endif