EcalTrigPrimCompactColl.h

Go to the documentation of this file.

#ifndef ECALTRIGPRIMCOMPACTCOLL_H
#define ECALTRIGPRIMCOMPACTCOLL_H

#include <vector>
#include <inttypes.h>
#include "FWCore/Utilities/interface/Exception.h"

#include "DataFormats/EcalDigi/interface/EcalTriggerPrimitiveDigi.h"
#include "DataFormats/Common/interface/SortedCollection.h"
typedef edm::SortedCollection<EcalTriggerPrimitiveDigi> EcalTrigPrimDigiCollection;

class EcalTrigPrimCompactColl {

private:
  static const int nPhiBins = 72;
  static const int nEtaBins = 56;
  static const int nBins = nPhiBins*nEtaBins;

private:
  static size_t index(int ieta, int iphi){
    size_t r = unsigned(((ieta<0) ? (ieta+28) : (ieta+27))*nPhiBins + (iphi -1));
    if(r >= (unsigned)nBins) throw cms::Exception("Invalid argument") << "Trigger tower index (" << ieta << "," << iphi << ") are out of range";
    return r;
  }
  
public:
  EcalTrigPrimCompactColl(): formatVers_(0), data_(nBins){};
  
  void setValue(int ieta, int iphi, uint16_t sample){ data_[index(ieta, iphi)] = sample;}

  uint16_t raw(int ieta, int iphi) const { return data_[index(ieta, iphi)]; }
  uint16_t raw(const EcalTrigTowerDetId& ttId) const { return raw(ttId.ieta(), ttId.iphi()); }

  int compressedEt(int ieta, int iphi) const { return raw(ieta, iphi)&0xFF; }
  int compressedEt(const EcalTrigTowerDetId& ttId) const { return compressedEt(ttId.ieta(), ttId.iphi()); }

  bool fineGrain(int ieta, int iphi) const { return (raw(ieta, iphi)&0x100)!=0; }
  bool fineGrain(const EcalTrigTowerDetId& ttId) const { return fineGrain(ttId.ieta(), ttId.iphi()); }

  int ttFlag(int ieta, int iphi) const { return (raw(ieta, iphi)>>9)&0x7; }
  int ttFlag(const EcalTrigTowerDetId& ttId) const { return ttFlag(ttId.ieta(), ttId.iphi()); }

  int sFGVB(int ieta, int iphi) const { return (raw(ieta, iphi) >>12) & 0x1; }
  int sFGVB(const EcalTrigTowerDetId& ttId) const { return sFGVB(ttId.ieta(), ttId.iphi()); }
  
  int l1aSpike(int ieta, int iphi) const { return sFGVB(ieta, iphi); }
  int l1aSpike(const EcalTrigTowerDetId& ttId) const { return sFGVB(ttId); }
  
  void toEcalTrigPrimDigiCollection(EcalTrigPrimDigiCollection& dest) const;
  
 private:
  int16_t formatVers_;
  std::vector<uint16_t> data_;
};

#endif //ECALTRIGPRIMCOMPACTCOLL_H not defined