CSCConditions.h

Go to the documentation of this file.

#ifndef CSCCalibration_CSCConditions_h
#define CSCCalibration_CSCConditions_h

#include "CondFormats/DataRecord/interface/CSCChamberTimeCorrectionsRcd.h"
#include "CondFormats/DataRecord/interface/CSCDBChipSpeedCorrectionRcd.h"
#include "CondFormats/DataRecord/interface/CSCDBCrosstalkRcd.h"
#include "CondFormats/DataRecord/interface/CSCDBGasGainCorrectionRcd.h"
#include "CondFormats/DataRecord/interface/CSCDBNoiseMatrixRcd.h"
#include "CondFormats/DataRecord/interface/CSCDBPedestalsRcd.h"
#include "CondFormats/CSCObjects/interface/CSCDBNoiseMatrix.h"
#include "CondFormats/DataRecord/interface/CSCBadStripsRcd.h"
#include "CondFormats/DataRecord/interface/CSCBadWiresRcd.h"
#include "CondFormats/DataRecord/interface/CSCDBGainsRcd.h"
#include "CalibMuon/CSCCalibration/interface/CSCChannelMapperRecord.h"
#include "CalibMuon/CSCCalibration/interface/CSCIndexerRecord.h"
#include "CondFormats/DataRecord/interface/CSCBadChambersRcd.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESWatcher.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Utilities/interface/ESGetToken.h"
#include <bitset>
#include <vector>

class CSCDBGains;
class CSCDBPedestals;
class CSCDBCrosstalk;
class CSCBadStrips;
class CSCBadWires;
class CSCBadChambers;
class CSCDBChipSpeedCorrection;
class CSCChamberTimeCorrections;
class CSCDBGasGainCorrection;
class CSCIndexerBase;
class CSCChannelMapperBase;

class CSCConditions {
public:
  explicit CSCConditions(const edm::ParameterSet &ps, edm::ConsumesCollector);
  ~CSCConditions();

  void initializeEvent(const edm::EventSetup &es);

  float gain(const CSCDetId &detId, int channel) const;
  float gainSigma(const CSCDetId &detId, int channel) const { return 0.005; }

  float pedestal(const CSCDetId &detId, int channel) const;
  float pedestalSigma(const CSCDetId &detId, int channel) const;

  float crosstalkSlope(const CSCDetId &detId, int channel, bool leftRight) const;
  float crosstalkIntercept(const CSCDetId &detId, int channel, bool leftRight) const;

  const CSCDBNoiseMatrix::Item &noiseMatrix(const CSCDetId &detId, int channel) const;

  void noiseMatrixElements(const CSCDetId &id, int channel, std::vector<float> &me) const;

  void crossTalk(const CSCDetId &id, int channel, std::vector<float> &ct) const;

  float chipCorrection(const CSCDetId &detId, int channel) const;

  float chamberTimingCorrection(const CSCDetId &detId) const;

  float anodeBXoffset(const CSCDetId &detId) const;

  const std::bitset<112> &badStripWord() const { return badStripWord_; }

  const std::bitset<112> &badWireWord() const { return badWireWord_; }

  const CSCDetId &idOfBadChannelWords() const { return idOfBadChannelWords_; }

  void print() const;

  bool isInBadChamber(const CSCDetId &id) const;

  bool readBadChannels() const { return readBadChannels_; }

  bool readBadChambers() const { return readBadChambers_; }

  bool useTimingCorrections() const { return useTimingCorrections_; }

  void fillBadChannelWords(const CSCDetId &id);

  float averageGain() const;

  float gasGainCorrection(const CSCDetId &detId, int strip, int wire) const;

  bool useGasGainCorrections() const { return useGasGainCorrections_; }

  int channelFromStrip(const CSCDetId &id, int geomStrip) const;
  int rawStripChannel(const CSCDetId &id, int geomChannel) const;

private:
  void fillBadStripWord(const CSCDetId &id);
  void fillBadWireWord(const CSCDetId &id);
  void setIdOfBadChannelWords(const CSCDetId &id) { idOfBadChannelWords_ = id; }

  // handles to conditions data

  edm::ESHandle<CSCDBGains> theGains;
  edm::ESHandle<CSCDBCrosstalk> theCrosstalk;
  edm::ESHandle<CSCDBPedestals> thePedestals;
  edm::ESHandle<CSCDBNoiseMatrix> theNoiseMatrix;
  edm::ESHandle<CSCBadStrips> theBadStrips;
  edm::ESHandle<CSCBadWires> theBadWires;
  edm::ESHandle<CSCBadChambers> theBadChambers;
  edm::ESHandle<CSCDBChipSpeedCorrection> theChipCorrections;
  edm::ESHandle<CSCChamberTimeCorrections> theChamberTimingCorrections;
  edm::ESHandle<CSCDBGasGainCorrection> theGasGainCorrections;

  // handles to algorithm versions

  edm::ESHandle<CSCIndexerBase> indexer_;
  edm::ESHandle<CSCChannelMapperBase> mapper_;

  //EventSetup Tokens for Handles
  edm::ESGetToken<CSCDBGains, CSCDBGainsRcd> gainsToken_;
  edm::ESGetToken<CSCDBCrosstalk, CSCDBCrosstalkRcd> crosstalkToken_;
  edm::ESGetToken<CSCDBPedestals, CSCDBPedestalsRcd> pedestalsToken_;
  edm::ESGetToken<CSCDBNoiseMatrix, CSCDBNoiseMatrixRcd> noiseMatrixToken_;
  edm::ESGetToken<CSCBadStrips, CSCBadStripsRcd> badStripsToken_;
  edm::ESGetToken<CSCBadWires, CSCBadWiresRcd> badWiresToken_;
  edm::ESGetToken<CSCBadChambers, CSCBadChambersRcd> badChambersToken_;
  edm::ESGetToken<CSCDBChipSpeedCorrection, CSCDBChipSpeedCorrectionRcd> chipCorrectionsToken_;
  edm::ESGetToken<CSCChamberTimeCorrections, CSCChamberTimeCorrectionsRcd> chamberTimingCorrectionsToken_;
  edm::ESGetToken<CSCDBGasGainCorrection, CSCDBGasGainCorrectionRcd> gasGainCorrectionsToken_;
  edm::ESGetToken<CSCIndexerBase, CSCIndexerRecord> indexerToken_;
  edm::ESGetToken<CSCChannelMapperBase, CSCChannelMapperRecord> mapperToken_;

  // logical flags controlling some conditions data usage

  bool readBadChannels_;        // flag whether or not to even attempt reading bad
                                // channel info from db
  bool readBadChambers_;        // flag whether or not to even attempt reading bad
                                // chamber info from db
  bool useTimingCorrections_;   // flag whether or not to even attempt reading
                                // timing correction info from db
  bool useGasGainCorrections_;  // flag whether or not to even attempt reading
                                // gas-gain correction info from db

  // Cache bad channel content for current CSC layer
  CSCDetId idOfBadChannelWords_;
  std::bitset<112> badStripWord_;
  std::bitset<112> badWireWord_;

  mutable float theAverageGain;  // average over entire system, subject to some
                                 // constraints!

  edm::ESWatcher<CSCDBGainsRcd> gainsWatcher_;
  //@@ remove until we have real information to use
  //  edm::ESWatcher<CSCBadStripsRcd> badStripsWatcher_;
  //  edm::ESWatcher<CSCBadWiresRcd> badWiresWatcher_;

  // Total number of CSC layers in the system, with full ME42 installed.
  enum elayers { MAX_LAYERS = 3240 };
};

#endif