#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

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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > 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 &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Types
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, const std::string &collName, 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_

std::string	ebSrFlagCollection_

std::string	eeDigiCollection_

std::string	eeSrFlagCollection_

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_

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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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

◆ 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 185 of file EcalSimRawData.h.

185 { littleEndian, bigEndian, ascii };

Constructor & Destructor Documentation

◆ EcalSimRawData()

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

explicit

Constructor

Parameters

pset	CMSSW configuration

Definition at line 42 of file EcalSimRawData.cc.

                                                             {
   // sets up parameters:
   digiProducer_ = params.getParameter<string>("unsuppressedDigiProducer");
   ebDigiCollection_ = params.getParameter<std::string>("EBdigiCollection");
   eeDigiCollection_ = params.getParameter<std::string>("EEdigiCollection");
   srDigiProducer_ = params.getParameter<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<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;
  
   string writeMode = params.getParameter<string>("writeMode");
  
   if (writeMode == string("littleEndian")) {
     writeMode_ = littleEndian;
   } else if (writeMode == string("bigEndian")) {
     writeMode_ = bigEndian;
   } else {
     writeMode_ = ascii;
   }
 }

References iEvent, CalibrationSummaryClient_cfi::params, AlCaHLTBitMon_QueryRunRegistry::string, and EcalSimRawData_cfi::writeMode.

◆ ~EcalSimRawData()

EcalSimRawData::~EcalSimRawData ( )

inlineoverride

Destructor

Definition at line 51 of file EcalSimRawData.h.

51 {};

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 77 of file EcalSimRawData.cc.

                                                                            {
   // Event counter:
   ++iEvent;
  
   if (xtalVerbose_ | tpVerbose_) {
     cout << "=================================================================="
             "====\n"
          << " Event " << iEvent << "\n"
          << "------------------------------------------------------------------"
             "----\n";
   }
  
   if (fe2dcc_) {
     vector<uint16_t> adc[nEbEta][nEbPhi];
     getEbDigi(event, adc);
     genFeData(basename_, iEvent, adc);
   }
  
   if (fe2tcc_) {
     int tcp[nTtEta][nTtPhi] = {{0}};
     getTp(event, tcpDigiCollection_, tcp);
     genTccIn(basename_, iEvent, tcp);
   }
  
   if (tcc2dcc_) {
     int tp[nTtEta][nTtPhi] = {{0}};
     getTp(event, tpDigiCollection_, 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);
   }
 }

References gpuClustering::adc, gather_cfg::cout, iEvent, nEndcaps, and cmsswSequenceInfo::tp.

◆ 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 235 of file EcalSimRawData.h.

235 { 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 241 of file EcalSimRawData.h.

241 { 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 117 of file EcalSimRawData.cc.

                                                                                                              {
   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);
 }

References cms::cuda::assert().

◆ 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 129 of file EcalSimRawData.cc.

                                                                                    {
   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" << setfill('0') << setw(4) << hex << data << "\t" << dec
         << setfill(' ');
       break;
   }
   ++iWord;
 }

References c, data, TauDecayModes::dec, and f.

