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

Calculates OFCs on the fly using pulse shapes and A/C matrix from database. More...

#include <TileCondToolOfc.h>

Inheritance diagram for TileCondToolOfc:
Collaboration diagram for TileCondToolOfc:

Public Member Functions

 TileCondToolOfc (const std::string &type, const std::string &name, const IInterface *parent)
virtual ~TileCondToolOfc ()
virtual StatusCode initialize () override
virtual StatusCode finalize () override
virtual StatusCode getOfcWeights (unsigned int drawerIdx, unsigned int channel, unsigned int adc, float &phase, bool of2, TileOfcWeightsStruct &weights, const EventContext &ctx) const override

Private Attributes

ToolHandle< TileCondToolPulseShape > m_tileToolPulseShape
ToolHandle< TileCondToolAutoCr > m_tileToolAutoCr
const TileInfo * m_tileInfo
int m_nSamples
int m_t0Sample
bool m_deltaCorrelation

Detailed Description

Calculates OFCs on the fly using pulse shapes and A/C matrix from database.

Optionally, it can create cache table of OFCs with 1-ns step to minimize CPU time. Also, by request, unity A/C matrix can be used.

Definition at line 40 of file TileCondToolOfc.h.

Constructor & Destructor Documentation

◆ TileCondToolOfc()

TileCondToolOfc::TileCondToolOfc ( const std::string & type,
const std::string & name,
const IInterface * parent )

Definition at line 16 of file TileCondToolOfc.cxx.

18 : base_class(type, name, parent)
19 , m_tileInfo(nullptr)
20{
21 declareProperty("nSamples", m_nSamples = 7, "number of samples used in the run");
22 declareProperty("OptFilterDeltaCorrelation", m_deltaCorrelation = false
23 , "true=> use delta correlation; false=>use calculation obtained from data");
24
25 m_t0Sample = (m_nSamples - 1) / 2;
26}
TileCondToolOfc::m_tileInfo
const TileInfo * m_tileInfo
Definition TileCondToolOfc.h:69

◆ ~TileCondToolOfc()

TileCondToolOfc::~TileCondToolOfc ( )
virtual

Definition at line 30 of file TileCondToolOfc.cxx.

30 {
31
32}

Member Function Documentation

◆ finalize()

StatusCode TileCondToolOfc::finalize ( )
overridevirtual

Definition at line 90 of file TileCondToolOfc.cxx.

90 {
91 ATH_MSG_DEBUG( "TileCondToolOfc finalize called" );
92
93 return StatusCode::SUCCESS;
94}
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29

◆ getOfcWeights()

StatusCode TileCondToolOfc::getOfcWeights ( unsigned int drawerIdx,
unsigned int channel,
unsigned int adc,
float & phase,
bool of2,
TileOfcWeightsStruct & weights,
const EventContext & ctx ) const
overridevirtual

Definition at line 100 of file TileCondToolOfc.cxx.

