#include <DTTFFEDSim.h>

Inheritance diagram for DTTFFEDSim:

Public Member Functions
	DTTFFEDSim (const edm::ParameterSet &pset)
	Constructor. More...

bool	fillRawData (edm::Event &e, FEDRawDataCollection &data)
	Generate and fill FED raw data for a full event. More...

void	produce (edm::Event &e, const edm::EventSetup &c)
	Produce digis out of raw data. More...

virtual	~DTTFFEDSim ()
	Destructor. More...

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
int	bxNr (int channel)

void	calcCRC (int myD1, int myD2, int &myC)

int	channel (int wheel, int sector, int bx)

edm::InputTag	getDTDigiInputTag ()

edm::InputTag	getDTPHTFInputTag ()

int	sector (int channel)

int	wheel (int channel)

Private Attributes
edm::EDGetTokenT < L1MuDTChambPhContainer >	ChPh_tok

edm::EDGetTokenT < L1MuDTChambThContainer >	ChTh_tok

edm::InputTag	DTDigiInputTag

edm::InputTag	DTPHTFInputTag

unsigned int	eventNum

edm::EDGetTokenT < L1MuDTTrackContainer >	Trk_tok

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

void	consumesMany (const TypeToGet &id)

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

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)

Detailed Description

L1 DT Track Finder Digi-to-Raw

J. Troconiz UAM Madrid

Definition at line 29 of file DTTFFEDSim.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 25 of file DTTFFEDSim.cc.

References ChPh_tok, ChTh_tok, DTDigiInputTag, DTPHTFInputTag, edm::ParameterSet::getParameter(), and Trk_tok.

                                                   : eventNum(0) {
 
   produces<FEDRawDataCollection>();
 
   DTDigiInputTag = pset.getParameter<edm::InputTag>("DTDigi_Source");
   DTPHTFInputTag = pset.getParameter<edm::InputTag>("DTTracks_Source");
 
   ChPh_tok = consumes<L1MuDTChambPhContainer>(DTDigiInputTag);
   ChTh_tok = consumes<L1MuDTChambThContainer>(DTDigiInputTag);
   Trk_tok = consumes<L1MuDTTrackContainer>(DTPHTFInputTag);
 
 }

DTTFFEDSim::~DTTFFEDSim ( )

virtual

Destructor.

Definition at line 38 of file DTTFFEDSim.cc.

38 {}

Member Function Documentation

int DTTFFEDSim::bxNr ( int channel )

private

Definition at line 291 of file DTTFFEDSim.cc.

References channel().

                                  {
 
   int myChannel = channel;
 
   if (myChannel > 127) myChannel -= 2;
 
   if (myChannel < 0 || myChannel > 251 ){ return -999; }
 
   int myBx = 1-(myChannel%3);
 
   return myBx;
 }

void DTTFFEDSim::calcCRC	(	int	myD1,
		int	myD2,
		int &	myC
	)

private

Definition at line 328 of file DTTFFEDSim.cc.

References funct::C, and i.

