]>
Commit | Line | Data |
---|---|---|
1b41ab20 | 1 | // $Id$ |
2 | ||
b85c748e | 3 | /************************************************************************** |
4 | * Copyright(c) 2006, ALICE Experiment at CERN, All rights reserved. * | |
5 | * * | |
6 | * Authors: Per Thomas Hille for the ALICE * | |
7 | * offline/HLT Project. Contributors are mentioned in the code where * | |
8 | * appropriate. * | |
9 | * * | |
10 | * Permission to use, copy, modify and distribute this software and its * | |
11 | * documentation strictly for non-commercial purposes is hereby granted * | |
12 | * without fee, provided that the above copyright notice appears in all * | |
13 | * copies and that both the copyright notice and this permission notice * | |
14 | * appear in the supporting documentation. The authors make no claims * | |
15 | * about the suitability of this software for any purpose. It is * | |
16 | * provided "as is" without express or implied warranty. * | |
17 | **************************************************************************/ | |
18 | ||
19 | ||
20 | ||
2589c3a3 | 21 | |
b85c748e | 22 | #include "AliHLTPHOSRcuAltroPatternTestComponent.h" |
23 | #include "AliHLTPHOSSharedMemoryInterface.h" | |
24 | #include "AliHLTPHOSValidCellDataStruct.h" | |
25 | #include "AliHLTPHOSRcuAltroPatternTest.h" | |
26 | ||
27 | AliHLTPHOSRcuAltroPatternTestComponent gAliHLTPHOSRcuAltroPatternTestComponent; | |
28 | ||
29 | AliHLTPHOSRcuAltroPatternTestComponent:: AliHLTPHOSRcuAltroPatternTestComponent() : AliHLTPHOSRcuProcessor(), | |
30 | fPatternTestPtr(0), | |
31 | fShmPtr(0), | |
32 | fNTotalPatterns(0), | |
33 | fNWrongPatterns(0), | |
34 | fNTotalSamples(0), | |
503dd508 | 35 | fNWrongSamples(0), |
36 | fNSamples(0) | |
b85c748e | 37 | { |
38 | fShmPtr = new AliHLTPHOSSharedMemoryInterface(); | |
39 | } | |
40 | ||
41 | ||
42 | AliHLTPHOSRcuAltroPatternTestComponent::~ AliHLTPHOSRcuAltroPatternTestComponent() | |
43 | { | |
44 | //Destructor | |
45 | } | |
46 | ||
47 | ||
48 | int | |
49 | AliHLTPHOSRcuAltroPatternTestComponent::Deinit() | |
50 | { | |
51 | //See html documentation of base class | |
52 | return 0; | |
53 | } | |
54 | ||
55 | ||
56 | const char* | |
57 | AliHLTPHOSRcuAltroPatternTestComponent::GetComponentID() | |
58 | { | |
59 | //See html documentation of base class | |
60 | return "AltroPatternTester"; | |
61 | } | |
62 | ||
63 | ||
64 | void | |
65 | AliHLTPHOSRcuAltroPatternTestComponent::GetInputDataTypes( vector<AliHLTComponentDataType>& list) | |
66 | { | |
67 | //See html documentation of base class | |
68 | const AliHLTComponentDataType* pType=fgkInputDataTypes; | |
69 | while (pType->fID!=0) | |
70 | { | |
71 | list.push_back(*pType); | |
72 | pType++; | |
73 | } | |
74 | } | |
75 | ||
76 | ||
77 | AliHLTComponentDataType | |
78 | AliHLTPHOSRcuAltroPatternTestComponent::GetOutputDataType() | |
79 | { | |
80 | //See html documentation of base class | |
81 | return AliHLTPHOSDefinitions::fgkCellEnergyDataType; | |
82 | } | |
83 | ||
84 | ||
85 | void | |
86 | AliHLTPHOSRcuAltroPatternTestComponent::GetOutputDataSize(unsigned long& constBase, double& inputMultiplier ) | |
87 | { | |
88 | //See html documentation of base class | |
89 | constBase = 30; | |
90 | inputMultiplier = 1; | |
91 | } | |
92 | ||
93 | ||
94 | ||
95 | int AliHLTPHOSRcuAltroPatternTestComponent::DoEvent( const AliHLTComponentEventData& evtData, const AliHLTComponentBlockData* blocks, | |
96 | AliHLTComponentTriggerData& /*trigData*/, AliHLTUInt8_t* outputPtr, | |
97 | AliHLTUInt32_t& size, vector<AliHLTComponentBlockData>& /*outputBlocks */) | |
98 | { | |
99 | //See html documentation of base class | |
100 | //cout << "AliHLTPHOSRcuAltroPatternTestComponent::DoEvent, processing event "<< fPhosEventCount << endl; | |
101 | ||
102 | AliHLTPHOSValidCellDataStruct *currentChannel =0; | |
103 | unsigned long ndx = 0; | |
104 | // UInt_t offset = 0; | |
105 | // UInt_t mysize = 0; | |
106 | UInt_t tSize = 0; | |
107 | const AliHLTComponentBlockData* iter = NULL; | |
108 | AliHLTPHOSRcuCellEnergyDataStruct *cellDataPtr; | |
109 | AliHLTUInt8_t* outBPtr; | |
110 | ||
b85c748e | 111 | for( ndx = 0; ndx < evtData.fBlockCnt; ndx++ ) |
112 | { | |
113 | iter = blocks+ndx; | |
114 | if(iter->fDataType != AliHLTPHOSDefinitions::fgkCellEnergyDataType) | |
115 | { | |
116 | continue; | |
117 | } | |
118 | ||
119 | cellDataPtr = (AliHLTPHOSRcuCellEnergyDataStruct*)( iter->fPtr); | |
120 | fShmPtr->SetMemory(cellDataPtr); | |
121 | currentChannel = fShmPtr->NextChannel(); | |
122 | ||
7ce40e5c | 123 | Int_t* tmpRawPtr = 0; |
124 | ||
b85c748e | 125 | int tmp =0; |
126 | while(currentChannel != 0) | |
127 | { | |
128 | ||
129 | int tmpZ = currentChannel->fZ; | |
130 | int tmpX = currentChannel->fX; | |
131 | int tmpGain = currentChannel->fGain; | |
7ce40e5c | 132 | int tmpSamples = 0; |
133 | tmpRawPtr = fShmPtr->GetRawData(tmpSamples); | |
b85c748e | 134 | |
135 | ||
27029341 | 136 | if( (tmpZ > NZROWSRCU) || (tmpX > NXCOLUMNSRCU) || (tmpGain > NGAINS)) |
b85c748e | 137 | { |
138 | cout <<" ERROR parameters out of range z = "<< tmpZ <<" x = "<< tmpX<< " gain = " << tmpGain<<" nSamples = " << tmpSamples <<endl; | |
139 | } | |
140 | else | |
141 | { | |
142 | // cout <<" all parameters in range" << endl; | |
143 | } | |
144 | ||
145 | // cout << "analyzing channelnr " << tmp <<" of event " << fPhosEventCount << endl; | |
146 | // cout << " z = " << currentChannel->fZ << " x = " << currentChannel->fX << " gain = " << currentChannel->fGain << " samples =" << currentChannel->fNSamples <<endl; | |
147 | tmp ++; | |
148 | ||
7ce40e5c | 149 | fPatternTestPtr->AddPattern(tmpRawPtr, currentChannel->fZ, currentChannel->fX, currentChannel->fGain, tmpSamples, fNPresamples); |
150 | int ret = fPatternTestPtr->ValidateAltroPattern(tmpRawPtr, tmpSamples); | |
b85c748e | 151 | |
152 | if(ret >= 0) | |
153 | { | |
154 | ||
7ce40e5c | 155 | fNTotalSamples += tmpSamples; |
b85c748e | 156 | fNTotalPatterns ++; |
157 | ||
158 | if(ret > 0) | |
159 | { | |
160 | fNWrongSamples += ret; | |
161 | fNWrongPatterns ++; | |
162 | ||
163 | if(fPhosEventCount%100 == 0) | |
164 | { | |
165 | /* | |
166 | cout << "warning: incorrect pattern found for event " << fPhosEventCount <<" X = " << currentChannel->fX <<" Z = " << currentChannel->fZ ; | |
167 | cout << "number of correct aptterns is " << fNTotalPatterns << "number of wrong patterns is " << fNWrongPatterns<<endl; | |
168 | float percent = 100*((float)fNWrongPatterns)/( float(fNTotalPatterns)); | |
169 | cout <<"The corrpution rate is " << percent << " percent " <<endl; | |
170 | cout <<" ERROR: incorrect parameters" << endl; | |
171 | */ | |
172 | } | |
173 | ||
174 | ||
175 | if(fPhosEventCount%100 == 0 && fPhosEventCount !=0) | |
176 | { | |
177 | fPatternTestPtr->countAllPatterns(fNSamples); | |
178 | fPatternTestPtr->PrintStatistics(); | |
179 | ||
180 | } | |
181 | } | |
182 | } | |
183 | else | |
184 | { | |
185 | cout <<" ERROR: incorrect parameters" << endl; | |
186 | } | |
187 | currentChannel = fShmPtr->NextChannel(); | |
188 | ||
189 | } | |
190 | } | |
191 | ||
192 | ||
193 | outBPtr = outputPtr; | |
194 | ||
195 | // mysize += sizeof(AliHLTPHOSRcuCellAccumulatedEnergyDataStruct); | |
196 | ||
197 | /* | |
198 | AliHLTComponentBlockData bd; | |
199 | FillBlockData( bd ); | |
200 | bd.fOffset = offset; | |
201 | bd.fSize = mysize; | |
202 | bd.fDataType = AliHLTPHOSDefinitions::fgkCellAccumulatedEnergyDataType; | |
203 | bd.fSpecification = 0xFFFFFFFF; | |
204 | outputBlocks.push_back( bd ); | |
205 | tSize += mysize; | |
206 | outBPtr += mysize; | |
207 | */ | |
208 | ||
209 | if( tSize > size ) | |
210 | { | |
211 | Logging( kHLTLogFatal, "HLT::AliHLTPHOSRcuAltroPatternTestComponent::DoEvent", "Too much data", | |
212 | "Data written over allowed buffer. Amount written: %lu, allowed amount: %lu." | |
213 | , tSize, size ); | |
214 | return EMSGSIZE; | |
215 | } | |
216 | ||
217 | fPhosEventCount++; | |
218 | ||
219 | ||
220 | return 0; | |
221 | }//end DoEvent | |
222 | ||
223 | ||
224 | ||
225 | ||
226 | ||
227 | int | |
228 | AliHLTPHOSRcuAltroPatternTestComponent::DoInit(int argc, const char** argv ) | |
229 | { | |
230 | //See html documentation of base class | |
231 | cout << "AliHLTPHOSRcuAltroPatternTestComponent::DoInit, argc =" << argc << endl; | |
232 | char patternFilename[2000]; | |
233 | bool isSetPartternArg = false; | |
234 | fPrintInfo = kFALSE; | |
235 | int iResult=0; | |
236 | TString argument=""; | |
237 | iResult = ScanArguments(argc, argv); | |
238 | // fPatternTestPtr = new AliHLTPHOSRcuAltroPatternTest(fModuleID, fRcuX, fRcuZ); | |
239 | ||
240 | if(iResult < 0) | |
241 | { | |
242 | return iResult; | |
243 | } | |
244 | ||
245 | if(fIsSetEquippmentID == kFALSE) | |
246 | { | |
247 | Logging( kHLTLogFatal, "HLT::AliHLTPHOSRcuAltroPatternTestComponent::DoInit( int argc, const char** argv )", "Missing argument", | |
248 | "The argument equippmentID is not set: set it with a component argumet like this: -equippmentID <number>"); | |
249 | // iResult = -2; | |
250 | return -5; | |
251 | } | |
252 | ||
253 | for(int i = 0; i < argc; i++) | |
254 | { | |
255 | if(!strcmp("-patternfile", argv[i])) | |
256 | { | |
257 | if(argc < (i+1)) | |
258 | { | |
259 | Logging( kHLTLogFatal, "HLT::AliHLTPHOSRcuAltroPatternTestComponent::DoInit( int argc, const char** argv )", "Missing argument", | |
260 | "The argument -patternfile requires that a filename follows: set it with a component argumet like this: -patternfile <filename>"); | |
261 | return -6; | |
262 | } | |
263 | else | |
264 | { | |
e304ea31 | 265 | sprintf(patternFilename, "%s", argv[i+1]); |
b85c748e | 266 | isSetPartternArg = true; |
267 | } | |
268 | } | |
503dd508 | 269 | if(!strcmp("-nsamples", argv[i])) |
270 | { | |
271 | fNSamples = atoi(argv[i+1]); | |
272 | } | |
b85c748e | 273 | } |
274 | ||
275 | if(isSetPartternArg == false) | |
276 | { | |
277 | Logging( kHLTLogFatal, "HLT::AliHLTPHOSRcuAltroPatternTestComponent::DoInit( int argc, const char** argv )", "Missing argument", | |
278 | "The argument patternfile is not set: set it with a component argumet like this: -patternfile <filename>"); | |
279 | return -7; | |
280 | ||
281 | } | |
2589c3a3 | 282 | else if(fUtilitiesPtr->CheckFile(patternFilename, "r") == false) |
b85c748e | 283 | { |
284 | char tmpMessage[1024]; | |
285 | sprintf(tmpMessage, "the file %s could not be found, please check that file exist and that you have read access to it", patternFilename); | |
286 | Logging( kHLTLogFatal, "HLT::AliHLTPHOSRcuAltroPatternTestComponent::DoInit( int argc, const char** argv )", "File not found",tmpMessage); | |
287 | ||
288 | return -8; | |
289 | ||
290 | } | |
291 | else | |
292 | { | |
27029341 | 293 | int tmpPattern[ALTROMAXSAMPLES]; |
294 | ScanPatternFromFile(patternFilename, tmpPattern, ALTROMAXSAMPLES) ; | |
295 | fPatternTestPtr = new AliHLTPHOSRcuAltroPatternTest(fModuleID, fRcuX, fRcuZ, tmpPattern, ALTROMAXSAMPLES); | |
b85c748e | 296 | } |
297 | ||
298 | // return iResult; | |
299 | return 0; | |
300 | } | |
301 | ||
302 | ||
303 | AliHLTComponent* | |
304 | AliHLTPHOSRcuAltroPatternTestComponent::Spawn() | |
305 | { | |
306 | //See html documentation of base class | |
307 | return new AliHLTPHOSRcuAltroPatternTestComponent; | |
308 | } | |
309 | ||
310 | ||
311 | void | |
6e587e00 | 312 | AliHLTPHOSRcuAltroPatternTestComponent::ScanPatternFromFile(const char *filename, int *pattern, const int /*length*/) const |
b85c748e | 313 | { |
314 | FILE *fp = fopen(filename, "r"); | |
315 | ||
316 | // int tmpPattern[ALTRO_MAX_SAMPLES]; | |
317 | int dummy = 0; | |
e304ea31 | 318 | int res = 0; |
27029341 | 319 | for(int i=0; i<ALTROMAXSAMPLES; i++) |
b85c748e | 320 | { |
e304ea31 | 321 | res = fscanf(fp,"w 0x%X 0x%X\n", &dummy, &pattern[i]); |
b85c748e | 322 | // cout << tmpPattern[i] << endl; |
323 | } | |
6e587e00 | 324 | |
b85c748e | 325 | |
326 | // fPatternTestPtr->SetAltroPattern(tmpPattern); | |
327 | } | |
328 | ||
329 | //w 0x6801 0x1 |