◆ 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 169 of file EcalSimRawData.cc.

                                                                                                                  {
   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;
  
       stringstream s;
       s.str("");
       const string &ext = getExt();
       s << basename << "_fe2dcc" << setfill('0') << setw(2) << iDcc1 << setfill(' ') << ext;
       ofstream f(s.str().c_str(), (iEvent == 1 ? ios::ate : 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_) {
                 cout << 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" << hex << adcCount[iEta0][iPhi0][5] << dec << "\n";
               }
  
               const 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
 }

References gpuClustering::adc, gather_cfg::cout, data, TauDecayModes::dec, f, iEvent, alignCSCRings::s, and testProducerWithPsetDescEmpty_cfi::x1.

◆ 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 244 of file EcalSimRawData.cc.

                                                                                            {
   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;
       stringstream s;
       s.str("");
       s << basename << "_ab2dcc" << setfill('0') << setw(2) << iDcc1 << setfill(' ') << getExt();
       ofstream f(s.str().c_str(), (iEvent == 1 ? ios::ate : ios::app));
  
       if (f.fail())
         throw cms::Exception(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;
           // cout << "TtEta0: " << iTtEta0 << "\tTtPhi0: " << iTtPhi0 << "\n";
           // cout << "#" << oct << (int)srf[iTtEta0][iTtPhi0] << "o ****> #" <<
           // oct << (srf[iTtEta0][iTtPhi0] << (field*3)) << "o\n" << dec;
  
           data |= srf[iTtEta0][iTtPhi0] << (field * 3);
  
           if (field == 3) {
             // cout <<  srf[iTtEta0][iTtPhi0] << "----> 0x" << hex << data <<
             // "\n";
             fwrite(f, data, iWord, true);
             data = 0;
           }
           ++iFlag;
         }  // next TT along phi
       }    // next TT along eta
     }      // next DCC
   }        // next half-barrel
 }

References data, f, iEvent, alignCSCRings::s, and testProducerWithPsetDescEmpty_cfi::x1.

◆ 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 307 of file EcalSimRawData.cc.

                                                                                               {
   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;
  
       stringstream s;
       s.str("");
       const char *ext = ".txt";  // only ascii mode supported for TCP
  
       s << basename << "_tcc" << setfill('0') << setw(2) << iTcc1 << setfill(' ') << ext;
       ofstream fe2tcc(s.str().c_str(), (iEvent == 1 ? ios::ate : ios::app));
  
       if (fe2tcc.fail())
         throw cms::Exception(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) {
             cout << dec << "iTcc1 = " << iTcc1 << "\t"
                  << "iTtEta0 = " << iTtEta0 << "\t"
                  << "iTtPhi0 = " << iTtPhi0 << "\t"
                  << "iCh1 = " << iCh1 << "\t"
                  << "memPos = " << memPos << "\t"
                  << "tp = 0x" << setfill('0') << hex << setw(3) << tp_fe2tcc << dec << setfill(' ') << "\n";
           }
           fe2tcc << iCh1 << "\t" << memPos << "\t" << setfill('0') << hex << "0x" << setw(4) << tp_fe2tcc << "\t"
                  << "0" << dec << setfill(' ') << "\n";
           ++iCh1;
         }  // next TT along phi
       }    // next TT along eta
       fe2tcc << std::flush;
       fe2tcc.close();
     }  // next TCC
   }    // next half-barrel
 }

References gather_cfg::cout, TauDecayModes::dec, iEvent, and alignCSCRings::s.

◆ 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 356 of file EcalSimRawData.cc.

                                                                                                {
   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;
  
       stringstream s;
       s.str("");
       const char *ext = ".txt";  // only ascii mode supported for TCP
  
       s << basename << "_tcc" << setfill('0') << setw(2) << iTcc1 << setfill(' ') << ext;
  
       s.str("");
       s << basename << "_tcc2dcc" << setfill('0') << setw(2) << iTcc1 << setfill(' ') << getExt();
       ofstream dccF(s.str().c_str(), (iEvent == 1 ? ios::ate : ios::app));
  
       if (dccF.fail()) {
         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_) {
             cout << dec << "iTcc1 = " << iTcc1 << "\t"
                  << "iTtEta0 = " << iTtEta0 << "\t"
                  << "iTtPhi0 = " << iTtPhi0 << "\t"
                  << "iCh1 = " << iCh1 << "\t"
                  << "memPos = " << memPos << "\t"
                  << "tp = 0x" << hex << tps[iTtEta0][iTtPhi0] << dec << "\n";
           }
           fwrite(dccF, tps[iTtEta0][iTtPhi0], iDccWord, false);
           ++iCh1;
         }  // next TT along phi
       }    // next TT along eta
     }      // next TCC
   }        // next half-barrel
 }

References gather_cfg::cout, data, TauDecayModes::dec, iEvent, PresampleTask_cfi::nSamples, alignCSCRings::s, and testProducerWithPsetDescEmpty_cfi::x1.