Referenced by fillRawData().

                                                     {
 
   int myCRC[16],D[64],C[16];
 
   for( int i=0; i < 32; i++ ){ D[i]=(myD2>>i)&0x1; }
   for( int i=0; i < 32; i++ ){ D[i+32]=(myD1>>i)&0x1; }
   for( int i=0; i < 16; i++ ){ C[i]=(myC>>i)&0x1; }
 
   myCRC[0] = ( D[63] + D[62] + D[61] + D[60] + D[55] + D[54] +
                D[53] + D[52] + D[51] + D[50] + D[49] + D[48] +
                D[47] + D[46] + D[45] + D[43] + D[41] + D[40] +
                D[39] + D[38] + D[37] + D[36] + D[35] + D[34] +
                D[33] + D[32] + D[31] + D[30] + D[27] + D[26] +
                D[25] + D[24] + D[23] + D[22] + D[21] + D[20] +
                D[19] + D[18] + D[17] + D[16] + D[15] + D[13] +
                D[12] + D[11] + D[10] + D[9]  + D[8]  + D[7]  +
                D[6]  + D[5]  + D[4]  + D[3]  + D[2]  + D[1]  +
                D[0]  + C[0]  + C[1]  + C[2]  + C[3]  + C[4]  +
                C[5]  + C[6]  + C[7]  + C[12] + C[13] + C[14] +
                C[15] )%2;
 
   myCRC[1] = ( D[63] + D[62] + D[61] + D[56] + D[55] + D[54] +
                D[53] + D[52] + D[51] + D[50] + D[49] + D[48] +
                D[47] + D[46] + D[44] + D[42] + D[41] + D[40] +
                D[39] + D[38] + D[37] + D[36] + D[35] + D[34] +
                D[33] + D[32] + D[31] + D[28] + D[27] + D[26] +
            D[25] + D[24] + D[23] + D[22] + D[21] + D[20] +
            D[19] + D[18] + D[17] + D[16] + D[14] + D[13] +
            D[12] + D[11] + D[10] + D[9]  + D[8]  + D[7]  +
            D[6]  + D[5]  + D[4]  + D[3]  + D[2]  + D[1]  +
            C[0]  + C[1]  + C[2]  + C[3]  + C[4]  + C[5]  +
            C[6]  + C[7]  + C[8]  + C[13] + C[14] + C[15] )%2;
 
   myCRC[2] = ( D[61] + D[60] + D[57] + D[56] + D[46] + D[42] +
            D[31] + D[30] + D[29] + D[28] + D[16] + D[14] +
            D[1]  + D[0]  + C[8]  + C[9]  + C[12] + C[13] )%2;
 
   myCRC[3] = ( D[62] + D[61] + D[58] + D[57] + D[47] + D[43] +
            D[32] + D[31] + D[30] + D[29] + D[17] + D[15] +
            D[2]  + D[1]  + C[9]  + C[10] + C[13] + C[14] )%2;
 
   myCRC[4] = ( D[63] + D[62] + D[59] + D[58] + D[48] + D[44] +
            D[33] + D[32] + D[31] + D[30] + D[18] + D[16] + 
            D[3]  + D[2]  + C[0]  + C[10] + C[11] + C[14] +
            C[15] )%2;
 
   myCRC[5] = ( D[63] + D[60] + D[59] + D[49] + D[45] + D[34] +
            D[33] + D[32] + D[31] + D[19] + D[17] + D[4]  +
            D[3]  + C[1]  + C[11] + C[12] + C[15] )%2;
 
   myCRC[6] = ( D[61] + D[60] + D[50] + D[46] + D[35] + D[34] +
            D[33] + D[32] + D[20] + D[18] + D[5]  + D[4]  +
            C[2]  + C[12] + C[13] )%2;
 
   myCRC[7] = ( D[62] + D[61] + D[51] + D[47] + D[36] + D[35] +
            D[34] + D[33] + D[21] + D[19] + D[6]  + D[5]  +
            C[3]  + C[13] + C[14] )%2;
 
   myCRC[8] = ( D[63] + D[62] + D[52] + D[48] + D[37] + D[36] +
            D[35] + D[34] + D[22] + D[20] + D[7]  + D[6]  +
            C[0]  + C[4]  + C[14] + C[15] )%2;
 
   myCRC[9] = ( D[63] + D[53] + D[49] + D[38] + D[37] + D[36] +
            D[35] + D[23] + D[21] + D[8]  + D[7]  + C[1]  +
            C[5]  + C[15] )%2;
 
   myCRC[10] = ( D[54] + D[50] + D[39] + D[38] + D[37] + D[36] + 
         D[24] + D[22] + D[9]  + D[8]  + C[2]  + C[6] )%2;
 
   myCRC[11] = ( D[55] + D[51] + D[40] + D[39] + D[38] + D[37] +
         D[25] + D[23] + D[10] + D[9]  + C[3]  + C[7] )%2;
 
   myCRC[12] = ( D[56] + D[52] + D[41] + D[40] + D[39] + D[38] +
         D[26] + D[24] + D[11] + D[10] + C[4]  + C[8] )%2;
 
   myCRC[13] = ( D[57] + D[53] + D[42] + D[41] + D[40] + D[39] +
         D[27] + D[25] + D[12] + D[11] + C[5]  + C[9] )%2;
 
   myCRC[14] = ( D[58] + D[54] + D[43] + D[42] + D[41] + D[40] +
         D[28] + D[26] + D[13] + D[12] + C[6]  + C[10] )%2;
 
   myCRC[15] = ( D[63] + D[62] + D[61] + D[60] + D[59] + D[54] +
         D[53] + D[52] + D[51] + D[50] + D[49] + D[48] + 
         D[47] + D[46] + D[45] + D[44] + D[42] + D[40] +
         D[39] + D[38] + D[37] + D[36] + D[35] + D[34] + 
         D[33] + D[32] + D[31] + D[30] + D[29] + D[26] +
         D[25] + D[24] + D[23] + D[22] + D[21] + D[20] + 
         D[19] + D[18] + D[17] + D[16] + D[15] + D[14] +
         D[12] + D[11] + D[10] + D[9]  + D[8]  + D[7]  + 
         D[6]  + D[5]  + D[4]  + D[3]  + D[2]  + D[1]  +
         D[0]  + C[0]  + C[1]  + C[2]  + C[3]  + C[4]  + 
         C[5]  + C[6]  + C[11] + C[12] + C[13] + C[14] +
         C[15] )%2;
 
   int tempC = 0x0;  
   for(int i=0; i<16 ; i++){ tempC = tempC + (myCRC[i]<<i); }
   myC = tempC ;
   return;
 }

