CSCGEMData Class Reference

#include <CSCGEMData.h>

Public Member Functions
void	addEtaPadCluster (const GEMPadDigiCluster &digi, int gem_chamber, int eta_roll)
	Add and pack GEMPadDigiCluster digis. More...
bool	check () const
	CSCGEMData (int ntbins=12, int gem_fibers_mask=0xf)
	default constructor More...
	CSCGEMData (const unsigned short *c04buf, int length, int gem_fibers_mask=0xf)
unsigned short *	data ()
std::vector< GEMPadDigiCluster >	digis (int gem_chamber) const
std::vector< GEMPadDigiCluster >	etaDigis (int gem_chamber, int eta, int correctionToALCTbx) const
int	gemFibersMask () const
int	nTbins () const
int	numGEMEnabledFibers () const
int	numGEMs () const
void	print () const
int	sizeInWords () const

Static Public Member Functions
static void	setDebug (bool debugValue)
	turns on the debug flag for this class More...

Private Member Functions
int	getPartitionNumber (int address, int nPads) const
int	getPartitionStripNumber (int address, int nPads, int etaPart) const

Private Attributes
int	gems_enabled_
int	ngems_
int	ntbins_
int	size_
unsigned short	theData [8 *2 *32+2]

Static Private Attributes
static std::atomic< bool >	debug {false}

Detailed Description

Definition at line 11 of file CSCGEMData.h.

Constructor & Destructor Documentation

CSCGEMData() [1/2]

CSCGEMData::CSCGEMData	(	int	ntbins = 12,
		int	gem_fibers_mask = 0xf
	)

default constructor

initialize GEM data

Definition at line 18 of file CSCGEMData.cc.

References gems_enabled_, mps_fire::i, ngems_, ntbins_, size_, and theData.

                                                  : ntbins_(ntbins), size_(0) {
  theData[0] = 0x6C04;
  ntbins_ = ntbins;
  gems_enabled_ = gems_fibers;
  ngems_ = 0;
  /* Not implemented in the firmware yet */
  /*
  for (int i=0; i<4; i++)
    ngems_ += (gems_fibers>>i) & 0x1;
  */
  ngems_ = 4;
  size_ = 2 + ntbins_ * ngems_ * 4;

  for (int i = 0; i < (size_ - 2); i++) {
    int gem_chamber = (i % 16) / 8;
    theData[i + 1] = 0x3FFF | (gem_chamber << 14);
  }

  theData[size_ - 1] = 0x6D04;
}

CSCGEMData() [2/2]

CSCGEMData::CSCGEMData	(	const unsigned short *	c04buf,
		int	length,
		int	gem_fibers_mask = 0xf
	)

Definition at line 40 of file CSCGEMData.cc.

References visDQMUpload::buf, gems_enabled_, ngems_, ntbins_, size_, and theData.

                                                                             : size_(length) {
  size_ = length;
  ngems_ = 0;
  /* Not implemented in the firmware yet */
  /*
  for (int i=0; i<4; i++)
    ngems_ += (gems_fibers>>i) & 0x1;
  */
  ngems_ = 4;
  ntbins_ = (size_ - 2) / (4 * ngems_);
  gems_enabled_ = gems_fibers;
  memcpy(theData, buf, size_ * 2);
}

Member Function Documentation

addEtaPadCluster()

void CSCGEMData::addEtaPadCluster	(	const GEMPadDigiCluster &	digi,
		int	gem_chamber,
		int	eta_roll
	)

Add and pack GEMPadDigiCluster digis.

Add/pack GEMPadDigiCluster digi trigger objects per eta/roll gem_chamber - GEM GE11 layer gemA/B [0,1] eta_roll - GEM eta/roll 8 rolls per GEM layer [0-7]

Check that bx < GEM data max allocated tbins and that gem_chamber/layer is in 0,1 range

search for first free/empty cluster word

fill free cluster word if it was found

Definition at line 191 of file CSCGEMData.cc.

