EcalSimRawData Class Reference

#include <EcalSimRawData.h>

Inheritance diagram for EcalSimRawData:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	EcalSimRawData (const edm::ParameterSet &pset)
	~EcalSimRawData () override=default
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Types
enum	tokenType { tcp = 0, tp = 1 }
enum	writeMode_t { littleEndian, bigEndian, ascii }

Private Member Functions
void	checkGeometry (const edm::EventSetup &eventSetup)
void	checkTriggerMap (const edm::EventSetup &eventSetup)
int	cIndex2iTtEta (int iTtEta0) const
int	cIndex2TtPhi (int iTtPhi0) const
void	elec2GeomNum (int ittEta0, int ittPhi0, int strip1, int ch1, int &iEta0, int &iPhi0) const
void	fwrite (std::ofstream &f, uint16_t data, int &iword, bool hpar=true) const
void	genFeData (std::string &basename, int iEvent, const std::vector< uint16_t > adcCount[nEbEta][nEbPhi]) const
void	genSrData (std::string &basename, int iEvent, int ttf[nEbTtEta][nTtPhi]) const
void	genTccIn (std::string &basename, int iEvent, const int tps[nTtEta][nTtPhi]) const
void	genTccOut (std::string &basename, int iEvent, const int tps[nTtEta][nTtPhi]) const
void	getEbDigi (const edm::Event &event, std::vector< uint16_t > adc[nEbEta][nEbPhi]) const
std::string	getExt () const
void	getSrfs (const edm::Event &event, int ebSrf[nTtEta][nTtPhi], int eeSrf[nEndcaps][nScX][nScY]) const
void	getTp (const edm::Event &event, tokenType type, int tp[nTtEta][nTtPhi]) const
int	iEta2cIndex (int iEta) const
int	iPhi2cIndex (int iPhi) const
int	iTtEta2cIndex (int iTtEta) const
int	iTtPhi2cIndex (int iTtPhi) const
void	setHParity (uint16_t &a) const

Private Attributes
std::string	basename_
int	dccNum_
std::string	digiProducer_
std::string	ebDigiCollection_
edm::EDGetTokenT< EBDigiCollection >	ebDigisToken_
std::string	ebSrFlagCollection_
edm::EDGetTokenT< EBSrFlagCollection >	ebSrFlagToken_
std::string	eeDigiCollection_
std::string	eeSrFlagCollection_
edm::EDGetTokenT< EESrFlagCollection >	eeSrFlagToken_
std::unique_ptr< EcalSelectiveReadout >	esr_
bool	fe2dcc_
bool	fe2tcc_
int	iEvent
std::string	srDigiProducer_
std::ofstream	srfFile
bool	srp2dcc_
bool	tcc2dcc_
int	tccInDefaultVal_
int	tccNum_
std::string	tcpDigiCollection_
const CaloGeometry *	theGeometry
const EcalTrigTowerConstituentsMap *	theTriggerTowerMap
std::string	tpDigiCollection_
std::string	tpProducer_
bool	tpVerbose_
edm::EDGetTokenT< EcalTrigPrimDigiCollection >	trigPrimDigisToken_ [2]
std::ofstream	ttfFile
writeMode_t	writeMode_
bool	xtalVerbose_

Static Private Attributes
static const int	ebDccPhiEdge = 20
static const int	ebTccPhiEdge = 20
static const int	nDccEndcap = 9
static const int	nDccInPhi = 18
static const int	nEbEta = 170
static const int	nEbPhi = 360
static const int	nEbTtEta = nEbEta / ttEdge
static const int	nEeTtEta = 11
static const int	nEeX = 100
static const int	nEeY = 100
static const int	nEndcaps = 2
static const int	nScX = 20
static const int	nScY = 20
static const int	nTccEndcap = 36
static const int	nTccInPhi = 18
static const int	nTtEta = nEbTtEta + 2 * nEeTtEta
static const int	nTtPhi = nEbPhi / ttEdge
static const int	nTtPhisPerEbDcc = 4
static const int	nTtPhisPerEbTcc = 4
static const int	nTtSmEta = 17
static const int	nTtSmPhi = 4
static const int	nTtTypes = 2
static const int	scEdge = 5
static const int	strip2Eta [nTtTypes][ttEdge]
static const int	stripCh2Phi [nTtTypes][ttEdge][ttEdge]
static const int	ttEdge = 5
static const int	ttType [nEbTtEta]

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

The EcalSimRawData CMSSW module produces raw data from digis. The raw data are written into files which can be loaded into the TCC DCC and SRP boards in order to emulate the front-end. Only barrel is fully supported. The produced files for TCC assumes a special LUT in the TCC forged for FE emulation mode.

Module Parameters:

• string digiProducer: digi label
• string EBDigiCollection: EB crystal digi product instance name
• string EEDigiCollection: EE crystal digi product instance name
• string srProducer: sr flags digi label"
• string EBSrFlagCollection: EB SR flags product instance name
• string EESrFlagCollection: EE SR flags product instance name
• string trigPrimProducer: trigger primitive digi label"
• string trigPrimDigiCollection: trigger primitive digi product instance name
• string tcpPrimDigiCollection: TCP FENIX output trigger primitive digi product instance name
• string writeMode: output format. "write", "littleEndian", "bigEndian"
• untracked bool tpVerbose: make verbose the trigger primitive processing
• untracked bool xtalVerbose: make verbose the crystal digi processing
• untracked int32 dccNum: Id of the dcc raw data must be produced for. -1 means every DCC
• untracked int32 tccNum: Id of the tcc raw data must be produced for. -1 means every TCC
• untracked bool tcc2dccData: switch for TCC->DCC data stream production
• untracked bool srp2dccData: switch for SRP->DCC data stream production
• untracked int32 tccInDefaultVal: default TriggerPrimitive values if the trigger tower is abscent
• untracked string outputBaseName: basename for output files

Definition at line 42 of file EcalSimRawData.h.

Member Enumeration Documentation

tokenType

enum EcalSimRawData::tokenType

private

Enumerator
tcp
tp

Definition at line 59 of file EcalSimRawData.h.

{ tcp = 0, tp = 1 };

writeMode_t

enum EcalSimRawData::writeMode_t

private

Output format mode

littleEndian: little endian binary bigEndian: big endian binary ascii: ascii mode. The one accepted by the TCC, DCC and SRP board control software.

Enumerator
littleEndian
bigEndian
ascii

Definition at line 187 of file EcalSimRawData.h.

{ littleEndian, bigEndian, ascii };

Constructor & Destructor Documentation

EcalSimRawData()

EcalSimRawData::EcalSimRawData ( const edm::ParameterSet &  pset )

explicit

Constructor

Parameters

pset	CMSSW configuration

Definition at line 36 of file EcalSimRawData.cc.

