|
|
|
#include <MuonTestSummary.h>
* DQM Client for global summary
Definition at line 34 of file MuonTestSummary.h.
| MuonTestSummary::MuonTestSummary | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 30 of file MuonTestSummary.cc.
References ExpressReco_HICollisions_FallBack::chi2Fraction, ExpressReco_HICollisions_FallBack::chi2Spread, ExpressReco_HICollisions_FallBack::etaExpected, ExpressReco_HICollisions_FallBack::expMultiplicityGlb_max, ExpressReco_HICollisions_FallBack::expMultiplicityGlb_min, ExpressReco_HICollisions_FallBack::expMultiplicitySta_max, ExpressReco_HICollisions_FallBack::expMultiplicitySta_min, ExpressReco_HICollisions_FallBack::expMultiplicityTk_max, ExpressReco_HICollisions_FallBack::expMultiplicityTk_min, ExpressReco_HICollisions_FallBack::expPeakEcalS9_max, ExpressReco_HICollisions_FallBack::expPeakEcalS9_min, ExpressReco_HICollisions_FallBack::expPeakHadS9_max, ExpressReco_HICollisions_FallBack::expPeakHadS9_min, edm::ParameterSet::getParameter(), ExpressReco_HICollisions_FallBack::matchesFractionCsc_max, ExpressReco_HICollisions_FallBack::matchesFractionCsc_min, ExpressReco_HICollisions_FallBack::matchesFractionDt_max, ExpressReco_HICollisions_FallBack::matchesFractionDt_min, ExpressReco_HICollisions_FallBack::numMatchedExpected_max, ExpressReco_HICollisions_FallBack::numMatchedExpected_min, cmsCodeRules::cppFunctionSkipper::operator, ExpressReco_HICollisions_FallBack::phiExpected, ExpressReco_HICollisions_FallBack::pullEtaSpread, ExpressReco_HICollisions_FallBack::pullOneOvPSpread, ExpressReco_HICollisions_FallBack::pullPhiSpread, ExpressReco_HICollisions_FallBack::resChargeLimit_glbSta, ExpressReco_HICollisions_FallBack::resChargeLimit_tkGlb, ExpressReco_HICollisions_FallBack::resChargeLimit_tkSta, ExpressReco_HICollisions_FallBack::resEtaSpread_glbSta, ExpressReco_HICollisions_FallBack::resEtaSpread_tkGlb, ExpressReco_HICollisions_FallBack::resOneOvPSpread_glbSta, ExpressReco_HICollisions_FallBack::resOneOvPSpread_tkGlb, ExpressReco_HICollisions_FallBack::resPhiSpread_glbSta, ExpressReco_HICollisions_FallBack::resPhiSpread_tkGlb, ExpressReco_HICollisions_FallBack::resSegmTrack_mean_max, ExpressReco_HICollisions_FallBack::resSegmTrack_mean_min, ExpressReco_HICollisions_FallBack::resSegmTrack_rms_max, and ExpressReco_HICollisions_FallBack::resSegmTrack_rms_min.
{
dbe = Service<DQMStore>().operator->();
// parameter initialization for kinematics test
etaExpected = ps.getParameter<double>("etaExpected");
phiExpected = ps.getParameter<double>("phiExpected");
chi2Fraction = ps.getParameter<double>("chi2Fraction");
chi2Spread = ps.getParameter<double>("chi2Spread");
resEtaSpread_tkGlb = ps.getParameter<double>("resEtaSpread_tkGlb");
resEtaSpread_glbSta = ps.getParameter<double>("resEtaSpread_glbSta");
resPhiSpread_tkGlb = ps.getParameter<double>("resPhiSpread_tkGlb");
resPhiSpread_glbSta = ps.getParameter<double>("resPhiSpread_glbSta");
resOneOvPSpread_tkGlb = ps.getParameter<double>("resOneOvPSpread_tkGlb");
resOneOvPSpread_glbSta = ps.getParameter<double>("resOneOvPSpread_glbSta");
pullEtaSpread = ps.getParameter<double>("pullEtaSpread");
pullPhiSpread = ps.getParameter<double>("pullPhiSpread");
pullOneOvPSpread = ps.getParameter<double>("pullOneOvPSpread");
resChargeLimit_tkGlb = ps.getParameter<double>("resChargeLimit_tkGlb");
resChargeLimit_glbSta = ps.getParameter<double>("resChargeLimit_glbSta");
resChargeLimit_tkSta = ps.getParameter<double>("resChargeLimit_tkSta");
numMatchedExpected_min = ps.getParameter<double>("numMatchedExpected_min");
numMatchedExpected_max = ps.getParameter<double>("numMatchedExpected_max");
matchesFractionDt_min = ps.getParameter<double>("matchesFractionDt_min");
matchesFractionDt_max = ps.getParameter<double>("matchesFractionDt_max");
matchesFractionCsc_min = ps.getParameter<double>("matchesFractionCsc_min");
matchesFractionCsc_max = ps.getParameter<double>("matchesFractionCsc_max");
resSegmTrack_rms_min = ps.getParameter<double>("resSegmTrack_rms_min");
resSegmTrack_rms_max = ps.getParameter<double>("resSegmTrack_rms_max");
resSegmTrack_mean_min = ps.getParameter<double>("resSegmTrack_mean_min");
resSegmTrack_mean_max = ps.getParameter<double>("resSegmTrack_mean_max");
expPeakEcalS9_min = ps.getParameter<double>("expPeakEcalS9_min");
expPeakEcalS9_max = ps.getParameter<double>("expPeakEcalS9_max");
expPeakHadS9_min = ps.getParameter<double>("expPeakHadS9_min");
expPeakHadS9_max = ps.getParameter<double>("expPeakHadS9_max");
expMultiplicityGlb_max = ps.getParameter<double>("expMultiplicityGlb_max");
expMultiplicityTk_max = ps.getParameter<double>("expMultiplicityTk_max");
expMultiplicitySta_max = ps.getParameter<double>("expMultiplicitySta_max");
expMultiplicityGlb_min = ps.getParameter<double>("expMultiplicityGlb_min");
expMultiplicityTk_min = ps.getParameter<double>("expMultiplicityTk_min");
expMultiplicitySta_min = ps.getParameter<double>("expMultiplicitySta_min");
}
| MuonTestSummary::~MuonTestSummary | ( | ) | [virtual] |
| void MuonTestSummary::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [inline, protected, virtual] |
| void MuonTestSummary::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 76 of file MuonTestSummary.cc.
References LogTrace, and metname.
{
metname = "muonTestSummary";
LogTrace(metname)<<"[MuonTestSummary] beginJob: Histo booking";
// book the summary histos
dbe->setCurrentFolder("Muons/TestSummary");
// kinematics test report
kinematicsSummaryMap = dbe->book2D("kinematicsSummaryMap","Kinematics test summary",5,1,6,3,1,4);
kinematicsSummaryMap->setAxisTitle("track monitored",1);
kinematicsSummaryMap->setBinLabel(1,"GLB",1);
kinematicsSummaryMap->setBinLabel(2,"TKfromGLB",1);
kinematicsSummaryMap->setBinLabel(3,"STAfromGLB",1);
kinematicsSummaryMap->setBinLabel(4,"TK",1);
kinematicsSummaryMap->setBinLabel(5,"STA",1);
kinematicsSummaryMap->setAxisTitle("parameter tested",2);
kinematicsSummaryMap->setBinLabel(1,"#chi^{2}",2);
kinematicsSummaryMap->setBinLabel(2,"#eta",2);
kinematicsSummaryMap->setBinLabel(3,"#phi",2);
//chi2 kinematics quality test report
chi2TestSummaryMap = dbe->book2D("chi2TestSummaryMap","#chi2 quality test summary",5,1,6,5,1,6);
chi2TestSummaryMap->setAxisTitle("track monitored",1);
chi2TestSummaryMap->setBinLabel(1,"GLB",1);
chi2TestSummaryMap->setBinLabel(2,"TKfromGLB",1);
chi2TestSummaryMap->setBinLabel(3,"STAfromGLB",1);
chi2TestSummaryMap->setBinLabel(4,"TK",1);
chi2TestSummaryMap->setBinLabel(5,"STA",1);
chi2TestSummaryMap->setAxisTitle("parameter tested",2);
chi2TestSummaryMap->setBinLabel(1,"#chi^{2}",2);
chi2TestSummaryMap->setBinLabel(2,"#eta",2);
chi2TestSummaryMap->setBinLabel(3,"#phi",2);
chi2TestSummaryMap->setBinLabel(4,"#pt",2);
chi2TestSummaryMap->setBinLabel(5,"#q",2);
//Kolmogorov kinematics quality test report
KolmogorovTestSummaryMap = dbe->book2D("KolmogorovTestSummaryMap","Kolmogorov quality test summary",5,1,6,5,1,6);
KolmogorovTestSummaryMap->setAxisTitle("track monitored",1);
KolmogorovTestSummaryMap->setBinLabel(1,"GLB",1);
KolmogorovTestSummaryMap->setBinLabel(2,"TKfromGLB",1);
KolmogorovTestSummaryMap->setBinLabel(3,"STAfromGLB",1);
KolmogorovTestSummaryMap->setBinLabel(4,"TK",1);
KolmogorovTestSummaryMap->setBinLabel(5,"STA",1);
KolmogorovTestSummaryMap->setAxisTitle("parameter tested",2);
KolmogorovTestSummaryMap->setBinLabel(1,"#chi^{2}",2);
KolmogorovTestSummaryMap->setBinLabel(2,"#eta",2);
KolmogorovTestSummaryMap->setBinLabel(3,"#phi",2);
KolmogorovTestSummaryMap->setBinLabel(4,"#pt",2);
KolmogorovTestSummaryMap->setBinLabel(5,"#q",2);
// residuals test report
residualsSummaryMap = dbe->book2D("residualsSummaryMap","Residuals test summary",4,1,5,4,1,5);
residualsSummaryMap->setAxisTitle("residuals",1);
residualsSummaryMap->setBinLabel(1,"TK-GLB",1);
residualsSummaryMap->setBinLabel(2,"GLB-STA",1);
residualsSummaryMap->setBinLabel(3,"TK-STA",1);
residualsSummaryMap->setBinLabel(4,"TK-STA Pull",1);
residualsSummaryMap->setAxisTitle("parameter tested",2);
residualsSummaryMap->setBinLabel(1,"#eta",2);
residualsSummaryMap->setBinLabel(2,"#phi",2);
residualsSummaryMap->setBinLabel(3,"1/p",2);
residualsSummaryMap->setBinLabel(4,"q",2);
// muonId test report
muonIdSummaryMap = dbe->book2D("muonIdSummaryMap","muonId test summary",4,1,5, 5,1,6);
muonIdSummaryMap->setAxisTitle("muons",1);
muonIdSummaryMap->setBinLabel(1,"GLB DT",1);
muonIdSummaryMap->setBinLabel(2,"GLB CSC",1);
muonIdSummaryMap->setBinLabel(3,"TK DT",1);
muonIdSummaryMap->setBinLabel(4,"TK CSC",1);
muonIdSummaryMap->setAxisTitle("tests",2);
muonIdSummaryMap->setBinLabel(1,"#assSeg",2);
muonIdSummaryMap->setBinLabel(2,"x mean",2);
muonIdSummaryMap->setBinLabel(3,"x rms",2);
muonIdSummaryMap->setBinLabel(4,"y mean",2);
muonIdSummaryMap->setBinLabel(5,"y rms",2);
// energy test report
energySummaryMap = dbe->book2D("energySummaryMap","Energy deposits test summary",3,1,4,3,1,4);
energySummaryMap->setAxisTitle("muons",1);
energySummaryMap->setBinLabel(1,"GLB",1);
energySummaryMap->setBinLabel(2,"TK",1);
energySummaryMap->setBinLabel(3,"STA",1);
energySummaryMap->setAxisTitle("calorimeter tested",2);
energySummaryMap->setBinLabel(1,"ECAL",2);
energySummaryMap->setBinLabel(2,"HAD",2);
energySummaryMap->setBinLabel(3,"H0",2);
// multiplicity tests report
multiplicitySummaryMap = dbe->book1D("multiplicitySummaryMap","muon multiplicity test summary",3,1,4);
multiplicitySummaryMap->setAxisTitle("muon");
multiplicitySummaryMap->setBinLabel(1,"GLB");
multiplicitySummaryMap->setBinLabel(2,"TK");
multiplicitySummaryMap->setBinLabel(3,"STA");
// summary test report
dbe->setCurrentFolder("Muons/EventInfo");
summaryReport = dbe->bookFloat("reportSummary");
summaryReportMap = dbe->book2D("reportSummaryMap","Muon Report Summary Map",3,1,4,7,1,8);
summaryReportMap->setAxisTitle("muons",1);
summaryReportMap->setBinLabel(1,"GLB",1);
summaryReportMap->setBinLabel(2,"TK",1);
summaryReportMap->setBinLabel(3,"STA",1);
summaryReportMap->setAxisTitle("test",2);
summaryReportMap->setBinLabel(1,"#chi^{2}/Df",2);
summaryReportMap->setBinLabel(2,"#eta",2);
summaryReportMap->setBinLabel(3,"#phi",2);
summaryReportMap->setBinLabel(4,"residuals",2);
summaryReportMap->setBinLabel(5,"muonId",2);
summaryReportMap->setBinLabel(6,"energyDeposits",2);
summaryReportMap->setBinLabel(7,"multiplicity",2);
dbe->setCurrentFolder("Muons/EventInfo/reportSummaryContents");
theSummaryContents.push_back(dbe->bookFloat("kinematics_GLB"));
theSummaryContents.push_back(dbe->bookFloat("muonId_GLB"));
theSummaryContents.push_back(dbe->bookFloat("residuals_GLB"));
theSummaryContents.push_back(dbe->bookFloat("GLB"));
theSummaryContents.push_back(dbe->bookFloat("kinematics_TK"));
theSummaryContents.push_back(dbe->bookFloat("muonId_TK"));
theSummaryContents.push_back(dbe->bookFloat("residuals_TK"));
theSummaryContents.push_back(dbe->bookFloat("TK"));
theSummaryContents.push_back(dbe->bookFloat("kinematics_STA"));
theSummaryContents.push_back(dbe->bookFloat("residuals_STA"));
theSummaryContents.push_back(dbe->bookFloat("STA"));
theSummaryContents.push_back(dbe->bookFloat("energyDeposits"));
theSummaryContents.push_back(dbe->bookFloat("multiplicity"));
// certification report
dbe->setCurrentFolder("Muons/EventInfo");
summaryCertification = dbe->bookFloat("CertificationSummary");
summaryCertification->Fill(-1);
summaryCertificationMap = dbe->book2D("CertificationSummaryMap","Muon Certification Summary Map",9,1,10,7,1,8);
summaryCertificationMap->setAxisTitle("muons",1);
summaryCertificationMap->setBinLabel(1,"GLB_Tot",1);
summaryCertificationMap->setBinLabel(2,"TK_Tot",1);
summaryCertificationMap->setBinLabel(3,"STA_tot",1);
summaryCertificationMap->setBinLabel(4,"GLB_B",1);
summaryCertificationMap->setBinLabel(5,"TK_B",1);
summaryCertificationMap->setBinLabel(6,"STA_B",1);
summaryCertificationMap->setBinLabel(7,"GLB_EC",1);
summaryCertificationMap->setBinLabel(8,"TK_EC",1);
summaryCertificationMap->setBinLabel(9,"STA_EC",1);
summaryCertificationMap->setAxisTitle("test",2);
summaryCertificationMap->setBinLabel(1,"#chi^{2}/Df",2);
summaryCertificationMap->setBinLabel(2,"#eta",2);
summaryCertificationMap->setBinLabel(3,"#phi",2);
summaryCertificationMap->setBinLabel(4,"residuals",2);
summaryCertificationMap->setBinLabel(5,"muonId",2);
summaryCertificationMap->setBinLabel(6,"energyDeposits",2);
summaryCertificationMap->setBinLabel(7,"multiplicity",2);
dbe->setCurrentFolder("Muons/EventInfo/CertificationContents");
theCertificationContents.push_back(dbe->bookFloat("GLB_Tot"));
theCertificationContents.push_back(dbe->bookFloat("STA_Tot"));
theCertificationContents.push_back(dbe->bookFloat("TK_Tot"));
theCertificationContents.push_back(dbe->bookFloat("GLB_B"));
theCertificationContents.push_back(dbe->bookFloat("STA_B"));
theCertificationContents.push_back(dbe->bookFloat("TK_B"));
theCertificationContents.push_back(dbe->bookFloat("GLB_EC"));
theCertificationContents.push_back(dbe->bookFloat("STA_EC"));
theCertificationContents.push_back(dbe->bookFloat("TK_EC"));
for (unsigned int icert=0;icert <theCertificationContents.size();icert++){
theCertificationContents[icert]->Fill(-1);
}
}
| void MuonTestSummary::beginLuminosityBlock | ( | edm::LuminosityBlock const & | lumiSeg, |
| edm::EventSetup const & | c | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 296 of file MuonTestSummary.cc.
{
// LogTrace(metname)<<"[MuonTestSummary]: beginLuminosityBlock";
}
| void MuonTestSummary::beginRun | ( | edm::Run const & | run, |
| edm::EventSetup const & | eSetup | ||
| ) | [protected, virtual] |
Histograms initialisation.
Reimplemented from edm::EDAnalyzer.
Definition at line 250 of file MuonTestSummary.cc.
References LogTrace, and metname.
{
LogTrace(metname)<<"[MuonTestSummary]: beginRun";
// initialisation of histo bins
for(int xBin=1; xBin<=5; xBin++){
for(int yBin=1; yBin<=3; yBin++){
kinematicsSummaryMap->Fill(xBin,yBin,1);
}
for(int yBin=1; yBin<=5; yBin++){
chi2TestSummaryMap->Fill(xBin,yBin,1);
KolmogorovTestSummaryMap->Fill(xBin,yBin,1);
}
}
for(int xBin=1; xBin<=residualsSummaryMap->getNbinsX(); xBin++){
for(int yBin=1; yBin<=residualsSummaryMap->getNbinsY(); yBin++){
residualsSummaryMap->Fill(xBin,yBin,1);
}
}
residualsSummaryMap->setBinContent(4, 4, 1); //not used for now
for(int xBin=1; xBin<=muonIdSummaryMap->getNbinsX(); xBin++){
for(int yBin=1; yBin<=muonIdSummaryMap->getNbinsY(); yBin++){
muonIdSummaryMap->Fill(xBin,yBin,1);
}
}
for(int xBin=1; xBin<=3; xBin++){
for(int yBin=1; yBin<=3; yBin++){
energySummaryMap->Fill(xBin,yBin,1);
}
}
for(int xBin=1; xBin<=3; xBin++){
multiplicitySummaryMap->Fill(xBin,1);
}
for(int xBin=1; xBin<=3; xBin++){
for(int yBin=1; yBin<=7; yBin++){
summaryReportMap->Fill(xBin,yBin,1);
}
}
for(int xBin=1; xBin<=9; xBin++){
for(int yBin=1; yBin<=7; yBin++){
summaryCertificationMap->Fill(xBin,yBin,1);
}
}
}
| void MuonTestSummary::doEnergyTests | ( | std::string | nameHisto, |
| std::string | muonType, | ||
| int | bin | ||
| ) | [protected] |
Definition at line 1191 of file MuonTestSummary.cc.
References ExpressReco_HICollisions_FallBack::expPeakEcalS9_max, ExpressReco_HICollisions_FallBack::expPeakEcalS9_min, ExpressReco_HICollisions_FallBack::expPeakHadS9_max, ExpressReco_HICollisions_FallBack::expPeakHadS9_min, MonitorElement::getEntries(), MonitorElement::getTH1F(), langaupro(), LogTrace, metname, and path().
{
// num matches test
string path = "Muons/MuonEnergyDepositAnalyzer/"+histname+muonType;
MonitorElement * energyHisto = dbe->get(path);
Double_t hPeak=-1, hFWHM=-1;
if(energyHisto){
TH1F * energyHisto_root = energyHisto->getTH1F();
// Setting fit range and start values
Double_t fitRange[2];
Double_t startValues[4], parlimitslo[4], parlimitshi[4], fitPar[4], fitParErr[4];
if (energyHisto->getEntries()>20){
if(histname=="ecalS9PointingMuDepositedEnergy_"){
fitRange[0]=0.04;
fitRange[1]=3.0;
startValues[0]=0.036; startValues[1]=0.193; startValues[2]=110.0; startValues[3]=0.06;
parlimitslo[0]=0.0; parlimitslo[1]=0.; parlimitslo[2]=1.0; parlimitslo[3]=0.;
parlimitshi[0]=0.05; parlimitshi[1]=0.5; parlimitshi[2]=80000.0; parlimitshi[3]=0.1;
Double_t chisqr;
Int_t ndf;
TF1 *fit = langaufit(energyHisto_root,fitRange,startValues,parlimitslo,parlimitshi,fitPar,fitParErr,&chisqr,&ndf);
if(fit){
langaupro(fitPar,hPeak,hFWHM);
LogTrace(metname)<<"hPeak from langau fit: "<<hPeak<<" for: "<<histname+muonType<<endl;
LogTrace(metname)<<"hFWHM from langau fit: "<<hFWHM<<" for: "<<histname+muonType<<endl;
}
}
if(histname=="hadS9PointingMuDepositedEnergy_"){
fitRange[0]=0.0;
fitRange[1]=7.0;
startValues[0]=2.0; startValues[1]=2.4; startValues[2]=110.0; startValues[3]=4.0;
parlimitslo[0]=0.0; parlimitslo[1]=0.; parlimitslo[2]=1.0; parlimitslo[3]=0.;
parlimitshi[0]=4.0; parlimitshi[1]=4.0; parlimitshi[2]=80000.0; parlimitshi[3]=8.0;
Double_t chisqr;
Int_t ndf;
TF1 *fit = langaufit(energyHisto_root,fitRange,startValues,parlimitslo,parlimitshi,fitPar,fitParErr,&chisqr,&ndf);
if(fit){
langaupro(fitPar,hPeak,hFWHM);
LogTrace(metname)<<"hPeak from langau fit: "<<hPeak<<" for: "<<histname+muonType<<endl;
LogTrace(metname)<<"hFWHM from langau fit: "<<hFWHM<<" for: "<<histname+muonType<<endl;
}
}
}
else{
LogTrace(metname) << "[MuonTestSummary]: Test of Energy for "<< histname+muonType << " not performed because # entries < 20 ";
}
}
if(histname=="ecalS9PointingMuDepositedEnergy_" && hPeak>expPeakEcalS9_min && hPeak<expPeakEcalS9_max)
energySummaryMap->setBinContent(binNumber,1, 1);
if(histname=="ecalS9PointingMuDepositedEnergy_" && (hPeak!=-1) &&!(hPeak>expPeakEcalS9_min && hPeak<expPeakEcalS9_max))
energySummaryMap->setBinContent(binNumber,1, 0);
if(histname=="hadS9PointingMuDepositedEnergy_" && hPeak>expPeakHadS9_min && hPeak<expPeakHadS9_max)
energySummaryMap->setBinContent(binNumber,2, 1);
if(histname=="hadS9PointingMuDepositedEnergy_" && (hPeak!=-1) && !(hPeak>expPeakHadS9_min && hPeak<expPeakHadS9_max))
energySummaryMap->setBinContent(binNumber,2, 0);
//missing test on ho distributions
}
| void MuonTestSummary::doKinematicsTests | ( | std::string | muonType, |
| int | bin | ||
| ) | [protected] |
test operations
Definition at line 523 of file MuonTestSummary.cc.
References ExpressReco_HICollisions_FallBack::chi2Fraction, ExpressReco_HICollisions_FallBack::chi2Spread, ExpressReco_HICollisions_FallBack::etaExpected, QReport::getMessage(), MonitorElement::getQReport(), QReport::getQTresult(), QReport::getStatus(), MonitorElement::getTH1F(), LogTrace, metname, path(), ExpressReco_HICollisions_FallBack::phiExpected, and mathSSE::sqrt().
{
// chi2 test
string path = "Muons/MuonRecoAnalyzer/" + muonType + "chi2OverDf";
MonitorElement * chi2Histo = dbe->get(path);
if(chi2Histo){
TH1F * chi2Histo_root = chi2Histo->getTH1F();
if(chi2Histo_root->GetEntries()>20){
//Standard QT based on fraction of events below and above a cut
LogTrace(metname)<<"chi2 kin test based on fraction for "<<muonType<<endl;
int maxBin = chi2Histo_root->GetMaximumBin();
if(chi2Histo_root->Integral(maxBin+1,chi2Histo_root->GetNbinsX())!=0){
double fraction = double(chi2Histo_root->Integral(1,maxBin))/double(chi2Histo_root->Integral(maxBin+1,chi2Histo_root->GetNbinsX()));
LogTrace(metname)<<"chi2 fraction for "<<muonType<<" : "<<fraction<< " must be within "<< chi2Fraction-chi2Spread<<","<<chi2Fraction+chi2Spread<<endl;
if(fraction>(chi2Fraction-chi2Spread) && fraction<(chi2Fraction+chi2Spread))
kinematicsSummaryMap->setBinContent(bin,1,1);
else
kinematicsSummaryMap->setBinContent(bin,1,0);
}
//QT based on comp wrt reference based on chi2
LogTrace(metname)<<"chi2 kin test based on comp wrt reference on terms of chi2 for "<<muonType<<endl;
const QReport * myQReport = chi2Histo->getQReport("Mu_Comp2RefChi2"); //get QReport associated to your ME
if(myQReport) {
LogTrace(metname) << "Chi2Report exists!!";
float qtresult = myQReport->getQTresult(); // get QT result value
int qtstatus = myQReport->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReport->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult << " qtstatus "<<qtstatus<<endl;
chi2TestSummaryMap->setBinContent(bin,1,qtresult);
}
//QT based on comp wrt reference based on kolmogorov test
LogTrace(metname)<<"chi2 kin test based on comp wrt reference using kolmogorov for "<<muonType<<endl;
const QReport * myQReportKolmo = chi2Histo->getQReport("Mu_Comp2RefKolmogorov");
if(myQReportKolmo) {
LogTrace(metname) << "Chi2Report Kolmogorov exists!!";
float qtresult = myQReportKolmo->getQTresult(); // get QT result value
int qtstatus = myQReportKolmo->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReportKolmo->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult<< " qtstatus "<<qtstatus<<endl;
KolmogorovTestSummaryMap->setBinContent(bin,1,qtresult);
}
}
else{
LogTrace(metname) << "[MuonTestSummary]: Test of Chi2 Kin not performed for "<< muonType << " because # entries < 20 ";
}
}
// pseudorapidity test
path = "Muons/MuonRecoAnalyzer/" + muonType + "eta";
MonitorElement * etaHisto = dbe->get(path);
if(etaHisto){
TH1F * etaHisto_root = etaHisto->getTH1F();
if (etaHisto_root ->GetEntries()>20){
//Standard QT based on fraction of events below and above a cut
LogTrace(metname)<<"eta kin test based on fraction for "<<muonType<<endl;
double binSize = (etaHisto_root->GetXaxis()->GetXmax()-etaHisto_root->GetXaxis()->GetXmin())/etaHisto_root->GetNbinsX();
int binZero = int((0-etaHisto_root->GetXaxis()->GetXmin())/binSize);
if(etaHisto_root->Integral(1,binZero-1)!=0 && etaHisto_root->Integral(binZero,etaHisto_root->GetNbinsX())!=0){
double symmetryFactor =
double(etaHisto_root->Integral(1,binZero-1)) / double(etaHisto_root->Integral(binZero,etaHisto_root->GetNbinsX()));
double errSymmetryFactor =
symmetryFactor*sqrt(1.0/double(etaHisto_root->Integral(1,binZero-1)) + 1.0/double(etaHisto_root->Integral(binZero,etaHisto_root->GetNbinsX())));
LogTrace(metname)<<"eta symmetryFactor for "<<muonType<<" : "<<symmetryFactor<< " (expected :" <<etaExpected <<")"<<endl;
LogTrace(metname)<<"eta errSymmetryFactor for "<<muonType<<" : "<<errSymmetryFactor<<endl;
double tParameter;
if((symmetryFactor-etaExpected)>0) tParameter=double(symmetryFactor-etaExpected)/errSymmetryFactor;
else tParameter=double(-symmetryFactor+etaExpected)/errSymmetryFactor;
LogTrace(metname)<<"eta tParameter for "<<muonType<<" : "<<tParameter<<" (expected < 1.95)" <<endl;
if (tParameter<1.95) //2sigma rejection
kinematicsSummaryMap->setBinContent(bin,2,1);
else
kinematicsSummaryMap->setBinContent(bin,2,0);
}
//QT based on comp wrt reference based on chi2
LogTrace(metname)<<"eta kin test based on comp wrt reference on terms of chi2 for "<<muonType<<endl;
const QReport * myQReport = etaHisto->getQReport("Mu_Comp2RefChi2"); //get QReport associated to your ME
if(myQReport) {
LogTrace(metname) << "EtaReport exists!!";
float qtresult = myQReport->getQTresult(); // get QT result value
int qtstatus = myQReport->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReport->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult << " qtstatus "<<qtstatus<<endl;
chi2TestSummaryMap->setBinContent(bin,2,qtresult);
}
//QT based on comp wrt reference based on kolmogorov test
LogTrace(metname)<<"eta kin test based on comp wrt reference using kolmogorov for "<<muonType<<endl;
const QReport * myQReportKolmo = etaHisto->getQReport("Mu_Comp2RefKolmogorov");
if(myQReportKolmo) {
LogTrace(metname) << "EtaReport Kolmogorov exists!!";
float qtresult = myQReportKolmo->getQTresult(); // get QT result value
int qtstatus = myQReportKolmo->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReportKolmo->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult<< " qtstatus "<<qtstatus<<endl;
KolmogorovTestSummaryMap->setBinContent(bin,2,qtresult);
}
}
else{
LogTrace(metname) << "[MuonTestSummary]: Test of Eta Kin not performed for "<< muonType << " because # entries < 20 ";
}
}
// phi test
path = "Muons/MuonRecoAnalyzer/" + muonType + "phi";
MonitorElement * phiHisto = dbe->get(path);
if(phiHisto ){
TH1F * phiHisto_root = phiHisto->getTH1F();
if(phiHisto_root ->GetEntries()>20){
//Standard QT based on fraction of events below and above a cut
LogTrace(metname)<<"phi kin test based on fraction for "<<muonType<<endl;
double binSize = (phiHisto_root->GetXaxis()->GetXmax()-phiHisto_root->GetXaxis()->GetXmin())/phiHisto_root->GetNbinsX();
int binZero = int((0-phiHisto_root->GetXaxis()->GetXmin())/binSize);
if(phiHisto_root->Integral(binZero+1,phiHisto_root->GetNbinsX())!=0 && phiHisto_root->Integral(1,binZero)!=0){
double symmetryFactor =
double(phiHisto_root->Integral(binZero+1,phiHisto_root->GetNbinsX())) / double(phiHisto_root->Integral(1,binZero));
double errSymmetryFactor =
symmetryFactor*sqrt(1.0/double(phiHisto_root->Integral(binZero+1,phiHisto_root->GetNbinsX())) + 1.0/double(phiHisto_root->Integral(1,binZero)));
LogTrace(metname)<<"phi symmetryFactor for "<<muonType<<" : "<<symmetryFactor<< "(phi expected :" <<phiExpected <<")"<<endl;
LogTrace(metname)<<"phi errSymmetryFactor for "<<muonType<<" : "<<errSymmetryFactor<<endl;
double tParameter;
if((symmetryFactor-phiExpected)>0) tParameter=double(symmetryFactor-phiExpected)/errSymmetryFactor;
else tParameter=double(-symmetryFactor+phiExpected)/errSymmetryFactor;
LogTrace(metname)<<"phi tParameter for "<<muonType<<" : "<<tParameter<<" (expected < 1.95)" <<endl;
if (tParameter<1.95) //2sigma rejection
kinematicsSummaryMap->setBinContent(bin,3,1);
else
kinematicsSummaryMap->setBinContent(bin,3,0);
}
//QT based on comp wrt reference based on chi2
LogTrace(metname)<<"phi kin test based on comp wrt reference on terms of chi2 for "<<muonType<<endl;
const QReport * myQReport = phiHisto->getQReport("Mu_Comp2RefChi2"); //get QReport associated to your ME
if(myQReport) {
LogTrace(metname) << "PhiReport exists!!";
float qtresult = myQReport->getQTresult(); // get QT result value
int qtstatus = myQReport->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReport->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult<< " qtstatus "<<qtstatus<<endl;
chi2TestSummaryMap->setBinContent(bin,3,qtresult);
}
//QT based on comp wrt reference based on kolmogorov test
LogTrace(metname)<<"phi kin test based on comp wrt reference using kolmogorov for "<<muonType<<endl;
const QReport * myQReportKolmo = phiHisto->getQReport("Mu_Comp2RefKolmogorov");
if(myQReportKolmo) {
LogTrace(metname) << "PhiReport Kolmogorov exists!!";
float qtresult = myQReportKolmo->getQTresult(); // get QT result value
int qtstatus = myQReportKolmo->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReportKolmo->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult<< " qtstatus "<<qtstatus<<endl;
KolmogorovTestSummaryMap->setBinContent(bin,3,qtresult);
}
}
else{
LogTrace(metname) << "[MuonTestSummary]: Test of Phi Kin not performed for "<< muonType << " because # entries < 20 ";
}
}
// pt test
path = "Muons/MuonRecoAnalyzer/" + muonType + "pt";
MonitorElement * ptHisto = dbe->get(path);
if(ptHisto ){
TH1F * ptHisto_root = ptHisto->getTH1F();
if(ptHisto_root ->GetEntries()>20){
//QT based on comp wrt reference based on chi2
LogTrace(metname)<<"pt kin test based on comp wrt reference on terms of chi2 for "<<muonType<<endl;
const QReport * myQReport = ptHisto->getQReport("Mu_Comp2RefChi2"); //get QReport associated to your ME
if(myQReport) {
LogTrace(metname) << "PtReport exists!!";
float qtresult = myQReport->getQTresult(); // get QT result value
int qtstatus = myQReport->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReport->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult<< " qtstatus "<<qtstatus<<endl;
chi2TestSummaryMap->setBinContent(bin,4,qtresult);
}
//QT based on comp wrt reference based on kolmogorov test
LogTrace(metname)<<"pt kin test based on comp wrt reference using kolmogorov for "<<muonType<<endl;
const QReport * myQReportKolmo = ptHisto->getQReport("Mu_Comp2RefKolmogorov");
if(myQReportKolmo) {
LogTrace(metname) << "PtReport Kolmogorov exists!!";
float qtresult = myQReportKolmo->getQTresult(); // get QT result value
int qtstatus = myQReportKolmo->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReportKolmo->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult<< " qtstatus "<<qtstatus<<endl;
KolmogorovTestSummaryMap->setBinContent(bin,4,qtresult);
}
}
else{
LogTrace(metname) << "[MuonTestSummary]: Test of Pt Kin not performed for "<< muonType << " because # entries < 20 ";
}
}
// q test
path = "Muons/MuonRecoAnalyzer/" + muonType + "q";
MonitorElement * qHisto = dbe->get(path);
if(ptHisto ){
TH1F * qHisto_root = qHisto->getTH1F();
if(qHisto_root ->GetEntries()>20){
//QT based on comp wrt reference based on chi2
LogTrace(metname)<<"q kin test based on comp wrt reference on terms of chi2 for "<<muonType<<endl;
const QReport * myQReport = qHisto->getQReport("Mu_Comp2RefChi2"); //get QReport associated to your ME
if(myQReport) {
LogTrace(metname) << "qReport exists!!";
float qtresult = myQReport->getQTresult(); // get QT result value
int qtstatus = myQReport->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReport->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult<< " qtstatus "<<qtstatus<<endl;
chi2TestSummaryMap->setBinContent(bin,5,qtresult);
}
//QT based on comp wrt reference based on kolmogorov test
LogTrace(metname)<<"q kin test based on comp wrt reference using kolmogorov for "<<muonType<<endl;
const QReport * myQReportKolmo = qHisto->getQReport("Mu_Comp2RefKolmogorov");
if(myQReportKolmo) {
LogTrace(metname) << "qReport Kolmogorov exists!!";
float qtresult = myQReportKolmo->getQTresult(); // get QT result value
int qtstatus = myQReportKolmo->getStatus() ; // get QT status value (see table below)
std::string qtmessage = myQReportKolmo->getMessage() ; // get the whole QT result message
LogTrace(metname) << "qtresult " <<qtresult<< " qtstatus "<<qtstatus<<endl;
KolmogorovTestSummaryMap->setBinContent(bin,5,qtresult);
}
}
else{
LogTrace(metname) << "[MuonTestSummary]: Test of q Kin not performed for "<< muonType << " because # entries < 20 ";
}
}
}
| void MuonTestSummary::doMultiplicityTests | ( | ) | [protected] |
Definition at line 1260 of file MuonTestSummary.cc.
References ExpressReco_HICollisions_FallBack::expMultiplicityGlb_max, ExpressReco_HICollisions_FallBack::expMultiplicityGlb_min, ExpressReco_HICollisions_FallBack::expMultiplicitySta_max, ExpressReco_HICollisions_FallBack::expMultiplicitySta_min, ExpressReco_HICollisions_FallBack::expMultiplicityTk_max, ExpressReco_HICollisions_FallBack::expMultiplicityTk_min, MonitorElement::getBinContent(), MonitorElement::getEntries(), LogTrace, and metname.
{
MonitorElement* multiplicityHisto = dbe->get("Muons/MuonRecoAnalyzer/muReco");
if(multiplicityHisto){
if(multiplicityHisto->getEntries()>20){
double multiplicity_GLB = double(multiplicityHisto->getBinContent(1)+multiplicityHisto->getBinContent(2))/double(multiplicityHisto->getEntries());
LogTrace(metname)<<"multiplicity_GLB: "<<multiplicity_GLB<< " ExpMultiplicityGlb_min " <<expMultiplicityGlb_min <<" ExpMultiplicityGlb_max " <<expMultiplicityGlb_max<<endl;
double multiplicity_TK = double(multiplicityHisto->getBinContent(3)+multiplicityHisto->getBinContent(4))/double(multiplicityHisto->getEntries());
LogTrace(metname)<<"multiplicity_TK: "<<multiplicity_TK<< " ExpMultiplicityTk_min " <<expMultiplicityTk_min <<" ExpMultiplicityTk_max " <<expMultiplicityTk_max<<endl;
double multiplicity_STA = double(multiplicityHisto->getBinContent(3)+multiplicityHisto->getBinContent(5))/double(multiplicityHisto->getEntries());
LogTrace(metname)<<"multiplicity_STA: "<<multiplicity_STA<< " ExpMultiplicitySta_min " <<expMultiplicitySta_min <<" ExpMultiplicitySta_max " <<expMultiplicitySta_max<<endl;
if(multiplicity_GLB>expMultiplicityGlb_min && multiplicity_GLB<expMultiplicityGlb_max)
multiplicitySummaryMap->setBinContent(1,1);
else
multiplicitySummaryMap->setBinContent(1,0);
if(multiplicity_TK>expMultiplicityTk_min && multiplicity_TK<expMultiplicityTk_max)
multiplicitySummaryMap->setBinContent(2,1);
else
multiplicitySummaryMap->setBinContent(2,0);
if(multiplicity_STA>expMultiplicitySta_min && multiplicity_STA<expMultiplicitySta_max)
multiplicitySummaryMap->setBinContent(3,1);
else
multiplicitySummaryMap->setBinContent(3,0);
}
else{
LogTrace(metname) << "[MuonTestSummary]: Test of Multiplicity not performed because # entries < 20 ";
}
}
}
| void MuonTestSummary::doMuonIDTests | ( | ) | [protected] |
Definition at line 923 of file MuonTestSummary.cc.
References MonitorElement::getBinContent(), MonitorElement::getEntries(), MonitorElement::getTH1F(), i, LogTrace, ExpressReco_HICollisions_FallBack::matchesFractionCsc_max, ExpressReco_HICollisions_FallBack::matchesFractionCsc_min, ExpressReco_HICollisions_FallBack::matchesFractionDt_max, ExpressReco_HICollisions_FallBack::matchesFractionDt_min, metname, ExpressReco_HICollisions_FallBack::numMatchedExpected_max, ExpressReco_HICollisions_FallBack::numMatchedExpected_min, path(), ResidualCheck(), ExpressReco_HICollisions_FallBack::resSegmTrack_mean_max, ExpressReco_HICollisions_FallBack::resSegmTrack_mean_min, ExpressReco_HICollisions_FallBack::resSegmTrack_rms_max, and ExpressReco_HICollisions_FallBack::resSegmTrack_rms_min.
{
vector<string> muType;
muType.push_back("GlobalMuons");
muType.push_back("TrackerMuons");
for(int i=0; i<=1; i++){
// num matches test
string path = "Muons/MuonIdDQM/" + muType[i] + "/hNumMatches";
MonitorElement * matchesHisto = dbe->get(path);
if(matchesHisto){
TH1F * matchesHisto_root = matchesHisto->getTH1F();
if(matchesHisto_root->GetMaximumBin() >= numMatchedExpected_min && matchesHisto_root->GetMaximumBin() <= numMatchedExpected_max)
muonIdSummaryMap->setBinContent(i+1,1,1);
else
muonIdSummaryMap->setBinContent(i+1,1,0);
}
// num of 0 associated segments
double numOneSegm_dt = 0;
int numHistos_dt=0;
int numHistos_csc=0;
MonitorElement * DT1Histo = dbe->get("Muons/MuonIdDQM/" + muType[i] + "/hDT1NumSegments");
if(DT1Histo)
{numHistos_dt++;
if(DT1Histo->getEntries()!=0) numOneSegm_dt+=double(DT1Histo->getBinContent(2))/double(DT1Histo->getEntries());}
MonitorElement * DT2Histo = dbe->get("Muons/MuonIdDQM/" + muType[i] + "/hDT2NumSegments");
if(DT2Histo)
{numHistos_dt++;
if(DT2Histo->getEntries()!=0) numOneSegm_dt+=double(DT2Histo->getBinContent(2))/double(DT2Histo->getEntries());}
MonitorElement * DT3Histo = dbe->get("Muons/MuonIdDQM/" + muType[i] + "/hDT3NumSegments");
if(DT3Histo)
{numHistos_dt++;
if(DT3Histo->getEntries()!=0) numOneSegm_dt+=double(DT3Histo->getBinContent(2))/double(DT3Histo->getEntries());}
MonitorElement * DT4Histo = dbe->get("Muons/MuonIdDQM/" + muType[i] + "/hDT4NumSegments");
if(DT4Histo)
{numHistos_dt++;
if(DT4Histo->getEntries()!=0) numOneSegm_dt+=double(DT4Histo->getBinContent(2))/double(DT4Histo->getEntries());}
double fraction_dt=0;
if(numOneSegm_dt!=0){
fraction_dt = numOneSegm_dt/double(numHistos_dt);
LogTrace(metname)<<"fraction_dt: "<<fraction_dt<<" for "<<muType[i]<<endl;
}
double numOneSegm_csc = 0;
MonitorElement * CSC1Histo = dbe->get("Muons/MuonIdDQM/" + muType[i] + "/hCSC1NumSegments");
if(CSC1Histo)
{numHistos_csc++;
if(CSC1Histo->getEntries()!=0) numOneSegm_csc+=double(CSC1Histo->getBinContent(2))/double(CSC1Histo->getEntries());}
MonitorElement * CSC2Histo = dbe->get("Muons/MuonIdDQM/" + muType[i] + "/hCSC2NumSegments");
if(CSC2Histo)
{numHistos_csc++;
if(CSC2Histo->getEntries()!=0) numOneSegm_csc+=double(CSC2Histo->getBinContent(2))/double(CSC2Histo->getEntries());}
MonitorElement * CSC3Histo = dbe->get("Muons/MuonIdDQM/" + muType[i] + "/hCSC3NumSegments");
if(CSC3Histo)
{numHistos_csc++;
if(CSC3Histo->getEntries()!=0) numOneSegm_csc+=double(CSC3Histo->getBinContent(2))/double(CSC3Histo->getEntries());}
MonitorElement * CSC4Histo = dbe->get("Muons/MuonIdDQM/" + muType[i] + "/hCSC4NumSegments");
if(CSC4Histo)
{numHistos_csc++;
if(CSC4Histo->getEntries()!=0) numOneSegm_csc+=double(CSC4Histo->getBinContent(2))/double(CSC4Histo->getEntries());}
double fraction_csc=0;
if(numOneSegm_csc!=0){
fraction_csc = numOneSegm_csc/double(numHistos_csc);
LogTrace(metname)<<"fraction_csc: "<<fraction_csc<<" for "<<muType[i]<<endl;
}
//--GH modified
if(fraction_dt>matchesFractionDt_min && fraction_dt<matchesFractionDt_max)
muonIdSummaryMap->setBinContent(2*i+1,1,1);
else
muonIdSummaryMap->setBinContent(2*i+1,1,0);
if(fraction_csc>matchesFractionCsc_min && fraction_csc<matchesFractionCsc_max)
muonIdSummaryMap->setBinContent(2*i+2,1,1);
else
muonIdSummaryMap->setBinContent(2*i+2,1,0);
//--GH modified
// residuals test
vector<string> DTXresHistos, DTYresHistos, CSCXresHistos, CSCYresHistos;
DTXresHistos.push_back("hDT1Pullx");
DTXresHistos.push_back("hDT2Pullx");
DTXresHistos.push_back("hDT3Pullx");
DTXresHistos.push_back("hDT4Pullx");
DTYresHistos.push_back("hDT1Pully");
DTYresHistos.push_back("hDT2Pully");
DTYresHistos.push_back("hDT3Pully");
CSCXresHistos.push_back("hCSC1Pullx");
CSCXresHistos.push_back("hCSC2Pullx");
CSCXresHistos.push_back("hCSC3Pullx");
CSCXresHistos.push_back("hCSC4Pullx");
CSCYresHistos.push_back("hCSC1Pully");
CSCYresHistos.push_back("hCSC2Pully");
CSCYresHistos.push_back("hCSC3Pully");
CSCYresHistos.push_back("hCSC4Pully");
int numPlot_dtX, numPlot_dtY, numPlot_cscX, numPlot_cscY;
double dtSigmaX, dtSigmaY, cscSigmaX, cscSigmaY;
double dtSigmaX_err, dtSigmaY_err, cscSigmaX_err, cscSigmaY_err;
double dtMeanX, dtMeanY, cscMeanX, cscMeanY;
double dtMeanX_err, dtMeanY_err, cscMeanX_err, cscMeanY_err;
MuonTestSummary::ResidualCheck(muType[i], DTXresHistos, numPlot_dtX, dtMeanX, dtMeanX_err, dtSigmaX, dtSigmaX_err);
MuonTestSummary::ResidualCheck(muType[i], DTYresHistos, numPlot_dtY, dtMeanY, dtMeanY_err, dtSigmaY, dtSigmaY_err);
MuonTestSummary::ResidualCheck(muType[i], CSCXresHistos, numPlot_cscX, cscMeanX, cscMeanX_err, cscSigmaX, cscSigmaX_err);
MuonTestSummary::ResidualCheck(muType[i], CSCYresHistos, numPlot_cscY, cscMeanY, cscMeanY_err, cscSigmaY, cscSigmaY_err);
LogTrace(metname)<<"DT mean must be between: "<<resSegmTrack_mean_min <<" and "<<resSegmTrack_mean_max<<endl;
LogTrace(metname)<<"DT rms must be between: "<<resSegmTrack_rms_min <<" and "<<resSegmTrack_rms_max<<endl;
LogTrace(metname)<<"DT X residual "<< muType[i]<<" mean: " << dtMeanX<<" +- "<< dtMeanX_err
<<", sigma: "<< dtSigmaX <<" +- "<<dtSigmaX_err<< endl;
LogTrace(metname)<<"DT Y residual "<< muType[i]<<" mean: " << dtMeanY<<" +- "<< dtMeanY_err
<<", sigma: "<< dtSigmaY <<" +- "<<dtSigmaY_err<< endl;
LogTrace(metname)<<"CSC X residual "<< muType[i]<<" mean: " << cscMeanX<<" +- "<< cscMeanX_err
<<", sigma: "<< cscSigmaX <<" +- "<<cscSigmaX_err<< endl;
LogTrace(metname)<<"CSC Y residual "<< muType[i]<<" mean: " << cscMeanY<<" +- "<< cscMeanY_err
<<", sigma: "<< cscSigmaY <<" +- "<<cscSigmaY_err<< endl;
//require the mean and rms to be within nsig sigma of preferred range;
const int nsig=2;
if(numPlot_dtX > 0 ) {
if( dtMeanX + nsig*dtMeanX_err>resSegmTrack_mean_min &&
dtMeanX - nsig*dtMeanX_err<resSegmTrack_mean_max)
muonIdSummaryMap->setBinContent(2*i+1,2,1);
else
muonIdSummaryMap->setBinContent(2*i+1,2,0);
if( dtSigmaX + nsig*dtSigmaX_err>resSegmTrack_rms_min &&
dtSigmaX - nsig*dtSigmaX_err<resSegmTrack_rms_max)
muonIdSummaryMap->setBinContent(2*i+1,3,1);
else
muonIdSummaryMap->setBinContent(2*i+1,3,0);
}
if(numPlot_dtY > 0 ){
if( dtMeanY + nsig*dtMeanY_err>resSegmTrack_mean_min &&
dtMeanY - nsig*dtMeanY_err<resSegmTrack_mean_max)
muonIdSummaryMap->setBinContent(2*i+1,4,1);
else
muonIdSummaryMap->setBinContent(2*i+1,4,0);
if( dtSigmaY + nsig*dtSigmaY_err>resSegmTrack_rms_min &&
dtSigmaY - nsig*dtSigmaX_err<resSegmTrack_rms_max)
muonIdSummaryMap->setBinContent(2*i+1,5,1);
else
muonIdSummaryMap->setBinContent(2*i+1,5,0);
}
if(numPlot_cscX > 0 ) {
if( cscMeanX + nsig*cscMeanX_err>resSegmTrack_mean_min &&
cscMeanX - nsig*cscMeanX_err<resSegmTrack_mean_max)
muonIdSummaryMap->setBinContent(2*i+2,2,1);
else
muonIdSummaryMap->setBinContent(2*i+2,2,0);
if( cscSigmaX + nsig*cscSigmaX_err>resSegmTrack_rms_min &&
cscSigmaX - nsig*cscSigmaX_err<resSegmTrack_rms_max)
muonIdSummaryMap->setBinContent(2*i+2,3,1);
else
muonIdSummaryMap->setBinContent(2*i+2,3,0);
}
if(numPlot_cscY > 0 ){
if( cscMeanY + nsig*cscMeanY_err>resSegmTrack_mean_min &&
cscMeanY - nsig*cscMeanY_err<resSegmTrack_mean_max)
muonIdSummaryMap->setBinContent(2*i+2,4,1);
else
muonIdSummaryMap->setBinContent(2*i+2,4,0);
if( cscSigmaY + nsig*cscSigmaY_err>resSegmTrack_rms_min &&
cscSigmaY - nsig*cscSigmaY_err<resSegmTrack_rms_max)
muonIdSummaryMap->setBinContent(2*i+2,5,1);
else
muonIdSummaryMap->setBinContent(2*i+2,5,0);
}
//---- end of modification
}
}
| void MuonTestSummary::doResidualsTests | ( | std::string | type, |
| std::string | parameter, | ||
| int | bin | ||
| ) | [protected] |
Definition at line 817 of file MuonTestSummary.cc.
References MonitorElement::getBinContent(), LogTrace, plotscripts::mean(), metname, path(), ExpressReco_HICollisions_FallBack::pullEtaSpread, ExpressReco_HICollisions_FallBack::pullOneOvPSpread, ExpressReco_HICollisions_FallBack::pullPhiSpread, ExpressReco_HICollisions_FallBack::resChargeLimit_glbSta, ExpressReco_HICollisions_FallBack::resChargeLimit_tkGlb, ExpressReco_HICollisions_FallBack::resChargeLimit_tkSta, ExpressReco_HICollisions_FallBack::resEtaSpread_glbSta, ExpressReco_HICollisions_FallBack::resEtaSpread_tkGlb, ExpressReco_HICollisions_FallBack::resOneOvPSpread_glbSta, ExpressReco_HICollisions_FallBack::resOneOvPSpread_tkGlb, ExpressReco_HICollisions_FallBack::resPhiSpread_glbSta, ExpressReco_HICollisions_FallBack::resPhiSpread_tkGlb, and ExpressReco_HICollisions_FallBack::sigma.
{
// residuals test
if(type!="GlbMuon"){
string path = "Muons/MuonRecoAnalyzer/Res_" + type + "_" + parameter;
MonitorElement * residualsHisto = dbe->get(path);
float mean = -1;
float mean_err = -1;
float sigma = -1;
float sigma_err = -1;
if(residualsHisto){
LogTrace(metname) << "[MuonTestSummary]: Starting Gaussian fit for Test of Res_"<<type<<"_"<<parameter<< endl;
GaussFit( type, parameter, residualsHisto, mean, mean_err, sigma, sigma_err);
if(sigma!=-1 && parameter=="eta" && type=="TkGlb"){
if(sigma-sigma_err<resEtaSpread_tkGlb) residualsSummaryMap->setBinContent(bin, 1, 1);
else residualsSummaryMap->setBinContent(bin, 1, 0);
}
if(sigma!=-1 && parameter=="eta" && (type=="GlbSta" || type=="TkSta")) {
if(sigma-sigma_err<resEtaSpread_glbSta) residualsSummaryMap->setBinContent(bin, 1, 1);
else residualsSummaryMap->setBinContent(bin, 1, 0);
}
if(sigma!=-1 && parameter=="phi" && type=="TkGlb"){
if(sigma-sigma_err<resPhiSpread_tkGlb) residualsSummaryMap->setBinContent(bin, 2, 1);
else residualsSummaryMap->setBinContent(bin, 2, 0);
}
if(sigma!=-1 && parameter=="phi" && (type=="GlbSta" || type=="TkSta")){
if(sigma-sigma_err<resPhiSpread_glbSta) residualsSummaryMap->setBinContent(bin, 2, 1);
else residualsSummaryMap->setBinContent(bin, 2, 0);
}
if(sigma!=-1 && parameter=="oneOverp" && type=="TkGlb"){
if(sigma-sigma_err<resOneOvPSpread_tkGlb) residualsSummaryMap->setBinContent(bin, 3, 1);
else residualsSummaryMap->setBinContent(bin, 3, 0);
}
if(sigma!=-1 && parameter=="oneOverp" && (type=="GlbSta" || type=="TkSta")) {
if(sigma-sigma_err<resOneOvPSpread_glbSta) residualsSummaryMap->setBinContent(bin, 3, 1);
else residualsSummaryMap->setBinContent(bin, 3, 0);
}
}
//--GH modified
if(type=="TkSta"){
//look at the pull:
string path = "Muons/MuonRecoAnalyzer/Pull_" + type + "_" + parameter;
MonitorElement * pullHisto = dbe->get(path);
if(pullHisto){
LogTrace(metname) << "[MuonTestSummary]: Starting Gaussian fit for Test of Pull_"<<type<<"_"<<parameter<< endl;
GaussFit( type, parameter, pullHisto, mean, mean_err, sigma, sigma_err);
if(sigma!=-1 && parameter=="eta" ){
if(sigma-sigma_err<pullEtaSpread) residualsSummaryMap->setBinContent(4, 1, 1);
else residualsSummaryMap->setBinContent(4, 1, 0);
}
if(sigma!=-1 && parameter=="phi"){
if(sigma-sigma_err<pullPhiSpread) residualsSummaryMap->setBinContent(4, 2, 1);
else residualsSummaryMap->setBinContent(4, 2, 0);
}
if(sigma!=-1 && parameter=="oneOverp"){
if(sigma-sigma_err<pullOneOvPSpread) residualsSummaryMap->setBinContent(4, 3, 1);
else residualsSummaryMap->setBinContent(4, 3, 0);
}
}//have pull histo
} //TkSta muons
}
//this part for Global Muons:
else{
string path = "Muons/MuonRecoAnalyzer/" + type + "_" + parameter;
MonitorElement * residualsHisto = dbe->get(path);
if(residualsHisto){
LogTrace(metname) << "[MuonTestSummary]: Test of Charge Comparison "<<type<<"_"<<parameter<< endl;
if((residualsHisto->getBinContent(3)+residualsHisto->getBinContent(4))!=0){
LogTrace(metname)<<"Charge comparison TkGlb: "<<residualsHisto->getBinContent(4)/double(residualsHisto->getBinContent(3)+residualsHisto->getBinContent(4))<<endl;
if(residualsHisto->getBinContent(4)/double(residualsHisto->getBinContent(3)+residualsHisto->getBinContent(4)) < resChargeLimit_tkGlb)
residualsSummaryMap->setBinContent(1, 4, 1);
else
residualsSummaryMap->setBinContent(1, 4, 0);
}
if((residualsHisto->getBinContent(1)+residualsHisto->getBinContent(2))!=0){
LogTrace(metname)<<"charge comparison GlbSta: "<<residualsHisto->getBinContent(2)/double(residualsHisto->getBinContent(1)+residualsHisto->getBinContent(2))<<endl;
if(residualsHisto->getBinContent(2)/double(residualsHisto->getBinContent(1)+residualsHisto->getBinContent(2)) < resChargeLimit_glbSta)
residualsSummaryMap->setBinContent(2, 4, 1);
else
residualsSummaryMap->setBinContent(2, 4, 0);
}
if(residualsHisto->getBinContent(5)+residualsHisto->getBinContent(6)!=0){
LogTrace(metname)<<"charge comparison TkSta: "<<residualsHisto->getBinContent(6)/double(residualsHisto->getBinContent(5)+residualsHisto->getBinContent(6))<<endl;
if(residualsHisto->getBinContent(6)/double(residualsHisto->getBinContent(5)+residualsHisto->getBinContent(6)) < resChargeLimit_tkSta)
residualsSummaryMap->setBinContent(3, 4, 1);
else
residualsSummaryMap->setBinContent(3, 4, 0);
}
}
}
}
| void MuonTestSummary::endLuminosityBlock | ( | edm::LuminosityBlock const & | lumiSeg, |
| edm::EventSetup const & | c | ||
| ) | [protected, virtual] |
Diagnostic.
Reimplemented from edm::EDAnalyzer.
Definition at line 300 of file MuonTestSummary.cc.
{
// LogTrace(metname)<<"[MuonTestSummary]: endLuminosityBlock, performing the DQM LS client operation";
}
| void MuonTestSummary::endRun | ( | edm::Run const & | run, |
| edm::EventSetup const & | eSetup | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 306 of file MuonTestSummary.cc.
References i, j, LogTrace, and metname.
{
LogTrace(metname)<<"[MuonTestSummary]: endRun, performing the DQM end of run client operation";
// fill the kinematics report summary
doKinematicsTests("GlbMuon_Glb_", 1);
doKinematicsTests("GlbMuon_Tk_", 2);
doKinematicsTests("GlbMuon_Sta_",3);
doKinematicsTests("TkMuon_", 4);
doKinematicsTests("StaMuon_", 5);
// fill the residuals report summary
doResidualsTests("TkGlb", "eta", 1);
doResidualsTests("GlbSta", "eta", 2);
doResidualsTests("TkSta", "eta", 3);
doResidualsTests("TkGlb", "phi", 1);
doResidualsTests("GlbSta", "phi", 2);
doResidualsTests("TkSta", "phi", 3);
doResidualsTests("TkGlb", "oneOverp", 1);
doResidualsTests("GlbSta", "oneOverp", 2);
doResidualsTests("TkSta", "oneOverp", 3);
doResidualsTests("GlbMuon", "qComparison", -1);
// fill the muonID report summary
doMuonIDTests();
// fill the energy report summary
doEnergyTests("ecalS9PointingMuDepositedEnergy_","Glb_muons", 1);
doEnergyTests("hadS9PointingMuDepositedEnergy_", "Glb_muons", 1);
doEnergyTests("hoS9PointingMuDepositedEnergy_", "Glb_muons", 1);
doEnergyTests("ecalS9PointingMuDepositedEnergy_", "Tk_muons", 2);
doEnergyTests("hadS9PointingMuDepositedEnergy_", "Tk_muons", 2);
doEnergyTests("hoS9PointingMuDepositedEnergy_", "Tk_muons", 2);
doEnergyTests("ecalS9PointingMuDepositedEnergy_", "Sta_muons", 3);
doEnergyTests("hadS9PointingMuDepositedEnergy_", "Sta_muons", 3);
doEnergyTests("hoS9PointingMuDepositedEnergy_", "Sta_muons", 3);
// fill the multiplicity test summary
doMultiplicityTests();
// fill the final report summary
/*
summaryReportMap->setBinContent(1,1,double(kinematicsSummaryMap->getBinContent(1,1)+kinematicsSummaryMap->getBinContent(2,1)+kinematicsSummaryMap->getBinContent(3,1))/3.0);
summaryReportMap->setBinContent(2,1,kinematicsSummaryMap->getBinContent(4,1));
summaryReportMap->setBinContent(3,1,kinematicsSummaryMap->getBinContent(5,1));
summaryReportMap->setBinContent(1,2,double(kinematicsSummaryMap->getBinContent(1,2)+kinematicsSummaryMap->getBinContent(2,2)+kinematicsSummaryMap->getBinContent(3,2))/3.0);
summaryReportMap->setBinContent(2,2,kinematicsSummaryMap->getBinContent(4,2));
summaryReportMap->setBinContent(3,2,kinematicsSummaryMap->getBinContent(5,2));
summaryReportMap->setBinContent(1,3,double(kinematicsSummaryMap->getBinContent(1,3)+kinematicsSummaryMap->getBinContent(2,3)+kinematicsSummaryMap->getBinContent(3,3))/3.0);
summaryReportMap->setBinContent(2,3,kinematicsSummaryMap->getBinContent(4,3));
summaryReportMap->setBinContent(3,3,kinematicsSummaryMap->getBinContent(5,3));
*/
//Changed to KolmogorovQuality test--------------------------
summaryReportMap->setBinContent(1,1,double(KolmogorovTestSummaryMap->getBinContent(1,1)+KolmogorovTestSummaryMap->getBinContent(2,1)+KolmogorovTestSummaryMap->getBinContent(3,1))/3.0);
summaryReportMap->setBinContent(2,1,KolmogorovTestSummaryMap->getBinContent(4,1));
summaryReportMap->setBinContent(3,1,KolmogorovTestSummaryMap->getBinContent(5,1));
summaryReportMap->setBinContent(1,2,double(KolmogorovTestSummaryMap->getBinContent(1,2)+KolmogorovTestSummaryMap->getBinContent(2,2)+KolmogorovTestSummaryMap->getBinContent(3,2))/3.0);
summaryReportMap->setBinContent(2,2,KolmogorovTestSummaryMap->getBinContent(4,2));
summaryReportMap->setBinContent(3,2,KolmogorovTestSummaryMap->getBinContent(5,2));
summaryReportMap->setBinContent(1,3,double(KolmogorovTestSummaryMap->getBinContent(1,3)+KolmogorovTestSummaryMap->getBinContent(2,3)+KolmogorovTestSummaryMap->getBinContent(3,3))/3.0);
summaryReportMap->setBinContent(2,3,KolmogorovTestSummaryMap->getBinContent(4,3));
summaryReportMap->setBinContent(3,3,KolmogorovTestSummaryMap->getBinContent(5,3));
//-- modified GH
double residualsSummary = 0;
//put the TRK-STA resid & pulls in the first bin ("GLB")
//then the GLB-TRK and GLB-STA residuals in the 2nd and 3rd
for(int i=3; i<=residualsSummaryMap->getNbinsX(); i++)
for (int j=1; j<=residualsSummaryMap->getNbinsY(); j++)
residualsSummary += residualsSummaryMap->getBinContent(i,j);
residualsSummary /=2*residualsSummaryMap->getNbinsY();
summaryReportMap->setBinContent(1,4, residualsSummary);
residualsSummary=0;
for(int i=1; i<=1; i++)
for (int j=1; j<=residualsSummaryMap->getNbinsY(); j++)
residualsSummary += residualsSummaryMap->getBinContent(i,j);
residualsSummary /=1*residualsSummaryMap->getNbinsY();
summaryReportMap->setBinContent(2,4,residualsSummary);
residualsSummary=0;
for(int i=2; i<=2; i++)
for (int j=1; j<=residualsSummaryMap->getNbinsY(); j++)
residualsSummary += residualsSummaryMap->getBinContent(i,j);
residualsSummary /=1*residualsSummaryMap->getNbinsY();
summaryReportMap->setBinContent(3,4,residualsSummary);
//--
//-- modified GH
float idtest=0;
for(int i=1; i<=2; i++)
for(int j=1; j<=5; j++) {
if(j==3 || j==5) continue; //ignore pull widths for now
idtest+=muonIdSummaryMap->getBinContent(i,j);
}
// idtest/=10.;
idtest/=6.;
summaryReportMap->setBinContent(1,5, idtest);
idtest=0;
for(int i=3; i<=4; i++)
for(int j=1; j<=5; j++) {
if(j==3 || j==5) continue; //ignore pull widths for now
idtest+=muonIdSummaryMap->getBinContent(i,j);
}
// idtest/=10.;
idtest/=6.;
summaryReportMap->setBinContent(2,5,idtest);
summaryReportMap->setBinContent(3,5,-1.0/6.0);
//--
//summaryReportMap->setBinContent(1,6,double(energySummaryMap->getBinContent(1,1)+energySummaryMap->getBinContent(1,2)+energySummaryMap->getBinContent(1,3))/3.0);
//summaryReportMap->setBinContent(2,6,double(energySummaryMap->getBinContent(2,1)+energySummaryMap->getBinContent(2,2)+energySummaryMap->getBinContent(2,3))/3.0);
//summaryReportMap->setBinContent(3,6,double(energySummaryMap->getBinContent(3,1)+energySummaryMap->getBinContent(3,2)+energySummaryMap->getBinContent(3,3))/3.0);
summaryReportMap->setBinContent(1,6,double(energySummaryMap->getBinContent(1,1)+energySummaryMap->getBinContent(1,2))/2.0);
summaryReportMap->setBinContent(2,6,double(energySummaryMap->getBinContent(2,1)+energySummaryMap->getBinContent(2,2))/2.0);
summaryReportMap->setBinContent(3,6,double(energySummaryMap->getBinContent(3,1)+energySummaryMap->getBinContent(3,2))/2.0);
summaryReportMap->setBinContent(1,7,multiplicitySummaryMap->getBinContent(1));
summaryReportMap->setBinContent(2,7,multiplicitySummaryMap->getBinContent(2));
summaryReportMap->setBinContent(3,7,multiplicitySummaryMap->getBinContent(3));
double kinematics_GLB = double(summaryReportMap->getBinContent(1,1)+summaryReportMap->getBinContent(1,2)+summaryReportMap->getBinContent(1,3))/3.0;
theSummaryContents[0]->Fill(kinematics_GLB);
double muonId_GLB = double(summaryReportMap->getBinContent(1,5));
theSummaryContents[1]->Fill(muonId_GLB);
double residuals_GLB = double(summaryReportMap->getBinContent(1,4));
theSummaryContents[2]->Fill(residuals_GLB);
double GLB = (kinematics_GLB+muonId_GLB+residuals_GLB)/3.0;
theSummaryContents[3]->Fill(GLB);
double kinematics_TK = double(summaryReportMap->getBinContent(2,1)+summaryReportMap->getBinContent(2,2)+summaryReportMap->getBinContent(2,3))/3.0;
theSummaryContents[4]->Fill(kinematics_TK);
double muonId_TK = double(summaryReportMap->getBinContent(2,5));
theSummaryContents[5]->Fill(muonId_TK);
double residuals_TK = double(summaryReportMap->getBinContent(2,4));
theSummaryContents[6]->Fill(residuals_TK);
double TK = double(kinematics_TK+muonId_TK+residuals_TK)/3.0;
theSummaryContents[7]->Fill(TK);
double kinematics_STA = double(summaryReportMap->getBinContent(3,1)+summaryReportMap->getBinContent(3,2)+summaryReportMap->getBinContent(3,3))/3.0;
theSummaryContents[8]->Fill(kinematics_STA);
double residuals_STA = double(summaryReportMap->getBinContent(3,4));
theSummaryContents[9]->Fill(residuals_STA);
double STA = double(kinematics_STA+residuals_STA)/2.0;
theSummaryContents[10]->Fill(STA);
double energyDeposits = double(summaryReportMap->getBinContent(1,6)+summaryReportMap->getBinContent(2,6)+summaryReportMap->getBinContent(3,6))/3.0;
theSummaryContents[11]->Fill(energyDeposits);
double multiplicity = double(summaryReportMap->getBinContent(1,7)+summaryReportMap->getBinContent(2,7)+summaryReportMap->getBinContent(3,7))/3.0;
theSummaryContents[12]->Fill(multiplicity);
summaryReport->Fill((GLB+TK+STA+energyDeposits+multiplicity)/5.0);
//global barrel:
float muonIDsummary=0;
// for(int i=2; i<=5; i++)
// muonIDsummary += muonIdSummaryMap->getBinContent(2, i);
// summaryCertificationMap->setBinContent(4, 5, muonIDsummary/4.);
//for now, just report the mean:
muonIDsummary += muonIdSummaryMap->getBinContent(1, 2);
muonIDsummary += muonIdSummaryMap->getBinContent(1, 4);
summaryCertificationMap->setBinContent(4, 5, muonIDsummary/2.);
//global EC:
muonIDsummary=0;
// for(int i=2; i<=5; i++)
// muonIDsummary += muonIdSummaryMap->getBinContent(2, i);
// summaryCertificationMap->setBinContent(7, 5, muonIDsummary/4.);
muonIDsummary += muonIdSummaryMap->getBinContent(2, 2);
muonIDsummary += muonIdSummaryMap->getBinContent(2, 4);
summaryCertificationMap->setBinContent(7, 5, muonIDsummary/2.);
//tracker barrel:
muonIDsummary=0;
// for(int i=2; i<=5; i++)
// muonIDsummary += muonIdSummaryMap->getBinContent(3, i);
// summaryCertificationMap->setBinContent(5, 5, muonIDsummary/4.);
muonIDsummary += muonIdSummaryMap->getBinContent(3, 2);
muonIDsummary += muonIdSummaryMap->getBinContent(3, 4);
summaryCertificationMap->setBinContent(5, 5, muonIDsummary/2.);
//tracker EC:
muonIDsummary=0;
// for(int i=2; i<=5; i++)
// muonIDsummary += muonIdSummaryMap->getBinContent(4, i);
// summaryCertificationMap->setBinContent(8, 5, muonIDsummary/4.);
muonIDsummary += muonIdSummaryMap->getBinContent(4, 2);
muonIDsummary += muonIdSummaryMap->getBinContent(4, 4);
summaryCertificationMap->setBinContent(8, 5, muonIDsummary/2.);
double muonId_GLB_B = double(summaryCertificationMap->getBinContent(4,5));
theCertificationContents[3]->Fill(muonId_GLB_B);
double muonId_GLB_EC = double(summaryCertificationMap->getBinContent(7,5));
theCertificationContents[6]->Fill(muonId_GLB_EC);
double muonId_TK_B = double(summaryCertificationMap->getBinContent(5,5));
theCertificationContents[5]->Fill(muonId_TK_B);
double muonId_TK_EC = double(summaryCertificationMap->getBinContent(8,5));
theCertificationContents[8]->Fill(muonId_TK_EC);
}
| void MuonTestSummary::GaussFit | ( | std::string | type, |
| std::string | parameter, | ||
| MonitorElement * | Histo, | ||
| float & | mean, | ||
| float & | mean_err, | ||
| float & | sigma, | ||
| float & | sigma_err | ||
| ) | [protected] |
Definition at line 780 of file MuonTestSummary.cc.
References MonitorElement::getMean(), MonitorElement::getRMS(), MonitorElement::getTH1F(), LogTrace, and metname.
{
// Gaussian Fit
float statMean = Histo->getMean(1);
float statSigma = Histo->getRMS(1);
TH1F * histo_root = Histo->getTH1F();
if(histo_root->GetEntries()>20){
TF1 *gfit = new TF1("Gaussian","gaus",(statMean-(2*statSigma)),(statMean+(2*statSigma)));
try {
histo_root->Fit(gfit, "Q0");
} catch (...) {
edm::LogError (metname)<< "[MuonTestSummary]: Exception when fitting Res_"<<type<<"_"<<parameter;
}
if(gfit){
mean = gfit->GetParameter(1);
mean_err = gfit->GetParErrors()[2];
sigma = gfit->GetParameter(2);
sigma_err = gfit->GetParErrors()[2];
LogTrace(metname)<<"Gaussian fit mean: "<<mean<<" +- "<<mean_err<<" for "<<type<<"_"<<parameter<<endl;
LogTrace(metname)<<"Gaussina fit sigma: "<<sigma<<" +- "<<sigma_err<<" for "<<type<<"_"<<parameter<<endl;
}
}
else{
LogTrace(metname) << "[MuonTestSummary]: Test of Res_"<<type<<"_"<<parameter<< " not performed because # entries < 20 ";
//auto-pass if not enough events.
mean=1;
mean_err=1;
sigma=1;
sigma_err=1;
}
}
| void MuonTestSummary::ResidualCheck | ( | std::string | muType, |
| std::vector< std::string > | resHistos, | ||
| int & | numPlot, | ||
| double & | Mean, | ||
| double & | Mean_err, | ||
| double & | Sigma, | ||
| double & | Sigma_err | ||
| ) | [protected] |
Definition at line 1119 of file MuonTestSummary.cc.
References MonitorElement::getTH1F(), LogTrace, plotscripts::mean(), metname, AlCaRecoCosmics_cfg::name, ExpressReco_HICollisions_FallBack::sigma, and mathSSE::sqrt().
Referenced by doMuonIDTests().
{
numPlot=0;
Mean=0;
Mean_err=0;
Sigma=0;
Sigma_err=0;
for(uint name=0; name<resHistos.size(); name++){
MonitorElement * resHisto = dbe->get("Muons/MuonIdDQM/" + muType + "/"+resHistos[name]);
if(resHisto){
TH1F * resHisto_root = resHisto->getTH1F();
if(resHisto_root->GetEntries() < 20) {
LogTrace(metname) << "[MuonTestSummary]: Test of "<< muType<<" for " <<resHistos[name]<< " not performed because # entries < 20 ";
continue;
}
//we also want to check if the peak is away from zero.
//so, try fitting in 3 sigma around the histogram mean.
//alternatively, could use the maximum bin (less sensitive to 1-sided tails).
// float mean = resHisto_root->GetMean();
float mean = resHisto_root->GetBinLowEdge(resHisto_root->GetMaximumBin());
TF1 *gfit = new TF1("Gaussian","gaus",mean-3,mean+3);
try {
resHisto_root->Fit(gfit, "Q0");
} catch (...) {
edm::LogError (metname)<< "[MuonTestSummary]: Exception when fitting "<<resHistos[name];
}
if(gfit){
double mean = gfit->GetParameter(1);
double mean_err = gfit->GetParError(1);
double sigma = gfit->GetParameter(2);
double sigma_err = gfit->GetParError(2);
LogTrace(metname)<<"meanRes: "<<mean<<" +- "<<mean_err<<" for "<<resHistos[name]<<endl;
LogTrace(metname)<<"sigmaRes: "<<sigma<<" +- "<<sigma_err<<" for "<<resHistos[name]<<endl;
Mean+=mean;
Mean_err +=mean_err * mean_err;
Sigma+=sigma;
Sigma_err +=sigma_err * sigma_err;
numPlot++;
} //if gfit? why would we not have gfit?
}//histogram exists...
} // loop over residuals histos
if(numPlot==0){ //eg no stats
Mean_err = 1;
Mean=1;
Sigma_err =1;
Sigma = 1;
}
else{
Mean_err = sqrt(Mean_err);
Mean_err/=numPlot;
Mean/=numPlot;
Sigma_err = sqrt(Sigma_err);
Sigma_err/=numPlot;
Sigma/=numPlot;
}
return;
}
double MuonTestSummary::chi2Fraction [private] |
Definition at line 79 of file MuonTestSummary.h.
double MuonTestSummary::chi2Spread [private] |
Definition at line 80 of file MuonTestSummary.h.
Definition at line 118 of file MuonTestSummary.h.
DQMStore* MuonTestSummary::dbe [private] |
Definition at line 72 of file MuonTestSummary.h.
MonitorElement* MuonTestSummary::energySummaryMap [private] |
Definition at line 123 of file MuonTestSummary.h.
double MuonTestSummary::etaExpected [private] |
Definition at line 77 of file MuonTestSummary.h.
double MuonTestSummary::expMultiplicityGlb_max [private] |
Definition at line 111 of file MuonTestSummary.h.
double MuonTestSummary::expMultiplicityGlb_min [private] |
Definition at line 108 of file MuonTestSummary.h.
double MuonTestSummary::expMultiplicitySta_max [private] |
Definition at line 113 of file MuonTestSummary.h.
double MuonTestSummary::expMultiplicitySta_min [private] |
Definition at line 110 of file MuonTestSummary.h.
double MuonTestSummary::expMultiplicityTk_max [private] |
Definition at line 112 of file MuonTestSummary.h.
double MuonTestSummary::expMultiplicityTk_min [private] |
Definition at line 109 of file MuonTestSummary.h.
double MuonTestSummary::expPeakEcalS9_max [private] |
Definition at line 105 of file MuonTestSummary.h.
double MuonTestSummary::expPeakEcalS9_min [private] |
Definition at line 104 of file MuonTestSummary.h.
double MuonTestSummary::expPeakHadS9_max [private] |
Definition at line 107 of file MuonTestSummary.h.
double MuonTestSummary::expPeakHadS9_min [private] |
Definition at line 106 of file MuonTestSummary.h.
Definition at line 120 of file MuonTestSummary.h.
Definition at line 117 of file MuonTestSummary.h.
double MuonTestSummary::matchesFractionCsc_max [private] |
Definition at line 98 of file MuonTestSummary.h.
double MuonTestSummary::matchesFractionCsc_min [private] |
Definition at line 97 of file MuonTestSummary.h.
double MuonTestSummary::matchesFractionDt_max [private] |
Definition at line 96 of file MuonTestSummary.h.
double MuonTestSummary::matchesFractionDt_min [private] |
Definition at line 95 of file MuonTestSummary.h.
std::string MuonTestSummary::metname [private] |
Definition at line 74 of file MuonTestSummary.h.
Definition at line 124 of file MuonTestSummary.h.
MonitorElement* MuonTestSummary::muonIdSummaryMap [private] |
Definition at line 122 of file MuonTestSummary.h.
double MuonTestSummary::numMatchedExpected_max [private] |
Definition at line 94 of file MuonTestSummary.h.
double MuonTestSummary::numMatchedExpected_min [private] |
Definition at line 93 of file MuonTestSummary.h.
double MuonTestSummary::phiExpected [private] |
Definition at line 78 of file MuonTestSummary.h.
double MuonTestSummary::pullEtaSpread [private] |
Definition at line 87 of file MuonTestSummary.h.
double MuonTestSummary::pullOneOvPSpread [private] |
Definition at line 89 of file MuonTestSummary.h.
double MuonTestSummary::pullPhiSpread [private] |
Definition at line 88 of file MuonTestSummary.h.
double MuonTestSummary::resChargeLimit_glbSta [private] |
Definition at line 91 of file MuonTestSummary.h.
double MuonTestSummary::resChargeLimit_tkGlb [private] |
Definition at line 90 of file MuonTestSummary.h.
double MuonTestSummary::resChargeLimit_tkSta [private] |
Definition at line 92 of file MuonTestSummary.h.
double MuonTestSummary::resEtaSpread_glbSta [private] |
Definition at line 82 of file MuonTestSummary.h.
double MuonTestSummary::resEtaSpread_tkGlb [private] |
Definition at line 81 of file MuonTestSummary.h.
Definition at line 121 of file MuonTestSummary.h.
double MuonTestSummary::resOneOvPSpread_glbSta [private] |
Definition at line 86 of file MuonTestSummary.h.
double MuonTestSummary::resOneOvPSpread_tkGlb [private] |
Definition at line 85 of file MuonTestSummary.h.
double MuonTestSummary::resPhiSpread_glbSta [private] |
Definition at line 84 of file MuonTestSummary.h.
double MuonTestSummary::resPhiSpread_tkGlb [private] |
Definition at line 83 of file MuonTestSummary.h.
double MuonTestSummary::resSegmTrack_mean_max [private] |
Definition at line 102 of file MuonTestSummary.h.
double MuonTestSummary::resSegmTrack_mean_min [private] |
Definition at line 101 of file MuonTestSummary.h.
double MuonTestSummary::resSegmTrack_rms_max [private] |
Definition at line 100 of file MuonTestSummary.h.
double MuonTestSummary::resSegmTrack_rms_min [private] |
Definition at line 99 of file MuonTestSummary.h.
double MuonTestSummary::sigmaResSegmTrackExp [private] |
Definition at line 103 of file MuonTestSummary.h.
Definition at line 128 of file MuonTestSummary.h.
Definition at line 129 of file MuonTestSummary.h.
MonitorElement* MuonTestSummary::summaryReport [private] |
Definition at line 125 of file MuonTestSummary.h.
MonitorElement* MuonTestSummary::summaryReportMap [private] |
Definition at line 126 of file MuonTestSummary.h.
std::vector<MonitorElement*> MuonTestSummary::theCertificationContents [private] |
Definition at line 130 of file MuonTestSummary.h.
std::vector<MonitorElement*> MuonTestSummary::theSummaryContents [private] |
Definition at line 127 of file MuonTestSummary.h.
1.7.3