int DTTFFEDSim::channel	(	int	wheel,
		int	sector,
		int	bx
	)

private

Definition at line 271 of file DTTFFEDSim.cc.

References funct::abs().

Referenced by bxNr(), fillRawData(), sector(), and wheel().

                                                        {
 
   // wheel  :  -3 -2 -1 +1 +2 +3 <=> PHTF's : N2, N1, N0, P0, P1, P2
   //                           0 <=> ETTF  
   // sector :  0 -> 11
   // bx     : -1 -> +1
 
   int myChannel = 255;
 
   if ( abs(bx) > 1)               { return myChannel; }
   if ( sector < 0 || sector > 11) { return myChannel; }
   if ( abs(wheel) > 3)            { return myChannel; }
 
   myChannel = sector*21 + wheel*3 - bx + 10 ; 
 
   if (myChannel > 125) myChannel += 2;
 
   return myChannel;
 }

bool DTTFFEDSim::fillRawData	(	edm::Event &	e,
		FEDRawDataCollection &	data
	)

Generate and fill FED raw data for a full event.

Definition at line 52 of file DTTFFEDSim.cc.

References calcCRC(), channel(), ChPh_tok, ChTh_tok, FEDRawData::data(), edm::EventID::event(), eventNum, FEDRawDataCollection::FEDData(), edm::Event::getByToken(), edm::EventBase::id(), groupFilesInBlocks::lines, FEDRawData::resize(), and Trk_tok.

