|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
|
#include "DQM/SiPixelMonitorClient/interface/SiPixelCertification.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "CondFormats/DataRecord/interface/RunSummaryRcd.h"
#include "CondFormats/RunInfo/interface/RunInfo.h"
#include "CondFormats/RunInfo/interface/RunSummary.h"
#include <cmath>
#include <cstdio>
#include <iomanip>
#include <iostream>
#include <sstream>
#include <string>
using namespace std;
using namespace edm;
SiPixelCertification::SiPixelCertification(const edm::ParameterSet &ps) {
edm::LogInfo("SiPixelCertification") << "SiPixelCertification::Creating SiPixelCertification ";
firstLumi = true;
}
SiPixelCertification::~SiPixelCertification() {
edm::LogInfo("SiPixelCertification") << "SiPixelCertification::Deleting SiPixelCertification ";
}
void SiPixelCertification::dqmEndLuminosityBlock(DQMStore::IBooker &iBooker,
DQMStore::IGetter &iGetter,
const edm::LuminosityBlock &lumiBlock,
const edm::EventSetup &iSetup) {
// cout<<"Entering SiPixelCertification::endLuminosityBlock: "<<endl;
// If first run, book some histograms
if (firstLumi) {
iBooker.setCurrentFolder("Pixel/EventInfo");
CertificationPixel = iBooker.bookFloat("CertificationSummary");
iBooker.setCurrentFolder("Pixel/EventInfo/CertificationContents");
CertificationBarrel = iBooker.bookFloat("PixelBarrelFraction");
CertificationEndcap = iBooker.bookFloat("PixelEndcapFraction");
CertificationPixel->Fill(1.);
CertificationBarrel->Fill(1.);
CertificationEndcap->Fill(1.);
firstLumi = false;
}
edm::LogInfo("SiPixelCertification") << "SiPixelCertification::endLuminosityBlock ";
// Compute and fill overall certification bits, for now use smallest single
// value:
float dcsFrac = (iGetter.getElement("Pixel/EventInfo/DCSSummary"))->getFloatValue();
float daqFrac = (iGetter.getElement("Pixel/EventInfo/DAQSummary"))->getFloatValue();
float dqmFrac = (iGetter.getElement("Pixel/EventInfo/reportSummary"))->getFloatValue();
float pixel_all = std::min(dcsFrac, daqFrac);
pixel_all = std::min(pixel_all, dqmFrac);
// std::cout<<"Pixel numbers: "<<dcsFrac<<" , "<<daqFrac<<" , "<<dqmFrac<<" ,
// "<<pixel_all<<std::endl;
CertificationPixel = iGetter.getElement("Pixel/EventInfo/CertificationSummary");
if (CertificationPixel)
CertificationPixel->Fill(pixel_all);
dcsFrac = (iGetter.getElement("Pixel/EventInfo/DCSContents/PixelBarrelFraction"))->getFloatValue();
daqFrac = (iGetter.getElement("Pixel/EventInfo/DAQContents/PixelBarrelFraction"))->getFloatValue();
dqmFrac = (iGetter.getElement("Pixel/EventInfo/reportSummaryContents/PixelBarrelFraction"))->getFloatValue();
float pixel_barrel = std::min(dcsFrac, daqFrac);
pixel_barrel = std::min(pixel_barrel, dqmFrac);
// std::cout<<"Barrel numbers: "<<dcsFrac<<" , "<<daqFrac<<" , "<<dqmFrac<<" ,
// "<<pixel_barrel<<std::endl;
CertificationBarrel = iGetter.getElement("Pixel/EventInfo/CertificationContents/PixelBarrelFraction");
if (CertificationBarrel)
CertificationBarrel->Fill(pixel_barrel);
dcsFrac = (iGetter.getElement("Pixel/EventInfo/DCSContents/PixelEndcapFraction"))->getFloatValue();
daqFrac = (iGetter.getElement("Pixel/EventInfo/DAQContents/PixelEndcapFraction"))->getFloatValue();
dqmFrac = (iGetter.getElement("Pixel/EventInfo/reportSummaryContents/PixelEndcapFraction"))->getFloatValue();
float pixel_endcap = std::min(dcsFrac, daqFrac);
pixel_endcap = std::min(pixel_endcap, dqmFrac);
// std::cout<<"Endcap numbers: "<<dcsFrac<<" , "<<daqFrac<<" , "<<dqmFrac<<" ,
// "<<pixel_endcap<<std::endl;
CertificationEndcap = iGetter.getElement("Pixel/EventInfo/CertificationContents/PixelEndcapFraction");
if (CertificationEndcap)
CertificationEndcap->Fill(pixel_endcap);
}
void SiPixelCertification::dqmEndJob(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter) {
// cout<<"Entering SiPixelCertification::endRun: "<<endl;
edm::LogInfo("SiPixelCertification") << "SiPixelCertification::endRun ";
// Compute and fill overall certification bits, for now use smallest single
// value:
float dcsFrac = (iGetter.getElement("Pixel/EventInfo/DCSSummary"))->getFloatValue();
float daqFrac = (iGetter.getElement("Pixel/EventInfo/DAQSummary"))->getFloatValue();
float dqmFrac = (iGetter.getElement("Pixel/EventInfo/reportSummary"))->getFloatValue();
float pixel_all = std::min(dcsFrac, daqFrac);
pixel_all = std::min(pixel_all, dqmFrac);
// std::cout<<"Pixel numbers: "<<dcsFrac<<" , "<<daqFrac<<" , "<<dqmFrac<<" ,
// "<<pixel_all<<std::endl;
if (CertificationPixel)
CertificationPixel->Fill(pixel_all);
dcsFrac = (iGetter.getElement("Pixel/EventInfo/DCSContents/PixelBarrelFraction"))->getFloatValue();
daqFrac = (iGetter.getElement("Pixel/EventInfo/DAQContents/PixelBarrelFraction"))->getFloatValue();
dqmFrac = (iGetter.getElement("Pixel/EventInfo/reportSummaryContents/PixelBarrelFraction"))->getFloatValue();
float pixel_barrel = std::min(dcsFrac, daqFrac);
pixel_barrel = std::min(pixel_barrel, dqmFrac);
// std::cout<<"Barrel numbers: "<<dcsFrac<<" , "<<daqFrac<<" , "<<dqmFrac<<" ,
// "<<pixel_barrel<<std::endl;
if (CertificationBarrel)
CertificationBarrel->Fill(pixel_barrel);
dcsFrac = (iGetter.getElement("Pixel/EventInfo/DCSContents/PixelEndcapFraction"))->getFloatValue();
daqFrac = (iGetter.getElement("Pixel/EventInfo/DAQContents/PixelEndcapFraction"))->getFloatValue();
dqmFrac = (iGetter.getElement("Pixel/EventInfo/reportSummaryContents/PixelEndcapFraction"))->getFloatValue();
float pixel_endcap = std::min(dcsFrac, daqFrac);
pixel_endcap = std::min(pixel_endcap, dqmFrac);
// std::cout<<"Endcap numbers: "<<dcsFrac<<" , "<<daqFrac<<" , "<<dqmFrac<<" ,
// "<<pixel_endcap<<std::endl;
if (CertificationEndcap)
CertificationEndcap->Fill(pixel_endcap);
}
|