References ascii, basename_, bigEndian, dccNum_, digiProducer_, ebDigiCollection_, ebDigisToken_, ebSrFlagCollection_, ebSrFlagToken_, eeDigiCollection_, eeSrFlagCollection_, eeSrFlagToken_, fe2dcc_, fe2tcc_, iEvent, ProducerED_cfi::InputTag, littleEndian, submitPVValidationJobs::params, srDigiProducer_, srp2dcc_, AlCaHLTBitMon_QueryRunRegistry::string, tcc2dcc_, tccInDefaultVal_, tccNum_, tcp, tcpDigiCollection_, tp, tpDigiCollection_, tpProducer_, tpVerbose_, trigPrimDigisToken_, EcalSimRawData_cfi::writeMode, writeMode_, and xtalVerbose_.

                                                            {
  // sets up parameters:
  digiProducer_ = params.getParameter<std::string>("unsuppressedDigiProducer");
  ebDigiCollection_ = params.getParameter<std::string>("EBdigiCollection");
  eeDigiCollection_ = params.getParameter<std::string>("EEdigiCollection");
  srDigiProducer_ = params.getParameter<std::string>("srProducer");
  ebSrFlagCollection_ = params.getParameter<std::string>("EBSrFlagCollection");
  eeSrFlagCollection_ = params.getParameter<std::string>("EESrFlagCollection");
  tpDigiCollection_ = params.getParameter<std::string>("trigPrimDigiCollection");
  tcpDigiCollection_ = params.getParameter<std::string>("tcpDigiCollection");
  tpProducer_ = params.getParameter<std::string>("trigPrimProducer");
  xtalVerbose_ = params.getUntrackedParameter<bool>("xtalVerbose", false);
  tpVerbose_ = params.getUntrackedParameter<bool>("tpVerbose", false);
  tcc2dcc_ = params.getUntrackedParameter<bool>("tcc2dccData", true);
  srp2dcc_ = params.getUntrackedParameter<bool>("srp2dccData", true);
  fe2dcc_ = params.getUntrackedParameter<bool>("fe2dccData", true);
  fe2tcc_ = params.getUntrackedParameter<bool>("fe2tccData", true);
  dccNum_ = params.getUntrackedParameter<int>("dccNum", -1);
  tccNum_ = params.getUntrackedParameter<int>("tccNum", -1);
  tccInDefaultVal_ = params.getUntrackedParameter<int>("tccInDefaultVal", 0xffff);
  basename_ = params.getUntrackedParameter<std::string>("outputBaseName");

  iEvent = 0;

  std::string writeMode = params.getParameter<std::string>("writeMode");

  if (writeMode == std::string("littleEndian")) {
    writeMode_ = littleEndian;
  } else if (writeMode == std::string("bigEndian")) {
    writeMode_ = bigEndian;
  } else {
    writeMode_ = ascii;
  }

  eeSrFlagToken_ = consumes<EESrFlagCollection>(edm::InputTag(srDigiProducer_, eeSrFlagCollection_));
  ebSrFlagToken_ = consumes<EBSrFlagCollection>(edm::InputTag(srDigiProducer_, ebSrFlagCollection_));
  ebDigisToken_ = consumes<EBDigiCollection>(edm::InputTag(digiProducer_, ebDigiCollection_));
  trigPrimDigisToken_[EcalSimRawData::tcp] =
      consumes<EcalTrigPrimDigiCollection>(edm::InputTag(tpProducer_, tcpDigiCollection_));
  trigPrimDigisToken_[EcalSimRawData::tp] =
      consumes<EcalTrigPrimDigiCollection>(edm::InputTag(tpProducer_, tpDigiCollection_));
}

~EcalSimRawData()

EcalSimRawData::~EcalSimRawData ( )

overridedefault

Destructor

Member Function Documentation

analyze()

void EcalSimRawData::analyze ( const edm::Event &  event, const edm::EventSetup &  es  )

overridevirtual

Main method. Called back for each event. This method produced the raw data and write them to disk.

Implements edm::one::EDAnalyzerBase.

Definition at line 79 of file EcalSimRawData.cc.

References gpuClustering::adc, basename_, gather_cfg::cout, fe2dcc_, fe2tcc_, genFeData(), genSrData(), genTccIn(), genTccOut(), getEbDigi(), getSrfs(), getTp(), iEvent, nEbEta, nEbPhi, nEndcaps, nScX, nScY, nTtEta, nTtPhi, srp2dcc_, tcc2dcc_, tcp, tp, tpVerbose_, and xtalVerbose_.

                                                                           {
  // Event counter:
  ++iEvent;

  if (xtalVerbose_ | tpVerbose_) {
    std::cout << "=================================================================="
                 "====\n"
              << " Event " << iEvent << "\n"
              << "------------------------------------------------------------------"
              << "----\n";
  }

  if (fe2dcc_) {
    std::vector<uint16_t> adc[nEbEta][nEbPhi];
    getEbDigi(event, adc);
    genFeData(basename_, iEvent, adc);
  }

  if (fe2tcc_) {
    int tcp[nTtEta][nTtPhi] = {{0}};
    getTp(event, EcalSimRawData::tcp, tcp);
    genTccIn(basename_, iEvent, tcp);
  }

  if (tcc2dcc_) {
    int tp[nTtEta][nTtPhi] = {{0}};
    getTp(event, EcalSimRawData::tp, tp);
    genTccOut(basename_, iEvent, tp);
  }

  // SR flags:
  int ebSrf[nTtEta][nTtPhi];
  int eeSrf[nEndcaps][nScX][nScY];

  if (srp2dcc_) {
    getSrfs(event, ebSrf, eeSrf);
    genSrData(basename_, iEvent, ebSrf);
  }
}

checkGeometry()

void EcalSimRawData::checkGeometry ( const edm::EventSetup &  eventSetup )

private

call these once an event, to make sure everything is up-to-date

checkTriggerMap()

void EcalSimRawData::checkTriggerMap ( const edm::EventSetup &  eventSetup )

private

cIndex2iTtEta()

int EcalSimRawData::cIndex2iTtEta ( int  iTtEta0 ) const

inlineprivate

Converse of iTtEta2cIndex

Parameters

iTtEta0	c eta index of TT
std	CMSSW TT eta index

Definition at line 237 of file EcalSimRawData.h.

{ return (iTtEta0 < 28) ? (iTtEta0 - 28) : (iTtEta0 - 27); }

cIndex2TtPhi()

int EcalSimRawData::cIndex2TtPhi ( int  iTtPhi0 ) const

inlineprivate

Converse of iTtPhi2cIndex

Parameters

iTtPhi0

phi index of TT

Returns

std CMSS TT index

Definition at line 243 of file EcalSimRawData.h.

{ return iTtPhi0 + 1; }

elec2GeomNum()

void EcalSimRawData::elec2GeomNum ( int  ittEta0, int  ittPhi0, int  strip1, int  ch1, int &  iEta0, int &  iPhi0  ) const

private

Converts electronic number of an ECAL barrel channel to geometrical indices

Parameters

	ittEta0	trigger tower c index
	ittPhi0	trigger tower c index
	strip1	strip index within the TT. Runs from 1 to 5.
	ch1	channel electronics number within the VFE. Runs from 1 to 5.
[out]	iEta0	eta c index of the channel
[out]	iPhi0	eta c index of the channel

Definition at line 119 of file EcalSimRawData.cc.

References cms::cuda::assert(), nEbEta, nEbPhi, nEbTtEta, nTtPhi, strip2Eta, stripCh2Phi, ttEdge, and ttType.

Referenced by genFeData().

                                                                                                             {
  assert(0 <= ittEta0 && ittEta0 < nEbTtEta);
  assert(0 <= ittPhi0 && ittPhi0 < nTtPhi);
  assert(1 <= strip1 && strip1 <= ttEdge);
  assert(1 <= ch1 && ch1 <= ttEdge);
  const int type = ttType[ittEta0];
  iEta0 = ittEta0 * ttEdge + strip2Eta[type][strip1 - 1];
  iPhi0 = ittPhi0 * ttEdge + stripCh2Phi[type][strip1 - 1][ch1 - 1];
  assert(0 <= iEta0 && iEta0 < nEbEta);
  assert(0 <= iPhi0 && iPhi0 < nEbPhi);
}

