|
|
|
#include <TowerBlockFormatter.h>
Public Member Functions | |
| void | DigiToRaw (const EBDataFrame &dataframe, FEDRawData &rawdata, const EcalElectronicsMapping *TheMapping) |
| void | DigiToRaw (const EEDataFrame &dataframe, FEDRawData &rawdata, const EcalElectronicsMapping *TheMapping) |
| void | EndEvent (FEDRawDataCollection *productRawData) |
| std::map< int, std::map< int, int > > * | GetFEDorder () |
| void | StartEvent () |
| TowerBlockFormatter () | |
| ~TowerBlockFormatter () | |
Static Public Attributes | |
| static const int | kCardsPerTower = 5 |
Private Attributes | |
| std::map< int, std::map< int, int > > * | FEDmap |
| std::map< int, std::map< int, int > > * | FEDorder |
Definition at line 22 of file TowerBlockFormatter.h.
| TowerBlockFormatter::TowerBlockFormatter | ( | ) |
Definition at line 15 of file TowerBlockFormatter.cc.
{
}
| TowerBlockFormatter::~TowerBlockFormatter | ( | ) |
Definition at line 19 of file TowerBlockFormatter.cc.
{
}
| void TowerBlockFormatter::DigiToRaw | ( | const EBDataFrame & | dataframe, |
| FEDRawData & | rawdata, | ||
| const EcalElectronicsMapping * | TheMapping | ||
| ) |
Definition at line 25 of file TowerBlockFormatter.cc.
References gather_cfg::cout, FEDRawData::data(), end, Exception, ExpressReco_HICollisions_FallBack::FEDid, createTree::ff, spr::find(), i, EBDataFrame::id(), edm::eventsetup::heterocontainer::insert(), FEDNumbering::MINECALFEDID, reco::print(), FEDRawData::resize(), EcalDataFrame::sample(), FEDRawData::size(), EcalDataFrame::size(), EcalElectronicsId::stripId(), EcalElectronicsId::towerId(), and EcalElectronicsId::xtalId().
{
int bx = *pbx_;
int lv1 = *plv1_ - 1;
int rdsize = rawdata.size() / 8; // size in Word64
bool newFE = false;
const EBDetId& ebdetid = dataframe.id();
int DCCid = TheMapping -> DCCid(ebdetid);
int FEDid = FEDNumbering::MINECALFEDID + DCCid ;
int nsamples = dataframe.size();
// -- FE number
const EcalElectronicsId& elid = TheMapping -> getElectronicsId(ebdetid);
int iFE = elid.towerId();
if (iFE <= 0 || iFE > 68) throw cms::Exception("InvalidFEid") <<
"TowerBlockFormatter::DigiToRaw : Invalid iFE " << iFE << endl;
map<int, map<int,int> >::iterator fen = FEDorder -> find(FEDid);
map<int, map<int,int> >::iterator fed = FEDmap -> find(FEDid);
if (fen == FEDorder -> end()) {
if (debug_) cout << "New FED in TowerBlockFormatter " << dec << FEDid << " 0x" << hex << FEDid << endl;
map<int,int> FEorder;
pair<map<int, map<int,int> >::iterator, bool> t1 = FEDorder -> insert(map<int, map<int,int> >::value_type(FEDid,FEorder));
map<int,int> FEmap;
pair<map<int, map<int,int> >::iterator, bool> t2 = FEDmap -> insert(map<int, map<int,int> >::value_type(FEDid,FEmap));
fen = t1.first;
fed = t2.first;
}
map<int, int>& FEorder = (*fen).second;
map<int, int>& FEmap = (*fed).second;
map<int,int>::iterator fe = FEorder.find(iFE);
int FE_order;
int FE_index;
if (fe != FEorder.end()) {
FE_order = (*fe).second;
map<int,int>::iterator ff = FEmap.find(FE_order);
if (ff == FEmap.end()) cout << "Error with maps... " << endl;
FE_index = (*ff).second;
if (debug_) cout << "FE already there, FE_index = " << dec << FE_index << " FEorder " << FE_order << endl;
}
else {
if (debug_) cout << "New FE in TowerBlockFormatter FE " << dec << iFE << " 0x" << hex << iFE << " in FED id " << dec << FEDid << endl;
newFE = true;
int inser = rdsize;
int number_FEs = FEorder.size() -1;
FE_order = number_FEs+1;
pair<map<int,int>::iterator, bool> t2 = FEorder.insert(map<int,int>::value_type(iFE,FE_order));
if (! t2.second) cout << " FE insertion failed...";
pair<map<int,int>::iterator, bool> tt = FEmap.insert(map<int,int>::value_type(FE_order,inser));
fe = tt.first;
FE_index = (*fe).second;
if (debug_) cout << "Build the Tower Block header for FE id " << iFE << " start at line " << rdsize << endl;
if (debug_) cout << "This is the Fe number (order) " << number_FEs+1 << endl;
rawdata.resize( 8*rdsize + 8);
unsigned char* pData = rawdata.data();
pData[8*FE_index] = iFE & 0xFF;
pData[8*FE_index+1] = (nsamples & 0x7F);
pData[8*FE_index+2] = bx & 0xFF;
pData[8*FE_index+3] = (bx >>8) & 0x0F;
pData[8*FE_index+3] |= 0xa0;
pData[8*FE_index+4] = lv1 & 0xFF;
pData[8*FE_index+5] = (lv1 >>8) & 0x0F;
pData[8*FE_index+6] = 1;
pData[8*FE_index+7] = 0xc0;
if (debug_) print(rawdata);
}
// -- Crystal number inside the SM :
int istrip = elid.stripId();
int ichannel = elid.xtalId();
if (debug_) cout << "Now add crystal : strip channel " << dec << istrip << " " << ichannel << endl;
unsigned char* pData = rawdata.data();
vector<unsigned char> vv(&pData[0],&pData[rawdata.size()]);
int n_add = 2 + 2*nsamples; // 2 bytes per sample, plus 2 bytes before sample 0
if (n_add % 8 != 0) n_add = n_add/8 +1;
else
n_add = n_add/8;
if (debug_) cout << "will add " << n_add << " lines of 64 bits at line " << (FE_index+1) << endl;
rawdata.resize( rawdata.size() + 8*n_add );
unsigned char* ppData = rawdata.data();
vector<unsigned char>::iterator iter = vv.begin() + 8*(FE_index+1);
vector<unsigned char> toadd(n_add*8);
int tzs=0;
toadd[0] = (istrip & 0x7) + ((ichannel & 0x7)<<4);
toadd[1] = (tzs & 0x1) <<12;
for (int isample=0; isample < (n_add*8-2)/2; isample++) {
if (isample < nsamples) {
uint16_t word = (dataframe.sample(isample)).raw(); // 16 bits word corresponding to this sample
toadd[2 + isample*2] = word & 0x00FF;
toadd[2 + isample*2 +1] = (word & 0xFF00)>>8;
}
else {
toadd[2 + isample*2] = 0;
toadd[2 + isample*2 +1] = 0;
}
if (isample % 2 == 0) toadd[2 + isample*2 +1] |= 0xc0; // need to add the B11 header...
}
vv.insert(iter,toadd.begin(),toadd.end());
// update the pData for this FED :
for (int i=0; i < (int)vv.size(); i++) {
ppData[i] = vv[i];
}
if (debug_) {
cout << "pData for this FED is now " << endl;
print(rawdata);
}
// and update the FEmap for this FED :
for (int i=FE_order+1; i < (int)FEorder.size(); i++) {
FEmap[i] += n_add;
if (debug_) cout << "FEmap updated for fe number " << dec << i << endl;
if (debug_) cout << " FEmap[" << i << "] = " << FEmap[i] << endl;
}
// update the block length
int blocklength = ppData[8*FE_index+6]
+ ((ppData[8*FE_index+7] & 0x1)<<8);
blocklength += n_add;
ppData[8*FE_index+6] = blocklength & 0xFF;
ppData[8*FE_index+7] |= (blocklength & 0x100)>>8;
}
| void TowerBlockFormatter::DigiToRaw | ( | const EEDataFrame & | dataframe, |
| FEDRawData & | rawdata, | ||
| const EcalElectronicsMapping * | TheMapping | ||
| ) |
Definition at line 406 of file TowerBlockFormatter.cc.
References gather_cfg::cout, FEDRawData::data(), EcalElectronicsId::dccId(), end, ExpressReco_HICollisions_FallBack::FEDid, createTree::ff, spr::find(), i, EEDataFrame::id(), edm::eventsetup::heterocontainer::insert(), FEDNumbering::MINECALFEDID, reco::print(), FEDRawData::resize(), EcalDataFrame::sample(), EcalDataFrame::size(), FEDRawData::size(), EcalElectronicsId::stripId(), EcalElectronicsId::towerId(), and EcalElectronicsId::xtalId().
{
// debug_ = false;
int bx = *pbx_;
int lv1 = *plv1_;
int rdsize = rawdata.size() / 8; // size in Word64
bool newFE = false;
const EEDetId& eedetid = dataframe.id();
EcalElectronicsId elid = TheMapping -> getElectronicsId(eedetid);
int DCCid = elid.dccId();
int FEDid = FEDNumbering::MINECALFEDID + DCCid ;
int iFE = elid.towerId();
if (debug_) cout << "enter in TowerBlockFormatter::DigiToRaw DCCid FEDid iFE " <<
dec << DCCid << " " << FEDid << " " << iFE << endl;
int nsamples = dataframe.size();
if (iFE <= 0 || iFE > 68) {
cout << "invalid iFE for EndCap DCCid iFE " << DCCid << " " << iFE << endl;
return;
}
map<int, map<int,int> >::iterator fen = FEDorder -> find(FEDid);
map<int, map<int,int> >::iterator fed = FEDmap -> find(FEDid);
if (fen == FEDorder -> end()) {
if (debug_) cout << "New FED in TowerBlockFormatter " << dec << FEDid << " 0x" << hex << FEDid << endl;
map<int,int> FEorder;
pair<map<int, map<int,int> >::iterator, bool> t1 = FEDorder -> insert(map<int, map<int,int> >::value_type(FEDid,FEorder));
map<int,int> FEmap;
pair<map<int, map<int,int> >::iterator, bool> t2 = FEDmap -> insert(map<int, map<int,int> >::value_type(FEDid,FEmap));
fen = t1.first;
fed = t2.first;
}
map<int, int>& FEorder = (*fen).second;
map<int, int>& FEmap = (*fed).second;
map<int,int>::iterator fe = FEorder.find(iFE);
int FE_order;
int FE_index;
if (fe != FEorder.end()) {
FE_order = (*fe).second;
map<int,int>::iterator ff = FEmap.find(FE_order);
if (ff == FEmap.end()) cout << "Error with maps... " << endl;
FE_index = (*ff).second;
if (debug_) cout << "FE already there, FE_index = " << dec << FE_index << " FEorder " << FE_order << endl;
}
else {
if (debug_) cout << "New FE in TowerBlockFormatter FE " << dec << iFE << " 0x" << hex << iFE << " in FED id " << dec << FEDid << endl;
newFE = true;
int inser = rdsize;
int number_FEs = FEorder.size() -1;
FE_order = number_FEs+1;
pair<map<int,int>::iterator, bool> t2 = FEorder.insert(map<int,int>::value_type(iFE,FE_order));
if (! t2.second) cout << " FE insertion failed...";
pair<map<int,int>::iterator, bool> tt = FEmap.insert(map<int,int>::value_type(FE_order,inser));
fe = tt.first;
FE_index = (*fe).second;
if (debug_) cout << "Build the Tower Block header for FE id " << iFE << " start at line " << rdsize << endl;
if (debug_) cout << "This is the Fe number (order) " << number_FEs+1 << endl;
rawdata.resize( 8*rdsize + 8);
unsigned char* pData = rawdata.data();
pData[8*FE_index] = iFE & 0xFF;
pData[8*FE_index+1] = (nsamples & 0x7F);
pData[8*FE_index+2] = bx & 0xFF;
pData[8*FE_index+3] = (bx >>8) & 0x0F;
pData[8*FE_index+3] |= 0xa0;
pData[8*FE_index+4] = lv1 & 0xFF;
pData[8*FE_index+5] = (lv1 >>8) & 0x0F;
pData[8*FE_index+6] = 1;
pData[8*FE_index+7] = 0xc0;
if (debug_) print(rawdata);
}
// -- Crystal number inside the SM :
int istrip = elid.stripId();
int ichannel = elid.xtalId();
if (debug_) cout << "Now add crystal strip channel " << dec << istrip << " " << ichannel << endl;
unsigned char* pData = rawdata.data();
vector<unsigned char> vv(&pData[0],&pData[rawdata.size()]);
int n_add = 2 + 2*nsamples; // 2 bytes per sample, plus 2 bytes before sample 0
if (n_add % 8 != 0) n_add = n_add/8 +1;
else
n_add = n_add/8;
if (debug_) cout << "nsamples = " << dec << nsamples << endl;
if (debug_) cout << "will add " << n_add << " lines of 64 bits at line " << (FE_index+1) << endl;
rawdata.resize( rawdata.size() + 8*n_add );
unsigned char* ppData = rawdata.data();
vector<unsigned char>::iterator iter = vv.begin() + 8*(FE_index+1);
vector<unsigned char> toadd(n_add*8);
int tzs=0;
toadd[0] = (istrip & 0x7) + ((ichannel & 0x7)<<4);
toadd[1] = (tzs & 0x1) <<12;
for (int isample=0; isample < (n_add*8-2)/2; isample++) {
if (isample < nsamples) {
uint16_t word = (dataframe.sample(isample)).raw(); // 16 bits word corresponding to this sample
toadd[2 + isample*2] = word & 0x00FF;
toadd[2 + isample*2 +1] = (word & 0xFF00)>>8;
}
else {
toadd[2 + isample*2] = 0;
toadd[2 + isample*2 +1] = 0;
}
if (isample % 2 == 0) toadd[2 + isample*2 +1] |= 0xc0; // need to add the B11 header...
}
vv.insert(iter,toadd.begin(),toadd.end());
// update the pData for this FED :
for (int i=0; i < (int)vv.size(); i++) {
ppData[i] = vv[i];
}
if (debug_) {
cout << "pData for this FED is now " << endl;
print(rawdata);
}
// and update the FEmap for this FED :
for (int i=FE_order+1; i < (int)FEorder.size(); i++) {
FEmap[i] += n_add;
if (debug_) cout << "FEmap updated for fe number " << dec << i << endl;
if (debug_) cout << " FEmap[" << i << "] = " << FEmap[i] << endl;
}
// update the block length
int blocklength = ppData[8*FE_index+6]
+ ((ppData[8*FE_index+7] & 0x1)<<8);
blocklength += n_add;
ppData[8*FE_index+6] = blocklength & 0xFF;
ppData[8*FE_index+7] |= (blocklength & 0x100)>>8;
}
| void TowerBlockFormatter::EndEvent | ( | FEDRawDataCollection * | productRawData | ) |
Definition at line 183 of file TowerBlockFormatter.cc.
References gather_cfg::cout, FEDRawData::data(), ExpressReco_HICollisions_FallBack::FEDid, first, i, j, prof2calltree::l, geometryCSVtoXML::line, FEDNumbering::MINECALFEDID, dbtoconf::nlines, edm::second(), FEDRawData::size(), and strip().
{
// -- Need to reorder the FE's in teh TowerBlock. They come in the right
// order when reading the unsuppressed digis, but ganz durcheinander
// when reading the SelectiveReadout_Suppressed digis.
if (debug_) cout << "enter in TowerBlockFormatter::EndEvent. First reorder the FE's. " << endl;
for (int idcc=1; idcc <= 54; idcc++) {
// debug_ = (idcc == 34);
//if (idcc != 34) continue;
int FEDid = FEDNumbering::MINECALFEDID + idcc;
// cout << "Process FED " << FEDid << endl;
FEDRawData& fedData = productRawData -> FEDData(FEDid);
if (fedData.size() <= 16) continue;
if (debug_) cout << "This is FEDid = " << FEDid << endl;
unsigned char * pData = fedData.data();
// Word64* words = reinterpret_cast<Word64*>(const_cast<unsigned char*>(pData));
Word64* words = reinterpret_cast<Word64*>(pData);
int length = fedData.size() / 8;
int iDAQ_header(-1), iDCC_header(-1), iTCCBlock_header(-1),
iSRBlock_header(-1), iTowerBlock_header(-1), iDAQ_trailer(-1);
for (int i=length-1; i > -1; i--) {
if ( ( (words[i] >> 60) & 0xF) == 0x5 ) iDAQ_header=i;
if ( ( (words[i] >> 62) & 0x3) == 0x0 ) iDCC_header=i;
if ( ( (words[i] >> 61) & 0x7) == 0x3 ) iTCCBlock_header=i;
if ( ( (words[i] >> 61) & 0x7) == 0x4 ) iSRBlock_header=i;
if ( ( (words[i] >> 62) & 0x3) == 0x3 ) iTowerBlock_header=i;
if ( ( (words[i] >> 60) & 0xF) == 0xA ) iDAQ_trailer=i;
}
if (iTowerBlock_header < 0) iTowerBlock_header = iDAQ_trailer;
if (iSRBlock_header < 0) iSRBlock_header = iTowerBlock_header;
if (iTCCBlock_header < 0) iTCCBlock_header = iSRBlock_header;
if (debug_) {
cout << "iDAQ_header = " << iDAQ_header << endl;
cout << " iDCC_header = " << iDCC_header << endl;
cout << " iTCCBlock_header = " << iTCCBlock_header << endl;
cout << " iSRBlock_header = " << iSRBlock_header << endl;
cout << " iTowerBlock_header = " << iTowerBlock_header << endl;
cout << " iDAQ_trailer = " << iDAQ_trailer << endl;
}
std::map<int, int> FrontEnd;
std::map<int, std::vector<Word64> > Map_xtal_data;
int iTowerBlock_header_keep = iTowerBlock_header;
while (iTowerBlock_header < iDAQ_trailer) {
int fe = words[iTowerBlock_header] & 0xFF;
int nlines = (words[iTowerBlock_header] >> 48) & 0x1FF;
if (debug_) cout << "This is FE number " << fe << "needs nlines = " << nlines << endl;
FrontEnd[fe] = nlines;
std::vector<Word64> xtal_data;
for (int j=0; j < nlines; j++) {
Word64 ww = words[iTowerBlock_header+j];
xtal_data.push_back(ww);
}
Map_xtal_data[fe] = xtal_data;
iTowerBlock_header += nlines;
}
if (debug_) {
cout << "vector of FrontEnd : " << FrontEnd.size() << endl;
for (std::map<int, int>::const_iterator it=FrontEnd.begin();
it != FrontEnd.end(); it++) {
int fe = it -> first;
int l = it -> second;
cout << "FE line " << fe << " " << l << endl;
}
}
iTowerBlock_header = iTowerBlock_header_keep;
for (std::map<int, int>::const_iterator it=FrontEnd.begin();
it != FrontEnd.end(); it++) {
int fe = it -> first;
int nlines = it -> second;
if (debug_) cout << "iTowerBlock_header = " << iTowerBlock_header << endl;
vector<Word64> xtal_data = Map_xtal_data[fe];
for (int j=0; j < nlines; j++) {
words[iTowerBlock_header+j] = xtal_data[j];
if (debug_) cout << "update line " << iTowerBlock_header+j << endl;
}
if (debug_) {
int jFE = pData[8*(iTowerBlock_header)];
cout << "Front End on RD : " << jFE << endl;
}
iTowerBlock_header += nlines;
}
// -- now the FEs are ordered. STill need to order the xtals within FEs;
// need : xtal 1,2,3,4, 5 in strip 1, xtal 1,2,3,4,5 in strip 2 etc..
// with possibly missing ones.
if (debug_) cout << "now reorder the xtals within the FEs" << endl;
iTowerBlock_header = iTowerBlock_header_keep;
for (std::map<int, int>::const_iterator it=FrontEnd.begin();
it != FrontEnd.end(); it++) {
int fe = it -> first;
if (fe > 68) cout << "Problem... fe = " << fe << " in FEDid = " << FEDid << endl;
if (debug_) cout << " This is for FE = " << fe << endl;
int nlines = it -> second;
int timesamples = pData[8*iTowerBlock_header+1] & 0x7F;
int n4=timesamples-3;
int n_lines4 = n4/4;
if ( n4 % 4 != 0) n_lines4 ++;
if (n_lines4<0) n_lines4=0;
int Nxtal_max = (nlines-1)/(1+n_lines4);
int Nxtal = 0;
map< int, map<int, vector<Word64> > > Strip_Map;
while (Nxtal < Nxtal_max) {
int i_xtal = iTowerBlock_header+1 + Nxtal*(1+n_lines4);
int strip = words[i_xtal] & 0x7;
int xtal = ( words[i_xtal] >>4) & 0x7;
map< int, map<int, vector<Word64> > >::iterator iit = Strip_Map.find(strip);
map<int, vector<Word64> > NewMap;
map<int, vector<Word64> > Xtal_Map;
if (iit == Strip_Map.end()) { // new strip
Xtal_Map = NewMap;
}
else {
Xtal_Map = iit -> second;
}
std::vector<Word64> xtal_data;
for (int j=0; j < n_lines4 +1; j++) {
Word64 ww = words[i_xtal +j];
xtal_data.push_back(ww);
}
Xtal_Map[xtal] = xtal_data;
Strip_Map[strip] = Xtal_Map;
Nxtal ++;
}
// now, update the xtals for this FE :
int idx = 0;
for (map< int, map<int, vector<Word64> > >::const_iterator jt = Strip_Map.begin();
jt != Strip_Map.end(); jt++) {
int strip = jt -> first;
if (debug_) cout << " this is strip number " << strip << endl;
map<int, vector<Word64> > Xtal_Map = jt -> second;
for (map<int, vector<Word64> >::const_iterator kt = Xtal_Map.begin();
kt != Xtal_Map.end(); kt++) {
int xtal = kt -> first;
if (debug_) cout << " this is xtal number " << xtal << endl;
vector<Word64> xtal_data = kt -> second;
int mlines = (int)xtal_data.size();
if (debug_) cout << " mlines = " << mlines << endl;
for (int j=0; j < mlines; j++) {
int line = iTowerBlock_header+1+idx+j;
if (line >= iDAQ_trailer) cout << "smth wrong... line " << line << " trailer " << iDAQ_trailer << endl;
words[line] = xtal_data[j] ;
if (debug_) cout << " updated line " << iTowerBlock_header+idx+j << endl;
}
idx += mlines;
} // end loop on xtals
Xtal_Map.clear();
} // end loop on strips
Strip_Map.clear();
iTowerBlock_header += nlines;
} // end loop on FEs
if (debug_) cout << " DONE FOR FED " << FEDid << endl;
FrontEnd.clear();
Map_xtal_data.clear();
} // end loop on DCC
// cout << " finished reorder, now clean up " << endl;
// -- clean up
// FEDmap -> empty();
// FEDorder -> empty();
FEDmap -> clear();
FEDorder -> clear();
delete FEDmap;
delete FEDorder;
FEDmap = 0;
FEDorder = 0;
debug_ = false;
// cout << "end of EndEvent " << endl;
}
| std::map<int, std::map<int,int> >* TowerBlockFormatter::GetFEDorder | ( | ) | [inline] |
| void TowerBlockFormatter::StartEvent | ( | ) |
Definition at line 396 of file TowerBlockFormatter.cc.
{
FEDmap = new map<int, map<int,int> >;
FEDorder = new map<int, map<int,int> >;
}
std::map<int, std::map<int,int> >* TowerBlockFormatter::FEDmap [private] |
Definition at line 36 of file TowerBlockFormatter.h.
std::map<int, std::map<int,int> >* TowerBlockFormatter::FEDorder [private] |
Definition at line 37 of file TowerBlockFormatter.h.
Referenced by GetFEDorder().
const int TowerBlockFormatter::kCardsPerTower = 5 [static] |
Reimplemented from BlockFormatter.
Definition at line 27 of file TowerBlockFormatter.h.
1.7.3