◆ 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 478 of file EcalSimRawData.cc.

                                                                                                 {
   edm::Handle<EBDigiCollection> hEbDigis;
   event.getByLabel(digiProducer_, ebDigiCollection_, hEbDigis);
  
   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] = 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_)
       cout << 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());
  
       //     cout << "xtl indices conv: (" << frame.id().ieta() << ","
       //     << frame.id().iphi() << ") -> ("
       //     << iEta0 << "," << iPhi0 << ")\n";
  
       if (iEta0 < 0 || iEta0 >= nEbEta) {
         cout << "iEta0 (= " << iEta0 << ") is out of range ("
              << "[0," << nEbEta - 1 << "])\n";
       }
       if (iPhi0 < 0 || iPhi0 >= nEbPhi) {
         cout << "iPhi0 (= " << iPhi0 << ") is out of range ("
              << "[0," << nEbPhi - 1 << "])\n";
       }
  
       if (xtalVerbose_) {
         cout << iEta0 << "\t" << iPhi0 << ":\t";
         cout << 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_) {
           cout << (iSample > 0 ? " " : "") << "0x" << setw(4) << encodedAdc;
         }
       }
       if (xtalVerbose_)
         cout << "\n" << dec;
     }
     if (xtalVerbose_)
       cout << setfill(' ');
   }
 }

References gpuClustering::adc, edm::DataFrameContainer::begin(), gather_cfg::cout, TauDecayModes::dec, edm::DataFrameContainer::empty(), edm::DataFrameContainer::end(), Exception, amptDefault_cfi::frame, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, edm::HandleBase::isValid(), PresampleTask_cfi::nSamples, and simplePhotonAnalyzer_cfi::sample.

◆ getExt()

string EcalSimRawData::getExt ( ) const

private

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

Definition at line 156 of file EcalSimRawData.cc.

                                     {
   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 428 of file EcalSimRawData.cc.

                                                                     {
   // EE
   edm::Handle<EESrFlagCollection> hEeSrFlags;
   event.getByLabel(srDigiProducer_, eeSrFlagCollection_, hEeSrFlags);
   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 {
     LogWarning("EcalSimRawData") << "EE SR flag not found ("
                                  << "Product label: " << srDigiProducer_ << "Producet instance: " << eeSrFlagCollection_
                                  << ")";
   }
  
   // EB
   edm::Handle<EBSrFlagCollection> hEbSrFlags;
   event.getByLabel(srDigiProducer_, ebSrFlagCollection_, hEbSrFlags);
   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 {
     LogWarning("EcalSimRawData") << "EB SR flag not found ("
                                  << "Product label: " << srDigiProducer_ << "Producet instance: " << ebSrFlagCollection_
                                  << ")";
   }
 }

References cms::cuda::assert(), edm::SortedCollection< T, SORT >::begin(), edm::SortedCollection< T, SORT >::end(), RemoveAddSevLevel::flag, mps_fire::i, L1TowerCalibrationProducer_cfi::iEta, edm::HandleBase::isValid(), and nEndcaps.

◆ getTp()

void EcalSimRawData::getTp	(	const edm::Event &	event,
		const std::string &	collName,
		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 545 of file EcalSimRawData.cc.

                                                                                                           {
   edm::Handle<EcalTrigPrimDigiCollection> hTpDigis;
   event.getByLabel(tpProducer_, collName, hTpDigis);
   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_) {
       cout << 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) {
         cout << "iTtEta0 (= " << iTtEta0 << ") is out of range ("
              << "[0," << nEbTtEta - 1 << "])\n";
       }
       if (iTtPhi0 < 0 || iTtPhi0 >= nTtPhi) {
         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
           cout << collName << (collName.empty() ? "" : " ") << "TP(" << setw(2) << iTtEta0 << "," << iTtPhi0 << ") = "
                << "0x" << setw(4) << tcp[iTtEta0][iTtPhi0] << "\tcmssw indices: " << tp.id().ieta() << " "
                << tp.id().iphi() << "\n";
       }
     }  // next TP
     if (tpVerbose_)
       cout << setfill(' ');
   }
 }