fwrite()

void EcalSimRawData::fwrite ( std::ofstream &  f, uint16_t  data, int &  iword, bool  hpar = true  ) const

private

Write a data 16-bit word into file according to selected format.

Parameters

	f	the file stream to write to
	data	the peace of data to write
[in,out]	iword	pass zero when writing for the first time in a file, then the value returned by the previous call. Counts the number of words written into the file.
	hpar	if true the horizontal odd word parity is set before writing the word into the file.

Definition at line 131 of file EcalSimRawData.cc.

References ascii, bigEndian, DummyCfis::c, data, TauDecayModes::dec, f, littleEndian, setHParity(), and writeMode_.

Referenced by genFeData(), genSrData(), and genTccOut().

                                                                                      {
  if (hpar) {
    // set horizontal odd parity bit:
    setHParity(data);
  }

  switch (writeMode_) {
    case littleEndian: {
      char c = data & 0x00FF;
      f.write(&c, sizeof(c));
      c = (data >> 8) & 0x00FF;
      f.write(&c, sizeof(c));
    } break;
    case bigEndian: {
      char c = (data >> 8) & 0x00FF;
      f.write(&c, sizeof(c));
      c = data & 0x00FF;
      f.write(&c, sizeof(c));
    } break;
    case ascii:
      f << ((iWord % 8 == 0 && iWord != 0) ? "\n" : "") << "0x" << std::setfill('0') << std::setw(4) << std::hex << data
        << "\t" << std::dec << std::setfill(' ');
      break;
  }
  ++iWord;
}

genFeData()

void EcalSimRawData::genFeData ( std::string &  basename, int  iEvent, const std::vector< uint16_t >  adcCount[nEbEta][nEbPhi]  ) const

private

Generates FE crystal data

Parameters

basename	base for the output file name. DCC number is appended to the name
iEvent	event index
adcCount	the payload, the ADC count of the channels.

stripHeader |= parity(stripHeader) << 15;

Definition at line 171 of file EcalSimRawData.cc.

References gpuClustering::adc, ascii, gather_cfg::cout, data, dccNum_, TauDecayModes::dec, elec2GeomNum(), f, fwrite(), getExt(), iEvent, nDccEndcap, nDccInPhi, nTtPhi, nTtPhisPerEbDcc, nTtSmEta, nTtSmPhi, alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, ttEdge, writeMode_, testProducerWithPsetDescEmpty_cfi::x1, and xtalVerbose_.

Referenced by analyze().

                                                                                         {
  int smf = 0;
  int gmf = 0;
  int nPendingEvt = 0;
  int monitorFlag = 0;
  int chFrameLen = adcCount[0][0].size() + 1;

  int iWord = 0;

  for (int iZ0 = 0; iZ0 < 2; ++iZ0) {
    for (int iDccPhi0 = 0; iDccPhi0 < nDccInPhi; ++iDccPhi0) {
      int iDcc1 = iDccPhi0 + iZ0 * nDccInPhi + nDccEndcap + 1;

      if (dccNum_ != -1 && dccNum_ != iDcc1)
        continue;

      std::stringstream s;
      s.str("");
      const std::string &ext = getExt();
      s << basename << "_fe2dcc" << std::setfill('0') << std::setw(2) << iDcc1 << std::setfill(' ') << ext;
      std::ofstream f(s.str().c_str(), (iEvent == 1 ? std::ios::ate : std::ios::app));

      if (f.fail())
        return;

      if (writeMode_ == ascii) {
        f << (iEvent == 1 ? "" : "\n") << "[Event:" << iEvent << "]\n";
      }

      for (int iTtEtaInSm0 = 0; iTtEtaInSm0 < nTtSmEta; ++iTtEtaInSm0) {
        int iTtEta0 = iZ0 * nTtSmEta + iTtEtaInSm0;
        for (int iTtPhiInSm0 = 0; iTtPhiInSm0 < nTtSmPhi; ++iTtPhiInSm0) {
          // phi=0deg at middle of 1st barrel DCC:
          int iTtPhi0 = -nTtPhisPerEbDcc / 2 + iDccPhi0 * nTtPhisPerEbDcc + iTtPhiInSm0;
          if (iTtPhi0 < 0)
            iTtPhi0 += nTtPhi;
          for (int stripId1 = 1; stripId1 <= ttEdge; ++stripId1) {
            uint16_t stripHeader =
                0xF << 11 | (nPendingEvt & 0x3F) << 5 | (gmf & 0x1) << 4 | (smf & 0x1) << 3 | (stripId1 & 0x7);
            fwrite(f, stripHeader, iWord);

            for (int xtalId1 = 1; xtalId1 <= ttEdge; ++xtalId1) {
              uint16_t crystalHeader = 1 << 14 | (chFrameLen & 0xFF) << 4 | (monitorFlag & 0x1) << 3 | (xtalId1 & 0x7);
              //          crystalHeader |=parity(crystalHeader) << 15;
              fwrite(f, crystalHeader, iWord);

              int iEta0;
              int iPhi0;
              elec2GeomNum(iTtEta0, iTtPhi0, stripId1, xtalId1, iEta0, iPhi0);
              if (xtalVerbose_) {
                std::cout << std::dec << "iDcc1 = " << iDcc1 << "\t"
                          << "iEbTtEta0 = " << iTtEta0 << "\t"
                          << "iEbTtPhi0 = " << iTtPhi0 << "\t"
                          << "stripId1 = " << stripId1 << "\t"
                          << "xtalId1 = " << xtalId1 << "\t"
                          << "iEta0 = " << iEta0 << "\t"
                          << "iPhi0 = " << iPhi0 << "\t"
                          << "adc[5] = 0x" << std::hex << adcCount[iEta0][iPhi0][5] << std::dec << "\n";
              }

              const std::vector<uint16_t> &adc = adcCount[iEta0][iPhi0];
              for (unsigned iSample = 0; iSample < adc.size(); ++iSample) {
                uint16_t data = adc[iSample] & 0x3FFF;
                //      data |= parity(data);
                fwrite(f, data, iWord);
              }  // next time sample
            }    // next crystal in strip
          }      // next strip in TT
        }        // next TT along phi
      }          // next TT along eta
    }            // next DCC
  }              // next half-barrel
}

genSrData()

void EcalSimRawData::genSrData ( std::string &  basename, int  iEvent, int  ttf[nEbTtEta][nTtPhi]  ) const

private

Generates SR flags

Parameters

basename	base for the output file name. DCC number is appended to the name
iEvent	event index
the	trigger tower flags

Definition at line 248 of file EcalSimRawData.cc.

