d9/d40/DTTFFEDSim_8cc_source.html

 //-------------------------------------------------

 //

 //   Class: DTTFFEDSim

 //

 //   L1 DT Track Finder Digi-to-Raw

 //

 //

 //

 //   Author :

 //   J. Troconiz  UAM Madrid

 //

 //--------------------------------------------------


 #include "EventFilter/DTTFRawToDigi/interface/DTTFFEDSim.h"


 #include "DataFormats/Common/interface/Handle.h"

 #include "FWCore/Framework/interface/Event.h"


 #include <DataFormats/FEDRawData/interface/FEDRawData.h>


 #include <iostream>


 using namespace std;


 DTTFFEDSim::DTTFFEDSim(const edm::ParameterSet& pset) : 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(){}


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


   FEDRawDataCollection data;


   if (!fillRawData(e, data)) return;


   auto_ptr<FEDRawDataCollection> fed_product(new FEDRawDataCollection(data));


   e.put(fed_product);


 }


 bool DTTFFEDSim::fillRawData(edm::Event& e,

                              FEDRawDataCollection& data) {


   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;

 }


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


   // 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;

 }


 int DTTFFEDSim::bxNr( int channel ){


   int myChannel = channel;


   if (myChannel > 127) myChannel -= 2;


   if (myChannel < 0 || myChannel > 251 ){ return -999; }


   int myBx = 1-(myChannel%3);


   return myBx;

 }


 int DTTFFEDSim::sector( int channel ){


   int myChannel = channel;


   if (myChannel > 127) myChannel -= 2;


   if (myChannel < 0 || myChannel > 251 ){ return -999; }


   return myChannel/21;

 }


 int DTTFFEDSim::wheel( int channel ){


   int myChannel = channel;


   if (myChannel > 127) myChannel -= 2;


   if (myChannel < 0 || myChannel > 251 ){ return -999; }


   int myWheel = ((myChannel%21)/3)-3;


   return myWheel;

 }


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


   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;

 }

edm::ParameterSet::getParameter
T getParameter(std::string const &) const

edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:44

i
int i
Definition: DBlmapReader.cc:9

edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434

DTTFFEDSim::~DTTFFEDSim
virtual ~DTTFFEDSim()
Destructor.
Definition: DTTFFEDSim.cc:38

Event.h

DTTFFEDSim::ChTh_tok
edm::EDGetTokenT< L1MuDTChambThContainer > ChTh_tok
Definition: DTTFFEDSim.h:68

funct::C
C
Definition: Factorize.h:141

Handle.h

edm::Handle
Definition: AssociativeIterator.h:47

L1MuDTChambThContainer::The_iterator
The_Container::const_iterator The_iterator
Definition: L1MuDTChambThContainer.h:39

DTTFFEDSim::sector
int sector(int channel)
Definition: DTTFFEDSim.cc:304

DTTFFEDSim::eventNum
unsigned int eventNum
Definition: DTTFFEDSim.h:48

DTTFFEDSim.h

FEDRawDataCollection::FEDData
const FEDRawData & FEDData(int fedid) const
retrieve data for fed
Definition: FEDRawDataCollection.cc:28

DTTFFEDSim::DTDigiInputTag
edm::InputTag DTDigiInputTag
Definition: DTTFFEDSim.h:50

FEDRawData::resize
void resize(size_t newsize)
Definition: FEDRawData.cc:32

edm::Event::put
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:116

FEDRawData
Definition: FEDRawData.h:20

DTTFFEDSim::wheel
int wheel(int channel)
Definition: DTTFFEDSim.cc:315

DTTFFEDSim::calcCRC
void calcCRC(int myD1, int myD2, int &myC)
Definition: DTTFFEDSim.cc:328

funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22

edm::EventSetup
Definition: EventSetup.h:44

groupFilesInBlocks.lines
tuple lines
Definition: groupFilesInBlocks.py:93

L1MuDTChambPhContainer::Phi_iterator
Phi_Container::const_iterator Phi_iterator
Definition: L1MuDTChambPhContainer.h:39

DTTFFEDSim::channel
int channel(int wheel, int sector, int bx)
Definition: DTTFFEDSim.cc:271

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

L1MuDTTrackContainer::Trackiterator
TrackContainer::const_iterator Trackiterator
Definition: L1MuDTTrackContainer.h:39

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

trackerHits.c
tuple c
Definition: trackerHits.py:26

FEDRawDataCollection
Definition: FEDRawDataCollection.h:19

alignCSCRings.e
list e
Definition: alignCSCRings.py:90

DTTFFEDSim::bxNr
int bxNr(int channel)
Definition: DTTFFEDSim.cc:291

edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56

data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82

edm::InputTag
Definition: InputTag.h:17

FEDRawData.h

DTTFFEDSim::Trk_tok
edm::EDGetTokenT< L1MuDTTrackContainer > Trk_tok
Definition: DTTFFEDSim.h:69

DTTFFEDSim::ChPh_tok
edm::EDGetTokenT< L1MuDTChambPhContainer > ChPh_tok
Definition: DTTFFEDSim.h:67

FEDRawData::data
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:28

DTTFFEDSim::DTPHTFInputTag
edm::InputTag DTPHTFInputTag
Definition: DTTFFEDSim.h:51

edm::ParameterSet
Definition: ParameterSet.h:35

edm::Event
Definition: Event.h:62

funct::D
DecomposeProduct< arg, typename Div::arg > D
Definition: Factorize.h:150

DTTFFEDSim::DTTFFEDSim
DTTFFEDSim(const edm::ParameterSet &pset)
Constructor.
Definition: DTTFFEDSim.cc:25