|
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
|
/** \file
*
* \author N. Amapane - Torino
*/
#include "CondFormats/MFObjects/interface/MagFieldConfig.h"
#include <FWCore/MessageLogger/interface/MessageLogger.h>
#include <iostream>
#include <vector>
#include <memory>
#include <boost/algorithm/string/split.hpp>
#include <boost/algorithm/string/classification.hpp>
using namespace std;
using namespace magneticfield;
MagFieldConfig::MagFieldConfig(const edm::ParameterSet& pset, bool debug) {
version = pset.getParameter<std::string>("version");
geometryVersion = pset.getParameter<int>("geometryVersion");
// Get specification for the grid tables to be used.
typedef vector<edm::ParameterSet> VPSet;
VPSet fileSpec = pset.getParameter<VPSet>("gridFiles");
if (!fileSpec.empty()) {
for (VPSet::const_iterator rule = fileSpec.begin(); rule != fileSpec.end(); ++rule) {
string s_volumes = rule->getParameter<string>("volumes");
string s_sectors = rule->getParameter<string>("sectors"); // 0 means all volumes
int master = rule->getParameter<int>("master");
string path = rule->getParameter<string>("path");
vector<unsigned> volumes = expandList(s_volumes);
vector<unsigned> sectors = expandList(s_sectors);
if (debug) {
cout << "Volumes: " << s_volumes << " Sectors: " << s_sectors << " Master: " << master << " Path: " << path
<< endl;
cout << " Expanded volumes: ";
copy(volumes.begin(), volumes.end(), ostream_iterator<unsigned>(cout, " "));
cout << endl;
cout << " Expanded sectors: ";
copy(sectors.begin(), sectors.end(), ostream_iterator<unsigned>(cout, " "));
cout << endl;
}
for (vector<unsigned>::iterator i = volumes.begin(); i != volumes.end(); ++i) {
for (vector<unsigned>::iterator j = sectors.begin(); j != sectors.end(); ++j) {
unsigned vpacked = (*i) * 100 + (*j);
if (gridFiles.find(vpacked) == gridFiles.end()) {
gridFiles[vpacked] = make_pair(path, master);
} else {
throw cms::Exception("ConfigurationError")
<< "VolumeBasedMagneticFieldESProducer: malformed gridFiles config parameter" << endl;
}
}
}
}
}
// Get scaling factors
keys = pset.getParameter<vector<int> >("scalingVolumes");
values = pset.getParameter<vector<double> >("scalingFactors");
// Slave field label. Either a label of an existing map (legacy support), or the
// type of parametrization to be constructed with the "paramData" parameters.
slaveFieldVersion = pset.getParameter<string>("paramLabel");
// Check for compatibility with older configurations
if (pset.existsAs<vector<double> >("paramData")) {
slaveFieldParameters = pset.getParameter<vector<double> >("paramData");
}
}
vector<unsigned> MagFieldConfig::expandList(const string& list) {
typedef vector<string> vstring;
vector<unsigned> values;
vstring v1;
boost::split(v1, list, boost::is_any_of(","));
for (vstring::const_iterator i = v1.begin(); i != v1.end(); ++i) {
vstring v2;
boost::split(v2, *i, boost::is_any_of("-"));
unsigned start = std::stoul(v2.front());
unsigned end = std::stoul(v2.back());
if ((v2.size() > 2) || (start > end)) {
throw cms::Exception("ConfigurationError")
<< "VolumeBasedMagneticFieldESProducerFromDB: malformed configuration" << list << endl;
}
for (unsigned k = start; k <= end; ++k) {
values.push_back(k);
}
}
return values;
}
|