|
|
|
#include <DTLocalTriggerLutTest.h>
Public Member Functions | |
| DTLocalTriggerLutTest (const edm::ParameterSet &ps) | |
| Constructor. | |
| virtual | ~DTLocalTriggerLutTest () |
| Destructor. | |
Protected Member Functions | |
| void | beginJob () |
| BeginJob. | |
| void | beginRun (const edm::Run &r, const edm::EventSetup &c) |
| BeginRun. | |
| void | runClientDiagnostic () |
| Run client analysis. | |
Private Member Functions | |
| void | fillWhPlot (MonitorElement *plot, int sect, int stat, float value, bool lessIsBest=true) |
| Fill summary plots managing double MB4 chambers. | |
| int | performLutTest (double mean, double RMS, double thresholdMean, double thresholdRMS) |
| Perform Lut Test logical operations. | |
Private Attributes | |
| bool | doCorrStudy |
| double | thresholdPhibMean |
| double | thresholdPhibRMS |
| double | thresholdPhiMean |
| double | thresholdPhiRMS |
* DQM Test Client
Definition at line 20 of file DTLocalTriggerLutTest.h.
| DTLocalTriggerLutTest::DTLocalTriggerLutTest | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 36 of file DTLocalTriggerLutTest.cc.
References edm::ParameterSet::getUntrackedParameter().
{
setConfig(ps,"DTLocalTriggerLut");
baseFolderDCC = "DT/03-LocalTrigger-DCC/";
baseFolderDDU = "DT/04-LocalTrigger-DDU/";
thresholdPhiMean = ps.getUntrackedParameter<double>("thresholdPhiMean",1.5);
thresholdPhiRMS = ps.getUntrackedParameter<double>("thresholdPhiRMS",.5);
thresholdPhibMean = ps.getUntrackedParameter<double>("thresholdPhibMean",1.5);
thresholdPhibRMS = ps.getUntrackedParameter<double>("thresholdPhibRMS",.8);
doCorrStudy = ps.getUntrackedParameter<bool>("doCorrelationStudy",false);
}
| DTLocalTriggerLutTest::~DTLocalTriggerLutTest | ( | ) | [virtual] |
| void DTLocalTriggerLutTest::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from DTLocalTriggerBaseTest.
Definition at line 56 of file DTLocalTriggerLutTest.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){
bookWheelHistos(wh,"PhiResidualMean");
bookWheelHistos(wh,"PhiResidualRMS");
bookWheelHistos(wh,"PhibResidualMean");
bookWheelHistos(wh,"PhibResidualRMS");
if (doCorrStudy) {
bookWheelHistos(wh,"PhiTkvsTrigSlope");
bookWheelHistos(wh,"PhiTkvsTrigIntercept");
bookWheelHistos(wh,"PhiTkvsTrigCorr");
bookWheelHistos(wh,"PhibTkvsTrigSlope");
bookWheelHistos(wh,"PhibTkvsTrigIntercept");
bookWheelHistos(wh,"PhibTkvsTrigCorr");
}
}
}
}
}
// 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){
bookWheelHistos(wh,"PhiLutSummary");
bookWheelHistos(wh,"PhibLutSummary");
}
bookCmsHistos("PhiLutSummary");
bookCmsHistos("PhibLutSummary");
}
}
}
| void DTLocalTriggerLutTest::beginRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
BeginRun.
Reimplemented from DTLocalTriggerBaseTest.
Definition at line 108 of file DTLocalTriggerLutTest.cc.
References DTLocalTriggerBaseTest::beginRun().
{
DTLocalTriggerBaseTest::beginRun(r,c);
}
| void DTLocalTriggerLutTest::fillWhPlot | ( | MonitorElement * | plot, |
| int | sect, | ||
| int | stat, | ||
| float | value, | ||
| bool | lessIsBest = true |
||
| ) | [private] |
Fill summary plots managing double MB4 chambers.
Definition at line 333 of file DTLocalTriggerLutTest.cc.
References MonitorElement::getBinContent(), and MonitorElement::setBinContent().
{
if (sect>12) {
int scsect = sect==13 ? 4 : 10;
if ( (fabs(value)>fabs(plot->getBinContent(scsect,stat)))==lessIsBest) {
plot->setBinContent(scsect,stat,value);
}
}
else {
plot->setBinContent(sect,stat,value);
}
return;
}
| int DTLocalTriggerLutTest::performLutTest | ( | double | mean, |
| double | RMS, | ||
| double | thresholdMean, | ||
| double | thresholdRMS | ||
| ) | [private] |
Perform Lut Test logical operations.
Definition at line 324 of file DTLocalTriggerLutTest.cc.
{
bool meanErr = fabs(mean)>thresholdMean;
bool rmsErr = RMS>thresholdRMS;
return (meanErr || rmsErr) ? 2+(meanErr!=rmsErr) : 0 ;
}
| void DTLocalTriggerLutTest::runClientDiagnostic | ( | ) | [protected, virtual] |
Run client analysis.
Implements DTLocalTriggerBaseTest.
Definition at line 115 of file DTLocalTriggerLutTest.cc.
References python::rootplot::argparse::category, spr::find(), newFWLiteAna::fullName, edm::second(), and DTChamberId::wheel().
{
// 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);
vector<DTChamber*>::const_iterator chIt = muonGeom->chambers().begin();
vector<DTChamber*>::const_iterator chEnd = muonGeom->chambers().end();
for (; chIt != chEnd; ++chIt) {
DTChamberId chId((*chIt)->id());
int wh = chId.wheel();
int sect = chId.sector();
int stat = chId.station();
if (doCorrStudy) {
// Perform Correlation Plots analysis (DCC + segment Phi)
TH2F * TrackPhitkvsPhitrig = getHisto<TH2F>(dbe->get(getMEName("PhitkvsPhitrig","Segment", chId)));
if (TrackPhitkvsPhitrig && TrackPhitkvsPhitrig->GetEntries()>10) {
// Fill client histos
if( whME[wh].find(fullName("PhiTkvsTrigCorr")) == whME[wh].end() ){
bookWheelHistos(wh,"PhiTkvsTrigSlope");
bookWheelHistos(wh,"PhiTkvsTrigIntercept");
bookWheelHistos(wh,"PhiTkvsTrigCorr");
}
TProfile* PhitkvsPhitrigProf = TrackPhitkvsPhitrig->ProfileX();
double phiInt = 0;
double phiSlope = 0;
double phiCorr = 0;
try {
PhitkvsPhitrigProf->Fit("pol1","CQO");
TF1 *ffPhi= PhitkvsPhitrigProf->GetFunction("pol1");
if (ffPhi) {
phiInt = ffPhi->GetParameter(0);
phiSlope = ffPhi->GetParameter(1);
phiCorr = TrackPhitkvsPhitrig->GetCorrelationFactor();
}
} catch (...) {
edm::LogError(category()) << "[" << testName << "Test]: Error fitting PhitkvsPhitrig for Wheel " << wh
<<" Sector " << sect << " Station " << stat;
}
std::map<std::string,MonitorElement*> &innerME = whME[wh];
fillWhPlot(innerME.find(fullName("PhiTkvsTrigSlope"))->second,sect,stat,phiSlope-1);
fillWhPlot(innerME.find(fullName("PhiTkvsTrigIntercept"))->second,sect,stat,phiInt);
fillWhPlot(innerME.find(fullName("PhiTkvsTrigCorr"))->second,sect,stat,phiCorr,false);
}
// Perform Correlation Plots analysis (DCC + segment Phib)
TH2F * TrackPhibtkvsPhibtrig = getHisto<TH2F>(dbe->get(getMEName("PhibtkvsPhibtrig","Segment", chId)));
if (stat != 3 && TrackPhibtkvsPhibtrig && TrackPhibtkvsPhibtrig->GetEntries()>10) {// station 3 has no meaningful MB3 phi bending information
// Fill client histos
if( whME[wh].find(fullName("PhibTkvsTrigCorr")) == whME[wh].end() ){
bookWheelHistos(wh,"PhibTkvsTrigSlope");
bookWheelHistos(wh,"PhibTkvsTrigIntercept");
bookWheelHistos(wh,"PhibTkvsTrigCorr");
}
TProfile* PhibtkvsPhibtrigProf = TrackPhibtkvsPhibtrig->ProfileX();
double phibInt = 0;
double phibSlope = 0;
double phibCorr = 0;
try {
PhibtkvsPhibtrigProf->Fit("pol1","CQO");
TF1 *ffPhib= PhibtkvsPhibtrigProf->GetFunction("pol1");
if (ffPhib) {
phibInt = ffPhib->GetParameter(0);
phibSlope = ffPhib->GetParameter(1);
phibCorr = TrackPhibtkvsPhibtrig->GetCorrelationFactor();
}
} catch (...) {
edm::LogError(category()) << "[" << testName << "Test]: Error fitting PhibtkvsPhibtrig for Wheel " << wh
<<" Sector " << sect << " Station " << stat;
}
std::map<std::string,MonitorElement*> &innerME = whME[wh];
fillWhPlot(innerME.find(fullName("PhibTkvsTrigSlope"))->second,sect,stat,phibSlope-1);
fillWhPlot(innerME.find(fullName("PhibTkvsTrigIntercept"))->second,sect,stat,phibInt);
fillWhPlot(innerME.find(fullName("PhibTkvsTrigCorr"))->second,sect,stat,phibCorr,false);
}
}
// Make Phi Residual Summary
TH1F * PhiResidual = getHisto<TH1F>(dbe->get(getMEName("PhiResidual","Segment", chId)));
int phiSummary = 1;
if (PhiResidual && PhiResidual->GetEffectiveEntries()>10) {
// Fill client histos
if( whME[wh].find(fullName("PhiResidualMean")) == whME[wh].end() ){
bookWheelHistos(wh,"PhiResidualMean");
bookWheelHistos(wh,"PhiResidualRMS");
}
double peak = PhiResidual->GetBinCenter(PhiResidual->GetMaximumBin());
double phiMean = 0;
double phiRMS = 0;
try {
PhiResidual->Fit("gaus","CQO","",peak-5,peak+5);
TF1 *ffPhi = PhiResidual->GetFunction("gaus");
if ( ffPhi ) {
phiMean = ffPhi->GetParameter(1);
phiRMS = ffPhi->GetParameter(2);
}
} catch (...) {
edm::LogError(category()) << "[" << testName << "Test]: Error fitting PhiResidual for Wheel " << wh
<<" Sector " << sect << " Station " << stat;
}
std::map<std::string,MonitorElement*> &innerME = whME[wh];
fillWhPlot(innerME.find(fullName("PhiResidualMean"))->second,sect,stat,phiMean);
fillWhPlot(innerME.find(fullName("PhiResidualRMS"))->second,sect,stat,phiRMS);
phiSummary = performLutTest(phiMean,phiRMS,thresholdPhiMean,thresholdPhiRMS);
}
fillWhPlot(whME[wh].find(fullName("PhiLutSummary"))->second,sect,stat,phiSummary);
// Make Phib Residual Summary
TH1F * PhibResidual = getHisto<TH1F>(dbe->get(getMEName("PhibResidual","Segment", chId)));
int phibSummary = stat==3 ? 0 : 1; // station 3 has no meaningful MB3 phi bending information
if (stat != 3 && PhibResidual && PhibResidual->GetEffectiveEntries()>10) {// station 3 has no meaningful MB3 phi bending information
// Fill client histos
if( whME[wh].find(fullName("PhibResidualMean")) == whME[wh].end() ){
bookWheelHistos(wh,"PhibResidualMean");
bookWheelHistos(wh,"PhibResidualRMS");
}
double peak = PhibResidual->GetBinCenter(PhibResidual->GetMaximumBin());
double phibMean = 0;
double phibRMS = 0;
try {
PhibResidual->Fit("gaus","CQO","",peak-5,peak+5);
TF1 *ffPhib = PhibResidual->GetFunction("gaus");
if ( ffPhib ) {
phibMean = ffPhib->GetParameter(1);
phibRMS = ffPhib->GetParameter(2);
}
} catch (...) {
edm::LogError(category()) << "[" << testName << "Test]: Error fitting PhibResidual for Wheel " << wh
<<" Sector " << sect << " Station " << stat;
}
std::map<std::string,MonitorElement*> &innerME = whME[wh];
fillWhPlot(innerME.find(fullName("PhibResidualMean"))->second,sect,stat,phibMean);
fillWhPlot(innerME.find(fullName("PhibResidualRMS"))->second,sect,stat,phibRMS);
phibSummary = performLutTest(phibMean,phibRMS,thresholdPhibMean,thresholdPhibRMS);
}
fillWhPlot(whME[wh].find(fullName("PhibLutSummary"))->second,sect,stat,phibSummary);
}
}
}
// Barrel Summary Plots
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]);
TH2F* phiWhSummary = getHisto<TH2F>(innerME->find(fullName("PhiLutSummary"))->second);
TH2F* phibWhSummary = getHisto<TH2F>(innerME->find(fullName("PhibLutSummary"))->second);
for (int sect=1; sect<=12; ++sect){
int phiErr = 0;
int phibErr = 0;
int phiNoData = 0;
int phibNoData = 0;
for (int stat=1; stat<=4; ++stat){
switch (static_cast<int>(phiWhSummary->GetBinContent(sect,stat))) {
case 1:
phiNoData++;
case 2:
case 3:
phiErr++;
}
switch (static_cast<int>(phibWhSummary->GetBinContent(sect,stat))) {
case 1:
phibNoData++;
case 2:
case 3:
phibErr++;
}
}
if (phiNoData == 4) phiErr = 5;
if (phibNoData == 3) phibErr = 5; // MB3 has no phib information
cmsME.find(fullName("PhiLutSummary"))->second->setBinContent(sect,wh+3,phiErr);
cmsME.find(fullName("PhibLutSummary"))->second->setBinContent(sect,wh+3,phibErr);
}
}
}
}
}
bool DTLocalTriggerLutTest::doCorrStudy [private] |
Definition at line 52 of file DTLocalTriggerLutTest.h.
double DTLocalTriggerLutTest::thresholdPhibMean [private] |
Definition at line 50 of file DTLocalTriggerLutTest.h.
double DTLocalTriggerLutTest::thresholdPhibRMS [private] |
Definition at line 51 of file DTLocalTriggerLutTest.h.
double DTLocalTriggerLutTest::thresholdPhiMean [private] |
Definition at line 50 of file DTLocalTriggerLutTest.h.
double DTLocalTriggerLutTest::thresholdPhiRMS [private] |
Definition at line 51 of file DTLocalTriggerLutTest.h.
1.7.3