ATLAS Offline Software
Loading...
Searching...
No Matches
Public Member Functions | Private Attributes | List of all members
dqm_algorithms::MDTOverview Class Reference

#include <MDTOverview.h>

Inheritance diagram for dqm_algorithms::MDTOverview:
Collaboration diagram for dqm_algorithms::MDTOverview:

Public Member Functions

 MDTOverview (const std::string &name)
virtual ~MDTOverview ()
virtual dqm_core::Algorithm * clone ()
virtual dqm_core::Result * execute (const std::string &name, const TObject &object, const dqm_core::AlgorithmConfig &config)
virtual void printDescription (std::ostream &out)

Private Attributes

std::string m_name

Detailed Description

Definition at line 20 of file MDTOverview.h.

Constructor & Destructor Documentation

◆ MDTOverview()

dqm_algorithms::MDTOverview::MDTOverview ( const std::string & name)

Definition at line 36 of file MDTOverview.cxx.

37 : m_name ( name )
38{
39 dqm_core::AlgorithmManager::instance().registerAlgorithm( name, this );
40}
dqm_algorithms::MDTOverview::m_name
std::string m_name
Definition MDTOverview.h:33

◆ ~MDTOverview()

dqm_algorithms::MDTOverview::~MDTOverview ( )
virtual

Definition at line 43 of file MDTOverview.cxx.

44{
45}

Member Function Documentation

◆ clone()

dqm_core::Algorithm * dqm_algorithms::MDTOverview::clone ( )
virtual

Definition at line 49 of file MDTOverview.cxx.

50{
51 return new MDTOverview(m_name);
52}
dqm_algorithms::MDTOverview::MDTOverview
MDTOverview(const std::string &name)
Definition MDTOverview.cxx:36

◆ execute()

dqm_core::Result * dqm_algorithms::MDTOverview::execute ( const std::string & name,
const TObject & object,
const dqm_core::AlgorithmConfig & config )
virtual

Definition at line 56 of file MDTOverview.cxx.

