---

CSCConditions.h

Go to the documentation of this file.

#ifndef CSCObjects_CSCConditions_h
#define CSCObjects_CSCConditions_h

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESWatcher.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "CondFormats/CSCObjects/interface/CSCDBNoiseMatrix.h"
#include "CondFormats/DataRecord/interface/CSCDBGainsRcd.h"
#include "CondFormats/DataRecord/interface/CSCBadStripsRcd.h"
#include "CondFormats/DataRecord/interface/CSCBadWiresRcd.h"
#include "CondFormats/DataRecord/interface/CSCBadChambersRcd.h"
#include <vector>
#include <bitset>

class CSCDBGains;
class CSCDBPedestals;
class CSCDBCrosstalk;
class CSCBadStrips;
class CSCBadWires;
class CSCBadChambers;

class CSCConditions
{
public:
  CSCConditions( const edm::ParameterSet& ps );
  ~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;

  const std::bitset<80>& badStripWord( const CSCDetId& id ) const;
  const std::bitset<112>& badWireWord( const CSCDetId& id ) const;

  bool isInBadChamber( const CSCDetId& id ) const;

  void print() const;

  bool readBadChannels() const { return readBadChannels_; }

  bool readBadChambers() const { return readBadChambers_; }

  void fillBadStripWords();
  void fillBadWireWords();

  float averageGain() const;

private:

  const CSCDBNoiseMatrix * theNoiseMatrix;
  const CSCDBGains * theGains;
  const CSCDBPedestals * thePedestals;
  const CSCDBCrosstalk * theCrosstalk;

  const CSCBadStrips* theBadStrips;
  const CSCBadWires* theBadWires;
  const CSCBadChambers* theBadChambers;
  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

  // cache bad channel words once created
  std::vector< std::bitset<80> > badStripWords;
  std::vector< std::bitset<112> > badWireWords;

  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_; 
};

#endif

---