References ascii, data, dccNum_, f, fwrite(), getExt(), iEvent, nDccEndcap, nDccInPhi, nEeTtEta, nTtPhi, nTtPhisPerEbDcc, nTtSmEta, nTtSmPhi, alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, writeMode_, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by analyze().

                                                                                               {
  for (int iZ0 = 0; iZ0 < 2; ++iZ0) {
    for (int iDccPhi0 = 0; iDccPhi0 < nDccInPhi; ++iDccPhi0) {
      int iDcc1 = iDccPhi0 + iZ0 * nDccInPhi + nDccEndcap + 1;
      if (dccNum_ != -1 && dccNum_ != iDcc1)
        continue;
      std::stringstream s;
      s.str("");
      s << basename << "_ab2dcc" << std::setfill('0') << std::setw(2) << iDcc1 << std::setfill(' ') << getExt();
      std::ofstream f(s.str().c_str(), (iEvent == 1 ? std::ios::ate : std::ios::app));

      if (f.fail())
        throw cms::Exception(std::string("Cannot create/open file ") + s.str() + ".");

      int iWord = 0;

      if (writeMode_ == ascii) {
        f << (iEvent == 1 ? "" : "\n") << "[Event:" << iEvent << "]\n";
      }

      const uint16_t le1 = 0;
      const uint16_t le0 = 0;
      const uint16_t h1 = 1;
      const uint16_t nFlags = 68;
      uint16_t data = (h1 & 0x1) << 14 | (le1 & 0x1) << 12 | (le0 & 0x1) << 11 | (nFlags & 0x7F);

      fwrite(f, data, iWord, true);

      int iFlag = 0;
      data = 0;

      for (int iTtEtaInSm0 = 0; iTtEtaInSm0 < nTtSmEta; ++iTtEtaInSm0) {
        //  int iTtEbEta0 = iZ0*nTtSmEta + iTtEtaInSm0;
        int iTtEta0 = nEeTtEta + iZ0 * nTtSmEta + iTtEtaInSm0;
        for (int iTtPhiInSm0 = 0; iTtPhiInSm0 < nTtSmPhi; ++iTtPhiInSm0) {
          // phi=0deg at middle of 1st barrel DCC:
          int iTtPhi0 = -nTtPhisPerEbDcc / 2 + iDccPhi0 * nTtPhisPerEbDcc + iTtPhiInSm0;
          if (iTtPhi0 < 0)
            iTtPhi0 += nTtPhi;
          // flags are packed by four:
          //|15 |14 |13-12 |11      9|8      6|5      3|2      0|
          //| P | 0 | X  X |  srf i+3| srf i+2| srf i+1| srf i  |
          //|   |   |      | field 3 |field 2 | field 1| field 0|
          const int field = iFlag % 4;
          // std::cout << "TtEta0: " << iTtEta0 << "\tTtPhi0: " << iTtPhi0 << "\n";
          // std::cout << "#" << oct << (int)srf[iTtEta0][iTtPhi0] << "o ****> #" <<
          // oct << (srf[iTtEta0][iTtPhi0] << (field*3)) << "o\n" << std::dec;

          data |= srf[iTtEta0][iTtPhi0] << (field * 3);

          if (field == 3) {
            // std::cout <<  srf[iTtEta0][iTtPhi0] << "----> 0x" << std::hex << data <<
            // "\n";
            fwrite(f, data, iWord, true);
            data = 0;
          }
          ++iFlag;
        }  // next TT along phi
      }    // next TT along eta
    }      // next DCC
  }        // next half-barrel
}

genTccIn()

void EcalSimRawData::genTccIn ( std::string &  basename, int  iEvent, const int  tps[nTtEta][nTtPhi]  ) const

private

Generates FE trigger primitives data

Parameters

basename	base for the output file name. DCC number is appended to the name
iEvent	event index
tps	the payload, the trigger primitives

Definition at line 311 of file EcalSimRawData.cc.

References gather_cfg::cout, TauDecayModes::dec, iEvent, nTccEndcap, nTccInPhi, nTtPhi, nTtPhisPerEbTcc, nTtSmEta, nTtSmPhi, alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, tccNum_, tcp, and tpVerbose_.

Referenced by analyze().

                                                                                                  {
  for (int iZ0 = 0; iZ0 < 2; ++iZ0) {
    for (int iTccPhi0 = 0; iTccPhi0 < nTccInPhi; ++iTccPhi0) {
      int iTcc1 = iTccPhi0 + iZ0 * nTccInPhi + nTccEndcap + 1;

      if (tccNum_ != -1 && tccNum_ != iTcc1)
        continue;

      std::stringstream s;
      s.str("");
      const char *ext = ".txt";  // only ascii mode supported for TCP

      s << basename << "_tcc" << std::setfill('0') << std::setw(2) << iTcc1 << std::setfill(' ') << ext;
      std::ofstream fe2tcc(s.str().c_str(), (iEvent == 1 ? std::ios::ate : std::ios::app));

      if (fe2tcc.fail())
        throw cms::Exception(std::string("Failed to create file ") + s.str() + ".");

      int memPos = iEvent - 1;
      int iCh1 = 1;
      for (int iTtEtaInSm0 = 0; iTtEtaInSm0 < nTtSmEta; ++iTtEtaInSm0) {
        int iTtEta0 = (iZ0 == 0) ? 27 - iTtEtaInSm0 : 28 + iTtEtaInSm0;
        for (int iTtPhiInSm0 = 0; iTtPhiInSm0 < nTtSmPhi; ++iTtPhiInSm0) {
          // phi=0deg at middle of 1st barrel DCC:
          int iTtPhi0 = -nTtPhisPerEbTcc / 2 + iTccPhi0 * nTtPhisPerEbTcc + iTtPhiInSm0;
          iTtPhi0 += nTtPhisPerEbTcc * iTccPhi0;
          if (iTtPhi0 < 0)
            iTtPhi0 += nTtPhi;
          uint16_t tp_fe2tcc = (tcp[iTtEta0][iTtPhi0] & 0x7ff);  // keep only Et (9:0) and FineGrain (10)

          if (tpVerbose_ && tp_fe2tcc != 0) {
            std::cout << std::dec << "iTcc1 = " << iTcc1 << "\t"
                      << "iTtEta0 = " << iTtEta0 << "\t"
                      << "iTtPhi0 = " << iTtPhi0 << "\t"
                      << "iCh1 = " << iCh1 << "\t"
                      << "memPos = " << memPos << "\t"
                      << "tp = 0x" << std::setfill('0') << std::hex << std::setw(3) << tp_fe2tcc << std::dec
                      << std::setfill(' ') << "\n";
          }
          fe2tcc << iCh1 << "\t" << memPos << "\t" << std::setfill('0') << std::hex << "0x" << std::setw(4) << tp_fe2tcc
                 << "\t"
                 << "0" << std::dec << std::setfill(' ') << "\n";
          ++iCh1;
        }  // next TT along phi
      }    // next TT along eta
      fe2tcc << std::flush;
      fe2tcc.close();
    }  // next TCC
  }    // next half-barrel
}

genTccOut()

void EcalSimRawData::genTccOut ( std::string &  basename, int  iEvent, const int  tps[nTtEta][nTtPhi]  ) const

private

Generates TCC->DCC data

Parameters

basename	base for the output file name. DCC number is appended to the name
iEvent	event index
tps	the payload, the trigger primitives

Definition at line 362 of file EcalSimRawData.cc.