References nano_mu_digi_cff::bx, GEMPadDigiCluster::bx(), PVValHelper::eta, ntbins_, GEMPadDigiCluster::pads(), and theData.

Referenced by CSCEventData::add().

                                                                                              {
  int bx = digi.bx();
  if ((gem_chamber < 2) && (bx < ntbins_)) {
    int gem_layer = gem_chamber & 0x1;
    int cluster_size = digi.pads().size() - 1;
    int pad = digi.pads()[0];
    int eta = eta_roll;
    int cl_word = (gem_layer << 14) + (pad & 0xff) + ((eta & 0x7) << 8) + ((cluster_size & 0x7) << 11);
    int dataAddr = 1 + bx * 16 + 8 * gem_layer;
    int cluster_num = 0;
    while (((theData[dataAddr + cluster_num] & 0x3fff) != 0x3fff) && (cluster_num < 8)) {
      cluster_num++;
    }
    if (((theData[dataAddr + cluster_num] & 0x3fff) == 0x3fff) && (cluster_num < 8)) {
      theData[dataAddr + cluster_num] = cl_word;
    }
  }
}

check()

bool CSCGEMData::check ( ) const

inline

Definition at line 29 of file CSCGEMData.h.

References size_, and theData.

{ return ((theData[0] == 0x6C04) && (theData[size_ - 1] == 0x6D04)); }

data()

unsigned short* CSCGEMData::data ( )

inline

Definition at line 37 of file CSCGEMData.h.

References theData.

Referenced by CSCTMBData::pack().

{ return theData; }

digis()

std::vector< GEMPadDigiCluster > CSCGEMData::digis

(

int

gem_chamber

)

const

GEM data format v2

GEM data format v1 It is not used in the production Keeping this commented code just for v1 unpacking algo reference

Definition at line 67 of file CSCGEMData.cc.

References debug, TauDecayModes::dec, PVValHelper::eta, FCDTask_cfi::fiber, mps_fire::i, LogTrace, beam_dqm_sourceclient-live_cfg::maxClusters, ntbins_, mps_fire::result, theData, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by cscdqm::EventProcessor::processCSC().

                                                                    {
  std::vector<GEMPadDigiCluster> result;
  result.clear();
  int nPads = 192;  // From geometry
  int maxClusters = 4;
  int nGEMs = 4;
  // nGEMs = ngems_; // based on enabled fibers. not implemented in the firmware yet
  for (int i = 0; i < ntbins_; i++) {
    for (int fiber = 0; fiber < nGEMs; fiber++) {
      for (int cluster = 0; cluster < maxClusters; cluster++) {
        int dataAddr = 1 + (i * nGEMs + fiber) * maxClusters + cluster;
        int gem_layer = (theData[dataAddr] >> 14) & 0x1;  // gemA=0 or gemB=1
        if (gem_layer == gem_chamber) {
          int cl_word = theData[dataAddr] & 0x3fff;
          int pad = theData[dataAddr] & 0xff;
          int eta = (theData[dataAddr] >> 8) & 0x7;
          int cluster_size = (theData[dataAddr] >> 11) & 0x7;
          if (pad < nPads) {
            int padInPart = eta * nPads + pad;
            if (debug)
              LogTrace("CSCGEMData|CSCRawToDigi")
                  << "GEMlayer" << gem_layer << " cl_word" << dataAddr << ": 0x" << std::hex << cl_word << std::dec
                  << " tbin: " << i << " fiber#: " << (fiber + 1) << " cluster#: " << (cluster + 1)
                  << " padInPart: " << padInPart << " pad: " << pad << " eta: " << eta
                  << " cluster_size: " << cluster_size << std::endl;
            std::vector<short unsigned int> pads;
            for (int iP = 0; iP <= cluster_size; ++iP)
              pads.push_back(padInPart + iP);
            GEMPadDigiCluster pad_cluster(pads, i);
            result.push_back(pad_cluster);
          }
        }
      }
    }
  }

  /*
  std::vector<GEMPadDigiCluster> result;
  result.clear();
  int nPads = 192; // From geometry
  int maxAddr = 1536;
  int nGEMs = 2;
  int maxClusters = 8;
  // std::cout << std::hex << "markers " << theData[0] << ": " << theData[size_-1] << std::dec << " size: " << size_ << std::endl;
  for (int i=0; i<ntbins_; i++)
    {
      for (int gem=0; gem<nGEMs; gem++)
        {
          if (gem==gem_chamber) // Return only digis for specified GEM chamber
            {
              for (int TMBCluster=0; TMBCluster<maxClusters; TMBCluster++)
                {
                  int dataAddr = 1 + (i*nGEMs+gem)*maxClusters + TMBCluster;
                  int address = theData[dataAddr] & 0x7ff;
                  // std::cout << dataAddr << ": " << address <<std::endl;
                  int nExtraPads = (theData[dataAddr] >>11)&0x7;
                  if (address<maxAddr)
                    {
                      int  etaPart   = getPartitionNumber(address,nPads);
                      int  padInPart = getPartitionStripNumber(address,nPads,etaPart);
                      vector<short unsigned int> pads;
                      for(int iP = 0; iP <= nExtraPads; ++iP)
                        pads.push_back(padInPart + iP );
                      GEMPadDigiCluster cluster ( pads, i);
                      result.push_back(cluster);
                    }
                }
            }
        }
    }
  */
  return result;
}

