CSCEventData.h

Go to the documentation of this file.

#ifndef EventFilter_CSCRawToDigi_CSCEventData_h
#define EventFilter_CSCRawToDigi_CSCEventData_h

class CSCCFEBData;
class CSCTMBHeader;
class CSCComparatorData;
class CSCWireDigi;
class CSCStripDigi;
#include <map>
#include <vector>
#ifndef LOCAL_UNPACK
#include <atomic>
#endif
#include "EventFilter/CSCRawToDigi/interface/CSCDMBHeader.h"
#include "EventFilter/CSCRawToDigi/interface/CSCDMBTrailer.h"
#include "EventFilter/CSCRawToDigi/interface/CSCALCTHeader.h"
#include "EventFilter/CSCRawToDigi/interface/CSCAnodeData.h"
#include "EventFilter/CSCRawToDigi/interface/CSCALCTTrailer.h"
#include "EventFilter/CSCRawToDigi/interface/CSCTMBData.h"
#include "EventFilter/CSCRawToDigi/interface/CSCDMBTrailer.h"
#include "DataFormats/CSCDigi/interface/CSCRPCDigi.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/MuonDetId/interface/GEMDetId.h"
#include <boost/dynamic_bitset.hpp>

class CSCEventData {
public:
  explicit CSCEventData(int chamberType, uint16_t format_version = 2005);
  CSCEventData(const uint16_t *buf, uint16_t format_version = 2005);
  CSCEventData() {}
  ~CSCEventData();
  CSCEventData(const CSCEventData &data);
  CSCEventData operator=(const CSCEventData &data);

  unsigned short size() const { return size_; }

  static void setDebug(const bool value) { debug = value; }

  bool isALCT(const uint16_t *buf);

  bool isTMB(const uint16_t *buf);

  const CSCCFEBData *cfebData(unsigned icfeb) const;

  std::vector<CSCStripDigi> stripDigis(const CSCDetId &idlayer) const;

  std::vector<CSCStripDigi> stripDigis(unsigned idlayer, unsigned icfeb) const;

  std::vector<std::vector<CSCStripDigi> > stripDigis() const;

  std::vector<CSCWireDigi> wireDigis(unsigned ilayer) const;
  std::vector<std::vector<CSCWireDigi> > wireDigis() const;

  int nalct() const { return theDMBHeader.nalct(); }

  int nclct() const { return theDMBHeader.nclct(); }

  const CSCDMBHeader *dmbHeader() const { return &theDMBHeader; }
  CSCDMBHeader *dmbHeader() { return &theDMBHeader; }

  CSCALCTHeader *alctHeader() const;

  CSCALCTTrailer *alctTrailer() const;

  CSCAnodeData *alctData() const;

  CSCTMBData *tmbData() const;

  CSCTMBHeader *tmbHeader() const;

  CSCComparatorData *comparatorData() const;

  const CSCDMBTrailer *dmbTrailer() const { return &theDMBTrailer; }
  void add(const CSCStripDigi &, int layer);
  void add(const CSCWireDigi &, int layer);
  void add(const CSCComparatorDigi &, int layer);
  void add(const CSCComparatorDigi &, const CSCDetId &);
  void add(const std::vector<CSCALCTDigi> &);
  void add(const std::vector<CSCCLCTDigi> &);
  void add(const std::vector<CSCCorrelatedLCTDigi> &);
  void addShower(const std::vector<CSCShowerDigi> &);           // add LCT CSCShower
  void addAnodeShower(const std::vector<CSCShowerDigi> &);      // add Anode CSCShower
  void addCathodeShower(const std::vector<CSCShowerDigi> &);    // add Cathode CSCShower
  void addAnodeALCTShower(const std::vector<CSCShowerDigi> &);  // add Anode ALCT CSCShower
  void add(const std::vector<GEMPadDigiCluster> &, const GEMDetId &);

  void setEventInformation(int bxnum, int lvl1num);

  boost::dynamic_bitset<> pack();

  void addALCTStructures();

  int chamberType() const { return theChamberType; }

  uint16_t getFormatVersion() const { return theFormatVersion; }

  unsigned int calcALCTcrc(std::vector<std::pair<unsigned int, unsigned short *> > &vec);

#ifdef LOCAL_UNPACK
  static bool debug;
#else
  static std::atomic<bool> debug;
#endif
  //uint16_t dataPresent; // 7 bit word which will tell if alct, clct, and 5 cfebs are present
  static void selfTest();

private:
  void init();
  void unpack_data(const uint16_t *buf);
  void copy(const CSCEventData &);
  void destroy();

  void checkALCTClasses();
  void checkTMBClasses();

  void addComparatorInformation(std::vector<CSCStripDigi> &, int layer) const;

  CSCDMBHeader theDMBHeader;
  //these are empty data objects unless filled in CSCEventData.cc
  CSCALCTHeader *theALCTHeader;
  CSCAnodeData *theAnodeData;
  CSCALCTTrailer *theALCTTrailer;
  CSCTMBData *theTMBData;

  CSCCFEBData *theCFEBData[CSCConstants::MAX_CFEBS_RUN2];

  CSCDMBTrailer theDMBTrailer;

  int size_;
  int theChamberType;

  unsigned short *alctZSErecovered;
  int zseEnable;

  uint16_t theFormatVersion;
};

std::ostream &operator<<(std::ostream &os, const CSCEventData &evt);
#endif