CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
MagFieldConfig.cc
Go to the documentation of this file.
1 
8 
9 #include <iostream>
10 #include <vector>
11 #include <memory>
12 
13 #include <boost/algorithm/string/split.hpp>
14 #include <boost/algorithm/string/classification.hpp>
15 #include <boost/lexical_cast.hpp>
16 
17 using namespace std;
18 using namespace magneticfield;
19 
20 
22 
23  version = pset.getParameter<std::string>("version");
24  geometryVersion = pset.getParameter<int>("geometryVersion");
25 
26 
27  // Get specification for the grid tables to be used.
28  typedef vector<edm::ParameterSet> VPSet;
29 
30  VPSet fileSpec = pset.getParameter<VPSet>("gridFiles");
31  if (fileSpec.size()!=0) {
32  for(VPSet::const_iterator rule = fileSpec.begin(); rule != fileSpec.end(); ++rule){
33  string s_volumes = rule->getParameter<string>("volumes");
34  string s_sectors = rule->getParameter<string>("sectors"); // 0 means all volumes
35  int master = rule->getParameter<int>("master");
36  string path = rule->getParameter<string>("path");
37 
38  vector<unsigned> volumes = expandList(s_volumes);
39  vector<unsigned> sectors = expandList(s_sectors);
40 
41  if (debug) {
42  cout << "Volumes: " << s_volumes <<" Sectors: " << s_sectors
43  << " Master: " << master << " Path: " << path << endl;
44  cout << " Expanded volumes: ";
45  copy(volumes.begin(), volumes.end(), ostream_iterator<unsigned>(cout, " "));
46  cout << endl;
47  cout << " Expanded sectors: ";
48  copy(sectors.begin(), sectors.end(), ostream_iterator<unsigned>(cout, " "));
49  cout << endl;
50  }
51 
52  for (vector<unsigned>::iterator i = volumes.begin(); i!=volumes.end(); ++i){
53  for (vector<unsigned>::iterator j = sectors.begin(); j!=sectors.end(); ++j){
54  unsigned vpacked = (*i)*100+(*j);
55  if (gridFiles.find(vpacked)==gridFiles.end()) {
56  gridFiles[vpacked] = make_pair(path, master);
57  } else {
58  throw cms::Exception("ConfigurationError") << "VolumeBasedMagneticFieldESProducer: malformed gridFiles config parameter" << endl;
59  }
60  }
61  }
62  }
63  }
64 
65  // Get scaling factors
66  keys = pset.getParameter<vector<int> >("scalingVolumes");
67  values = pset.getParameter<vector<double> >("scalingFactors");
68 
69 
70  // Slave field label. Either a label of an existing map (legacy support), or the
71  // type of parametrization to be constructed with the "paramData" parameters.
72  slaveFieldVersion = pset.getParameter<string>("paramLabel");
73  // Check for compatibility with older configurations
74  if (pset.existsAs<vector<double> >("paramData")) {
75  slaveFieldParameters = pset.getParameter<vector<double> >("paramData");
76  }
77 
78 }
79 
80 vector<unsigned> MagFieldConfig::expandList(const string& list) {
81  typedef vector<string> vstring;
82  vector<unsigned> values;
83  vstring v1;
84  boost::split(v1, list, boost::is_any_of(","));
85  for (vstring::const_iterator i= v1.begin(); i!=v1.end(); ++i){
86  vstring v2;
87  boost::split(v2, *i, boost::is_any_of("-"));
88  unsigned start = boost::lexical_cast<unsigned>(v2.front());
89  unsigned end = boost::lexical_cast<unsigned>(v2.back());
90  if ((v2.size()>2) || (start>end)) {
91  throw cms::Exception("ConfigurationError") << "VolumeBasedMagneticFieldESProducerFromDB: malformed configuration" << list << endl;
92  }
93  for (unsigned k = start; k<=end; ++k){
94  values.push_back(k);
95  }
96  }
97  return values;
98 }
T getParameter(std::string const &) const
int i
Definition: DBlmapReader.cc:9
tuple start
Check for commandline option errors.
Definition: dqm_diff.py:58
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:186
vector< string > vstring
Definition: ExoticaDQM.cc:8
int j
Definition: DBlmapReader.cc:9
#define end
Definition: vmac.h:37
#define debug
Definition: HDRShower.cc:19
tuple cout
Definition: gather_cfg.py:145
double split
Definition: MVATrainer.cc:139
std::vector< unsigned > expandList(const std::string &list)
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger list("!*","!HLTx*"if it matches 2 triggers or more) will accept the event if all the matching triggers are FAIL.It will reject the event if any of the triggers are PASS or EXCEPTION(this matches the behavior of"!*"before the partial wildcard feature was incorporated).Triggers which are in the READY state are completely ignored.(READY should never be returned since the trigger paths have been run