107{
108 ATH_MSG_DEBUG( "TileCondToolOfc weights, drawerIdx:" << drawerIdx
109 << " channel: " << channel
110 << " gain: " << gain
111 << " phase: " << phase );
112
113 std::fill (std::begin(weights.g), std::end(weights.g), 0);
114 std::fill (std::begin(weights.dg), std::end(weights.dg), 0);
115 std::fill (std::begin(weights.w_a), std::end(weights.w_a), 0);
116 std::fill (std::begin(weights.w_b), std::end(weights.w_b), 0);
117 std::fill (std::begin(weights.w_c), std::end(weights.w_c), 0);
118
119
120 CLHEP::HepMatrix Correlation(m_nSamples, m_nSamples, 1);
121 //aa , Inverse(m_nSamples,m_nSamples,0), Zero(m_nSamples,m_nSamples,0);
122 CLHEP::HepMatrix PulseShape(m_nSamples, 1, 0), DPulseShape(m_nSamples, 1, 0);
123 CLHEP::HepMatrix a(m_nSamples, 1, 0), b(m_nSamples, 1, 0);
124
125
126 //int ierr=0;
127
128 if (msgLvl(MSG::DEBUG)) {
129 msg(MSG::DEBUG) << " Calculating " << (of2 ? "OF2" : "OF1") << " weights ";
130
131 if (m_deltaCorrelation)
132 msg(MSG::DEBUG) << "with Delta correlation matrix " << endmsg;
133 else
134 msg(MSG::DEBUG) << "with correlation matrix obtained from DB " << endmsg;
135
136 msg(MSG::DEBUG) << "for drawerIdx= " << drawerIdx
137 << " channel=" << channel
138 << " gain=" << gain
139 << " phase " << phase << endmsg;
140 }
141
142 int npr = 2;
143 if (of2) npr = 3;
144 weights.n_samples = m_nSamples;
145 weights.of2 = of2;
146
147 if (!m_deltaCorrelation) { //=== Retrieve autocorrelations from COOL DB
148 std::vector<float> vecAutoCr;
149 m_tileToolAutoCr->getAutoCorr(drawerIdx, channel, gain, vecAutoCr);
150
151 ATH_MSG_DEBUG( " vecAutoCr size: " << vecAutoCr.size() );
152 ATH_MSG_VERBOSE( " vecAutoCr " << vecAutoCr );
153 if (vecAutoCr[0] > -1233.){
154 if (vecAutoCr.size() == 28){
155 int AutoCrPosition = 0;
156 for (int i = 0; i < m_nSamples; ++i)
157 for (int j = i; j < m_nSamples; ++j){
158 Correlation[i][j] = vecAutoCr[AutoCrPosition];
159 Correlation[j][i] = vecAutoCr[AutoCrPosition];
160 ++AutoCrPosition;
161 }
162 } else {
163 for (int i = 0; i < m_nSamples; i++)
164 Correlation[i][i] = 1.;
165 for (int i = 0; i < m_nSamples - 1; i++) {// Fill non-diag. elements
166 for (int j = 0; j < m_nSamples - i - 1; j++) {
167 Correlation[i][j + i + 1] = vecAutoCr[j];
168 Correlation[j + i + 1][i] = vecAutoCr[j];
169 }
170 }
171 }
172 }
173 }
174
175 ATH_MSG_VERBOSE( " Correlation: " << Correlation );
176
177 float py = 0.;
178 float pdy = 0.;
179 //if (ierr==0)
180
181 {
182 for (int i = 0; i < m_nSamples; i++) {
183 m_tileToolPulseShape->getPulseShapeYDY(drawerIdx, channel, gain
184 , phase + 25 * (i - m_t0Sample), py, pdy, ctx);
185
186 // If overflow/underflow using boundary value of y and dy=0, OF affected when the number of samples is high. Set y=0 and dy=0 for OF weights computation
187 if( (phase + 25 * (i - m_t0Sample) < -150) || (phase + 25 * (i - m_t0Sample) > 150) ){
188 py=0;
189 pdy=0;
190 }
191
192
193 PulseShape[i][0] = py;
194 DPulseShape[i][0] = pdy;
195 weights.g[i] = PulseShape[i][0];
196 weights.dg[i] = DPulseShape[i][0];
197
198 ATH_MSG_VERBOSE( " Pulse shape: isamp " << i
199 << " phase " << phase
200 << " py " << py
201 << " pdy " << pdy );
202
203 }
204
205 // Build System Matrix with Correlations and pulse function points
206
207 CLHEP::HepMatrix SystemMatrix(m_nSamples + npr, m_nSamples + npr, 0);
208 CLHEP::HepVector Result(m_nSamples + npr, 0);
209 CLHEP::HepVector IndependTermsAmp(m_nSamples + npr, 0);
210 CLHEP::HepVector IndependTermsTime(m_nSamples + npr, 0);
211 CLHEP::HepVector IndependTermsPed(m_nSamples + npr, 0);
212
213 for (int i = 0; i < m_nSamples; i++) {
214 for (int j = 0; j < m_nSamples; j++) {
215 SystemMatrix[i][j] = Correlation[i][j];
216 }
217
218 SystemMatrix[m_nSamples][i] = PulseShape[i][0];
219 SystemMatrix[i][m_nSamples] = PulseShape[i][0];
220
221 SystemMatrix[m_nSamples + 1][i] = DPulseShape[i][0];
222 SystemMatrix[i][m_nSamples + 1] = DPulseShape[i][0];
223
224 if (of2) {
225 SystemMatrix[m_nSamples + 2][i] = 1.;
226 SystemMatrix[i][m_nSamples + 2] = 1.;
227 }
228 }
229
230 IndependTermsAmp[m_nSamples] = 1.;
231 IndependTermsTime[m_nSamples + 1] = -1.;
232
233 if (of2) IndependTermsPed[m_nSamples + 2] = 1.;
234
235 Result = solve(SystemMatrix, IndependTermsAmp);
236
237 for (int ismp = 0; ismp < m_nSamples; ismp++) {
238 weights.w_a[ismp] = (double) Result[ismp];
239// ATH_MSG_DEBUG( "w_a " << weights.w_a[ismp] << " ismp " << ismp );
240
241 }
242
243 Result = solve(SystemMatrix, IndependTermsTime);
244
245 for (int ismp = 0; ismp < m_nSamples; ismp++) {
246 weights.w_b[ismp] = (double) Result[ismp];
247// ATH_MSG_DEBUG( "w_b " << weights.w_b[ismp] << " ismp " << ismp );
248 }
249
250 if (of2) { // OF2
251 Result = solve(SystemMatrix, IndependTermsPed);
252
253 for (int ismp = 0; ismp < m_nSamples; ismp++) {
254 weights.w_c[ismp] = (double) Result[ismp];
255// ATH_MSG_DEBUG( "w_c " << weights.w_c[ismp] << " ismp " << ismp );
256 }
257 }
258
259 }
260
261#if 0
262 else {
263 for (int ismp=0; ismp<m_nSamples; ismp++) {
264 weights.w_a[ismp] = 0.;
265 weights.w_b[ismp] = 0.;
266 weights.w_c[ismp] = 0.;
267 weights.g[ismp] = 0.;
268 weights.dg[ismp] = 0.;
269 }
270 }
271#endif
272
273 ATH_MSG_DEBUG( "...OFC weights fixed-phase calculated" );
274
275 return StatusCode::SUCCESS;
276}
endmsg
#define endmsg
Definition AnalysisConfig_Ntuple.cxx:63
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition AthMsgStreamMacros.h:28
Result
ICscStripFitter::Result Result
Definition CalibCscStripFitter.cxx:13
a
static Double_t a
Definition LArPhysWaveHECTool.cxx:38
TileCondToolOfc::m_tileToolAutoCr
ToolHandle< TileCondToolAutoCr > m_tileToolAutoCr
Definition TileCondToolOfc.h:66
TileCondToolOfc::m_tileToolPulseShape
ToolHandle< TileCondToolPulseShape > m_tileToolPulseShape
Definition TileCondToolOfc.h:63
CaloCondBlobAlgs_fillNoiseFromASCII.gain
gain
Definition CaloCondBlobAlgs_fillNoiseFromASCII.py:109
plotting.yearwise_efficiency.channel
channel
Definition yearwise_efficiency.py:24
msg
MsgStream & msg
Definition testRead.cxx:32