Referenced by produce().

                                                          {
 
   eventNum = e.id().event();
 
   int lines = 2;
 
   edm::Handle<L1MuDTChambPhContainer> phtrig;
   e.getByToken(ChPh_tok,phtrig);
   lines += phtrig->bxSize(-1, 1);
 
   edm::Handle<L1MuDTChambThContainer> thtrig;
   e.getByToken(ChTh_tok,thtrig);
   lines += thtrig->bxSize(-1, 1);
 
   edm::Handle<L1MuDTTrackContainer>   trtrig;
   e.getByToken(Trk_tok,trtrig);
   lines += trtrig->bxSize(-1, 1)*3;
 
   FEDRawData& dttfdata = data.FEDData(0x30C);
   dttfdata.resize(lines*8); // size in bytes
   unsigned char* LineFED=dttfdata.data();
 
   int* dataWord1 = new int;
   int* dataWord2 = new int;
 
   //--> Header
 
   *dataWord1 = 0x50000000
              + (eventNum&0xFFFFFF);
   *dataWord2 = 0x00030C00;
 
   int newCRC =  0xFFFF;
   calcCRC(*dataWord1, *dataWord2, newCRC);
   *((int*)LineFED)=*dataWord2; 
   LineFED+=4;
   *((int*)LineFED)=*dataWord1; 
 
   //--> DTTF data 
 
   int TS1Id[4], TS2Id[4];   // word identifier for TS #1,#2 for stations 
   TS1Id[0] = 0x0E;      
   TS2Id[0] = 0x1E;
   TS1Id[1] = 0x2E;
   TS2Id[1] = 0x3E;
   TS1Id[3] = 0x4E;
   TS2Id[3] = 0x5E;
   TS1Id[2] = 0x8FFF8;
   TS2Id[2] = 0x9FFF8;
 
   //Input
   L1MuDTChambPhContainer::Phi_iterator tsphi;
 
   for ( tsphi =  phtrig->getContainer()->begin();
         tsphi != phtrig->getContainer()->end();
         tsphi++ ) {
     if ( tsphi->code() != 7 ) {
 
       int wheelID   = tsphi->whNum()+1;
       if ( wheelID <= 0 ) wheelID -= 2;
       int stationID = tsphi->stNum()-1;
       int is2nd     = tsphi->Ts2Tag();
 
       int channelNr = channel(wheelID, tsphi->scNum(), tsphi->bxNum()-is2nd);
       if ( channelNr == 255 ) continue;
       int TSId = ( is2nd == 0 ) ? TS1Id[stationID] : TS2Id[stationID];
 
       *dataWord1 = ((channelNr&0xFF)<<24)
                  + 0x00FFFFFF;
 
       if ( stationID != 2 ){
         *dataWord2 = (             (TSId&0x0FF)<<24)
                + (~(tsphi->code()+1)&0x007)
                + (   (~tsphi->phiB()&0x3FF)<<3) 
                + (    (~tsphi->phi()&0xFFF)<<13);
       }
       else {
         *dataWord2 = (             (TSId&0xFFFFF)<<12) 
                + (~(tsphi->code()+1)&0x00007)
                + (    (~tsphi->phi()&0x00FFF)<<3);
       }
 
       calcCRC(*dataWord1, *dataWord2, newCRC);
       LineFED+=4;
       *((int*)LineFED)=*dataWord2; 
       LineFED+=4;
       *((int*)LineFED)=*dataWord1; 
     }
   }
   //Input
 
   //Input
   L1MuDTChambThContainer::The_iterator tsthe;
 
   for ( tsthe =  thtrig->getContainer()->begin();
         tsthe != thtrig->getContainer()->end();
         tsthe++ ) {
 
     int wheelTh  = tsthe->whNum();
     int sectorID = tsthe->scNum();
 
     int channelNr = channel(0, sectorID, tsthe->bxNum());
     if ( channelNr == 255 ) continue;
     int TSId = wheelTh+2;
 
     *dataWord1 = ((channelNr&0xFF)<<24)
                + 0x00FFFFFF;
 
     *dataWord2 = ((TSId&0x07)<<28)
                + 0x0FFFFFFF;
 
     int stationID = tsthe->stNum()-1;
     for ( int bti = 0; bti < 7; bti++ )
       if ( wheelTh == -2 || wheelTh == -1 || 
        ( wheelTh == 0 && (sectorID == 0 || sectorID == 3 || sectorID == 4 || sectorID == 7 || sectorID == 8 || sectorID == 11) ) )
     *dataWord2 -= (tsthe->position(bti)&0x1)<<(stationID*7+bti);
       else
     *dataWord2 -= (tsthe->position(6-bti)&0x1)<<(stationID*7+bti);
 
     calcCRC(*dataWord1, *dataWord2, newCRC);
     LineFED+=4;
     *((int*)LineFED)=*dataWord2; 
     LineFED+=4;
     *((int*)LineFED)=*dataWord1; 
   }
   //Input
 
   //Output
   L1MuDTTrackContainer::Trackiterator tstrk;
 
   for ( tstrk =  trtrig->getContainer()->begin();
         tstrk != trtrig->getContainer()->end();
         tstrk++ ) {
 
     int channelNr = channel(tstrk->whNum(), tstrk->scNum(), tstrk->bx());
     if ( channelNr == 255 ) continue;
     int TSId = ( tstrk->TrkTag() == 0 ) ? 0xAFFF : 0xBFFF;
 
     *dataWord1 = ((channelNr&0xFF)<<24)
                + 0x00FFFFFF;
 
     *dataWord2 = (           (TSId&0xFFFF)<<16)
            + ( tstrk->stNum(4)&0x0000F)
                + ((tstrk->stNum(3)&0x0000F)<<4)
                + ((tstrk->stNum(2)&0x0000F)<<8)
            + ((tstrk->stNum(1)&0x00003)<<12);
 
     calcCRC(*dataWord1, *dataWord2, newCRC);
     LineFED+=4;
     *((int*)LineFED)=*dataWord2; 
     LineFED+=4;
     *((int*)LineFED)=*dataWord1; 
 
     TSId = ( tstrk->TrkTag() == 0 ) ? 0xCFFE : 0xDFFE;
 
     *dataWord1 = ((channelNr&0xFF)<<24)
                + 0x00FFFFFF;
 
     *dataWord2 = (                    (TSId&0xFFFE)<<16)
            + ( ~tstrk->quality_packed()&0x0007)
                + (     (tstrk->phi_packed()&0x00FF)<<3)
                + ( (~tstrk->charge_packed()&0x0001)<<11)
            + (     (~tstrk->pt_packed()&0x001F)<<12);
 
     calcCRC(*dataWord1, *dataWord2, newCRC);
     LineFED+=4;
     *((int*)LineFED)=*dataWord2; 
     LineFED+=4;
     *((int*)LineFED)=*dataWord1; 
 
     channelNr = channel(0, tstrk->scNum(), tstrk->bx());
     if ( channelNr == 255 ) continue;
     TSId = (tstrk->whNum()+3)<<16;
     TSId += ( tstrk->whNum() < 0 ) ? 0x8FFFC : 0x7FFFC;
 
     *dataWord1 = ((channelNr&0xFF)<<24)
                + 0x00FFFFFF;
 
     *dataWord2 = (TSId&0xFFFFC)<<12;
 
     if ( tstrk->TrkTag() == 0 ) {
       *dataWord2 +=                            0x3F80
                  + (       tstrk->eta_packed()&0x003F)
                  + ((~tstrk->finehalo_packed()&0x0001)<<6);
     }
     else {
       *dataWord2 +=                            0x007F
                  + (     ( tstrk->eta_packed()&0x003F)<<7)
                  + ((~tstrk->finehalo_packed()&0x0001)<<13);
     }
 
     calcCRC(*dataWord1, *dataWord2, newCRC);
     LineFED+=4;
     *((int*)LineFED)=*dataWord2; 
     LineFED+=4;
     *((int*)LineFED)=*dataWord1; 
   }
   //Output
 
   //--> Trailer
 
   *dataWord1 = 0xA0000000
              + (lines&0xFFFFFF);
   *dataWord2 = 0;
 
   calcCRC(*dataWord1, *dataWord2&0xFFFF, newCRC);
 
   *dataWord2 += (newCRC&0xFFFF)<<16;
 
   LineFED+=4;
   *((int*)LineFED)=*dataWord2; 
   LineFED+=4;
   *((int*)LineFED)=*dataWord1; 
 
   delete dataWord1;
   delete dataWord2;
   return true;
 }

