ATLAS Offline Software
Loading...
Searching...
No Matches
LArReadParamsFromFile.icc
Go to the documentation of this file.
1//Dear emacs, this is -*- c++ -*-
2
3/*
4 Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
5*/
6
7
8#include "LArRawConditions/LArConditionsContainerBase.h"
9
10#include <iostream>
11#include <sstream>
12#include <string>
13
14template <class DATA>
15LArReadParamsFromFile<DATA>::~LArReadParamsFromFile() = default;
16
17template <class DATA>
18StatusCode LArReadParamsFromFile<DATA>::initialize() {
19 ATH_MSG_DEBUG ( "in initialize" );
20
21 if ( m_groupingName == "Unknown" ) {
22 m_groupingType = LArConditionsContainerBase::Unknown ;
23 } else if ( m_groupingName == "SingleGroup" ) {
24 m_groupingType = LArConditionsContainerBase::SingleGroup ;
25 } else if ( m_groupingName == "SubDetectorGrouping" ) {
26 m_groupingType = LArConditionsContainerBase::SubDetectorGrouping ;
27 } else if ( m_groupingName == "FeedThroughGrouping" ) {
28 m_groupingType = LArConditionsContainerBase::FeedThroughGrouping ;
29 } else if ( m_groupingName == "ExtendedFThGrouping" ) {
30 m_groupingType = LArConditionsContainerBase::ExtendedFTGrouping ;
31 } else if ( m_groupingName == "ExtendedSubDetGrouping" ) {
32 m_groupingType = LArConditionsContainerBase::ExtendedSubDetGrouping ;
33 } else {
34 ATH_MSG_ERROR ( "Grouping type " << m_groupingName << " is not foreseen!" );
35 ATH_MSG_ERROR ( "Only \"Unknown\", \"SingleGroup\", \"SubDetectorGrouping\", \"FeedThroughGrouping\" are allowed" );
36 return StatusCode::FAILURE ;
37 }
38
39 const CaloCell_ID* idHelper = nullptr;
40 ATH_CHECK( detStore()->retrieve (idHelper, "CaloCell_ID") );
41 m_emId = idHelper->em_idHelper();
42 m_hecId = idHelper->hec_idHelper();
43 m_fcalId = idHelper->fcal_idHelper();
44 if (!m_emId) {
45 ATH_MSG_ERROR ( "Could not access lar EM ID helper" );
46 return StatusCode::FAILURE;
47 }
48 if (!m_hecId) {
49 ATH_MSG_ERROR ( "Could not access lar HEC ID helper" );
50 return StatusCode::FAILURE;
51 }
52 if (!m_fcalId) {
53 ATH_MSG_ERROR ( "Could not access lar FCAL ID helper" );
54 return StatusCode::FAILURE;
55 }
56
57 ATH_CHECK( detStore()->retrieve(m_onlineHelper, "LArOnlineID") );
58
59 ATH_CHECK( m_cablingKey.initialize() );
60
61 if ( m_chIdType == std::string("UNKNOWN") ) {
62 // use jO specification (default=false) for m_useOfflineIdentifier
63 // and default=false for m_useCalibLines
64 m_useCalibLines = false ;
65 } else if ( m_chIdType == std::string("OfflineID") ) {
66 m_useOfflineIdentifier = true ;
67 m_useCalibLines = false ;
68 } else if ( m_chIdType == std::string("OnlineID") ) {
69 m_useOfflineIdentifier = false ;
70 m_useCalibLines = false ;
71 } else if ( m_chIdType == std::string("OnlineCalibID") ) {
72 m_useOfflineIdentifier = false ;
73 m_useCalibLines = true ;
74 } else {
75 ATH_MSG_ERROR ( "ChannelIdType=" << m_chIdType << " is not valid!" );
76 ATH_MSG_ERROR ( "Allowed options are: OfflineID, OnlineID, OnlineCalibID" );
77 return StatusCode::FAILURE ;
78 }
79
80 if( m_useOfflineIdentifier ) {
81 ATH_MSG_INFO ( "Reading parameters using *** offline identifiers ***" );
82 } else if ( m_useCalibLines ) {
83 ATH_MSG_INFO ( "Reading parameters for *** calibration lines ***" );
84 } else {
85 ATH_MSG_INFO ( "Reading parameters using *** online identifiers ***" );
86 }
87
88 return StatusCode::SUCCESS;
89}
90
91
92template <class DATA>
93StatusCode LArReadParamsFromFile<DATA>::stop ATLAS_NOT_THREAD_SAFE ()
94{
95 ATH_MSG_DEBUG ( "in stop" );
96
97 if ( m_file == std::string("") ) {
98 ATH_MSG_ERROR ( "Input file name was not specified!" );
99 return StatusCode::FAILURE ;
100 }
101
102 m_dataclass = new DATA ;
103 if ( m_groupingType != LArConditionsContainerBase::Unknown ) {
104 m_dataclass->setGroupingType( static_cast<LArConditionsContainerBase::GroupingType>(m_groupingType) ) ;
105 }
106 ATH_CHECK( m_dataclass->initialize() );
107
108 std::string classname = LArParamsProperties::getClassName(m_dataclass) ;
109 std::string key = LArParamsProperties::keyword( classname ) ;
110 if ( m_customKey != "") key = m_customKey;
111
112 ATH_CHECK( readFile() );
113 ATH_MSG_INFO ( "Object " << classname << " successfully uploaded from file " << m_file );
114
115 ATH_CHECK( detStore()->record(m_dataclass,key) );
116 ATH_MSG_INFO ( "Object " << classname << " recorded into detStore with key " << key );
117
118 ATH_CHECK ( do_symLink(m_dataclass) );
119 ATH_MSG_INFO ( "Object " << classname << " symLinked to its abstract interface" );
120
121 ATH_MSG_INFO ( "finished!" );
122 return StatusCode::SUCCESS;
123}// end stop-method.
124
125
126template <class DATA>
127StatusCode LArReadParamsFromFile<DATA>::readFile() {
128 // File format:
129 // for online ID: barrel_ec posneg FT slot channel parameters...
130 // for offline ID: det subdet barrel_ec layer region eta phi parameters...
131 int det, subdet, barrel_ec, posneg, ft, slot, channel, region, layer, eta, phi ;
132 std::vector<float> pRead ;
133 HWIdentifier chID ;
134
135 SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey};
136 const LArOnOffIdMapping* cabling{*cablingHdl};
137 if(!cabling) {
138 ATH_MSG_ERROR( "Do not have cabling from key " << m_cablingKey.key() );
139 return StatusCode::FAILURE;
140 }
141
142 std::ifstream f(m_file) ; // open input file
143 if ( ! f.good() ) {
144 ATH_MSG_ERROR ( "Could not open file " << m_file << ", are you sure it exists??" );
145 return StatusCode::FAILURE ;
146 }
147 ATH_MSG_DEBUG ( "File " << m_file << " opened" );
148
149 unsigned nrec = 0 ;
150 unsigned npar = LArParamsProperties::numberOfParams(m_dataclass) ;
151 ATH_MSG_DEBUG ( "The file should contain " << npar << " parameters per record" );
152
153 while (1) { // read through the file records
154
155 if ( f.eof() ) break ;
156 std::string s ;
157 getline(f,s) ;
158 std::istringstream buffer(s) ;
159 ATH_MSG_VERBOSE ( "*** Record: '" << s << "'" );
160
161 if ( s == std::string("") ) continue ; // skip empty records
162
163 bool good_chID = true ;
164
165 if ( m_useOfflineIdentifier ) { // get HWIdentifier from offline "coordinates"
166
167 buffer >> det ;
168 if ( buffer.fail() ) {
169 ATH_MSG_WARNING ( "Could not read 'det' field for offline identifier" );
170 continue ; // skip and go to next record
171 }
172 buffer >> subdet ;
173 if ( buffer.fail() ) {
174 ATH_MSG_WARNING ( "Could not read 'subdet' field for offline identifier" );
175 continue ; // skip and go to next record
176 }
177 buffer >> barrel_ec ;
178 if ( buffer.fail() ) {
179 ATH_MSG_WARNING ( "Could not read 'barrel_ec' field for offline identifier" );
180 continue ; // skip and go to next record
181 }
182 buffer >> layer ;
183 if ( buffer.fail() ) {
184 ATH_MSG_WARNING ( "Could not read 'layer' field for offline identifier" );
185 continue ; // skip and go to next record
186 }
187 buffer >> region ;
188 if ( buffer.fail() ) {
189 ATH_MSG_WARNING ( "Could not read 'region' field for offline identifier" );
190 continue ; // skip and go to next record
191 }
192 buffer >> eta ;
193 if ( buffer.fail() ) {
194 ATH_MSG_WARNING ( "Could not read 'eta' field for offline identifier" );
195 continue ; // skip and go to next record
196 }
197 buffer >> phi ;
198 if ( buffer.fail() ) {
199 ATH_MSG_WARNING ( "Could not read 'phi' field for offline identifier" );
200 continue ; // skip and go to next record
201 }
202 ATH_MSG_VERBOSE ( "Det=" << det << " SubDet=" << subdet << " BarrelEC=" << barrel_ec
203 << " Sampling=" << layer << " Region=" << region << " Eta=" << eta << " Phi=" << phi );
204 good_chID = true ;
205 if ( det != 4 ) {
206 ATH_MSG_WARNING ( "Det=" << det << " is not LArCalorimeter -- cannot handle it!" );
207 good_chID = false ;
208 } else {
209 Identifier id ;
210 if ( subdet == 1 ) { // LArEM
211 id = m_emId->channel_id(barrel_ec,layer,region,eta,phi) ;
212 ATH_MSG_VERBOSE ( "LAr EM Identifier=" << id );
213 } else if ( subdet == 2 ) { // LArHEC
214 id = m_hecId->channel_id(barrel_ec,layer,region,eta,phi) ;
215 ATH_MSG_VERBOSE ( "LAr HEC Identifier=" << id );
216 } else if ( subdet == 3 ) { // LArFCAL
217 // id = m_fcalId->channel_id(det,layer,region,eta,phi) ;
218 // log << MSG::VERBOSE << "LAr FCAL Identifier=" << id << endmsg ;
219 ATH_MSG_WARNING ( "Subdet=" << subdet
220 << " is LArFCAL -- conversion to online channel not implemented yet!" );
221 good_chID = false ;
222 } else {
223 ATH_MSG_WARNING ( "Subdet=" << subdet << " is unknown -- cannot handle it!" );
224 good_chID = false ;
225 }
226 if ( good_chID ) {
227 try {
228 chID = cabling->createSignalChannelID(id) ;
229 } catch ( LArID_Exception & except ) {
230 ATH_MSG_WARNING ( "Could not get HWId for offline Id " << id );
231 good_chID = false ;
232 }
233 }
234 }
235 if ( ! good_chID ) continue ; // skip and go to next record
236
237 } else { // get HWIdentifier from online "coordinates"
238
239 buffer >> barrel_ec ;
240 if ( buffer.fail() ) {
241 ATH_MSG_WARNING ( "Could not read 'barrel_ec' field for online identifier" );
242 continue ; // skip and go to next record
243 }
244 buffer >> posneg ;
245 if ( buffer.fail() ) {
246 ATH_MSG_WARNING ( "Could not read 'posneg' field for online identifier" );
247 continue ; // skip and go to next record
248 }
249 buffer >> ft ;
250 if ( buffer.fail() ) {
251 ATH_MSG_WARNING ( "Could not read 'feedthrough' field for online identifier" );
252 continue ; // skip and go to next record
253 }
254 buffer >> slot ;
255 if ( buffer.fail() ) {
256 ATH_MSG_WARNING ( "Could not read 'slot' field for online identifier" );
257 continue ; // skip and go to next record
258 }
259 buffer >> channel ;
260 if ( buffer.fail() ) {
261 ATH_MSG_WARNING ( "Could not read 'channel' field for online identifier" );
262 continue ; // skip and go to next record
263 }
264 if ( m_useCalibLines ) {
265 chID = m_onlineHelper->calib_channel_Id(barrel_ec,posneg,ft,slot,channel) ;
266 } else {
267 chID = m_onlineHelper->channel_Id(barrel_ec,posneg,ft,slot,channel) ;
268 }
269
270 } // here the HWIdentifier is obtained!
271
272 if ( ! good_chID ) continue ; // skip and go to next record
273
274 //
275 // get parameters
276 //
277 pRead.resize(npar) ;
278 const float defaultValue = -999 ;
279 for ( unsigned i=0 ; i<npar ; i++ ) pRead[i] = defaultValue ;
280 for ( unsigned i=0 ; i<npar ; i++ ) {
281 //if ( buffer.eof() ) {
282 //log << MSG::DEBUG << "end-of-record reached when reading parameter " << i << " of channel " << chID.get_compact() << endmsg ;
283 //break ;
284 //}
285 buffer >> pRead[i] ;
286 if ( buffer.fail() ) {
287 ATH_MSG_WARNING ( "failure when reading parameter " << i << " of channel " << chID.get_compact() );
288 break ;
289 }
290 }
291 // store into parameters vector
292 msg() << MSG::VERBOSE << "read ch=" << chID.get_compact() ;
293 for ( unsigned i=0 ; i<npar ; i++ ) msg() << " par[" << i << "]=" << pRead[i] ;
294 msg() << endmsg ;
295 for ( unsigned gain=0 ; gain<3 ; gain++ ) {
296 try {
297 StatusCode sc = set(m_dataclass,chID,gain,pRead) ;
298 if ( sc.isFailure() ) {
299 ATH_MSG_WARNING ( "Could not set parameters for ch=" << chID.get_compact() << " gain=" << gain );
300 }
301 } catch ( LArID_Exception & except ) {
302 ATH_MSG_WARNING ( "Bad HWId " << chID );
303 good_chID = false ;
304 }
305 }
306
307 nrec ++ ;
308
309 } // end of while(1) i.e. the loop through all file records
310
311 ATH_MSG_DEBUG ( "end of file reached, read " << nrec << " records" );
312 return StatusCode::SUCCESS ;
313}
Generated on for ATLAS Offline Software by doxygen 1.16.1