◆ initialize()

StatusCode TileCondToolOfc::initialize ( )
overridevirtual

Definition at line 36 of file TileCondToolOfc.cxx.

36 {
37
38 ATH_MSG_INFO( "In initialize()" );
39
40 //=== Get TileCondToolPulseShape
41 ATH_CHECK( m_tileToolPulseShape.retrieve() );
42
43 if (!m_deltaCorrelation) {
44 ATH_CHECK(m_tileToolAutoCr.retrieve());
45 } else {
46 m_tileToolAutoCr.disable();
47 }
48
49 //==== TileInfo
50 ATH_CHECK( detStore()->retrieve(m_tileInfo, "TileInfo") );
51
52 if (m_nSamples != m_tileInfo->NdigitSamples()) {
53 ATH_MSG_WARNING( "Changing number of samples from " << m_nSamples
54 << " to " << m_tileInfo->NdigitSamples() );
55
56 m_nSamples = m_tileInfo->NdigitSamples();
57 } else {
58 ATH_MSG_INFO( "Number of samples is " << m_nSamples );
59 }
60
61 m_t0Sample = (m_nSamples - 1) / 2;
62 if (m_t0Sample != m_tileInfo->ItrigSample()) {
63 ATH_MSG_WARNING( "Changing T0 sample from " << m_t0Sample
64 << " to " << m_tileInfo->ItrigSample() );
65
66 m_t0Sample = m_tileInfo->ItrigSample();
67 } else {
68 ATH_MSG_INFO( "T0 sample is " << m_t0Sample );
69 }
70
71 //=== Initialize max values
72 ServiceHandle<TileCablingSvc> cablingSvc("TileCablingSvc", name());
73 ATH_CHECK( cablingSvc.retrieve());
74
75 const TileCablingService* cabling = cablingSvc->cablingService();
76 if (!cabling) {
77 ATH_MSG_ERROR( "Unable to retrieve TileCablingService" );
78 return StatusCode::FAILURE;
79 }
80
81
82 //--------------------------------------------------------
83 ATH_MSG_INFO( "TileCondToolOfc initialization completed. " );
84
85 return StatusCode::SUCCESS;
86}
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition AthMsgStreamMacros.h:33
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition AthMsgStreamMacros.h:31
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition AthMsgStreamMacros.h:32
python.PyKernel.retrieve
retrieve(aClass, aKey=None)
Definition PyKernel.py:110

Member Data Documentation

◆ m_deltaCorrelation

bool TileCondToolOfc::m_deltaCorrelation
private

Definition at line 74 of file TileCondToolOfc.h.

◆ m_nSamples

int TileCondToolOfc::m_nSamples
private

Definition at line 72 of file TileCondToolOfc.h.

◆ m_t0Sample

int TileCondToolOfc::m_t0Sample
private

Definition at line 73 of file TileCondToolOfc.h.

◆ m_tileInfo

const TileInfo* TileCondToolOfc::m_tileInfo
private

Definition at line 69 of file TileCondToolOfc.h.

◆ m_tileToolAutoCr

ToolHandle<TileCondToolAutoCr> TileCondToolOfc::m_tileToolAutoCr
private
Initial value:
{this,
"TileCondToolAutoCr", "TileCondToolAutoCr", "Tile auto correlation matrix tool"}

Definition at line 66 of file TileCondToolOfc.h.

66 {this,
67 "TileCondToolAutoCr", "TileCondToolAutoCr", "Tile auto correlation matrix tool"};

◆ m_tileToolPulseShape

ToolHandle<TileCondToolPulseShape> TileCondToolOfc::m_tileToolPulseShape
private
Initial value:
{this,
"TileCondToolPulseShape", "TileCondToolPulseShape", "Tile pulse shape tool"}

Definition at line 63 of file TileCondToolOfc.h.

63 {this,
64 "TileCondToolPulseShape", "TileCondToolPulseShape", "Tile pulse shape tool"};
The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.14.0