edm::InputTag DTTFFEDSim::getDTDigiInputTag ( )

inlineprivate

Definition at line 64 of file DTTFFEDSim.h.

References DTDigiInputTag.

64 { return DTDigiInputTag; }

DTTFFEDSim::DTDigiInputTag

edm::InputTag DTDigiInputTag

Definition: DTTFFEDSim.h:50

edm::InputTag DTTFFEDSim::getDTPHTFInputTag ( )

inlineprivate

Definition at line 65 of file DTTFFEDSim.h.

References DTPHTFInputTag.

65 { return DTPHTFInputTag; }

DTTFFEDSim::DTPHTFInputTag

edm::InputTag DTPHTFInputTag

Definition: DTTFFEDSim.h:51

void DTTFFEDSim::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

virtual

Produce digis out of raw data.

Implements edm::EDProducer.

Definition at line 40 of file DTTFFEDSim.cc.

References data, fillRawData(), and edm::Event::put().

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                             {
 
   FEDRawDataCollection data;
 
   if (!fillRawData(e, data)) return;
 
   auto_ptr<FEDRawDataCollection> fed_product(new FEDRawDataCollection(data));
 
   e.put(fed_product);
 
 }

int DTTFFEDSim::sector ( int channel )

private

Definition at line 304 of file DTTFFEDSim.cc.

References channel().

Referenced by geometryXMLparser.DTAlignable::index().

                                    {
 
   int myChannel = channel;
 
   if (myChannel > 127) myChannel -= 2;
 
   if (myChannel < 0 || myChannel > 251 ){ return -999; }
 
   return myChannel/21;
 }

int DTTFFEDSim::wheel ( int channel )

private

Definition at line 315 of file DTTFFEDSim.cc.

References channel().

Referenced by geometryXMLparser.DTAlignable::index().

                                   {
 
   int myChannel = channel;
 
   if (myChannel > 127) myChannel -= 2;
 
   if (myChannel < 0 || myChannel > 251 ){ return -999; }
 
   int myWheel = ((myChannel%21)/3)-3;
 
   return myWheel;
 }