References gather_cfg::cout, data, TauDecayModes::dec, fwrite(), getExt(), iEvent, nEeTtEta, PresampleTask_cfi::nSamples, nTccEndcap, nTccInPhi, nTtPhi, nTtPhisPerEbTcc, nTtSmEta, nTtSmPhi, alignCSCRings::s, tccNum_, tpVerbose_, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by analyze().

                                                                                                   {
  int iDccWord = 0;

  for (int iZ0 = 0; iZ0 < 2; ++iZ0) {
    for (int iTccPhi0 = 0; iTccPhi0 < nTccInPhi; ++iTccPhi0) {
      int iTcc1 = iTccPhi0 + iZ0 * nTccInPhi + nTccEndcap + 1;

      if (tccNum_ != -1 && tccNum_ != iTcc1)
        continue;

      std::stringstream s;
      s.str("");
      const char *ext = ".txt";  // only ascii mode supported for TCP

      s << basename << "_tcc" << std::setfill('0') << std::setw(2) << iTcc1 << std::setfill(' ') << ext;

      s.str("");
      s << basename << "_tcc2dcc" << std::setfill('0') << std::setw(2) << iTcc1 << std::setfill(' ') << getExt();
      std::ofstream dccF(s.str().c_str(), (iEvent == 1 ? std::ios::ate : std::ios::app));

      if (dccF.fail()) {
        std::cout << "Warning: failed to create or open file " << s.str() << ".\n";
        return;
      }

      const uint16_t h1 = 1;
      const uint16_t le1 = 0;
      const uint16_t le0 = 0;
      const uint16_t nSamples = 1;
      const uint16_t nTts = 68;
      const uint16_t data =
          (h1 & 0x1) << 14 | (le1 & 0x1) << 12 | (le0 & 0x1) << 11 | (nSamples & 0xF) << 7 | (nTts & 0x7F);
      dccF << (iEvent == 1 ? "" : "\n") << "[Event:" << iEvent << "]\n";
      fwrite(dccF, data, iDccWord, false);

      int memPos = iEvent - 1;
      int iCh1 = 1;
      for (int iTtEtaInSm0 = 0; iTtEtaInSm0 < nTtSmEta; ++iTtEtaInSm0) {
        int iTtEta0 = nEeTtEta + iZ0 * nTtSmEta + iTtEtaInSm0;
        for (int iTtPhiInSm0 = 0; iTtPhiInSm0 < nTtSmPhi; ++iTtPhiInSm0) {
          // phi=0deg at middle of 1st barrel DCC:
          int iTtPhi0 = -nTtPhisPerEbTcc / 2 + iTccPhi0 * nTtPhisPerEbTcc + iTtPhiInSm0;
          if (iTtPhi0 < 0)
            iTtPhi0 += nTtPhi;

          if (tpVerbose_) {
            std::cout << std::dec << "iTcc1 = " << iTcc1 << "\t"
                      << "iTtEta0 = " << iTtEta0 << "\t"
                      << "iTtPhi0 = " << iTtPhi0 << "\t"
                      << "iCh1 = " << iCh1 << "\t"
                      << "memPos = " << memPos << "\t"
                      << "tp = 0x" << std::hex << tps[iTtEta0][iTtPhi0] << std::dec << "\n";
          }
          fwrite(dccF, tps[iTtEta0][iTtPhi0], iDccWord, false);
          ++iCh1;
        }  // next TT along phi
      }    // next TT along eta
    }      // next TCC
  }        // next half-barrel
}

getEbDigi()

void EcalSimRawData::getEbDigi ( const edm::Event &  event, std::vector< uint16_t >  adc[nEbEta][nEbPhi]  ) const

private

Retrieves barrel digis (APD ADC count).

Parameters

event	CMS event
adc	[out] the adc counts: adc[iEta0][iPhi0][iTimeSample0]

Definition at line 480 of file EcalSimRawData.cc.

