|
|
|
#include <DTTFFEDSim.h>
Public Member Functions | |
| DTTFFEDSim (const edm::ParameterSet &pset) | |
| Constructor. | |
| bool | fillRawData (edm::Event &e, FEDRawDataCollection &data) |
| Generate and fill FED raw data for a full event. | |
| void | produce (edm::Event &e, const edm::EventSetup &c) |
| Produce digis out of raw data. | |
| virtual | ~DTTFFEDSim () |
| Destructor. | |
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::InputTag | DTDigiInputTag |
| edm::InputTag | DTPHTFInputTag |
| unsigned int | eventNum |
L1 DT Track Finder Digi-to-Raw
J. Troconiz UAM Madrid
Definition at line 26 of file DTTFFEDSim.h.
| DTTFFEDSim::DTTFFEDSim | ( | const edm::ParameterSet & | pset | ) |
Constructor.
Definition at line 30 of file DTTFFEDSim.cc.
References DTDigiInputTag, DTPHTFInputTag, and edm::ParameterSet::getParameter().
: eventNum(0) { produces<FEDRawDataCollection>(); DTDigiInputTag = pset.getParameter<edm::InputTag>("DTDigi_Source"); DTPHTFInputTag = pset.getParameter<edm::InputTag>("DTTracks_Source"); }
| DTTFFEDSim::~DTTFFEDSim | ( | ) | [virtual] |
| int DTTFFEDSim::bxNr | ( | int | channel | ) | [private] |
Definition at line 292 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 329 of file DTTFFEDSim.cc.
References funct::C, funct::D, 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 272 of file DTTFFEDSim.cc.
References 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 53 of file DTTFFEDSim.cc.
References calcCRC(), channel(), FEDRawData::data(), edm::EventID::event(), eventNum, FEDRawDataCollection::FEDData(), edm::Event::getByLabel(), getDTDigiInputTag(), getDTPHTFInputTag(), edm::EventBase::id(), groupFilesInBlocks::lines, and FEDRawData::resize().
Referenced by produce().
{
eventNum = e.id().event();
int lines = 2;
edm::Handle<L1MuDTChambPhContainer> phtrig;
e.getByLabel(getDTDigiInputTag(),phtrig);
lines += phtrig->bxSize(-1, 1);
edm::Handle<L1MuDTChambThContainer> thtrig;
e.getByLabel(getDTDigiInputTag(),thtrig);
lines += thtrig->bxSize(-1, 1);
edm::Handle<L1MuDTTrackContainer> trtrig;
e.getByLabel(getDTPHTFInputTag(),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 | ( | ) | [inline, private] |
Definition at line 61 of file DTTFFEDSim.h.
References DTDigiInputTag.
Referenced by fillRawData().
{ return DTDigiInputTag; }
| edm::InputTag DTTFFEDSim::getDTPHTFInputTag | ( | ) | [inline, private] |
Definition at line 62 of file DTTFFEDSim.h.
References DTPHTFInputTag.
Referenced by fillRawData().
{ return DTPHTFInputTag; }
| void DTTFFEDSim::produce | ( | edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [virtual] |
Produce digis out of raw data.
Implements edm::EDProducer.
Definition at line 41 of file DTTFFEDSim.cc.
References data, fillRawData(), and edm::Event::put().
{
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 305 of file DTTFFEDSim.cc.
References channel().
{
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 316 of file DTTFFEDSim.cc.
References channel().
{
int myChannel = channel;
if (myChannel > 127) myChannel -= 2;
if (myChannel < 0 || myChannel > 251 ){ return -999; }
int myWheel = ((myChannel%21)/3)-3;
return myWheel;
}
edm::InputTag DTTFFEDSim::DTDigiInputTag [private] |
Definition at line 47 of file DTTFFEDSim.h.
Referenced by DTTFFEDSim(), and getDTDigiInputTag().
edm::InputTag DTTFFEDSim::DTPHTFInputTag [private] |
Definition at line 48 of file DTTFFEDSim.h.
Referenced by DTTFFEDSim(), and getDTPHTFInputTag().
unsigned int DTTFFEDSim::eventNum [private] |
Definition at line 45 of file DTTFFEDSim.h.
Referenced by fillRawData().
1.7.3