57{
58 const TH2 * hist;
59 const TH2 * ref;
60
61 if( object.IsA()->InheritsFrom( "TH2" ) ) {
62 hist = static_cast<const TH2*>(&object);
63 if (hist->GetDimension() != 2 ){
64 throw dqm_core::BadConfig( ERS_HERE, name, "dimension != 2 " );
65 }
66 } else {
67 throw dqm_core::BadConfig( ERS_HERE, name, "does not inherit from TH2" );
68 }
69
70 //Get Parameters and Thresholds
71
72 double greenTh;
73 double redTh;
74 double minstat;
75 double thresh;
76 double using_ref_for_flag;
77 std::string thresholdname;
78 if (m_name == "MDTOverview_Global") {
79 thresholdname="Percentage";
80 } else if (m_name == "MDTOverview_Station") {
81 thresholdname="Limits";
82 } else {
83 throw dqm_core::BadConfig( ERS_HERE, "None", m_name );
84 }
85
86 try {
87 minstat = dqm_algorithms::tools::GetFirstFromMap("MinStat", config.getParameters(), 5000);
88 thresh = dqm_algorithms::tools::GetFirstFromMap("thresh", config.getParameters(), 10);
89 using_ref_for_flag = dqm_algorithms::tools::GetFirstFromMap("using_ref_for_flag", config.getParameters(), 0);
90 redTh = dqm_algorithms::tools::GetFromMap( thresholdname, config.getRedThresholds());
91 greenTh = dqm_algorithms::tools::GetFromMap( thresholdname, config.getGreenThresholds() );
92 }
93 catch ( dqm_core::Exception & ex ) {
94 throw dqm_core::BadConfig( ERS_HERE, name, ex.what(), ex );
95 }
96
97 //Get Reference Histo
98 try {
99 ref = static_cast<const TH2 *>( config.getReference() );
100 }
101 catch ( dqm_core::Exception & ex ) {
102 throw dqm_core::BadRefHist(ERS_HERE,name," Could not retreive reference");
103 }
104 if (hist->GetDimension() != ref->GetDimension() ) {
105 throw dqm_core::BadRefHist( ERS_HERE, name, "Reference VS histo: Different dimension!" );
106 }
107
108
109 //Check of statistics
110 if (hist->GetEntries() < minstat ) {
111 ERS_INFO("Histogram does not satisfy MinStat requirement " <<hist->GetName());
112 dqm_core::Result *result = new dqm_core::Result(dqm_core::Result::Undefined);
113 result->tags_["InsufficientEntries"] = hist->GetEntries();
114 return result;
115 }
116 ERS_DEBUG(1,"Statistics: "<< hist->GetEntries()<< " entries ");
117
118
119 //Algo
120 int binX = hist->GetNbinsX();
121 int binY = hist->GetNbinsY();
122 double nML=0;
123
124 std::vector<int> hist_buffer;
125 std::vector<int> ref_buffer;
126
127 for(int x_index=1;x_index<=binX;x_index++){
128 for(int y_index=1;y_index<=binY;y_index++){
129 if( hist->GetBinContent(x_index,y_index) != 0 ) hist_buffer.push_back((int)hist->GetBinContent(x_index,y_index));
130 if( ref->GetBinContent(x_index,y_index) != 0 ) ref_buffer.push_back((int)ref->GetBinContent(x_index,y_index));
131 };
132 };
133
134 const char *htitle = hist->GetTitle();
135 ERS_DEBUG(1,"Histogram title: "<< htitle);
136
137 if (m_name == "MDTOverview_Global") {
138 //GLOBAL htitle="Global HitsPerML Endcap"
139 //if (strstr(htitle, "Global")) {
140 if (strstr(htitle, "BA")) {
141 nML = 594; //4 BOG0 included
142 } else if (strstr(htitle, "BC")) {
143 nML = 590;
144 } else if (strstr(htitle, "EA")) {
145 nML = 502; // 16 BIS7 + 8 BIS8 +16 BEE + 10 EEL included
146 } else if (strstr(htitle, "EC")) {
147 nML = 502; // BIS7/8 + BEE + EEL included
148 }
149
150 } else if (m_name == "MDTOverview_Station") {
151 // Barrel htitle="HitsPerML_BX_Station"
152 if (strstr(htitle, "B")) {
153 if (strstr(htitle,"Inner")) {
154 ERS_DEBUG(1,"Station Inner ");
155 nML = 212;
156 } else if (strstr(htitle,"Middle")) {
157 ERS_DEBUG(1,"Station Middle ");
158 nML = 178;
159 } else if (strstr(htitle,"Outer")) {
160 ERS_DEBUG(1,"Station Outer ");
161 if (strstr(htitle, "BA")) {
162 nML = 204; // BOG0 included
163 } else if (strstr(htitle, "BC")) {
164 nML = 200;
165 }
166 }
167 // Endcap htitle="HitsPerML_EX_Station"
168 } else if (strstr(htitle, "E")) {
169
170 if (strstr(htitle,"Inner")) {
171 ERS_DEBUG(1,"Station Inner ");
172 nML = 124;// 16 BIS7 + 8 BIS8 included
173 } else if (strstr(htitle,"Middle")) {
174 ERS_DEBUG(1,"Station Middle ");
175 nML = 160;
176 } else if (strstr(htitle,"Outer")) {
177 ERS_DEBUG(1,"Station Outer ");
178 nML = 192;
179 } else if (strstr(htitle,"extra")) {
180 ERS_DEBUG(1,"Extra chambers (BEE, EEL, BIS7, BIS8)");
181 nML = 26;// 16 BEE + 10 EEL included
182 }
183 }
184 }
185
186 auto median=[](std::vector<int> & v)->int {
187 const auto midPoint = v.begin()+v.size()/2;
188 std::nth_element(v.begin(), midPoint, v.end());
189 return v[v.size()/2];
190 };
191
192
193
194 int mediana_hist = median(hist_buffer);
195 int mediana_ref = median(ref_buffer);
196 const int size_hist = hist_buffer.size();
197
198 int count=0;
199 int count_non_0=0;
200 int sec_count=0;
201 int ML_count=0;
202 int sector=1;
203 std::pair<int,int> empty_bin;
204 std::vector< std::pair<int,int> > new_empty_bins;
205
206 for(int y_index=1;y_index<=binY;y_index++){
207 if(ML_count==2){
208 ML_count=0;
209 sec_count=0;
210 sector++;
211 };
212 ML_count++;
213
214 for(int x_index=1;x_index<=binX;x_index++){
215 if(using_ref_for_flag) {
216 if( hist->GetBinContent(x_index,y_index)< mediana_hist*thresh/100 && ref->GetBinContent(x_index,y_index) >= mediana_ref*thresh/100){
217 sec_count++;
218 count++;
219 if( hist->GetBinContent(x_index,y_index)>0) count_non_0++;
220 }
221 } else {
222 if( hist->GetBinContent(x_index,y_index)< mediana_hist*thresh/100){
223 sec_count++;
224 count++;
225 };
226 }
227 };
228
229 if(ML_count==2 && sec_count>0){
230 empty_bin.first=sector;
231 empty_bin.second=sec_count;
232 new_empty_bins.push_back(empty_bin);
233 }
234 };
235
236 if(nML<=0){
237 dqm_core::Result* result = new dqm_core::Result();
238 result->tags_["Configuration_ERROR_Setted_N_tot_ML"] = nML;
239 result->status_ = dqm_core::Result::Undefined;
240 return result;
241 };
242
243 //Result
244 dqm_core::Result* result = new dqm_core::Result();
245 result->tags_["00-%_ML_ON"] = 100*((double)(size_hist-count_non_0)/nML);
246
247 if(thresholdname == "Limits" ) {
248
249 if(count>0) {
250 std::string tag;
251 std::string tag1="Sector";
252 std::string tag2="_number_of_ML_off";
253 result->tags_["01-NEW_ML_OFF"] = count;
254 for(int j=0; j<(int)new_empty_bins.size(); j++){
255 std::string formatted = std::format("{:02d}", new_empty_bins[j].first);
256 tag = "01-" + tag1 + formatted + tag2;
257 result->tags_[tag] = new_empty_bins[j].second;
258 };
259 }
260 if (nML-count-size_hist > 0) result->tags_["02-KNOWN_ML_OFF"] = nML-count-size_hist;
261
262 // Flag assignment
263 if (count > redTh) {
264 result->status_ = dqm_core::Result::Red;
265 ERS_DEBUG(1,"Red");
266 } else if ( count <= redTh && count > greenTh ) {
267 result->status_ = dqm_core::Result::Yellow;
268 ERS_DEBUG(1,"Yellow");
269 } else if (count <= greenTh) {
270 result->status_ = dqm_core::Result::Green;
271 ERS_DEBUG(1,"Green");
272 } else {
273 result->status_ = dqm_core::Result::Undefined;
274 ERS_DEBUG(1,"Undefined");
275 }
276
277 } else if (thresholdname == "Percentage") {
278 // Flag assignment
279 if (100*((double)size_hist/nML)<redTh) {
280 result->status_ = dqm_core::Result::Red;
281 ERS_DEBUG(1,"Red");
282 } else if ( 100*((double)size_hist/nML)>=redTh && 100*((double)size_hist/nML)<greenTh ) {
283 result->status_ = dqm_core::Result::Yellow;
284 ERS_DEBUG(1,"Yellow");
285 } else if (100*((double)size_hist/nML)>=greenTh) {
286 result->status_ = dqm_core::Result::Green;
287 ERS_DEBUG(1,"Green");
288 } else {
289 result->status_ = dqm_core::Result::Undefined;
290 ERS_DEBUG(1,"Undefined");
291 }
292 }
293 return result;
294}
ref
const boost::regex ref(r_ef)
count
int count(std::string s, const std::string &regx)
count how many occurances of a regx are in a string
Definition hcg.cxx:146
CaloCondBlobAlgs_fillNoiseFromASCII.tag
str tag
Definition CaloCondBlobAlgs_fillNoiseFromASCII.py:23
InDet::median
float median(std::vector< float > &Vec)
Definition BTagVrtSec.cxx:35
Trk::binX
@ binX
Definition BinningType.h:47
Trk::binY
@ binY
Definition BinningType.h:48
dqm_algorithms::tools::GetFirstFromMap
double GetFirstFromMap(const std::string &paramName, const std::map< std::string, double > &params)
Definition AlgorithmHelper.cxx:339
dqm_algorithms::tools::GetFromMap
const T & GetFromMap(const std::string &pname, const std::map< std::string, T > &params)
Definition AlgorithmHelper.h:114
get_generator_info.result
result
Definition get_generator_info.py:21
physval_make_web_display.thresh
thresh
Definition physval_make_web_display.py:35
python.base_data.config
config
Definition base_data.py:20
IsA
#define IsA
Declare the TObject style functions.
Definition xAODTEventBranch.h:59

◆ printDescription()

void dqm_algorithms::MDTOverview::printDescription ( std::ostream & out)
virtual

Definition at line 298 of file MDTOverview.cxx.

299{
300 std::string message;
301 message += "\n";
302 message += "Algorithm: \"" + m_name + "\"\n";
303 message += "Description: check if there are any bins with less then a thresh% of the entries of the median bin\n";
304 message += " and perform a comparison with reference.\n";
305 message += "MDTOverview_Global: Mandatory Green/Red Threshold: Percentage: of ML OFF\n";
306 message += "MDTOverview_Station: Mandatory Green/Red Threshold: Limits: number of ML OFF\n";
307 message += "Optional Parameters: thresh = min % of hits from the median chamber accepted\n";
308 message += " MinStat = Minimum histogram statistics needed to perform Algorithm\n";
309 message += " using_ref_for_flag = if enabled (1) it compares the OFF MLs with reference, if not enabled (0)\n";
310 message += " the flag is setted checking the absolute % of ML on in the detector\n";
311 message += "\n";
312
313 out << message;
314}

Member Data Documentation

◆ m_name

std::string dqm_algorithms::MDTOverview::m_name
private

Definition at line 33 of file MDTOverview.h.

The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.14.0