ATLAS Offline Software
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Private Attributes | List of all members
MuonCalib::TubeDataFromFile Class Reference

Manage the I/O of MdtTubeFitContainer objects. More...

#include <TubeDataFromFile.h>

Collaboration diagram for MuonCalib::TubeDataFromFile:

Public Types

using TubeData = std::vector<const MdtTubeFitContainer*>

Public Member Functions

 TubeDataFromFile ()=default
 ~TubeDataFromFile ()=default
unsigned int nRegions () const
 return number of regions
const MdtTubeFitContainer * getTubes (unsigned int regionId) const
 retrieve MdtTubeFitContainer for a give regionId
void setNRegions (unsigned int n)
 set total number of regions
bool addTubes (int regionId, const MdtTubeFitContainer *tubes)
 TubeDataFromFile takes ownership of the MdtTubeFitContainer.
std::ostream & write (std::ostream &os) const
std::ostream & write_forDB (std::ostream &ftube, int mdt_head, int lowrun, int uprun) const

Private Attributes

unsigned int m_regions {0}
 total number of regions
TubeData m_tubeData
 data

Detailed Description

Manage the I/O of MdtTubeFitContainer objects.

Author
domiz.nosp@m.ia.o.nosp@m.resta.nosp@m.no@c.nosp@m.ern.c.nosp@m.h
Date
july 2005

Definition at line 24 of file TubeDataFromFile.h.

Member Typedef Documentation

◆ TubeData

using MuonCalib::TubeDataFromFile::TubeData = std::vector<const MdtTubeFitContainer*>

Definition at line 26 of file TubeDataFromFile.h.

Constructor & Destructor Documentation

◆ TubeDataFromFile()

MuonCalib::TubeDataFromFile::TubeDataFromFile ( )
default

◆ ~TubeDataFromFile()

MuonCalib::TubeDataFromFile::~TubeDataFromFile ( )
default

Member Function Documentation

◆ addTubes()

bool MuonCalib::TubeDataFromFile::addTubes ( int regionId,
const MdtTubeFitContainer * tubes )
inline

TubeDataFromFile takes ownership of the MdtTubeFitContainer.

Definition at line 53 of file TubeDataFromFile.h.

53 {
54 if (regionId < 0 || regionId >= (int)m_regions) {
55 MsgStream log(Athena::getMessageSvc(), "MdtTubeFitContainer");
56 log << MSG::WARNING << "TubeDataFromFile::addTubes: <regionId out of range> " << regionId << " size " << m_regions
57 << endmsg;
58 return false;
59 }
60 if (m_tubeData[regionId] != 0) {
61 MsgStream log(Athena::getMessageSvc(), "MdtTubeFitContainer");
62 log << MSG::WARNING << "TubeDataFromFile::addTubes: <tubes already set>" << endmsg;
63 return false;
64 }
65
66 m_tubeData[regionId] = tubes;
67
68 return true;
69 }
endmsg
#define endmsg
Definition AnalysisConfig_Ntuple.cxx:63
MuonCalib::TubeDataFromFile::m_regions
unsigned int m_regions
total number of regions
Definition TubeDataFromFile.h:77
Athena::getMessageSvc
IMessageSvc * getMessageSvc(bool quiet=false)
Definition getMessageSvc.cxx:20

◆ getTubes()

const MdtTubeFitContainer * MuonCalib::TubeDataFromFile::getTubes ( unsigned int regionId) const
inline

retrieve MdtTubeFitContainer for a give regionId

Definition at line 36 of file TubeDataFromFile.h.

36 {
37 if (regionId >= (unsigned int)m_regions) {
38 MsgStream log(Athena::getMessageSvc(), "MdtTubeFitContainer");
39 log << MSG::WARNING << "TubeDataFromFile::getTubes: <regionId out of range> " << regionId << " size " << m_regions
40 << endmsg;
41 return nullptr;
42 }
43 return m_tubeData[regionId];
44 }

◆ nRegions()

unsigned int MuonCalib::TubeDataFromFile::nRegions ( ) const
inline

return number of regions

Definition at line 33 of file TubeDataFromFile.h.

33{ return m_regions; }

◆ setNRegions()

void MuonCalib::TubeDataFromFile::setNRegions ( unsigned int n)
inline

set total number of regions

Definition at line 47 of file TubeDataFromFile.h.

47 {
48 m_regions = n;
49 m_tubeData.resize(n);
50 }

