140 std::set<RpcChamber> to_ret{};
141 std::cout<<
"Read the Rpc geometry tree dump from "<<
inputFile<<std::endl;
144 std::cerr<<__FILE__<<
":"<<__LINE__<<
" Failed to open "<<
inputFile<<std::endl;
149 std::cerr<<__FILE__<<
":"<<__LINE__<<
" The file "<<
inputFile<<
" does not contain the 'RpcGeoModelTree'"<<std::endl;
153 TTreeReaderValue<unsigned short> stationIndex{
treeReader,
"stationIndex"};
156 TTreeReaderValue<uint8_t> stationDoubletR{
treeReader,
"stationDoubletR"};
157 TTreeReaderValue<uint8_t> stationDoubletZ{
treeReader,
"stationDoubletZ"};
158 TTreeReaderValue<uint8_t> stationDoubletPhi{
treeReader,
"stationDoubletPhi"};
159 TTreeReaderValue<std::string> chamberDesign{
treeReader,
"chamberDesign"};
163 TTreeReaderValue<uint8_t> numStripsEta{
treeReader,
"numEtaStrips"};
164 TTreeReaderValue<uint8_t> numStripsPhi{
treeReader,
"numPhiStrips"};
167 TTreeReaderValue<uint8_t> numPhiPanels{
treeReader,
"numPhiPanels"};
171 TTreeReaderValue<float> stripEtaPitch{
treeReader,
"stripEtaPitch"};
172 TTreeReaderValue<float> stripPhiPitch{
treeReader,
"stripPhiPitch"};
173 TTreeReaderValue<float> stripEtaWidth{
treeReader,
"stripEtaWidth"};
174 TTreeReaderValue<float> stripPhiWidth{
treeReader,
"stripPhiWidth"};
175 TTreeReaderValue<float> stripEtaLength{
treeReader,
"stripEtaLength"};
176 TTreeReaderValue<float> stripPhiLength{
treeReader,
"stripPhiLength"};
179 TTreeReaderValue<std::vector<float>> geoModelTransformX{
treeReader,
"GeoModelTransformX"};
180 TTreeReaderValue<std::vector<float>> geoModelTransformY{
treeReader,
"GeoModelTransformY"};
181 TTreeReaderValue<std::vector<float>> geoModelTransformZ{
treeReader,
"GeoModelTransformZ"};
184 TTreeReaderValue<std::vector<float>> stripRotTranslationX{
treeReader,
"stripRotTranslationX"};
185 TTreeReaderValue<std::vector<float>> stripRotTranslationY{
treeReader,
"stripRotTranslationY"};
186 TTreeReaderValue<std::vector<float>> stripRotTranslationZ{
treeReader,
"stripRotTranslationZ"};
189 TTreeReaderValue<std::vector<float>> stripRotCol1X{
treeReader,
"stripRotLinearCol1X"};
190 TTreeReaderValue<std::vector<float>> stripRotCol1Y{
treeReader,
"stripRotLinearCol1Y"};
191 TTreeReaderValue<std::vector<float>> stripRotCol1Z{
treeReader,
"stripRotLinearCol1Z"};
193 TTreeReaderValue<std::vector<float>> stripRotCol2X{
treeReader,
"stripRotLinearCol2X"};
194 TTreeReaderValue<std::vector<float>> stripRotCol2Y{
treeReader,
"stripRotLinearCol2Y"};
195 TTreeReaderValue<std::vector<float>> stripRotCol2Z{
treeReader,
"stripRotLinearCol2Z"};
197 TTreeReaderValue<std::vector<float>> stripRotCol3X{
treeReader,
"stripRotLinearCol3X"};
198 TTreeReaderValue<std::vector<float>> stripRotCol3Y{
treeReader,
"stripRotLinearCol3Y"};
199 TTreeReaderValue<std::vector<float>> stripRotCol3Z{
treeReader,
"stripRotLinearCol3Z"};
201 TTreeReaderValue<std::vector<uint8_t>> stripRotGasGap{
treeReader,
"stripRotGasGap"};
202 TTreeReaderValue<std::vector<uint8_t>> stripRotDblPhi{
treeReader,
"stripRotDoubletPhi"};
203 TTreeReaderValue<std::vector<bool>> stripRotMeasPhi{
treeReader,
"stripRotMeasPhi"};
205 TTreeReaderValue<std::vector<float>> stripPosX{
treeReader,
"stripPosX"};
206 TTreeReaderValue<std::vector<float>> stripPosY{
treeReader,
"stripPosY"};
207 TTreeReaderValue<std::vector<float>> stripPosZ{
treeReader,
"stripPosZ"};
209 TTreeReaderValue<std::vector<float>> stripLocPosX{
treeReader,
"stripLocPosX"};
210 TTreeReaderValue<std::vector<float>> stripLocPosY{
treeReader,
"stripLocPosY"};
213 TTreeReaderValue<std::vector<bool>> stripPosMeasPhi{
treeReader,
"stripPosMeasPhi"};
214 TTreeReaderValue<std::vector<uint8_t>> stripPosGasGap{
treeReader,
"stripPosGasGap"};
215 TTreeReaderValue<std::vector<uint8_t>> stripPosNum{
treeReader,
"stripPosNum"};
216 TTreeReaderValue<std::vector<uint8_t>> stripDblPhi{
treeReader,
"stripPosDoubletPhi"};
222 newchamber.design = (*chamberDesign);
223 newchamber.id.eta = (*stationEta);
224 newchamber.id.phi = (*stationPhi);
225 newchamber.id.doubletPhi = (*stationDoubletPhi);
226 newchamber.id.doubletR = (*stationDoubletR);
227 newchamber.id.doubletZ = (*stationDoubletZ);
229 newchamber.stripPitchEta = (*stripEtaPitch);
230 newchamber.stripPitchPhi = (*stripPhiPitch);
231 newchamber.stripWidthEta = (*stripEtaWidth);
232 newchamber.stripWidthPhi = (*stripPhiWidth);
233 newchamber.stripLengthEta = (*stripEtaLength);
234 newchamber.stripLengthPhi = (*stripPhiLength);
236 newchamber.numStripsEta = (*numStripsEta);
237 newchamber.numStripsPhi = (*numStripsPhi);
239 newchamber.numPhiPanels = (*numPhiPanels);
240 newchamber.numLayers = (*numLayers);
243 Amg::Vector3D geoTrans{(*geoModelTransformX)[0], (*geoModelTransformY)[0], (*geoModelTransformZ)[0]};
245 geoRot.col(0) =
Amg::Vector3D((*geoModelTransformX)[1], (*geoModelTransformY)[1], (*geoModelTransformZ)[1]);
246 geoRot.col(1) =
Amg::Vector3D((*geoModelTransformX)[2], (*geoModelTransformY)[2], (*geoModelTransformZ)[2]);
247 geoRot.col(2) =
Amg::Vector3D((*geoModelTransformX)[3], (*geoModelTransformY)[3], (*geoModelTransformZ)[3]);
251 for (
size_t s = 0;
s < stripPosX->size(); ++
s){
254 newStrip.locPos =
Amg::Vector2D{(*stripLocPosX)[
s], (*stripLocPosY)[
s]};
256 newStrip.gasGap = (*stripPosGasGap)[
s];
257 newStrip.doubletPhi = (*stripDblPhi)[
s];
258 newStrip.measPhi = (*stripPosMeasPhi)[
s];
259 newStrip.strip = (*stripPosNum)[
s];
260 newchamber.strips.insert(std::move(newStrip));
262 for (
size_t l = 0;
l < stripRotMeasPhi->size(); ++
l){
264 newLayer.
measPhi = (*stripRotMeasPhi)[
l];
265 newLayer.gasGap = (*stripRotGasGap)[
l];
266 newLayer.doubletPhi = (*stripRotDblPhi)[
l];
268 Amg::Vector3D translation{(*stripRotTranslationX)[
l],(*stripRotTranslationY)[
l],(*stripRotTranslationZ)[
l]};
269 stripRot.col(0) =
Amg::Vector3D((*stripRotCol1X)[
l],(*stripRotCol1Y)[
l], (*stripRotCol1Z)[
l]);
270 stripRot.col(1) =
Amg::Vector3D((*stripRotCol2X)[
l],(*stripRotCol2Y)[
l], (*stripRotCol2Z)[
l]);
271 stripRot.col(2) =
Amg::Vector3D((*stripRotCol3X)[
l],(*stripRotCol3Y)[
l], (*stripRotCol3Z)[
l]);
273 newchamber.layers.insert(std::move(newLayer));
276 auto insert_itr = to_ret.insert(std::move(newchamber));
277 if (!insert_itr.second) {
278 std::stringstream
err{};
279 err<<__FILE__<<
":"<<__LINE__<<
" The chamber "<<(*insert_itr.first).
id
280 <<
" has already been inserted. "<<std::endl;
281 throw std::runtime_error(
err.str());
284 std::cout<<
"File parsing is finished. Found in total "<<to_ret.size()<<
" readout element dumps "<<std::endl;