|
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
|
#include "Geometry/CaloTopology/interface/HcalTopologyRestrictionParser.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <boost/tokenizer.hpp>
#include <sstream>
#include <iostream>
HcalTopologyRestrictionParser::HcalTopologyRestrictionParser(HcalTopology& target) : target_(target) {}
static HcalSubdetector determineSubdet(const std::string& item) {
if (item == "HB")
return HcalBarrel;
if (item == "HE")
return HcalEndcap;
if (item == "HF")
return HcalForward;
if (item == "HO")
return HcalOuter;
return (HcalSubdetector)0;
}
std::string HcalTopologyRestrictionParser::parse(const std::string& line) {
std::ostringstream errors;
boost::char_separator<char> sep(" \t", ";");
typedef boost::tokenizer<boost::char_separator<char> > myTokType;
std::string totaline(line);
totaline += ';'; // terminate
myTokType tok(totaline, sep);
int ieta1 = 0, ieta2 = 0, iphi1 = -1, iphi2 = -1, depth1 = 1, depth2 = 4;
HcalSubdetector subdet = (HcalSubdetector)0;
int phase = 0;
for (myTokType::iterator beg = tok.begin(); beg != tok.end() && phase >= 0; ++beg) {
edm::LogVerbatim("HCalGeom") << phase << " : <" << *beg << ">";
if (*beg == ";") {
if (phase == 0)
continue; // empty
if (phase != 1 && phase != 5 && phase != 7) {
errors << "Expect 1, 5, or 7 arguments, got " << phase;
phase = -1;
} else {
if (phase == 1) { // reject whole subdetector...
target_.excludeSubdetector(subdet);
} else {
target_.exclude(subdet, ieta1, ieta2, iphi1, iphi2, depth1, depth2);
}
phase = 0;
}
} else {
switch (phase) {
case (0):
subdet = determineSubdet(*beg);
break;
case (1):
ieta1 = atoi(beg->c_str());
break;
case (2):
ieta2 = atoi(beg->c_str());
break;
case (3):
iphi1 = atoi(beg->c_str());
break;
case (4):
iphi2 = atoi(beg->c_str());
depth1 = 1;
depth2 = 4;
break; // also set defaults...
case (5):
depth1 = atoi(beg->c_str());
break;
case (6):
depth2 = atoi(beg->c_str());
break;
}
phase++;
}
}
return errors.str();
}
|