|
|
|
#include <CalibCalorimetry/HcalLuttoDB/src/HcalLuttoDB.cc>
Public Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
| virtual void | beginJob () |
| virtual void | endJob () |
| HcalLuttoDB (const edm::ParameterSet &) | |
| ~HcalLuttoDB () | |
Private Member Functions | |
| std::vector< unsigned char > | extractOutputLut (const CaloTPGTranscoder &coder, HcalTrigTowerDetId chan) |
| std::ostream * | openChecksums () |
| std::ostream * | openPerCrate (int crate) |
| std::ostream * | openPerLut1 (HcalElectronicsId eid) |
| std::ostream * | openPerLut2 (HcalElectronicsId eid) |
| void | writeoutlut1 (HcalDetId id, HcalElectronicsId eid, const std::vector< unsigned short > &lut, std::ostream &os) |
| void | writeoutlut2 (HcalTrigTowerDetId id, HcalElectronicsId eid, const std::vector< unsigned char > &lut, std::ostream &os) |
Private Attributes | |
| std::string | creationstamp_ |
| std::string | fileformat_ |
| bool | filePerCrate_ |
| std::ostream * | oc_ |
Description: <one line="" class="" summary>=""> R Implementation: <Notes on="" implementation>="">
Definition at line 57 of file HcalLuttoDB.cc.
| HcalLuttoDB::HcalLuttoDB | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 87 of file HcalLuttoDB.cc.
References creationtag_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and targetfirmware_.
{
creationtag_ = iConfig.getParameter<std::string>("creationtag");
targetfirmware_ = iConfig.getParameter<std::string>("targetfirmware");
filePerCrate_ = iConfig.getUntrackedParameter<bool>("filePerCrate",true);
fileformat_=iConfig.getParameter<std::string>("filePrefix");
}
| HcalLuttoDB::~HcalLuttoDB | ( | ) |
Definition at line 96 of file HcalLuttoDB.cc.
{
}
| void HcalLuttoDB::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 262 of file HcalLuttoDB.cc.
References HcalElectronicsMap::allElectronicsId(), HcalGenericDetId::genericSubdet(), edm::EventSetup::get(), HcalGenericDetId::HcalGenBarrel, HcalGenericDetId::HcalGenEndcap, HcalGenericDetId::HcalGenForward, CaloTPGTranscoder::HcalTPG, HcalElectronicsMap::lookup(), HcalElectronicsMap::lookupTrigger(), lumiPlot::lut, and DetId::null().
{
//using namespace edm;
//using namespace std;
edm::LogInfo("Hcal") << "Beginning dump of Hcal TPG LUTS (this may take a minute or two)";
const HcalElectronicsMap* Map_;
ESHandle<HcalDbService> pSetup;
iSetup.get<HcalDbRecord>().get( pSetup );
Map_=pSetup->getHcalMapping();
// get the conditions, for the decoding
edm::ESHandle<HcalTPGCoder> inputCoder;
iSetup.get<HcalTPGRecord>().get(inputCoder);
edm::ESHandle<CaloTPGTranscoder> outTranscoder;
iSetup.get<CaloTPGRecord>().get(outTranscoder);
outTranscoder->setup(iSetup,CaloTPGTranscoder::HcalTPG);
std::vector<HcalElectronicsId> allEID = Map_->allElectronicsId();
std::vector<HcalElectronicsId>::iterator itreid;
std::ostream* pfile=0;
oc_=openChecksums();
for (int crate=0; crate<20; crate++) {
edm::LogInfo("Hcal") << "Beginning crate " << crate;
for(itreid = allEID.begin(); itreid != allEID.end(); itreid++)
{
if (itreid->readoutVMECrateId()!=crate) continue;
if (itreid->isTriggerChainId()) { // lut2
HcalTrigTowerDetId tid=Map_->lookupTrigger(*itreid);
if (tid.null()) continue;
if (filePerCrate_ && pfile==0) pfile=openPerCrate(crate);
else if (pfile==0) pfile=openPerLut2(*itreid);
std::vector<unsigned char> lut=extractOutputLut(*outTranscoder,tid);
writeoutlut2(tid,*itreid,lut,*pfile);
if (!filePerCrate_) { delete pfile; pfile=0; }
} else { // lut1
HcalGenericDetId gid=Map_->lookup(*itreid);
if (gid.null() || !(gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel || gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap || gid.genericSubdet()==HcalGenericDetId::HcalGenForward)) continue;
if (filePerCrate_ && pfile==0) pfile=openPerCrate(crate);
else if (pfile==0) pfile=openPerLut1(*itreid);
std::vector<unsigned short> lut=inputCoder->getLinearizationLUT(HcalDetId(gid));
writeoutlut1(HcalDetId(gid),*itreid,lut,*pfile);
if (!filePerCrate_) { delete pfile; pfile=0; }
}
}
if (pfile!=0) {
if (filePerCrate_) *pfile << "</CFGBrickSet>\n";
delete pfile;
pfile=0;
}
}
*oc_ << "</CFGBrick>\n";
delete oc_;
outTranscoder->releaseSetup();
}
| void HcalLuttoDB::beginJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 329 of file HcalLuttoDB.cc.
References cmsPerfSuiteHarvest::now, and cond::rpcobgas::time.
{
char buffer[120];
time_t now=time(0);
struct tm* tm=localtime(&now);
strftime(buffer,120,"%F %T",tm);
creationstamp_ = buffer;
}
| void HcalLuttoDB::endJob | ( | void | ) | [virtual] |
| std::vector< unsigned char > HcalLuttoDB::extractOutputLut | ( | const CaloTPGTranscoder & | coder, |
| HcalTrigTowerDetId | chan | ||
| ) | [private] |
Definition at line 251 of file HcalLuttoDB.cc.
References HcalTriggerPrimitiveSample::compressedEt(), CaloTPGTranscoder::hcalCompress(), i, lumiPlot::lut, and alignCSCRings::s.
{
std::vector<unsigned char> lut;
for (int i=0; i<1024; i++) {
HcalTriggerPrimitiveSample s=coder.hcalCompress(chan,i,false);
lut.push_back(s.compressedEt());
}
return lut;
}
| std::ostream * HcalLuttoDB::openChecksums | ( | ) | [private] |
Definition at line 105 of file HcalLuttoDB.cc.
References alignmentValidation::fname.
{
char fname[1024];
snprintf(fname,1024,"%s_checksums.xml",fileformat_.c_str());
std::ostream* os=new std::ofstream(fname);
(*os) << "<?xml version=\"1.0\"?>\n<CFGBrick>\n";
return os;
}
| std::ostream * HcalLuttoDB::openPerCrate | ( | int | crate | ) | [private] |
Definition at line 113 of file HcalLuttoDB.cc.
References alignmentValidation::fname.
{
char fname[1024];
snprintf(fname,1024,"%s_%d.xml",fileformat_.c_str(),crate);
std::ostream* os=new std::ofstream(fname);
(*os) << "<?xml version=\"1.0\"?>\n<CFGBrickSet>\n";
return os;
}
| std::ostream * HcalLuttoDB::openPerLut1 | ( | HcalElectronicsId | eid | ) | [private] |
Definition at line 121 of file HcalLuttoDB.cc.
References HcalElectronicsId::fiberChanId(), HcalElectronicsId::fiberIndex(), alignmentValidation::fname, HcalElectronicsId::htrSlot(), HcalElectronicsId::htrTopBottom(), and HcalElectronicsId::readoutVMECrateId().
{
char fname[1024];
snprintf(fname,1024,"%s_%d_%d%c_%d_%d_1.xml",fileformat_.c_str(),eid.readoutVMECrateId(),
eid.htrSlot(),((eid.htrTopBottom())?('t'):('b')),eid.fiberIndex(),eid.fiberChanId());
std::ostream* os=new std::ofstream(fname);
(*os) << "<?xml version=\"1.0\"?>\n";
return os;
}
| std::ostream * HcalLuttoDB::openPerLut2 | ( | HcalElectronicsId | eid | ) | [private] |
Definition at line 130 of file HcalLuttoDB.cc.
References alignmentValidation::fname, HcalElectronicsId::htrSlot(), HcalElectronicsId::htrTopBottom(), HcalElectronicsId::readoutVMECrateId(), HcalElectronicsId::slbChannelIndex(), and HcalElectronicsId::slbSiteNumber().
{
char fname[1024];
snprintf(fname,1024,"%s_%d_%d%c_%d_%d_2.xml",fileformat_.c_str(),eid.readoutVMECrateId(),
eid.htrSlot(),((eid.htrTopBottom())?('t'):('b')),eid.slbSiteNumber(),eid.slbChannelIndex());
std::ostream* os=new std::ofstream(fname);
(*os) << "<?xml version=\"1.0\"?>\n";
return os;
}
| void HcalLuttoDB::writeoutlut1 | ( | HcalDetId | id, |
| HcalElectronicsId | eid, | ||
| const std::vector< unsigned short > & | lut, | ||
| std::ostream & | os | ||
| ) | [private] |
Definition at line 141 of file HcalLuttoDB.cc.
References creationtag_, HcalElectronicsId::fiberChanId(), HcalElectronicsId::fiberIndex(), formatRevision_, HcalForward, HcalElectronicsId::htrSlot(), HcalElectronicsId::htrTopBottom(), i, md5_append(), md5_finish(), md5_init(), HcalElectronicsId::readoutVMECrateId(), and targetfirmware_.
{
os <<"<CFGBrick> "<<std::endl;
os <<" <Parameter name='IETA' type='int'>"<<id.ieta()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='IPHI' type='int'>"<<id.iphi()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='DEPTH' type='int'>"<<id.depth()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='CRATE' type='int'>"<<eid.readoutVMECrateId()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='SLOT' type='int'>"<<eid.htrSlot()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='TOPBOTTOM' type='int'>"<<eid.htrTopBottom()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='FIBER' type='int'>"<<eid.fiberIndex()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='FIBERCHAN' type='int'>"<<eid.fiberChanId()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='LUT_TYPE' type='int'>1</Parameter>"<<std::endl;
os <<" <Parameter name='CREATIONTAG' type='string'>"<<creationtag_<<"</Parameter>"<<std::endl;
os <<" <Parameter name='CREATIONSTAMP' type='string'>"<<creationstamp_<<"</Parameter>"<<std::endl;
os <<" <Parameter name='FORMATREVISION' type='string'>"<<formatRevision_<<"</Parameter>"<<std::endl;
os <<" <Parameter name='TARGETFIRMWARE' type='string'>"<<targetfirmware_<<"</Parameter>"<<std::endl;
int generalizedIndex=id.ietaAbs()+1000*id.depth()+10000*id.iphi()+
((id.ieta()<0)?(0):(100))+((id.subdet()==HcalForward && id.ietaAbs()==29)?(4*10000):(0));
os <<" <Parameter name='GENERALIZEDINDEX' type='int'>"<<generalizedIndex<<"</Parameter>"<<std::endl;
// do checksum
md5_state_t md5er;
md5_byte_t digest[16];
unsigned char tool[2];
md5_init(&md5er);
for (int i=0; i<128; i++) {
tool[0]=lut[i]&0xFF;
tool[1]=(lut[i]>>8)&0xFF;
md5_append(&md5er,tool,2);
}
md5_finish(&md5er,digest);
os <<" <Parameter name='CHECKSUM' type='string'>";
for (int i=0; i<16; i++) os << std::hex << (((int)(digest[i]))&0xFF);
os << "</Parameter>\n";
*oc_ << " <Data crate='" << eid.readoutVMECrateId()
<< "' slot='" << eid.htrSlot()
<< "' fpga='" << eid.htrTopBottom()
<< "' fiber='" << eid.fiberIndex()
<< "' fiberchan='" << eid.fiberChanId()
<< "' luttype='1' elements='1' encoding='hex'>";
for (int i=0; i<16; i++) *oc_ << std::hex << (((int)(digest[i]))&0xFF);
*oc_ << "</Data>\n";
os <<" <Data elements='128' encoding='hex'> "<<std::endl;
os << std::hex;
for(int initr2 = 0; initr2 < 128; initr2++){
os<<lut[initr2]<<" ";
}
os << std::dec;
os<<std::endl;
os <<" </Data> "<<std::endl;
os <<"</CFGBrick> "<<std::endl;
}
| void HcalLuttoDB::writeoutlut2 | ( | HcalTrigTowerDetId | id, |
| HcalElectronicsId | eid, | ||
| const std::vector< unsigned char > & | lut, | ||
| std::ostream & | os | ||
| ) | [private] |
Definition at line 198 of file HcalLuttoDB.cc.
References creationtag_, formatRevision_, HcalElectronicsId::htrSlot(), HcalElectronicsId::htrTopBottom(), i, md5_append(), md5_finish(), md5_init(), HcalElectronicsId::readoutVMECrateId(), HcalElectronicsId::slbChannelIndex(), HcalElectronicsId::slbSiteNumber(), and targetfirmware_.
{
os <<"<CFGBrick> "<<std::endl;
os <<" <Parameter name='IETA' type='int'>"<<id.ieta()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='IPHI' type='int'>"<<id.iphi()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='CRATE' type='int'>"<<eid.readoutVMECrateId()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='SLOT' type='int'>"<<eid.htrSlot()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='TOPBOTTOM' type='int'>"<<eid.htrTopBottom()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='SLB' type='int'>"<<eid.slbSiteNumber()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='SLBCHAN' type='int'>"<<eid.slbChannelIndex()<<"</Parameter>"<<std::endl;
os <<" <Parameter name='LUT_TYPE' type='int'>2</Parameter>"<<std::endl;
os <<" <Parameter name='CREATIONTAG' type='string'>"<<creationtag_<<"</Parameter>"<<std::endl;
os <<" <Parameter name='CREATIONSTAMP' type='string'>"<<creationstamp_<<"</Parameter>"<<std::endl;
os <<" <Parameter name='FORMATREVISION' type='string'>"<<formatRevision_<<"</Parameter>"<<std::endl;
os <<" <Parameter name='TARGETFIRMWARE' type='string'>"<<targetfirmware_<<"</Parameter>"<<std::endl;
int generalizedIndex=id.ietaAbs()+10000*id.iphi()+
((id.ieta()<0)?(0):(100));
os <<" <Parameter name='GENERALIZEDINDEX' type='int'>"<<generalizedIndex<<"</Parameter>"<<std::endl;
// do checksum
md5_state_t md5er;
md5_byte_t digest[16];
md5_init(&md5er);
md5_append(&md5er,&(lut[0]),1024);
md5_finish(&md5er,digest);
os <<" <Parameter name='CHECKSUM' type='string'>";
for (int i=0; i<16; i++) os << std::hex << (((int)(digest[i]))&0xFF);
os << "</Parameter>\n";
*oc_ << " <Data crate='" << eid.readoutVMECrateId()
<< "' slot='" << eid.htrSlot()
<< "' fpga='" << eid.htrTopBottom()
<< "' slb='" << eid.slbSiteNumber()
<< "' slbchan='" << eid.slbChannelIndex()
<< "' luttype='2' elements='1' encoding='hex'>";
for (int i=0; i<16; i++) *oc_ << std::hex << (((int)(digest[i]))&0xFF);
*oc_ << "</Data>\n";
os <<" <Data elements='1024' encoding='hex'> "<<std::endl;
os << std::hex;
for(int initr2 = 0; initr2 < 1024; initr2++){
os<< (int(lut[initr2])&0xFF)<<" ";
}
os << std::dec;
os<<std::endl;
os <<" </Data> "<<std::endl;
os <<"</CFGBrick> "<<std::endl;
}
std::string HcalLuttoDB::creationstamp_ [private] |
Definition at line 71 of file HcalLuttoDB.cc.
std::string HcalLuttoDB::fileformat_ [private] |
Definition at line 72 of file HcalLuttoDB.cc.
bool HcalLuttoDB::filePerCrate_ [private] |
Definition at line 70 of file HcalLuttoDB.cc.
std::ostream* HcalLuttoDB::oc_ [private] |
Definition at line 77 of file HcalLuttoDB.cc.
1.7.3