References gpuClustering::adc, gather_cfg::cout, TauDecayModes::dec, ebDigisToken_, Exception, hcalRecHitTable_cff::ieta, iEta2cIndex(), hcalRecHitTable_cff::iphi, iPhi2cIndex(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, nEbEta, nEbPhi, PresampleTask_cfi::nSamples, ecalGpuTask_cfi::sample, and xtalVerbose_.

Referenced by analyze().

                                                                                                   {
  const auto &hEbDigis = event.getHandle(ebDigisToken_);

  int nSamples = 0;
  if (hEbDigis.isValid() && !hEbDigis->empty()) {  // there is at least one digi
    nSamples = hEbDigis->begin()->size();          // gets the sample count from 1st digi
  }

  const uint16_t suppressed = 0xFFFF;

  adc[0][0] = std::vector<uint16_t>(nSamples, suppressed);

  for (int iEbEta = 0; iEbEta < nEbEta; ++iEbEta) {
    for (int iEbPhi = 0; iEbPhi < nEbPhi; ++iEbPhi) {
      adc[iEbEta][iEbPhi] = adc[0][0];
    }
  }
  if (hEbDigis.isValid()) {
    if (xtalVerbose_)
      std::cout << std::setfill('0');
    for (EBDigiCollection::const_iterator it = hEbDigis->begin(); it != hEbDigis->end(); ++it) {
      const EBDataFrame &frame = *it;

      int iEta0 = iEta2cIndex((frame.id()).ieta());
      int iPhi0 = iPhi2cIndex((frame.id()).iphi());

      //     std::cout << "xtl indices conv: (" << frame.id().ieta() << ","
      //     << frame.id().iphi() << ") -> ("
      //     << iEta0 << "," << iPhi0 << ")\n";

      if (iEta0 < 0 || iEta0 >= nEbEta) {
        std::cout << "iEta0 (= " << iEta0 << ") is out of range ("
                  << "[0," << nEbEta - 1 << "])\n";
      }
      if (iPhi0 < 0 || iPhi0 >= nEbPhi) {
        std::cout << "iPhi0 (= " << iPhi0 << ") is out of range ("
                  << "[0," << nEbPhi - 1 << "])\n";
      }

      if (xtalVerbose_) {
        std::cout << iEta0 << "\t" << iPhi0 << ":\t";
        std::cout << std::hex;
      }

      if (nSamples != frame.size()) {
        throw cms::Exception("EcalSimRawData",
                             "Found EB digis with different sample count! This "
                             "is not supported by EcalSimRawData.");
      }

      for (int iSample = 0; iSample < nSamples; ++iSample) {
        const EcalMGPASample &sample = frame.sample(iSample);
        uint16_t encodedAdc = sample.raw();
        adc[iEta0][iPhi0][iSample] = encodedAdc;
        if (xtalVerbose_) {
          std::cout << (iSample > 0 ? " " : "") << "0x" << std::setw(4) << encodedAdc;
        }
      }
      if (xtalVerbose_)
        std::cout << "\n" << std::dec;
    }
    if (xtalVerbose_)
      std::cout << std::setfill(' ');
  }
}

getExt()

std::string EcalSimRawData::getExt ( ) const

private

Help function to get the file extension which depends on the output formats.

Definition at line 158 of file EcalSimRawData.cc.

References ascii, bigEndian, littleEndian, and writeMode_.

Referenced by genFeData(), genSrData(), and genTccOut().

                                       {
  switch (writeMode_) {
    case littleEndian:
      return ".le";
    case bigEndian:
      return ".be";
    case ascii:
      return ".txt";
    default:
      return ".?";
  }
}

getSrfs()

void EcalSimRawData::getSrfs ( const edm::Event &  event, int  ebSrf[nTtEta][nTtPhi], int  eeSrf[nEndcaps][nScX][nScY]  ) const

private

Computes the selective readout flags.

Parameters

[in]	event	CMS event
[out]	ebSrf	the computed SR flags for barrel
[out]	eeSrf	the computed SR flags for endcaps

Definition at line 434 of file EcalSimRawData.cc.

References cms::cuda::assert(), ebSrFlagCollection_, ebSrFlagToken_, eeSrFlagCollection_, eeSrFlagToken_, RemoveAddSevLevel::flag, mps_fire::i, l1tPhase2CaloJetEmulator_cfi::iEta, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, nEbTtEta, nEeTtEta, nEndcaps, nScX, nScY, nTtEta, nTtPhi, and srDigiProducer_.

Referenced by analyze().

                                                                    {
  // EE
  const auto &hEeSrFlags = event.getHandle(eeSrFlagToken_);
  for (size_t i = 0; i < (nEndcaps * nScX * nScY); ((int *)eeSrf)[i++] = -1) {
  };
  if (hEeSrFlags.isValid()) {
    for (EESrFlagCollection::const_iterator it = hEeSrFlags->begin(); it != hEeSrFlags->end(); ++it) {
      const EESrFlag &flag = *it;
      int iZ0 = flag.id().zside() > 0 ? 1 : 0;
      int iX0 = flag.id().ix() - 1;
      int iY0 = flag.id().iy() - 1;
      assert(iZ0 >= 0 && iZ0 < nEndcaps);
      assert(iX0 >= 0 && iX0 < nScX);
      assert(iY0 >= 0 && iY0 < nScY);
      eeSrf[iZ0][iX0][iY0] = flag.value();
    }
  } else {
    edm::LogWarning("EcalSimRawData") << "EE SR flag not found (Product label: " << srDigiProducer_
                                      << "Producet instance: " << eeSrFlagCollection_ << ")";
  }

  // EB
  const auto &hEbSrFlags = event.getHandle(ebSrFlagToken_);
  for (size_t i = 0; i < (nTtEta * nTtPhi); ((int *)ebSrf)[i++] = -1) {
  };
  if (hEbSrFlags.isValid()) {
    for (EBSrFlagCollection::const_iterator it = hEbSrFlags->begin(); it != hEbSrFlags->end(); ++it) {
      const EBSrFlag &flag = *it;
      int iEta = flag.id().ieta();
      int iEta0 = iEta + nTtEta / 2 - (iEta >= 0 ? 1 : 0);  // 0->55 from eta=-3 to eta=3
      int iEbEta0 = iEta0 - nEeTtEta;                       // 0->33 from eta=-1.48 to eta=1.48
      int iPhi0 = flag.id().iphi() - 1;

      assert(iEbEta0 >= 0 && iEbEta0 < nEbTtEta);
      assert(iPhi0 >= 0 && iPhi0 < nTtPhi);

      ebSrf[iEbEta0][iPhi0] = flag.value();
    }
  } else {
    edm::LogWarning("EcalSimRawData") << "EB SR flag not found (Product label: " << srDigiProducer_
                                      << "Producet instance: " << ebSrFlagCollection_ << ")";
  }
}

getTp()

void EcalSimRawData::getTp ( const edm::Event &  event, EcalSimRawData::tokenType  type, int  tp[nTtEta][nTtPhi]  ) const

private

Extracts the trigger primitive (TP). The collection name parameter permits to select either the TCP Fenix output or the TCC output. Note: TCP output is only valid for the barrel.

Parameters

event	CMS event
collName	label of the EDM collection containing the TP.
tp	[out] the trigger primitives

Definition at line 546 of file EcalSimRawData.cc.

References edm::SortedCollection< T, SORT >::begin(), trigObjTnPSource_cfi::collName, gather_cfg::cout, edm::SortedCollection< T, SORT >::end(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, iTtEta2cIndex(), iTtPhi2cIndex(), nEbTtEta, nTtEta, nTtPhi, AlCaHLTBitMon_QueryRunRegistry::string, tccInDefaultVal_, tcp, tcpDigiCollection_, tp, tpDigiCollection_, tpVerbose_, and trigPrimDigisToken_.

Referenced by analyze().

                                                                                                             {
  const auto &hTpDigis = event.getHandle(trigPrimDigisToken_[type]);
  if (hTpDigis.isValid() && !hTpDigis->empty()) {
    const EcalTrigPrimDigiCollection &tpDigis = *hTpDigis.product();

    //    EcalSelectiveReadout::ttFlag_t ttf[nTtEta][nTtPhi];
    for (int iTtEta0 = 0; iTtEta0 < nTtEta; ++iTtEta0) {
      for (int iTtPhi0 = 0; iTtPhi0 < nTtPhi; ++iTtPhi0) {
        tcp[iTtEta0][iTtPhi0] = tccInDefaultVal_;
      }
    }
    if (tpVerbose_) {
      std::cout << std::setfill('0');
    }
    for (EcalTrigPrimDigiCollection::const_iterator it = tpDigis.begin(); it != tpDigis.end(); ++it) {
      const EcalTriggerPrimitiveDigi &tp = *it;
      int iTtEta0 = iTtEta2cIndex(tp.id().ieta());
      int iTtPhi0 = iTtPhi2cIndex(tp.id().iphi());
      if (iTtEta0 < 0 || iTtEta0 >= nTtEta) {
        std::cout << "iTtEta0 (= " << iTtEta0 << ") is out of range ("
                  << "[0," << nEbTtEta - 1 << "])\n";
      }
      if (iTtPhi0 < 0 || iTtPhi0 >= nTtPhi) {
        std::cout << "iTtPhi0 (= " << iTtPhi0 << ") is out of range ("
                  << "[0," << nTtPhi - 1 << "])\n";
      }

      tcp[iTtEta0][iTtPhi0] = tp[tp.sampleOfInterest()].raw();

      if (tpVerbose_) {
        if (tcp[iTtEta0][iTtPhi0] != 0) {  // print non-zero values only
          std::string collName = (type == 0) ? tcpDigiCollection_ : tpDigiCollection_;
          std::cout << collName << (collName.empty() ? "" : " ") << "TP(" << std::setw(2) << iTtEta0 << "," << iTtPhi0
                    << ") = "
                    << "0x" << std::setw(4) << tcp[iTtEta0][iTtPhi0] << "\tcmssw indices: " << tp.id().ieta() << " "
                    << tp.id().iphi() << "\n";
        }
      }
    }  // next TP
    if (tpVerbose_)
      std::cout << std::setfill(' ');
  }
}

iEta2cIndex()

int EcalSimRawData::iEta2cIndex ( int  iEta ) const

inlineprivate

Converts std CMSSW crystal eta index into a c-index (contiguous integer starting from 0 and increasing with pseudo-rapidity).

Parameters

iEta	std CMSSW crystal eta index

Returns

the c-array index

Definition at line 211 of file EcalSimRawData.h.

References l1tPhase2CaloJetEmulator_cfi::iEta.

Referenced by getEbDigi().

{ return (iEta < 0) ? iEta + 85 : iEta + 84; }

iPhi2cIndex()

int EcalSimRawData::iPhi2cIndex ( int  iPhi ) const

inlineprivate

Converts std CMSSW crystal phi index into a c-index (contiguous integer starting from 0 at phi=0deg and increasing with phi).

Parameters

iPhi	std CMSSW crystal phi index

Returns

the c-array index

Definition at line 218 of file EcalSimRawData.h.

References nEbPhi.

Referenced by getEbDigi().

                                  {
    int iPhi0 = iPhi - 11;
    if (iPhi0 < 0)
      iPhi0 += nEbPhi;
    return iPhi0;
  }

iTtEta2cIndex()

int EcalSimRawData::iTtEta2cIndex ( int  iTtEta ) const

inlineprivate

Converts std CMSSW ECAL trigger tower eta index into a c-index (contiguous integer starting from 0 and increasing with pseudo-rapidity).

Parameters

iEta	std CMSSW trigger tower eta index

Returns

the c-array index

Definition at line 231 of file EcalSimRawData.h.

Referenced by getTp().

{ return (iTtEta < 0) ? (iTtEta + 28) : (iTtEta + 27); }

iTtPhi2cIndex()

int EcalSimRawData::iTtPhi2cIndex ( int  iTtPhi ) const

inlineprivate

Converts std CMSSW ECAL trigger tower phi index into a c-index (contiguous integer starting from 0 at phi=0deg and increasing with phi).

Parameters

iPhi	std CMSSW ECAL trigger tower phi index

Returns

the c-array index

Definition at line 250 of file EcalSimRawData.h.

Referenced by getTp().

{ return iTtPhi - 1; }

setHParity()

void EcalSimRawData::setHParity ( uint16_t &  a ) const

private

Definition at line 423 of file EcalSimRawData.cc.

References a, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by fwrite().

                                                 {
  const int odd = 1 << 15;
  const int even = 0;
  // parity bit of numbers from 0x0 to 0xF:
  //                    0   1   2    3   4    5    6   7   8    9    A   B    C
  //                    D   E    F
  const int p[16] = {even, odd, odd, even, odd, even, even, odd, odd, even, even, odd, even, odd, odd, even};
  // inverts parity bit (LSB) of 'a' in case of even parity:
  a ^= p[a & 0xF] ^ p[(a >> 4) & 0xF] ^ p[(a >> 8) & 0xF] ^ p[a >> 12 & 0xF] ^ odd;
}

Member Data Documentation

basename_

std::string EcalSimRawData::basename_

private

basename for output files

Definition at line 444 of file EcalSimRawData.h.

Referenced by analyze(), and EcalSimRawData().

dccNum_

int EcalSimRawData::dccNum_

private

Index of the DCC, FE data must be produced for. -1 for all TTCs

Definition at line 436 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), genFeData(), and genSrData().

digiProducer_

std::string EcalSimRawData::digiProducer_

private

Name of module/plugin/producer making digis

Definition at line 348 of file EcalSimRawData.h.

Referenced by EcalSimRawData().

ebDccPhiEdge

const int EcalSimRawData::ebDccPhiEdge = 20

staticprivate

Number of barrel crystals along phi covered by a DCC

Definition at line 145 of file EcalSimRawData.h.

ebDigiCollection_

std::string EcalSimRawData::ebDigiCollection_

private

EB digi product instance name

Definition at line 352 of file EcalSimRawData.h.

Referenced by EcalSimRawData().

ebDigisToken_

edm::EDGetTokenT<EBDigiCollection> EcalSimRawData::ebDigisToken_

private

Definition at line 448 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getEbDigi().

ebSrFlagCollection_

std::string EcalSimRawData::ebSrFlagCollection_

private

EB SRP flag digi product instance name

Definition at line 364 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getSrfs().

ebSrFlagToken_

edm::EDGetTokenT<EBSrFlagCollection> EcalSimRawData::ebSrFlagToken_

private

Definition at line 447 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getSrfs().

ebTccPhiEdge

const int EcalSimRawData::ebTccPhiEdge = 20

staticprivate

number of TTs along phi of a TCC sector

Definition at line 133 of file EcalSimRawData.h.

eeDigiCollection_

std::string EcalSimRawData::eeDigiCollection_

private

EE digi product instance name

Definition at line 356 of file EcalSimRawData.h.

Referenced by EcalSimRawData().

eeSrFlagCollection_

std::string EcalSimRawData::eeSrFlagCollection_

private

EE SRP flag digi product instancename

Definition at line 368 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getSrfs().

eeSrFlagToken_

edm::EDGetTokenT<EESrFlagCollection> EcalSimRawData::eeSrFlagToken_

private

Definition at line 446 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getSrfs().

esr_

std::unique_ptr<EcalSelectiveReadout> EcalSimRawData::esr_

private

Selective readout simulator

Definition at line 420 of file EcalSimRawData.h.

fe2dcc_

bool EcalSimRawData::fe2dcc_

private

Switch for data of FE->DCC link

Definition at line 408 of file EcalSimRawData.h.

Referenced by analyze(), and EcalSimRawData().

fe2tcc_

bool EcalSimRawData::fe2tcc_

private

Switch for data of FE->TCC link

Definition at line 412 of file EcalSimRawData.h.

Referenced by analyze(), and EcalSimRawData().

iEvent

int EcalSimRawData::iEvent

private

Definition at line 61 of file EcalSimRawData.h.

Referenced by analyze(), EcalSimRawData(), genFeData(), genSrData(), genTccIn(), and genTccOut().

nDccEndcap

const int EcalSimRawData::nDccEndcap = 9

staticprivate

Number of DCCs for a single endcap

Definition at line 129 of file EcalSimRawData.h.

Referenced by genFeData(), and genSrData().

nDccInPhi

const int EcalSimRawData::nDccInPhi = 18

staticprivate

Number of barrel DCCs along Phi

Definition at line 125 of file EcalSimRawData.h.

Referenced by genFeData(), and genSrData().

nEbEta

const int EcalSimRawData::nEbEta = 170

staticprivate

Number of crystals in ECAL barrel along eta

Definition at line 65 of file EcalSimRawData.h.

Referenced by analyze(), elec2GeomNum(), and getEbDigi().

nEbPhi

const int EcalSimRawData::nEbPhi = 360

staticprivate

Number of crystals in ECAL barrel along phi

Definition at line 69 of file EcalSimRawData.h.

Referenced by analyze(), elec2GeomNum(), getEbDigi(), and iPhi2cIndex().

nEbTtEta

const int EcalSimRawData::nEbTtEta = nEbEta / ttEdge

staticprivate

Number of TTs along Ecal barrel eta

Definition at line 113 of file EcalSimRawData.h.

Referenced by elec2GeomNum(), getSrfs(), and getTp().

nEeTtEta

const int EcalSimRawData::nEeTtEta = 11

staticprivate

Number of TTs along eta for one endcap.

Definition at line 117 of file EcalSimRawData.h.

Referenced by genSrData(), genTccOut(), and getSrfs().

nEeX

const int EcalSimRawData::nEeX = 100

staticprivate

X-edge of endcap (x,y)- crystal grid

Definition at line 73 of file EcalSimRawData.h.

nEeY

const int EcalSimRawData::nEeY = 100

staticprivate

Y-edge of endcap (x,y)- crystal grid

Definition at line 77 of file EcalSimRawData.h.

nEndcaps

const int EcalSimRawData::nEndcaps = 2

staticprivate

Number of endcaps

Definition at line 81 of file EcalSimRawData.h.

Referenced by analyze(), and getSrfs().

nScX

const int EcalSimRawData::nScX = 20

staticprivate

Maximum number of supercrystal along x axis

Definition at line 89 of file EcalSimRawData.h.

Referenced by analyze(), and getSrfs().

nScY

const int EcalSimRawData::nScY = 20

staticprivate

Maximum number of supercrystal along y axis

Definition at line 93 of file EcalSimRawData.h.

Referenced by analyze(), and getSrfs().

nTccEndcap

const int EcalSimRawData::nTccEndcap = 36

staticprivate

Number of TCCs for a single endcap

Definition at line 141 of file EcalSimRawData.h.

Referenced by genTccIn(), and genTccOut().

nTccInPhi

const int EcalSimRawData::nTccInPhi = 18

staticprivate

Number of Barrel TTCs along phi

Definition at line 137 of file EcalSimRawData.h.

Referenced by genTccIn(), and genTccOut().

nTtEta

const int EcalSimRawData::nTtEta = nEbTtEta + 2 * nEeTtEta

staticprivate

Number of TTs along ECAL eta

Definition at line 121 of file EcalSimRawData.h.

Referenced by analyze(), getSrfs(), and getTp().

nTtPhi

const int EcalSimRawData::nTtPhi = nEbPhi / ttEdge

staticprivate

Number of TTs along Ecal Phi

Definition at line 109 of file EcalSimRawData.h.

Referenced by analyze(), elec2GeomNum(), genFeData(), genSrData(), genTccIn(), genTccOut(), getSrfs(), and getTp().

nTtPhisPerEbDcc

const int EcalSimRawData::nTtPhisPerEbDcc = 4

staticprivate

Number of trigger towers alng phi covered by a DCC

Definition at line 149 of file EcalSimRawData.h.

Referenced by genFeData(), and genSrData().

nTtPhisPerEbTcc

const int EcalSimRawData::nTtPhisPerEbTcc = 4

staticprivate

Number of trigger towers alng phi covered by a TCC

Definition at line 153 of file EcalSimRawData.h.

Referenced by genTccIn(), and genTccOut().

nTtSmEta

const int EcalSimRawData::nTtSmEta = 17

staticprivate

Number of TTs along SM eta

Definition at line 105 of file EcalSimRawData.h.

Referenced by genFeData(), genSrData(), genTccIn(), and genTccOut().

nTtSmPhi

const int EcalSimRawData::nTtSmPhi = 4

staticprivate

Number of TTs along SM phi

Definition at line 101 of file EcalSimRawData.h.

Referenced by genFeData(), genSrData(), genTccIn(), and genTccOut().

nTtTypes

const int EcalSimRawData::nTtTypes = 2

staticprivate

Number of barrel trigger tower types (in term of VFE card orientation)

Definition at line 157 of file EcalSimRawData.h.

scEdge

const int EcalSimRawData::scEdge = 5

staticprivate

Supercrystal edge in number of crystals

Definition at line 85 of file EcalSimRawData.h.

srDigiProducer_

std::string EcalSimRawData::srDigiProducer_

private

Label of SR flags and suppressed digis

Definition at line 360 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getSrfs().

srfFile

std::ofstream EcalSimRawData::srfFile

private

Output file for selective readout flags

Definition at line 428 of file EcalSimRawData.h.

srp2dcc_

bool EcalSimRawData::srp2dcc_

private

Switch for data of SRP->DCC link

Definition at line 400 of file EcalSimRawData.h.

Referenced by analyze(), and EcalSimRawData().

strip2Eta

const int EcalSimRawData::strip2Eta

staticprivate

Initial value:

= {
    {4, 3, 2, 1, 0},  
    {0, 1, 2, 3, 4}   
}

Maps strip_id to eta index within a TT. strip2Eta[iTtType][strip_id-1] will be the eta index of channels with strip_id 'strip_id'. The eta index runs from 0 to 4 and increasing with the eta of the std CMS coordinates.

Definition at line 178 of file EcalSimRawData.h.

Referenced by elec2GeomNum().

stripCh2Phi

const int EcalSimRawData::stripCh2Phi

staticprivate

Initial value:

= {
    
    
    {{4, 3, 2, 1, 0},  
     {0, 1, 2, 3, 4},  
     {4, 3, 2, 1, 0},  
     {0, 1, 2, 3, 4},  
     {4, 3, 2, 1, 0}}, 
    
    {{0, 1, 2, 3, 4}, {4, 3, 2, 1, 0}, {0, 1, 2, 3, 4}, {4, 3, 2, 1, 0}, {0, 1, 2, 3, 4}}}

Maps (strip_id, channel_id) to phi index within a TT. stripCh2Phi[iTtType][strip_id-1][ch_id-1] will be the phi index of the TT of type iTtType with strip ID 'strip_id' and channel VFE id 'ch_id'. The phi runs from 0 to 4 and increases with the phi of std CMS coordinates.

Definition at line 171 of file EcalSimRawData.h.

Referenced by elec2GeomNum().

tcc2dcc_

bool EcalSimRawData::tcc2dcc_

private

Switch for data of TCC->DCC link

Definition at line 404 of file EcalSimRawData.h.

Referenced by analyze(), and EcalSimRawData().

tccInDefaultVal_

int EcalSimRawData::tccInDefaultVal_

private

default TriggerPrimitive values if the trigger tower is abscent

Definition at line 440 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getTp().

tccNum_

int EcalSimRawData::tccNum_

private

Index of the TCC, FE data must be produced for. -1 for all TTCs

Definition at line 432 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), genTccIn(), and genTccOut().

