|
|
|
#include <DTLocalTriggerTest.h>
Public Member Functions | |
| DTLocalTriggerTest (const edm::ParameterSet &ps) | |
| Constructor. | |
| virtual | ~DTLocalTriggerTest () |
| Destructor. | |
Protected Member Functions | |
| void | beginJob () |
| BeginJob. | |
| void | beginRun (const edm::Run &r, const edm::EventSetup &c) |
| BeginRun. | |
| void | fillGlobalSummary () |
| void | runClientDiagnostic () |
| Run client analysis. | |
Private Attributes | |
| int | nMinEvts |
* DQM Test Client
Definition at line 20 of file DTLocalTriggerTest.h.
| DTLocalTriggerTest::DTLocalTriggerTest | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 36 of file DTLocalTriggerTest.cc.
References edm::ParameterSet::getUntrackedParameter().
{
setConfig(ps,"DTLocalTrigger");
baseFolderDCC = "DT/03-LocalTrigger-DCC/";
baseFolderDDU = "DT/04-LocalTrigger-DDU/";
nMinEvts = ps.getUntrackedParameter<int>("nEventsCert", 5000);
}
| DTLocalTriggerTest::~DTLocalTriggerTest | ( | ) | [virtual] |
| void DTLocalTriggerTest::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from DTLocalTriggerBaseTest.
Definition at line 50 of file DTLocalTriggerTest.cc.
References bk::beginJob(), and Parameters::parameters.
{
DTLocalTriggerBaseTest::beginJob();
vector<string>::const_iterator iTr = trigSources.begin();
vector<string>::const_iterator trEnd = trigSources.end();
vector<string>::const_iterator iHw = hwSources.begin();
vector<string>::const_iterator hwEnd = hwSources.end();
//Booking
if(parameters.getUntrackedParameter<bool>("staticBooking", true)){
for (; iTr != trEnd; ++iTr){
trigSource = (*iTr);
for (; iHw != hwEnd; ++iHw){
hwSource = (*iHw);
// Loop over the TriggerUnits
for (int wh=-2; wh<=2; ++wh){
if (hwSource=="COM") {
bookWheelHistos(wh,"MatchingPhi");
}
else {
for (int sect=1; sect<=12; ++sect){
bookSectorHistos(wh,sect,"BXDistribPhi");
bookSectorHistos(wh,sect,"QualDistribPhi");
}
bookWheelHistos(wh,"CorrectBXPhi");
bookWheelHistos(wh,"ResidualBXPhi");
bookWheelHistos(wh,"CorrFractionPhi");
bookWheelHistos(wh,"2ndFractionPhi");
bookWheelHistos(wh,"TriggerInclusivePhi");
bookWheelHistos(wh,"CorrectBXTheta");
if (hwSource=="DDU") {
bookWheelHistos(wh,"HFractionTheta");
}
}
}
}
}
}
// Summary test histo booking (only static)
for (iTr = trigSources.begin(); iTr != trEnd; ++iTr){
trigSource = (*iTr);
for (iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
hwSource = (*iHw);
// Loop over the TriggerUnits
for (int wh=-2; wh<=2; ++wh){
if (hwSource=="COM") {
bookWheelHistos(wh,"MatchingSummary","Summaries");
}
else {
bookWheelHistos(wh,"CorrFractionSummary","Summaries");
bookWheelHistos(wh,"2ndFractionSummary","Summaries");
}
}
if (hwSource=="COM") {
bookCmsHistos("MatchingSummary","Summaries");
}
else {
bookCmsHistos("CorrFractionSummary");
bookCmsHistos("2ndFractionSummary");
}
if (hwSource=="DCC") {
bookCmsHistos("TrigGlbSummary","",true);
}
}
}
}
| void DTLocalTriggerTest::beginRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
BeginRun.
Reimplemented from DTLocalTriggerBaseTest.
Definition at line 123 of file DTLocalTriggerTest.cc.
References DTLocalTriggerBaseTest::beginRun().
{
DTLocalTriggerBaseTest::beginRun(r,c);
}
| void DTLocalTriggerTest::fillGlobalSummary | ( | ) | [protected] |
Definition at line 392 of file DTLocalTriggerTest.cc.
References python::rootplot::argparse::category, corr, newFWLiteAna::fullName, MonitorElement::getBinContent(), MonitorElement::getFloatValue(), LogTrace, lumiPlot::lut, edm::min(), and edm::second().
{
float glbPerc[5] = { 1., 0.9, 0.6, 0.3, 0.01 };
trigSource = "";
hwSource = "DCC";
int nSecReadout = 0;
for (int wh=-2; wh<=2; ++wh) {
for (int sect=1; sect<=12; ++sect) {
float maxErr = 8.;
int corr = cmsME.find(fullName("CorrFractionSummary"))->second->getBinContent(sect,wh+3);
int second = cmsME.find(fullName("2ndFractionSummary"))->second->getBinContent(sect,wh+3);
int lut=0;
MonitorElement * lutsME = dbe->get(topFolder(hwSource=="DCC") + "Summaries/TrigLutSummary");
if (lutsME) {
lut = lutsME->getBinContent(sect,wh+3);
maxErr+=4;
} else {
LogTrace(category()) << "[" << testName
<< "Test]: DCC Lut test Summary histo not found." << endl;
}
(corr <5 || second<5) && nSecReadout++;
int errcode = ((corr<5 ? corr : 4) + (second<5 ? second : 4) + (lut<5 ? lut : 4) );
errcode = min(int((errcode/maxErr + 0.01)*5),5);
cmsME.find("TrigGlbSummary")->second->setBinContent(sect,wh+3,glbPerc[errcode]);
}
}
if (!nSecReadout)
cmsME.find("TrigGlbSummary")->second->Reset(); // white histo id DCC is not RO
string nEvtsName = "DT/EventInfo/Counters/nProcessedEventsTrigger";
MonitorElement * meProcEvts = dbe->get(nEvtsName);
if (meProcEvts) {
int nProcEvts = meProcEvts->getFloatValue();
cmsME.find("TrigGlbSummary")->second->setEntries(nProcEvts < nMinEvts ? 10. : nProcEvts);
} else {
cmsME.find("TrigGlbSummary")->second->setEntries(nMinEvts + 1);
LogVerbatim (category()) << "[" << testName
<< "Test]: ME: " << nEvtsName << " not found!" << endl;
}
}
| void DTLocalTriggerTest::runClientDiagnostic | ( | ) | [protected, virtual] |
Run client analysis.
Implements DTLocalTriggerBaseTest.
Definition at line 130 of file DTLocalTriggerTest.cc.
References end, python::tagInventory::entries, spr::find(), newFWLiteAna::fullName, Parameters::parameters, and testRegression::stat.
{
// Loop over Trig & Hw sources
for (vector<string>::const_iterator iTr = trigSources.begin(); iTr != trigSources.end(); ++iTr){
trigSource = (*iTr);
for (vector<string>::const_iterator iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
hwSource = (*iHw);
// Loop over the TriggerUnits
for (int stat=1; stat<=4; ++stat){
for (int wh=-2; wh<=2; ++wh){
for (int sect=1; sect<=12; ++sect){
DTChamberId chId(wh,stat,sect);
int sector_id = (wh+3)+(sect-1)*5;
// uint32_t indexCh = chId.rawId();
if (hwSource=="COM") {
// Perform DCC-DDU matching test and generates summaries (Phi view)
TH2F * DDUvsDCC = getHisto<TH2F>(dbe->get(getMEName("QualDDUvsQualDCC","LocalTriggerPhi", chId)));
if (DDUvsDCC) {
int matchSummary = 1;
if (DDUvsDCC->GetEntries()>1) {
double entries = DDUvsDCC->GetEntries();
double corrEntries = 0;
for (int ibin=2; ibin<=8; ++ibin) {
corrEntries += DDUvsDCC->GetBinContent(ibin,ibin);
}
double corrRatio = corrEntries/entries;
if (corrRatio < parameters.getUntrackedParameter<double>("matchingFracError",.65)){
matchSummary = 2;
}
else if (corrRatio < parameters.getUntrackedParameter<double>("matchingFracWarning",.85)){
matchSummary = 3;
}
else {
matchSummary = 0;
}
if( whME[wh].find(fullName("MatchingPhi")) == whME[wh].end() ){
bookWheelHistos(wh,"MatchingPhi");
}
whME[wh].find(fullName("MatchingPhi"))->second->setBinContent(sect,stat,corrRatio);
}
whME[wh].find(fullName("MatchingSummary"))->second->setBinContent(sect,stat,matchSummary);
}
}
else {
// Perform DCC/DDU common plot analysis (Phi ones)
TH2F * BXvsQual = getHisto<TH2F>(dbe->get(getMEName("BXvsQual","LocalTriggerPhi", chId)));
TH1F * BestQual = getHisto<TH1F>(dbe->get(getMEName("BestQual","LocalTriggerPhi", chId)));
TH2F * Flag1stvsQual = getHisto<TH2F>(dbe->get(getMEName("Flag1stvsQual","LocalTriggerPhi", chId)));
if (BXvsQual && Flag1stvsQual && BestQual) {
int corrSummary = 1;
int secondSummary = 1;
if (BestQual->GetEntries()>1) {
TH1D* BXHH = BXvsQual->ProjectionY("",6,7,"");
TH1D* Flag1st = Flag1stvsQual->ProjectionY();
int BXOK_bin = BXHH->GetEntries()>=1 ? BXHH->GetMaximumBin() : 51;
double BXMean = BXHH->GetEntries()>=1 ? BXHH->GetMean() : 51;
double BX_OK = BXvsQual->GetYaxis()->GetBinCenter(BXOK_bin);
double trigsFlag2nd = Flag1st->GetBinContent(2);
double trigs = Flag1st->GetEntries();
double besttrigs = BestQual->GetEntries();
double besttrigsCorr = BestQual->Integral(5,7,"");
delete BXHH;
delete Flag1st;
double corrFrac = besttrigsCorr/besttrigs;
double secondFrac = trigsFlag2nd/trigs;
if (corrFrac < parameters.getUntrackedParameter<double>("corrFracError",.5)){
corrSummary = 2;
}
else if (corrFrac < parameters.getUntrackedParameter<double>("corrFracWarning",.6)){
corrSummary = 3;
}
else {
corrSummary = 0;
}
if (secondFrac > parameters.getUntrackedParameter<double>("secondFracError",.2)){
secondSummary = 2;
}
else if (secondFrac > parameters.getUntrackedParameter<double>("secondFracWarning",.1)){
secondSummary = 3;
}
else {
secondSummary = 0;
}
if( secME[sector_id].find(fullName("BXDistribPhi")) == secME[sector_id].end() ){
bookSectorHistos(wh,sect,"QualDistribPhi");
bookSectorHistos(wh,sect,"BXDistribPhi");
}
TH1D* BXDistr = BXvsQual->ProjectionY();
TH1D* QualDistr = BXvsQual->ProjectionX();
std::map<std::string,MonitorElement*> *innerME = &(secME[sector_id]);
int nbinsBX = BXDistr->GetNbinsX();
int firstBinCenter = static_cast<int>(BXDistr->GetBinCenter(1));
int lastBinCenter = static_cast<int>(BXDistr->GetBinCenter(nbinsBX));
int iMin = firstBinCenter>-4 ? firstBinCenter : -4;
int iMax = lastBinCenter<20 ? lastBinCenter : 20;
for (int ibin=iMin+5;ibin<=iMax+5; ++ibin) {
innerME->find(fullName("BXDistribPhi"))->second->setBinContent(ibin,stat,BXDistr->GetBinContent(ibin-5-firstBinCenter+1));
}
for (int ibin=1;ibin<=7;++ibin) {
innerME->find(fullName("QualDistribPhi"))->second->setBinContent(ibin,stat,QualDistr->GetBinContent(ibin));
}
delete BXDistr;
delete QualDistr;
if( whME[wh].find(fullName("CorrectBXPhi")) == whME[wh].end() ){
bookWheelHistos(wh,"ResidualBXPhi");
bookWheelHistos(wh,"CorrectBXPhi");
bookWheelHistos(wh,"CorrFractionPhi");
bookWheelHistos(wh,"2ndFractionPhi");
bookWheelHistos(wh,"TriggerInclusivePhi");
}
innerME = &(whME[wh]);
innerME->find(fullName("CorrectBXPhi"))->second->setBinContent(sect,stat,BX_OK+0.00001);
innerME->find(fullName("ResidualBXPhi"))->second->setBinContent(sect,stat,round(25.*(BXMean-BX_OK))+0.00001);
innerME->find(fullName("CorrFractionPhi"))->second->setBinContent(sect,stat,corrFrac);
innerME->find(fullName("TriggerInclusivePhi"))->second->setBinContent(sect,stat,besttrigs);
innerME->find(fullName("2ndFractionPhi"))->second->setBinContent(sect,stat,secondFrac);
}
whME[wh].find(fullName("CorrFractionSummary"))->second->setBinContent(sect,stat,corrSummary);
whME[wh].find(fullName("2ndFractionSummary"))->second->setBinContent(sect,stat,secondSummary);
}
if (hwSource=="DDU") {
// Perform DDU plot analysis (Theta ones)
TH2F * ThetaBXvsQual = getHisto<TH2F>(dbe->get(getMEName("ThetaBXvsQual","LocalTriggerTheta", chId)));
TH1F * ThetaBestQual = getHisto<TH1F>(dbe->get(getMEName("ThetaBestQual","LocalTriggerTheta", chId)));
// no theta triggers in stat 4!
if (ThetaBXvsQual && ThetaBestQual && stat<4 && ThetaBestQual->GetEntries()>1) {
TH1D* BXH = ThetaBXvsQual->ProjectionY("",4,4,"");
int BXOK_bin = BXH->GetEffectiveEntries()>=1 ? BXH->GetMaximumBin(): 10;
double BX_OK = ThetaBXvsQual->GetYaxis()->GetBinCenter(BXOK_bin);
double trigs = ThetaBestQual->GetEntries();
double trigsH = ThetaBestQual->GetBinContent(4);
delete BXH;
// if( secME[sector_id].find(fullName("HFractionTheta")) == secME[sector_id].end() ){
// // bookSectorHistos(wh,sect,"CorrectBXTheta");
// bookSectorHistos(wh,sect,"HFractionTheta");
// }
//std::map<std::string,MonitorElement*> *innerME = &(secME.find(sector_id)->second);
// innerME->find(fullName("CorrectBXTheta"))->second->setBinContent(stat,BX_OK);
//innerME->find(fullName("HFractionTheta"))->second->setBinContent(stat,trigsH/trigs);
if( whME[wh].find(fullName("HFractionTheta")) == whME[wh].end() ){
bookWheelHistos(wh,"CorrectBXTheta");
bookWheelHistos(wh,"HFractionTheta");
}
std::map<std::string,MonitorElement*> *innerME = &(whME.find(wh)->second);
innerME->find(fullName("CorrectBXTheta"))->second->setBinContent(sect,stat,BX_OK+0.00001);
innerME->find(fullName("HFractionTheta"))->second->setBinContent(sect,stat,trigsH/trigs);
}
}
else if (hwSource=="DCC") {
// Perform DCC plot analysis (Theta ones)
TH2F * ThetaPosvsBX = getHisto<TH2F>(dbe->get(getMEName("PositionvsBX","LocalTriggerTheta", chId)));
// no theta triggers in stat 4!
if (ThetaPosvsBX && stat<4 && ThetaPosvsBX->GetEntries()>1) {
TH1D* BX = ThetaPosvsBX->ProjectionX();
int BXOK_bin = BX->GetEffectiveEntries()>=1 ? BX->GetMaximumBin(): 10;
double BX_OK = ThetaPosvsBX->GetXaxis()->GetBinCenter(BXOK_bin);
delete BX;
if( whME[wh].find(fullName("CorrectBXTheta")) == whME[wh].end() ){
bookWheelHistos(wh,"CorrectBXTheta");
}
std::map<std::string,MonitorElement*> *innerME = &(whME.find(wh)->second);
innerME->find(fullName("CorrectBXTheta"))->second->setBinContent(sect,stat,BX_OK+0.00001);
}
}
}
}
}
}
}
}
for (vector<string>::const_iterator iTr = trigSources.begin(); iTr != trigSources.end(); ++iTr){
trigSource = (*iTr);
for (vector<string>::const_iterator iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
hwSource = (*iHw);
for (int wh=-2; wh<=2; ++wh){
std::map<std::string,MonitorElement*> *innerME = &(whME[wh]);
if(hwSource=="COM") {
TH2F* matchWhSummary = getHisto<TH2F>(innerME->find(fullName("MatchingSummary"))->second);
for (int sect=1; sect<=12; ++sect){
int matchErr = 0;
int matchNoData = 0;
for (int stat=1; stat<=4; ++stat){
switch (static_cast<int>(matchWhSummary->GetBinContent(sect,stat))) {
case 1:
matchNoData++;
case 2:
matchErr++;
}
}
if (matchNoData == 4) matchErr = 5;
cmsME.find(fullName("MatchingSummary"))->second->setBinContent(sect,wh+3,matchErr);
}
}
else {
TH2F* corrWhSummary = getHisto<TH2F>(innerME->find(fullName("CorrFractionSummary"))->second);
TH2F* secondWhSummary = getHisto<TH2F>(innerME->find(fullName("2ndFractionSummary"))->second);
for (int sect=1; sect<=12; ++sect){
int corrErr = 0;
int secondErr = 0;
int corrNoData = 0;
int secondNoData = 0;
for (int stat=1; stat<=4; ++stat){
switch (static_cast<int>(corrWhSummary->GetBinContent(sect,stat))) {
case 1:
corrNoData++;
case 2:
corrErr++;
}
switch (static_cast<int>(secondWhSummary->GetBinContent(sect,stat))) {
case 1:
secondNoData++;
case 2:
secondErr++;
}
}
if (corrNoData == 4) corrErr = 5;
if (secondNoData == 4) secondErr = 5;
cmsME.find(fullName("CorrFractionSummary"))->second->setBinContent(sect,wh+3,corrErr);
cmsME.find(fullName("2ndFractionSummary"))->second->setBinContent(sect,wh+3,secondErr);
}
}
}
}
}
fillGlobalSummary();
}
int DTLocalTriggerTest::nMinEvts [private] |
Definition at line 45 of file DTLocalTriggerTest.h.
1.7.3