References edm::SortedCollection< T, SORT >::begin(), trigObjTnPSource_cfi::collName, gather_cfg::cout, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), edm::HandleBase::isValid(), edm::Handle< T >::product(), and cmsswSequenceInfo::tp.

◆ 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 209 of file EcalSimRawData.h.

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

References L1TowerCalibrationProducer_cfi::iEta.

◆ 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 216 of file EcalSimRawData.h.

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

References nEbPhi.

◆ 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 229 of file EcalSimRawData.h.

229 { 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 248 of file EcalSimRawData.h.

248 { return iTtPhi - 1; }

◆ setHParity()

void EcalSimRawData::setHParity ( uint16_t & a ) const

private

Definition at line 417 of file EcalSimRawData.cc.

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

References a, and AlCaHLTBitMon_ParallelJobs::p.

Member Data Documentation

◆ basename_

std::string EcalSimRawData::basename_

private

basename for output files

Definition at line 442 of file EcalSimRawData.h.

◆ dccNum_

int EcalSimRawData::dccNum_

private

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

Definition at line 434 of file EcalSimRawData.h.

◆ digiProducer_

std::string EcalSimRawData::digiProducer_

private

Name of module/plugin/producer making digis

Definition at line 346 of file EcalSimRawData.h.

◆ ebDccPhiEdge

const int EcalSimRawData::ebDccPhiEdge = 20

staticprivate

Number of barrel crystals along phi covered by a DCC

Definition at line 143 of file EcalSimRawData.h.

◆ ebDigiCollection_

std::string EcalSimRawData::ebDigiCollection_

private

EB digi product instance name

Definition at line 350 of file EcalSimRawData.h.

◆ ebSrFlagCollection_

std::string EcalSimRawData::ebSrFlagCollection_

private

EB SRP flag digi product instance name

Definition at line 362 of file EcalSimRawData.h.

◆ ebTccPhiEdge

const int EcalSimRawData::ebTccPhiEdge = 20

staticprivate

number of TTs along phi of a TCC sector

Definition at line 131 of file EcalSimRawData.h.

◆ eeDigiCollection_

std::string EcalSimRawData::eeDigiCollection_

private

EE digi product instance name

Definition at line 354 of file EcalSimRawData.h.

◆ eeSrFlagCollection_

std::string EcalSimRawData::eeSrFlagCollection_

private

EE SRP flag digi product instancename

Definition at line 366 of file EcalSimRawData.h.

◆ esr_

std::unique_ptr<EcalSelectiveReadout> EcalSimRawData::esr_

private

Selective readout simulator

Definition at line 418 of file EcalSimRawData.h.

◆ fe2dcc_

bool EcalSimRawData::fe2dcc_

private

Switch for data of FE->DCC link

Definition at line 406 of file EcalSimRawData.h.

◆ fe2tcc_

bool EcalSimRawData::fe2tcc_

private

Switch for data of FE->TCC link

Definition at line 410 of file EcalSimRawData.h.

◆ iEvent

int EcalSimRawData::iEvent

private

Definition at line 59 of file EcalSimRawData.h.

◆ nDccEndcap

const int EcalSimRawData::nDccEndcap = 9

staticprivate

Number of DCCs for a single endcap

Definition at line 127 of file EcalSimRawData.h.

◆ nDccInPhi

const int EcalSimRawData::nDccInPhi = 18

staticprivate

Number of barrel DCCs along Phi

Definition at line 123 of file EcalSimRawData.h.

◆ nEbEta

const int EcalSimRawData::nEbEta = 170

staticprivate

Number of crystals in ECAL barrel along eta

Definition at line 63 of file EcalSimRawData.h.

◆ nEbPhi

const int EcalSimRawData::nEbPhi = 360

staticprivate

Number of crystals in ECAL barrel along phi

Definition at line 67 of file EcalSimRawData.h.

Referenced by iPhi2cIndex().

◆ nEbTtEta

const int EcalSimRawData::nEbTtEta = nEbEta / ttEdge

staticprivate

Number of TTs along Ecal barrel eta

Definition at line 111 of file EcalSimRawData.h.

◆ nEeTtEta

const int EcalSimRawData::nEeTtEta = 11

staticprivate

Number of TTs along eta for one endcap.

Definition at line 115 of file EcalSimRawData.h.

◆ nEeX

const int EcalSimRawData::nEeX = 100

staticprivate

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

Definition at line 71 of file EcalSimRawData.h.

◆ nEeY

const int EcalSimRawData::nEeY = 100

staticprivate

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

Definition at line 75 of file EcalSimRawData.h.

◆ nEndcaps

const int EcalSimRawData::nEndcaps = 2

staticprivate

Number of endcaps

Definition at line 79 of file EcalSimRawData.h.

◆ nScX

const int EcalSimRawData::nScX = 20

staticprivate

Maximum number of supercrystal along x axis

Definition at line 87 of file EcalSimRawData.h.

◆ nScY

const int EcalSimRawData::nScY = 20

staticprivate

Maximum number of supercrystal along y axis

Definition at line 91 of file EcalSimRawData.h.

◆ nTccEndcap

const int EcalSimRawData::nTccEndcap = 36

staticprivate

Number of TCCs for a single endcap

Definition at line 139 of file EcalSimRawData.h.

◆ nTccInPhi

const int EcalSimRawData::nTccInPhi = 18

staticprivate

Number of Barrel TTCs along phi

Definition at line 135 of file EcalSimRawData.h.

◆ nTtEta

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

staticprivate

Number of TTs along ECAL eta

Definition at line 119 of file EcalSimRawData.h.

◆ nTtPhi

const int EcalSimRawData::nTtPhi = nEbPhi / ttEdge

staticprivate

Number of TTs along Ecal Phi

Definition at line 107 of file EcalSimRawData.h.

◆ nTtPhisPerEbDcc

const int EcalSimRawData::nTtPhisPerEbDcc = 4

staticprivate

Number of trigger towers alng phi covered by a DCC

Definition at line 147 of file EcalSimRawData.h.

◆ nTtPhisPerEbTcc

const int EcalSimRawData::nTtPhisPerEbTcc = 4

staticprivate

Number of trigger towers alng phi covered by a TCC

Definition at line 151 of file EcalSimRawData.h.

◆ nTtSmEta

const int EcalSimRawData::nTtSmEta = 17

staticprivate

Number of TTs along SM eta

Definition at line 103 of file EcalSimRawData.h.

◆ nTtSmPhi

const int EcalSimRawData::nTtSmPhi = 4

staticprivate

Number of TTs along SM phi

Definition at line 99 of file EcalSimRawData.h.

◆ nTtTypes

const int EcalSimRawData::nTtTypes = 2

staticprivate

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

Definition at line 155 of file EcalSimRawData.h.

◆ scEdge

const int EcalSimRawData::scEdge = 5

staticprivate

Supercrystal edge in number of crystals

Definition at line 83 of file EcalSimRawData.h.

◆ srDigiProducer_

std::string EcalSimRawData::srDigiProducer_

private

Label of SR flags and suppressed digis

Definition at line 358 of file EcalSimRawData.h.

◆ srfFile

std::ofstream EcalSimRawData::srfFile

private

Output file for selective readout flags

Definition at line 426 of file EcalSimRawData.h.

◆ srp2dcc_

bool EcalSimRawData::srp2dcc_

private

Switch for data of SRP->DCC link

Definition at line 398 of file EcalSimRawData.h.

◆ 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 176 of file EcalSimRawData.h.

◆ 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 169 of file EcalSimRawData.h.

◆ tcc2dcc_

bool EcalSimRawData::tcc2dcc_

private

Switch for data of TCC->DCC link

Definition at line 402 of file EcalSimRawData.h.

◆ tccInDefaultVal_

int EcalSimRawData::tccInDefaultVal_

private

default TriggerPrimitive values if the trigger tower is abscent

Definition at line 438 of file EcalSimRawData.h.

◆ tccNum_

int EcalSimRawData::tccNum_

private

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

Definition at line 430 of file EcalSimRawData.h.

◆ tcpDigiCollection_

std::string EcalSimRawData::tcpDigiCollection_

private

TCP Fenix output digi product instance name

Definition at line 374 of file EcalSimRawData.h.

◆ theGeometry

const CaloGeometry* EcalSimRawData::theGeometry

private

Calorimeter geometry

Definition at line 378 of file EcalSimRawData.h.

◆ theTriggerTowerMap

const EcalTrigTowerConstituentsMap* EcalSimRawData::theTriggerTowerMap

private

ECAL endcap trigger tower map

Definition at line 414 of file EcalSimRawData.h.

◆ tpDigiCollection_

std::string EcalSimRawData::tpDigiCollection_

private

Trigger primitive digi product instance name

Definition at line 370 of file EcalSimRawData.h.

◆ tpProducer_

std::string EcalSimRawData::tpProducer_

private

Name of the trigger primitive label

Definition at line 382 of file EcalSimRawData.h.

◆ tpVerbose_

bool EcalSimRawData::tpVerbose_

private

Verbosity switch for crystal data

Definition at line 394 of file EcalSimRawData.h.

◆ ttEdge

const int EcalSimRawData::ttEdge = 5

staticprivate

Definition at line 95 of file EcalSimRawData.h.

◆ ttfFile

std::ofstream EcalSimRawData::ttfFile

private

Output file for trigger tower flags

Definition at line 422 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 161 of file EcalSimRawData.h.

◆ writeMode_

writeMode_t EcalSimRawData::writeMode_

private

output format

Definition at line 386 of file EcalSimRawData.h.

◆ xtalVerbose_

bool EcalSimRawData::xtalVerbose_

private

Verbosity switch for crystal data

Definition at line 390 of file EcalSimRawData.h.

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ writeMode_t

Constructor & Destructor Documentation

◆ EcalSimRawData()

◆ ~EcalSimRawData()

Member Function Documentation

◆ analyze()

◆ checkGeometry()

◆ checkTriggerMap()

◆ cIndex2iTtEta()

◆ cIndex2TtPhi()

◆ elec2GeomNum()

◆ fwrite()

◆ genFeData()

◆ genSrData()

◆ genTccIn()

◆ genTccOut()

◆ getEbDigi()

◆ getExt()

◆ getSrfs()

◆ getTp()

◆ iEta2cIndex()

◆ iPhi2cIndex()

◆ iTtEta2cIndex()

◆ iTtPhi2cIndex()

◆ setHParity()

Member Data Documentation

◆ basename_

◆ dccNum_

◆ digiProducer_

◆ ebDccPhiEdge

◆ ebDigiCollection_

◆ ebSrFlagCollection_

◆ ebTccPhiEdge

◆ eeDigiCollection_

◆ eeSrFlagCollection_

◆ esr_

◆ fe2dcc_

◆ fe2tcc_

◆ iEvent

◆ nDccEndcap

◆ nDccInPhi

◆ nEbEta

◆ nEbPhi

◆ nEbTtEta

◆ nEeTtEta

◆ nEeX

◆ nEeY

◆ nEndcaps

◆ nScX

◆ nScY

◆ nTccEndcap

◆ nTccInPhi

◆ nTtEta

◆ nTtPhi

◆ nTtPhisPerEbDcc

◆ nTtPhisPerEbTcc

◆ nTtSmEta

◆ nTtSmPhi

◆ nTtTypes

◆ scEdge

◆ srDigiProducer_

◆ srfFile

◆ srp2dcc_

◆ strip2Eta

◆ stripCh2Phi

◆ tcc2dcc_

◆ tccInDefaultVal_

◆ tccNum_

◆ tcpDigiCollection_

◆ theGeometry

◆ theTriggerTowerMap

◆ tpDigiCollection_