etaDigis()

std::vector< GEMPadDigiCluster > CSCGEMData::etaDigis	(	int	gem_chamber,
		int	eta_roll,
		int	alctMatchTime
	)		const

Unpack GEMPadDigiCluster digi trigger objects per eta/roll gem_chamber - GEM GE11 layer gemA/B [0,1] eta_roll - GEM eta/roll 8 rolls per GEM layer [0-7]

GEM data format v2

Definition at line 148 of file CSCGEMData.cc.

References debug, TauDecayModes::dec, PVValHelper::eta, FCDTask_cfi::fiber, mps_fire::i, LogTrace, beam_dqm_sourceclient-live_cfg::maxClusters, ntbins_, mps_fire::result, GEMPadDigiCluster::setAlctMatchTime(), theData, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by cscdqm::EventProcessor::processCSC().

                                                                                                        {
  std::vector<GEMPadDigiCluster> result;
  result.clear();
  int nPads = 192;  // From geometry
  int maxClusters = 4;
  int nGEMs = 4;
  // nGEMs = ngems_; // based on enabled fibers. currently not implemented in the firmware
  for (int i = 0; i < ntbins_; i++) {
    for (int fiber = 0; fiber < nGEMs; fiber++) {
      for (int cluster = 0; cluster < maxClusters; cluster++) {
        int dataAddr = 1 + (i * nGEMs + fiber) * maxClusters + cluster;
        int gem_layer = (theData[dataAddr] >> 14) & 0x1;  // gemA=0 or gemB=1
        if (gem_layer == gem_chamber) {
          int cl_word = theData[dataAddr] & 0x3fff;
          int pad = theData[dataAddr] & 0xff;
          int eta = (theData[dataAddr] >> 8) & 0x7;
          int cluster_size = (theData[dataAddr] >> 11) & 0x7;
          if ((pad < nPads) && (eta == eta_roll)) {
            int padInPart = pad;
            if (debug)
              LogTrace("CSCGEMData|CSCRawToDigi")
                  << "GEMlayer" << gem_layer << " cl_word" << dataAddr << ": 0x" << std::hex << cl_word << std::dec
                  << " tbin: " << i << " fiber#: " << (fiber + 1) << " cluster#: " << (cluster + 1)
                  << " padInPart: " << padInPart << " pad: " << pad << " eta: " << eta
                  << " cluster_size: " << cluster_size << " alctMatchTime: " << alctMatchTime << std::endl;
            std::vector<short unsigned int> pads;
            for (int iP = 0; iP <= cluster_size; ++iP)
              pads.push_back(padInPart + iP);
            GEMPadDigiCluster pad_cluster(pads, i);
            pad_cluster.setAlctMatchTime(alctMatchTime);
            result.push_back(pad_cluster);
          }
        }
      }
    }
  }
  return result;
}

