|
|
|
#include <L1TdeGCT.h>
Definition at line 27 of file L1TdeGCT.h.
| L1TdeGCT::L1TdeGCT | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 6 of file L1TdeGCT.cc.
References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), NULL, cmsCodeRules::cppFunctionSkipper::operator, and validate_alignment_devdb10_cfg::verbose.
{
verbose_ = iConfig.getUntrackedParameter<int>("VerboseFlag",0);
if(verbose())
std::cout << "L1TdeGCT::L1TdeGCT()...\n" << std::flush;
DEsource_ = iConfig.getParameter<edm::InputTag>("DataEmulCompareSource");
histFolder_ = iConfig.getUntrackedParameter<std::string>("HistFolder", "L1TEMU/GCTexpert/");
dbe = NULL;
if (iConfig.getUntrackedParameter<bool>("DQMStore", false)) {
dbe = edm::Service<DQMStore>().operator->();
dbe->setVerbose(0);
}
histFile_ = iConfig.getUntrackedParameter<std::string>("HistFile", "");
if(iConfig.getUntrackedParameter<bool> ("disableROOToutput", true))
histFile_ = "";
if (histFile_.size()!=0) {
edm::LogInfo("OutputRootFile")
<< "L1TEmulator GCT specific histograms will be saved to "
<< histFile_.c_str()
<< std::endl;
}
if(dbe!=NULL)
dbe->setCurrentFolder(histFolder_);
hasRecord_=true;
if(verbose())
std::cout << "L1TdeGCT::L1TdeGCT()...done.\n" << std::flush;
}
| L1TdeGCT::~L1TdeGCT | ( | ) |
Definition at line 42 of file L1TdeGCT.cc.
{}
| void L1TdeGCT::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [protected, virtual] |
get the comparison results
get the de candidates
--- Fill histograms(me) ---
(note see L1Trigger/HardwareValidation/interface/DEtrait.h)
bitset loop
event based rate
object based rate
Implements edm::EDAnalyzer.
Definition at line 219 of file L1TdeGCT.cc.
References edm::HandleBase::clear(), gather_cfg::cout, eta(), fastmath::etaphi(), dedefs::GCT, dedefs::GCTisolaem, edm::Event::getByLabel(), i, edm::HandleBase::isValid(), j, LogDebug, phi, L1DataEmulDigi::reset(), evf::utils::sid, dedefs::SystLabel, and validate_alignment_devdb10_cfg::verbose.
{
if(!hasRecord_)
return;
if(verbose())
std::cout << "L1TdeGCT::analyze() start\n" << std::flush;
edm::Handle<L1DataEmulRecord> deRecord;
iEvent.getByLabel(DEsource_, deRecord);
if (!deRecord.isValid()) {
edm::LogInfo("DataNotFound")
<< "Cannot find L1DataEmulRecord with label "
<< DEsource_.label()
<< " Please verify that comparator was successfully executed."
<< " Emulator DQM for GCT will be skipped!"
<< std::endl;
hasRecord_=false;
return;
}
bool isComp = deRecord->get_isComp(GCT);
if(!isComp) {
if(verbose())
std::cout << "[L1TdeGCT] Gct information not generated in de-record."
<< " Skiping event!\n" << std::flush;
return;
}
int DEncand[2]={0};
for(int j=0; j<2; j++)
DEncand[j] = deRecord->getNCand(GCT,j);
if(verbose())
std::cout << "[L1TdeGCT] ncands"
<< " data: " << DEncand[0]
<< " emul: " << DEncand[1]
<< std::endl;
L1DEDigiCollection deColl;
deColl = deRecord->getColl();
// extract the GCT comparison digis
L1DEDigiCollection gctColl;
gctColl.reserve(20);
gctColl.clear();
for(L1DEDigiCollection::const_iterator it=deColl.begin();
it!=deColl.end(); it++)
if(!it->empty())
if(it->sid()==GCT)
gctColl.push_back(*it);
if(verbose()) {
std::cout << "[L1TdeGCT] record has " << gctColl.size()
<< " gct de digis\n" << std::flush;
for(L1DEDigiCollection::const_iterator it=gctColl.begin();
it!=gctColl.end(); it++)
std::cout << "\t" << *it << std::endl;
}
const int nullVal = L1DataEmulDigi().reset();
// d|e candidate loop
for(L1DEDigiCollection::const_iterator it=gctColl.begin();
it!=gctColl.end(); it++) {
// sid should be GCT
int sid = it->sid();
// cid: GCTisolaem, GCTnoisoem, GCTcenjets, GCTforjets, GCTtaujets
int cid = it->cid();
if(verbose())
std::cout << "[L1TdeGCT] processing digi "
<< " sys:" << sid
<< " type:" << cid
<< " \n\t"
<< *it << "\n"
<< std::flush;
//assert(cid==GCT);
if(sid!=GCT || it->empty()) {
LogDebug("L1TdeGCT") << "consistency check failure, non-gct digis!";
continue;
}
int type = it->type();
double phiv = it->x1();
double etav = it->x2();
float rankarr[2];
it->rank(rankarr);
float rnkv = rankarr[0];
double wei = 1.;
unsigned int mask = (~0x0);
// shift coll type for starting at zero
int ccid = cid - dedefs::GCTisolaem;
if(ccid<0 || ccid >= nGctColl_) {
LogDebug("L1TdeGCT") << "consistency check failure, col type outbounds:"
<< ccid << "\n";
ccid=0;
}
//type: 0:agree 1:loc.agree, 2:loc.disagree, 3:data.only, 4:emul.only
if(it->type()<4)
sysncand[0]->Fill(ccid);
if(it->type()<5&&it->type()!=3)
sysncand[1]->Fill(ccid);
errortype[ccid]->Fill(type);
wei=1.; if(!type) wei=0.;
if(etav!=nullVal && phiv!=nullVal)
etaphi[ccid]->Fill(etav,phiv,wei);
if(etav!=nullVal)
eta [ccid]->Fill(etav,wei);
if(phiv!=nullVal)
phi [ccid]->Fill(phiv,wei);
rnk[ccid]->Fill(rnkv,wei);
//exclude e-only cands (only data)
wei=1.;if(type==4) wei=0.;
if(etav!=nullVal)
etaData[ccid]->Fill(etav,wei);
if(phiv!=nullVal)
phiData[ccid]->Fill(phiv,wei);
rnkData[ccid]->Fill(rnkv,wei);
wei=1;
// GCT trigger bits
unsigned int word[2];
it->data(word);
std::bitset<32> dbits(word[0]);
std::bitset<32> ebits(word[1]);
unsigned int dexor = ( (word[0]) ^ (word[1]) );
//disagreeing bits
std::bitset<32> debits(dexor);
//disagreeing bits after masking
std::bitset<32> dembits( ( (dexor) & (mask) ) );
if(verbose())
std::cout << "l1degct"
<< " sid:" << sid << " cid:" << cid << "\n"
<< " data:0x" << std::hex << word[0] << std::dec
<< " bitset:" << dbits
<< "\n"
<< " emul:0x" << std::hex << word[1] << std::dec
<< " bitset:" << ebits
<< "\n"
<< " xor:0x" << std::hex << dexor << std::dec
<< " bitset:" << debits
<< " bitset:" << ( (dbits) ^ (ebits) )
<< "\n" << std::flush;
for(int ibit=0; ibit<32; ibit++) {
wei=1.;
//comparison gives no info if there's only 1 candidate
if(type==3 || type==4) wei=0.;
if(dbits [ibit]) dword[sid]->Fill(ibit,wei);
if(ebits [ibit]) eword[sid]->Fill(ibit,wei);
if(debits [ibit])deword[sid]->Fill(ibit,wei);
//if(dembits[ibit])masked[sid]->Fill(ibit,wei);
}
wei=1;
}
//error rates per GCT trigger object type
int hasCol[nGctColl_]={0};
int nagree[nGctColl_]={0};
for(L1DEDigiCollection::const_iterator it=gctColl.begin();
it!=gctColl.end(); it++) {
int ccid = it->cid()-dedefs::GCTisolaem;
ccid = (ccid<0 || ccid >= nGctColl_) ? 0:ccid;
hasCol[ccid]++;
if(!it->type())
nagree[ccid]++;
}
for(int i=0; i<nGctColl_; i++) {
if(!hasCol[i]) continue;
//nWithCol[i]++;
//if(nagree[i]<hasCol[i]) colCount[i]++;
nWithCol[i]+=hasCol[i];//#of objects
colCount[i]+=nagree[i];//#of agreements
}
for(int i=0; i<nGctColl_; i++) {
int ibin = i+1;
double rate = nWithCol[i] ? 1.-1.*colCount[i]/nWithCol[i]: 0.;
sysrates->setBinContent(ibin,rate);
if(verbose()) {
std::cout << "[L1TDEMON] analyze rate computation\t\n"
<< " colid:" << i
<< "(so far)"
<< " nWithCol: " << nWithCol[i]
<< " colCount: " << colCount[i]
<< "(this event)"
<< "hasCol: " << hasCol[i]
<< " nagree: " << nagree[i]
<< " rate:" << sysrates->getBinContent(ibin)
<< "\n" << std::flush;
if(rate>1. || rate<0.)
std::cout << "problem, error rate for " << SystLabel[i]
<<" is "<<sysrates->getBinContent(ibin)
<< "\n" << std::flush;
}
}
if(verbose())
std::cout << "L1TdeGCT::analyze() end.\n" << std::flush;
}
| void L1TdeGCT::beginJob | ( | void | ) | [protected, virtual] |
labeling
Reimplemented from edm::EDAnalyzer.
Definition at line 45 of file L1TdeGCT.cc.
References DQMStore::book1D(), DQMStore::book2D(), gather_cfg::cout, eta(), fastmath::etaphi(), i, j, cmsCodeRules::cppFunctionSkipper::operator, phi, DQMStore::rmdir(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), indexGen::title, and validate_alignment_devdb10_cfg::verbose.
{
if(verbose())
std::cout << "L1TdeGCT::beginJob() start\n" << std::flush;
DQMStore* dbe = 0;
dbe = edm::Service<DQMStore>().operator->();
if(dbe) {
dbe->setCurrentFolder(histFolder_);
dbe->rmdir(histFolder_);
}
// (em) iso, no-iso, (jets) cen, for, tau
std::string cLabel[nGctColl_]=
{"IsoEM", "NoisoEM", "CenJet", "ForJet", "TauJet", "HT", "MET", "ET", "MHT", "HFSums", "HFCnts"};
const int nerr = 5;
const int nbit = 32;
if(dbe) {
dbe->setCurrentFolder(histFolder_);
// book histograms here
const int phiNBins = 18 ;
const double phiMinim = -0.5;
const double phiMaxim = 17.5;
const int etaNBins = 22 ;
const double etaMinim = -0.5;
const double etaMaxim = 21.5;
const int rnkNBins = 63;
const double rnkMinim = 0.5;
const double rnkMaxim = 63.5;
sysrates = dbe->book1D("sysrates","RATE OF COMPARISON FAILURES",nGctColl_, 0, nGctColl_ );
for(int j=0; j<2; j++) {
std::string lbl("sysncand");
lbl += (j==0?"Data":"Emul");
std::string title("GCT OBJECT MULTIPLICITY ");
title += (j==0?"(DATA)":"(EMULATOR)");
sysncand[j] = dbe->book1D(lbl.data(),title.data(),nGctColl_, 0, nGctColl_ );
}
for(int j=0; j<nGctColl_; j++) {
dbe->setCurrentFolder(std::string(histFolder_+cLabel[j]));
std::string lbl("");
std::string title("");
lbl.clear();
title.clear();
lbl+=cLabel[j];lbl+="ErrorFlag";
title+=cLabel[j];title+=" ErrorFlag";
errortype[j] = dbe->book1D(lbl.data(),title.data(), nerr, 0, nerr);
lbl.clear();
title.clear();
lbl+=cLabel[j];lbl+="Eta";
title+=cLabel[j];title+=" ETA OF COMPARISON FAILURES";
eta[j] = dbe->book1D(lbl.data(),title.data(),
etaNBins, etaMinim, etaMaxim);
lbl.clear();
title.clear();
lbl+=cLabel[j];lbl+="Phi";
title+=cLabel[j];title+=" PHI OF COMPARISON FAILURES";
phi[j] = dbe->book1D(lbl.data(),title.data(),
phiNBins, phiMinim, phiMaxim);
lbl.clear();
title.clear();
lbl+=cLabel[j];lbl+="Etaphi";
title+=cLabel[j];title+=" ETA PHI OF COMPARISON FAILURES";
etaphi[j] = dbe->book2D(lbl.data(),title.data(),
etaNBins, etaMinim, etaMaxim,
phiNBins, phiMinim, phiMaxim
);
lbl.clear();
title.clear();
lbl+=cLabel[j];lbl+="Rank";
title+=cLabel[j];title+=" RANK OF COMPARISON FAILURES";
rnk[j] = dbe->book1D(lbl.data(),title.data(),
rnkNBins, rnkMinim, rnkMaxim);
//
lbl.clear();
title.clear();
lbl+=cLabel[j];lbl+="Eta"; lbl+="Data";
title+=cLabel[j];title+=" ETA (DATA)";
etaData[j] = dbe->book1D(lbl.data(),title.data(),
etaNBins, etaMinim, etaMaxim);
lbl.clear();
title.clear();
lbl+=cLabel[j];lbl+="Phi"; lbl+="Data";
title+=cLabel[j];title+=" PHI (DATA)";
phiData[j] = dbe->book1D(lbl.data(),title.data(),
phiNBins, phiMinim, phiMaxim);
lbl.clear();
title.clear();
lbl+=cLabel[j];lbl+="Rank"; lbl+="Data";
title+=cLabel[j];title+=" RANK (DATA)";
rnkData[j] = dbe->book1D(lbl.data(),title.data(),
rnkNBins, rnkMinim, rnkMaxim);
lbl.clear();
lbl+=cLabel[j];lbl+="Dword";
dword[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
lbl.clear();
lbl+=cLabel[j];lbl+="Eword";
eword[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
lbl.clear();
lbl+=cLabel[j];lbl+="DEword";
deword[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
//lbl.clear();
//lbl+=cLabel[j];lbl+="Masked";
//masked[j] = dbe->book1D(lbl.data(),lbl.data(),nbit,0,nbit);
}
}
std::string errLabel[nerr]= {
"Agree", "Loc. Agree", "L.Disagree", "Data only", "Emul only"
};
for(int i=0; i<nGctColl_; i++) {
sysrates ->setBinLabel(i+1,cLabel[i]);
sysncand[0]->setBinLabel(i+1,cLabel[i]);
sysncand[1]->setBinLabel(i+1,cLabel[i]);
}
for(int i=0; i<nGctColl_; i++) {
for(int j=0; j<nerr; j++) {
errortype[i]->setBinLabel(j+1,errLabel[j]);
}
}
for(int i=0; i<nGctColl_; i++) {
etaphi [i]->setAxisTitle("GCT #eta",1);
etaphi [i]->setAxisTitle("GCT #phi",2);
eta [i]->setAxisTitle("GCT #eta");
phi [i]->setAxisTitle("GCT #phi");
rnk [i]->setAxisTitle("Rank");
etaData[i]->setAxisTitle("GCT #eta");
phiData[i]->setAxisTitle("GCT #phi");
rnkData[i]->setAxisTitle("Rank");
dword [i]->setAxisTitle("trigger data word bit");
eword [i]->setAxisTitle("trigger data word bit");
deword [i]->setAxisTitle("trigger data word bit");
//masked [i]->setAxisTitle("trigger data word bit");
}
for(int i=0; i<nGctColl_; i++) {
colCount[i]=0;
nWithCol[i]=0;
}
if(verbose())
std::cout << "L1TdeGCT::beginJob() end.\n" << std::flush;
}
| void L1TdeGCT::endJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 205 of file L1TdeGCT.cc.
References gather_cfg::cout, and validate_alignment_devdb10_cfg::verbose.
| int L1TdeGCT::verbose | ( | ) | [inline, private] |
int L1TdeGCT::colCount[nGctColl_] [private] |
Definition at line 64 of file L1TdeGCT.h.
DQMStore* L1TdeGCT::dbe [private] |
Definition at line 57 of file L1TdeGCT.h.
edm::InputTag L1TdeGCT::DEsource_ [private] |
Definition at line 43 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::deword[nGctColl_] [private] |
Definition at line 83 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::dword[nGctColl_] [private] |
Definition at line 81 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::errortype[nGctColl_] [private] |
Definition at line 70 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::eta[nGctColl_] [private] |
Definition at line 73 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::etaData[nGctColl_] [private] |
Definition at line 76 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::etaphi[nGctColl_] [private] |
Definition at line 72 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::eword[nGctColl_] [private] |
Definition at line 82 of file L1TdeGCT.h.
bool L1TdeGCT::hasRecord_ [private] |
Definition at line 44 of file L1TdeGCT.h.
std::string L1TdeGCT::histFile_ [private] |
Definition at line 51 of file L1TdeGCT.h.
std::string L1TdeGCT::histFolder_ [private] |
Definition at line 54 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::masked[nGctColl_] [private] |
Definition at line 84 of file L1TdeGCT.h.
bool L1TdeGCT::monitorDaemon_ [private] |
Definition at line 58 of file L1TdeGCT.h.
const int L1TdeGCT::nGctColl_ = dedefs::GCThfbit-dedefs::GCTisolaem+1 [static, private] |
Definition at line 61 of file L1TdeGCT.h.
int L1TdeGCT::nWithCol[nGctColl_] [private] |
Definition at line 65 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::phi[nGctColl_] [private] |
Definition at line 74 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::phiData[nGctColl_] [private] |
Definition at line 77 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::rnk[nGctColl_] [private] |
Definition at line 75 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::rnkData[nGctColl_] [private] |
Definition at line 78 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::sysncand[2] [private] |
Definition at line 69 of file L1TdeGCT.h.
MonitorElement* L1TdeGCT::sysrates [private] |
Definition at line 68 of file L1TdeGCT.h.
int L1TdeGCT::verbose_ [private] |
Definition at line 47 of file L1TdeGCT.h.
Referenced by verbose().
1.7.3