◆ write()

std::ostream & MuonCalib::TubeDataFromFile::write ( std::ostream & os) const

Definition at line 17 of file TubeDataFromFile.cxx.

17 {
18 if (m_regions != m_tubeData.size()) {
19 MsgStream log(Athena::getMessageSvc(), "MdtTubeFitContainer");
20 log << MSG::WARNING << "TubeDataFromFile::write: <inconsistent count>" << endmsg;
21 }
22
23 os << "v0.0 " << m_regions << " ";
24 unsigned int counter{0};
25 for (const MdtTubeFitContainer* container : m_tubeData) {
26 if (!container) continue;
27 const unsigned int nml = container->numMultilayers();
28 const unsigned int nl = container->numLayers();
29 const unsigned int nt = container->numTubes();
30 const std::string datatype = container->name();
31 const std::string implementation = container->implementation();
32 const int stnEta = container->idHelperSvc()->stationEta(container->identify());
33 const int stnPhi = container->idHelperSvc()->stationPhi(container->identify());
34 const std::string stnName = container->idHelperSvc()->stationNameString(container->identify());
35 os << datatype << " " << implementation << " " << counter << " " << container->size() << std::endl;
36 ++counter;
37 const MdtIdHelper& id_helper{container->idHelperSvc()->mdtIdHelper()};
38 for (unsigned int km = 1; km <= nml; ++km) {
39 for (unsigned int kl = 1; kl <= nl; ++kl) {
40 for (unsigned int kt = 1; kt <= nt; ++kt) {
41 const Identifier chId{id_helper.channelID(container->identify(),km ,kl ,kt)};
42 const MdtTubeFitContainer::SingleTubeCalib* stc = container->getCalib(chId);
43 const MdtTubeFitContainer::SingleTubeFit* stf = container->getFit(chId);
44 double t0 = 999.;
45 if (stc) {
46 t0 = stc->t0;
47 if (std::isnan(t0) != 0) t0 = -99999.;
48 }
49 constexpr int technology = 0;
50 os << " " << stnName << " " << stnPhi << " " << stnEta
51 << " " << technology << " " << km << " " << kl << " " << kt;
52 if (stc) {
53 os << " " << t0 << " " << stc->adcCal << " " << stc->statusCode;
54 } else{
55 os << " " << 0 << " " << 0 << " " << 0;
56
57 }
58 if (stf)
59 os << " " << stf->statistics << " " << stf->chi2Tdc << " " << stf->cov[4] // stf->t0Err
60 << " " << stf->par[5] // stf->tMax
61 << " " << stf->cov[5] // stf->tMaxErr
62 << " " << stf->par[0] // stf->noise
63 << " " << stf->cov[0] // stf->noiseErr
64 << " " << stf->par[6] // stf->slope
65 << " " << stf->cov[6]; // stf->slopeErr;
66 else {
67 os << " " << -1 << " " << 0 << " " << 0 << " " << 0 << " "
68 << 0 << " " << 0 << " " << 0 << " " << 0 << " "
69 << 0;
70 }
71 os << std::endl;
72 }
73 }
74 }
75 }
76 return os;
77 }
t0
static Double_t t0
Definition LArPhysWaveHECTool.cxx:38
MdtIdHelper::channelID
Identifier channelID(int stationName, int stationEta, int stationPhi, int multilayer, int tubeLayer, int tube) const
Definition MdtIdHelper.cxx:664
keylayer_zslicemap.kl
dict kl
Definition keylayer_zslicemap.py:109
test_pyathena.counter
counter
Definition test_pyathena.py:15

◆ write_forDB()

std::ostream & MuonCalib::TubeDataFromFile::write_forDB ( std::ostream & ftube,
int mdt_head,
int lowrun,
int uprun ) const

Definition at line 80 of file TubeDataFromFile.cxx.