gemFibersMask()

int CSCGEMData::gemFibersMask ( ) const

inline

Definition at line 25 of file CSCGEMData.h.

References gems_enabled_.

{ return gems_enabled_; }

getPartitionNumber()

int CSCGEMData::getPartitionNumber ( int  address, int  nPads  ) const

private

Definition at line 61 of file CSCGEMData.cc.

References generateTowerEtThresholdLUT::addr.

{ return addr / (npads - 1); }

getPartitionStripNumber()

int CSCGEMData::getPartitionStripNumber ( int  address, int  nPads, int  etaPart  ) const

private

Definition at line 63 of file CSCGEMData.cc.

                                                                                 {
  return address - (nPads * etaPart);
}

nTbins()

int CSCGEMData::nTbins ( ) const

inline

Definition at line 27 of file CSCGEMData.h.

References ntbins_.

{ return ntbins_; }

numGEMEnabledFibers()

int CSCGEMData::numGEMEnabledFibers ( ) const

inline

Definition at line 26 of file CSCGEMData.h.

References ngems_.

{ return ngems_; }

numGEMs()

int CSCGEMData::numGEMs ( ) const

inline

Definition at line 22 of file CSCGEMData.h.

Referenced by cscdqm::EventProcessor::processCSC().

                      {
    return 2;  // !!! TODO actual number of GEM chambers in readout
  }

print()

void CSCGEMData::print

(

void

)

const

Definition at line 54 of file CSCGEMData.cc.

References mps_splice::line, LogTrace, size_, and theData.

                             {
  LogTrace("CSCGEMData|CSCRawToDigi") << "CSCGEMData.Print";
  for (int line = 0; line < ((size_)); ++line) {
    LogTrace("CSCGEMData|CSCRawToDigi") << std::hex << theData[line];
  }
}

setDebug()

static void CSCGEMData::setDebug ( bool  debugValue )

inlinestatic

turns on the debug flag for this class

Definition at line 32 of file CSCGEMData.h.

References debug.

{ debug = debugValue; }

sizeInWords()

int CSCGEMData::sizeInWords ( ) const

inline

Definition at line 21 of file CSCGEMData.h.

References size_.

Referenced by CSCTMBData::CSCTMBData(), CSCTMBData::pack(), and CSCTMBData::UnpackTMB().

{ return size_; }

Member Data Documentation

debug

std::atomic< bool > CSCGEMData::debug {false}

staticprivate

Definition at line 46 of file CSCGEMData.h.

Referenced by digis(), util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), etaDigis(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), runTauIdMVA.TauIDEmbedder::load_againstElectronMVA6(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), setDebug(), and pkg.AbstractPkg::write().

gems_enabled_

int CSCGEMData::gems_enabled_

private

Definition at line 50 of file CSCGEMData.h.

Referenced by CSCGEMData(), and gemFibersMask().

ngems_

int CSCGEMData::ngems_

private

Definition at line 51 of file CSCGEMData.h.

Referenced by CSCGEMData(), and numGEMEnabledFibers().

ntbins_

int CSCGEMData::ntbins_

private

Definition at line 49 of file CSCGEMData.h.

Referenced by addEtaPadCluster(), CSCGEMData(), digis(), etaDigis(), and nTbins().

size_

int CSCGEMData::size_

private

Definition at line 52 of file CSCGEMData.h.

Referenced by check(), CSCGEMData(), print(), and sizeInWords().

theData

unsigned short CSCGEMData::theData[8 *2 *32+2]

private

Definition at line 53 of file CSCGEMData.h.

Referenced by addEtaPadCluster(), check(), CSCGEMData(), data(), digis(), etaDigis(), and print().