tcpDigiCollection_

std::string EcalSimRawData::tcpDigiCollection_

private

TCP Fenix output digi product instance name

Definition at line 376 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getTp().

theGeometry

const CaloGeometry* EcalSimRawData::theGeometry

private

Calorimeter geometry

Definition at line 380 of file EcalSimRawData.h.

theTriggerTowerMap

const EcalTrigTowerConstituentsMap* EcalSimRawData::theTriggerTowerMap

private

ECAL endcap trigger tower map

Definition at line 416 of file EcalSimRawData.h.

tpDigiCollection_

std::string EcalSimRawData::tpDigiCollection_

private

Trigger primitive digi product instance name

Definition at line 372 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getTp().

tpProducer_

std::string EcalSimRawData::tpProducer_

private

Name of the trigger primitive label

Definition at line 384 of file EcalSimRawData.h.

Referenced by EcalSimRawData().

tpVerbose_

bool EcalSimRawData::tpVerbose_

private

Verbosity switch for crystal data

Definition at line 396 of file EcalSimRawData.h.

Referenced by analyze(), EcalSimRawData(), genTccIn(), genTccOut(), and getTp().

trigPrimDigisToken_

edm::EDGetTokenT<EcalTrigPrimDigiCollection> EcalSimRawData::trigPrimDigisToken_[2]

private

Definition at line 449 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), and getTp().

ttEdge

const int EcalSimRawData::ttEdge = 5

staticprivate

Definition at line 97 of file EcalSimRawData.h.

Referenced by elec2GeomNum(), and genFeData().

ttfFile

std::ofstream EcalSimRawData::ttfFile

private

Output file for trigger tower flags

Definition at line 424 of file EcalSimRawData.h.

ttType

const int EcalSimRawData::ttType

staticprivate

Initial value:

= {
    0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1,  
    0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1   
}

Map of trigger tower types (in term of VFE card orientation). ttType[iTtEta0]: trigger tower type of TTs with eta 'c-arrary' index iTtEta0

Definition at line 163 of file EcalSimRawData.h.

Referenced by elec2GeomNum().

writeMode_

writeMode_t EcalSimRawData::writeMode_

private

output format

Definition at line 388 of file EcalSimRawData.h.

Referenced by EcalSimRawData(), fwrite(), genFeData(), genSrData(), and getExt().

xtalVerbose_

bool EcalSimRawData::xtalVerbose_

private

Verbosity switch for crystal data

Definition at line 392 of file EcalSimRawData.h.

Referenced by analyze(), EcalSimRawData(), genFeData(), and getEbDigi().