80 {
82 // write out ascii files for Calibration DB //
84 if (m_regions != m_tubeData.size()) {
85 MsgStream log(Athena::getMessageSvc(), "MdtTubeFitContainer");
86 log << MSG::WARNING << "TubeDataFromFile::write_forDB: <inconsistent count>" << endmsg;
87 }
88
89 for (const MdtTubeFitContainer* container : m_tubeData) {
90 if (!container) continue;
91 const unsigned int nml = container->numMultilayers();
92 const unsigned int nl = container->numLayers();
93 const unsigned int nt = container->numTubes();
94 const std::string datatype = container->name();
95 const std::string implementation = container->implementation();
96 const int stnEta = container->idHelperSvc()->stationEta(container->identify());
97 const int stnPhi = container->idHelperSvc()->stationPhi(container->identify());
98 const std::string stnName = container->idHelperSvc()->stationNameString(container->identify());
99
100 const MdtIdHelper& id_helper{container->idHelperSvc()->mdtIdHelper()};
101 for (unsigned int km = 1; km <= nml; ++km) {
102 for (unsigned int kl = 1; kl <= nl; ++kl) {
103 for (unsigned int kt = 1; kt <= nt; ++kt) {
104 const Identifier chId{id_helper.channelID(container->identify(),km ,kl ,kt)};
105 const MdtTubeFitContainer::SingleTubeCalib* stc = container->getCalib(chId);
106 const MdtTubeFitContainer::SingleTubeFit* stf = container->getFit(chId);
107 constexpr int technology = 0;
108
109 MuonFixedId fixId{};
110 int sc{0}, stat{-1};
111 double t0{0.}, adcm{0.}, chi2tdc{0.}, t0err{0.},
112 tmax{0.}, tmaxerr{0.}, noise{0.}, noiseerr{0.},
113 slope{0.}, sloperr{0.};
114 int tube_id{0}, runflag{0}, validflag{0};
115
116 if (!fixId.setTechnology(technology) ||
117 !fixId.setStationName(fixId.stationStringToFixedStationNumber(stnName)) ||
118 !fixId.setStationEta(stnEta) ||
119 !fixId.setStationPhi(stnPhi) ||
120 !fixId.setMdtTube(kt) ||
121 !fixId.setMdtTubeLayer(kl) ||
122 !fixId.setMdtMultilayer(km)) {
123 std::stringstream except{};
124 except<<__FILE__<<":"<<__LINE__<<"nTubeDataFromFile::write_forDB() - Setting identifier failed!";
125 throw std::runtime_error(except.str());
126 }
127 tube_id = fixId.getIdInt();
128
129 if (stc) {
130 t0 = stc->t0;
131 adcm = stc->adcCal;
132 sc = stc->statusCode;
133 }
134 if (stf) {
135 stat = stf->statistics;
136 chi2tdc = stf->chi2Tdc;
137 t0err = stf->cov[4];
138 tmax = stf->par[5];
139 tmaxerr = stf->cov[5];
140 noise = stf->par[0];
141 noiseerr = stf->cov[0];
142 slope = stf->par[6];
143 sloperr = stf->cov[6];
144 }
145
146 if (std::isnan(t0) != 0) t0 = -99999.;
147
148 ftube<<mdt_head<<","<<tube_id<<","<<fixId.mdtChamberId().getIdInt()<<","
149 <<","<<lowrun<<","<<uprun<<","<<runflag<<","<<sc
150 <<validflag<<","<<stat<<",";
151
152 for (int ii = -2; ii < 35; ii++) {
153 if (ii == -2) {
154 ftube<<chi2tdc<<",";
155 } else if (ii == 0) {
156 ftube<<noise<<",";
157 } else if (ii == 4) {
158 ftube<<t0<<",";
159 } else if (ii == 5) {
160 ftube<<tmax<<",";
161 } else if (ii == 6) {
162 ftube<<slope<<",";
163 } else if (ii == 10) {
164 ftube<<t0err<<",";
165 } else if (ii == 11) {
166 ftube<<tmaxerr<<",";
167 } else if (ii == 12) {
168 ftube<<noiseerr<<",";
169 } else if (ii == 13) {
170 ftube<<sloperr<<",";
171 } else if (ii == 14) {
172 ftube<<adcm<<",";
173 } else {
174 ftube<<0.0<<",";
175 }
176 }
177 ftube<<"0,0,algoflag,tubegrpgr"<<std::endl;
178 }
179 }
180 }
181 }
182 return ftube;
183 }
sc
static Double_t sc
Definition LArPhysWaveHECTool.cxx:37

Member Data Documentation

◆ m_regions

unsigned int MuonCalib::TubeDataFromFile::m_regions {0}
private

total number of regions

Definition at line 77 of file TubeDataFromFile.h.

77{0};

◆ m_tubeData

TubeData MuonCalib::TubeDataFromFile::m_tubeData
private

data

Definition at line 80 of file TubeDataFromFile.h.

The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.14.0