04fa961a |
1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
3 | * * |
4 | * Author: The ALICE Off-line Project. * |
5 | * Contributors are mentioned in the code where appropriate. * |
6 | * * |
7 | * Permission to use, copy, modify and distribute this software and its * |
8 | * documentation strictly for non-commercial purposes is hereby granted * |
9 | * without fee, provided that the above copyright notice appears in all * |
10 | * copies and that both the copyright notice and this permission notice * |
11 | * appear in the supporting documentation. The authors make no claims * |
12 | * about the suitability of this software for any purpose. It is * |
13 | * provided "as is" without express or implied warranty. * |
14 | **************************************************************************/ |
15 | |
16 | /* $Id$ */ |
17 | |
18 | /////////////////////////////////////////////////////////////////////////////// |
19 | // // |
20 | // This is the class for perfoming the monitoring process. // |
21 | // It checks if a raw data file exists, loops over the events in the raw // |
22 | // data file, reconstructs TPC and ITS clusters and tracks, fills the // |
23 | // monitor histograms and sends the updated histograms to the clients. // |
24 | // Then the raw data file is deleted and it waits for a new file. // |
25 | // // |
26 | /////////////////////////////////////////////////////////////////////////////// |
27 | |
28 | |
29 | #include "AliMonitorProcess.h" |
30 | #include "AliMonitorTPC.h" |
31 | #include "AliMonitorITS.h" |
32 | #include "AliMonitorV0s.h" |
97d6eb66 |
33 | #include "AliMonitorHLT.h" |
04fa961a |
34 | #include "AliRawReaderRoot.h" |
35 | #include "AliLoader.h" |
36 | #include "AliRun.h" |
37 | #include "AliTPC.h" |
38 | #include "AliTPCclustererMI.h" |
39 | #include "AliTPCtrackerMI.h" |
40 | #include "AliITS.h" |
41 | #include "AliITSclustererV2.h" |
42 | #include "AliITStrackerV2.h" |
43 | #include "AliITSLoader.h" |
44 | #include "AliV0vertexer.h" |
c4bd737c |
45 | #include <TGrid.h> |
04fa961a |
46 | #include <TSocket.h> |
c4bd737c |
47 | #include <TServerSocket.h> |
04fa961a |
48 | #include <TMessage.h> |
49 | #include <TGridResult.h> |
50 | #include <TROOT.h> |
1899848d |
51 | #ifdef ALI_HLT |
c4bd737c |
52 | #include <AliLevel3.h> |
1899848d |
53 | #include <AliL3Transform.h> |
54 | #endif |
04fa961a |
55 | |
56 | ClassImp(AliMonitorProcess) |
57 | |
58 | |
c4bd737c |
59 | const Int_t AliMonitorProcess::fgkPort = 9327; |
04fa961a |
60 | |
61 | |
62 | //_____________________________________________________________________________ |
63 | AliMonitorProcess::AliMonitorProcess(const char* alienDir, |
64 | const char* fileNameGalice) |
65 | { |
66 | // initialize the monitoring process and the monitor histograms |
67 | |
68 | fGrid = TGrid::Connect("alien", gSystem->Getenv("USER")); |
69 | if (!fGrid || fGrid->IsZombie() || !fGrid->IsConnected()) { |
70 | delete fGrid; |
71 | Fatal("AliMonitorProcess", "could not connect to alien"); |
72 | } |
73 | fGrid->cd(alienDir); |
74 | fLogicalFileName = ""; |
75 | fFileName = ""; |
76 | |
77 | fRunLoader = AliRunLoader::Open(fileNameGalice); |
78 | if (!fRunLoader) Fatal("AliMonitorProcess", |
79 | "could not get run loader from file %s", |
80 | fileNameGalice); |
81 | |
82 | fRunLoader->CdGAFile(); |
83 | fTPCParam = AliTPC::LoadTPCParam(gFile); |
84 | if (!fTPCParam) Fatal("AliMonitorProcess", "could not load TPC parameters"); |
85 | |
86 | fRunLoader->LoadgAlice(); |
87 | gAlice = fRunLoader->GetAliRun(); |
88 | if (!gAlice) Fatal("AliMonitorProcess", "no gAlice object found"); |
c4bd737c |
89 | AliITS* its = (AliITS*) gAlice->GetModule("ITS"); |
90 | if (!its) Fatal("AliMonitorProcess", "no ITS detector found"); |
91 | fITSgeom = its->GetITSgeom(); |
04fa961a |
92 | if (!fITSgeom) Fatal("AliMonitorProcess", "could not load ITS geometry"); |
93 | |
1899848d |
94 | #ifdef ALI_HLT |
95 | // Init TPC parameters for HLT |
96 | Bool_t isinit=AliL3Transform::Init(const_cast<char*>(fileNameGalice),kTRUE); |
97 | if(!isinit){ |
98 | cerr << "Could not create transform settings, please check log for error messages!" << endl; |
99 | return; |
100 | } |
101 | #endif |
102 | |
04fa961a |
103 | fRunNumber = 0; |
104 | fSubRunNumber = 0; |
105 | fNEvents = 0; |
106 | fNEventsMin = 2; |
107 | fWriteHistoList = kFALSE; |
108 | |
109 | fTopFolder = new TFolder("Monitor", "monitor histograms"); |
110 | fTopFolder->SetOwner(kTRUE); |
111 | |
112 | fMonitors.Add(new AliMonitorTPC(fTPCParam)); |
113 | fMonitors.Add(new AliMonitorITS(fITSgeom)); |
114 | fMonitors.Add(new AliMonitorV0s); |
97d6eb66 |
115 | #ifdef ALI_HLT |
116 | fMonitors.Add(new AliMonitorHLT(fTPCParam)); |
117 | #endif |
04fa961a |
118 | |
119 | for (Int_t iMonitor = 0; iMonitor < fMonitors.GetEntriesFast(); iMonitor++) { |
120 | ((AliMonitor*) fMonitors[iMonitor])->CreateHistos(fTopFolder); |
121 | } |
122 | |
123 | TIterator* iFolder = fTopFolder->GetListOfFolders()->MakeIterator(); |
124 | while (TFolder* folder = (TFolder*) iFolder->Next()) folder->SetOwner(kTRUE); |
125 | delete iFolder; |
126 | |
127 | fFile = TFile::Open("monitor_tree.root", "RECREATE"); |
128 | if (!fFile || !fFile->IsOpen()) { |
129 | Fatal("AliMonitorProcess", "could not open file for tree"); |
130 | } |
131 | fTree = new TTree("MonitorTree", "tree for monitoring"); |
132 | for (Int_t iMonitor = 0; iMonitor < fMonitors.GetEntriesFast(); iMonitor++) { |
133 | ((AliMonitor*) fMonitors[iMonitor])->CreateBranches(fTree); |
134 | } |
135 | gROOT->cd(); |
136 | |
c4bd737c |
137 | fServerSocket = new TServerSocket(fgkPort, kTRUE); |
04fa961a |
138 | fServerSocket->SetOption(kNoBlock, 1); |
139 | fDisplaySocket = NULL; |
140 | CheckForConnections(); |
97d6eb66 |
141 | #ifdef ALI_HLT |
142 | fHLT = NULL; |
143 | #endif |
04fa961a |
144 | |
1899848d |
145 | SetStatus(kStopped); |
04fa961a |
146 | fStopping = kFALSE; |
147 | } |
148 | |
c4bd737c |
149 | //_____________________________________________________________________________ |
150 | AliMonitorProcess::AliMonitorProcess(const AliMonitorProcess& process) : |
151 | TObject(process) |
152 | { |
153 | Fatal("AliMonitorProcess", "copy constructor not implemented"); |
154 | } |
155 | |
156 | //_____________________________________________________________________________ |
157 | AliMonitorProcess& AliMonitorProcess::operator = (const AliMonitorProcess& |
158 | /*process*/) |
159 | { |
160 | Fatal("operator =", "assignment operator not implemented"); |
161 | return *this; |
162 | } |
163 | |
04fa961a |
164 | //_____________________________________________________________________________ |
165 | AliMonitorProcess::~AliMonitorProcess() |
166 | { |
167 | // clean up |
168 | |
169 | fMonitors.Delete(); |
170 | delete fTopFolder; |
171 | delete fRunLoader; |
172 | |
173 | delete fServerSocket; |
174 | fSockets.Delete(); |
175 | delete fDisplaySocket; |
176 | |
177 | fGrid->Close(); |
178 | delete fGrid; |
179 | |
180 | fFile->Close(); |
181 | delete fFile; |
182 | gSystem->Unlink("monitor_tree.root"); |
97d6eb66 |
183 | |
184 | #ifdef ALI_HLT |
185 | delete fHLT; |
186 | #endif |
04fa961a |
187 | } |
188 | |
189 | |
190 | //_____________________________________________________________________________ |
191 | const char* AliMonitorProcess::GetRevision() |
192 | { |
193 | return "$Revision$"; |
194 | } |
195 | |
196 | |
1899848d |
197 | //_____________________________________________________________________________ |
198 | void AliMonitorProcess::SetStatus(EStatus status) |
199 | { |
c4bd737c |
200 | // set the current status and process system events |
201 | |
1899848d |
202 | fStatus = status; |
203 | gSystem->ProcessEvents(); |
204 | } |
205 | |
206 | |
04fa961a |
207 | //_____________________________________________________________________________ |
208 | void AliMonitorProcess::Run() |
209 | { |
210 | // run the monitor process: |
211 | // check for a raw data file, process the raw data file and delete it |
212 | |
213 | fStopping = kFALSE; |
214 | |
215 | while (!fStopping) { |
1899848d |
216 | SetStatus(kWaiting); |
04fa961a |
217 | while (!CheckForNewFile()) { |
218 | CheckForConnections(); |
1899848d |
219 | SetStatus(kWaiting); |
04fa961a |
220 | if (fStopping) break; |
221 | gSystem->Sleep(10); |
222 | } |
223 | if (fStopping) break; |
224 | |
225 | ProcessFile(); |
226 | } |
227 | |
228 | WriteHistos(); |
229 | |
230 | fStopping = kFALSE; |
1899848d |
231 | SetStatus(kStopped); |
04fa961a |
232 | } |
233 | |
234 | |
235 | //_____________________________________________________________________________ |
236 | void AliMonitorProcess::Stop() |
237 | { |
238 | // set the fStopping flag to terminate the monitor process after the current |
239 | // event was processed |
240 | |
1899848d |
241 | if (GetStatus() != kStopped) fStopping = kTRUE; |
04fa961a |
242 | } |
243 | |
244 | |
245 | //_____________________________________________________________________________ |
246 | void AliMonitorProcess::ProcessFile(const char* fileName) |
247 | { |
248 | // create a file with monitor histograms for a single file |
249 | |
1899848d |
250 | if (GetStatus() != kStopped) { |
04fa961a |
251 | Error("ProcessFile", "ProcessFile can not be called" |
252 | " while the monitor process is running"); |
253 | return; |
254 | } |
255 | |
256 | fFileName = fileName; |
257 | Int_t nEventMin = fNEventsMin; |
258 | fNEventsMin = 1; |
259 | ProcessFile(); |
260 | WriteHistos(); |
261 | fNEventsMin = nEventMin; |
1899848d |
262 | SetStatus(kStopped); |
04fa961a |
263 | } |
264 | |
265 | |
266 | //_____________________________________________________________________________ |
267 | Bool_t AliMonitorProcess::CheckForNewFile() |
268 | { |
269 | // check whether a new file was registered in alien |
270 | |
271 | TGridResult* result = fGrid->Ls(); |
272 | Long_t maxDate = -1; |
273 | Long_t maxTime = -1; |
274 | TString fileName; |
275 | |
276 | while (const char* entry = result->Next()) { |
277 | // entry = host_date_time.root |
278 | TString entryCopy(entry); |
279 | char* p = const_cast<char*>(entryCopy.Data()); |
6ad9b9d5 |
280 | if (!strtok(p, "_") || !p) continue; // host name |
97d6eb66 |
281 | char* dateStr = strtok(NULL, "_"); |
04fa961a |
282 | if (!dateStr || !p) continue; |
97d6eb66 |
283 | char* timeStr = strtok(NULL, "."); |
04fa961a |
284 | if (!timeStr || !p) continue; |
285 | Long_t date = atoi(dateStr); |
286 | Long_t time = atoi(timeStr); |
287 | |
288 | if ((date > maxDate) || ((date == maxDate) && (time > maxTime))) { |
289 | maxDate = date; |
290 | maxTime = time; |
291 | fileName = entry; |
292 | } |
293 | } |
294 | |
295 | if (maxDate < 0) return kFALSE; // no files found |
296 | if (fLogicalFileName.CompareTo(fileName) == 0) return kFALSE; // no new file |
297 | |
298 | fLogicalFileName = fileName; |
1899848d |
299 | TGridResult* result2 = fGrid->GetPhysicalFileNames(fLogicalFileName.Data()); |
300 | fFileName = result2->Next(); |
301 | |
04fa961a |
302 | return kTRUE; |
303 | } |
304 | |
305 | //_____________________________________________________________________________ |
306 | Bool_t AliMonitorProcess::ProcessFile() |
307 | { |
308 | // loop over all events in the raw data file, run the reconstruction |
309 | // and fill the monitor histograms |
310 | |
311 | Int_t nEvents = GetNumberOfEvents(fFileName); |
312 | if (nEvents <= 0) return kFALSE; |
313 | Info("ProcessFile", "found %d event(s) in file %s", |
314 | nEvents, fFileName.Data()); |
97d6eb66 |
315 | #ifdef ALI_HLT |
316 | CreateHLT(fFileName); |
317 | #endif |
04fa961a |
318 | |
319 | // loop over the events |
320 | for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) { |
1899848d |
321 | SetStatus(kReading); |
04fa961a |
322 | fRunLoader->SetEventNumber(0); |
323 | AliRawReaderRoot rawReader(fFileName, iEvent); |
324 | if (fStopping) break; |
325 | if (rawReader.GetRunNumber() != fRunNumber) { |
326 | WriteHistos(); |
327 | StartNewRun(); |
328 | fRunNumber = rawReader.GetRunNumber(); |
329 | fEventNumber[0] = rawReader.GetEventId()[0]; |
330 | fEventNumber[1] = rawReader.GetEventId()[1]; |
331 | fSubRunNumber = 0; |
332 | if (fStopping) break; |
333 | } |
334 | |
335 | if (!ReconstructTPC(&rawReader)) return kFALSE; |
336 | if (fStopping) break; |
337 | if (!ReconstructITS(&rawReader)) return kFALSE; |
338 | if (fStopping) break; |
339 | if (!ReconstructV0s()) return kFALSE; |
340 | if (fStopping) break; |
97d6eb66 |
341 | if (!ReconstructHLT(iEvent)) return kFALSE; |
342 | if (fStopping) break; |
04fa961a |
343 | |
344 | if (fDisplaySocket) fDisplaySocket->Send("new event"); |
345 | |
346 | Info("ProcessFile", "filling histograms..."); |
1899848d |
347 | SetStatus(kFilling); |
04fa961a |
348 | for (Int_t iMonitor = 0; iMonitor < fMonitors.GetEntriesFast(); iMonitor++) { |
349 | ((AliMonitor*) fMonitors[iMonitor])->FillHistos(fRunLoader, &rawReader); |
350 | if (fStopping) break; |
351 | } |
352 | if (fStopping) break; |
353 | |
354 | Info("ProcessFile", "updating histograms..."); |
1899848d |
355 | SetStatus(kUpdating); |
04fa961a |
356 | TIterator* iFolder = fTopFolder->GetListOfFolders()->MakeIterator(); |
357 | while (TFolder* folder = (TFolder*) iFolder->Next()) { |
358 | TIterator* iHisto = folder->GetListOfFolders()->MakeIterator(); |
359 | while (AliMonitorPlot* histo = (AliMonitorPlot*) iHisto->Next()) { |
360 | histo->Update(); |
361 | } |
362 | delete iHisto; |
363 | } |
364 | delete iFolder; |
365 | if (fStopping) break; |
366 | |
367 | Info("ProcessFile", "filling the tree..."); |
368 | fTree->Fill(); |
369 | |
370 | Info("ProcessFile", "broadcasting histograms..."); |
371 | CheckForConnections(); |
372 | BroadcastHistos(); |
373 | |
374 | fNEvents++; |
375 | if (fStopping) break; |
376 | } |
377 | |
97d6eb66 |
378 | #ifdef ALI_HLT |
379 | delete fHLT; |
380 | fHLT = NULL; |
381 | #endif |
382 | |
04fa961a |
383 | return kTRUE; |
384 | } |
385 | |
386 | //_____________________________________________________________________________ |
387 | void AliMonitorProcess::Reset() |
388 | { |
389 | // write the current histograms to a file and reset them |
390 | |
391 | if (fSubRunNumber == 0) fSubRunNumber++; |
392 | WriteHistos(); |
393 | StartNewRun(); |
394 | fSubRunNumber++; |
395 | } |
396 | |
397 | |
398 | //_____________________________________________________________________________ |
c4bd737c |
399 | UInt_t AliMonitorProcess::GetEventPeriodNumber() const |
04fa961a |
400 | { |
401 | // get the period number from the event id |
402 | |
403 | return (fEventNumber[1] >> 4); |
404 | } |
405 | |
406 | //_____________________________________________________________________________ |
c4bd737c |
407 | UInt_t AliMonitorProcess::GetEventOrbitNumber() const |
04fa961a |
408 | { |
409 | // get the orbit number from the event id |
410 | |
411 | return ((fEventNumber[1] & 0x000F) << 20) + (fEventNumber[0] >> 12); |
412 | } |
413 | |
414 | //_____________________________________________________________________________ |
c4bd737c |
415 | UInt_t AliMonitorProcess::GetEventBunchNumber() const |
04fa961a |
416 | { |
417 | // get the bunch number from the event id |
418 | |
419 | return (fEventNumber[0] % 0x0FFF); |
420 | } |
421 | |
422 | //_____________________________________________________________________________ |
c4bd737c |
423 | Int_t AliMonitorProcess::GetNumberOfEvents(const char* fileName) const |
04fa961a |
424 | { |
425 | // determine the number of events in the given raw data file |
426 | |
427 | Int_t nEvents = -1; |
428 | |
429 | TFile* file = TFile::Open(fileName); |
430 | if (!file || !file->IsOpen()) { |
431 | Error("GetNumberOfEvents", "could not open file %s", fileName); |
432 | if (file) delete file; |
433 | return -1; |
434 | } |
435 | |
436 | TTree* tree = (TTree*) file->Get("RAW"); |
437 | if (!tree) { |
438 | Error("GetNumberOfEvents", "could not find tree with raw data"); |
439 | } else { |
440 | nEvents = (Int_t) tree->GetEntries(); |
441 | } |
442 | file->Close(); |
443 | delete file; |
444 | |
445 | return nEvents; |
446 | } |
447 | |
448 | //_____________________________________________________________________________ |
449 | Bool_t AliMonitorProcess::ReconstructTPC(AliRawReader* rawReader) |
450 | { |
451 | // find TPC clusters and tracks |
452 | |
1899848d |
453 | SetStatus(kRecTPC); |
04fa961a |
454 | |
455 | AliLoader* tpcLoader = fRunLoader->GetLoader("TPCLoader"); |
456 | if (!tpcLoader) { |
457 | Error("ReconstructTPC", "no TPC loader found"); |
458 | return kFALSE; |
459 | } |
460 | gSystem->Unlink("TPC.RecPoints.root"); |
461 | gSystem->Unlink("TPC.Tracks.root"); |
462 | |
463 | // cluster finder |
464 | Info("ReconstructTPC", "reconstructing clusters..."); |
465 | tpcLoader->LoadRecPoints("recreate"); |
466 | AliTPCclustererMI clusterer(fTPCParam); |
467 | tpcLoader->MakeRecPointsContainer(); |
468 | clusterer.SetOutput(tpcLoader->TreeR()); |
469 | clusterer.Digits2Clusters(rawReader); |
470 | tpcLoader->WriteRecPoints("OVERWRITE"); |
471 | |
472 | // track finder |
473 | Info("ReconstructTPC", "reconstructing tracks..."); |
474 | tpcLoader->LoadTracks("recreate"); |
475 | { |
476 | AliTPCtrackerMI tracker(fTPCParam); |
477 | tracker.Clusters2Tracks(); |
478 | } |
479 | |
480 | tpcLoader->UnloadRecPoints(); |
481 | tpcLoader->UnloadTracks(); |
482 | return kTRUE; |
483 | } |
484 | |
485 | //_____________________________________________________________________________ |
486 | Bool_t AliMonitorProcess::ReconstructITS(AliRawReader* rawReader) |
487 | { |
488 | // find ITS clusters and tracks |
489 | |
1899848d |
490 | SetStatus(kRecITS); |
04fa961a |
491 | |
492 | AliLoader* itsLoader = fRunLoader->GetLoader("ITSLoader"); |
493 | if (!itsLoader) { |
494 | Error("ReconstructITS", "no ITS loader found"); |
495 | return kFALSE; |
496 | } |
497 | AliLoader* tpcLoader = fRunLoader->GetLoader("TPCLoader"); |
498 | if (!tpcLoader) { |
499 | Error("ReconstructITS", "no TPC loader found"); |
500 | return kFALSE; |
501 | } |
502 | gSystem->Unlink("ITS.RecPoints.root"); |
503 | gSystem->Unlink("ITS.Tracks.root"); |
504 | |
505 | // cluster finder |
506 | Info("ReconstructITS", "reconstructing clusters..."); |
507 | itsLoader->LoadRecPoints("recreate"); |
508 | AliITSclustererV2 clusterer(fITSgeom); |
509 | itsLoader->MakeRecPointsContainer(); |
510 | clusterer.Digits2Clusters(rawReader); |
511 | |
512 | // track finder |
513 | Info("ReconstructITS", "reconstructing tracks..."); |
514 | itsLoader->LoadTracks("recreate"); |
515 | itsLoader->MakeTracksContainer(); |
516 | tpcLoader->LoadTracks(); |
517 | AliITStrackerV2 tracker(fITSgeom); |
518 | tracker.LoadClusters(itsLoader->TreeR()); |
519 | tracker.Clusters2Tracks(tpcLoader->TreeT(), itsLoader->TreeT()); |
520 | tracker.UnloadClusters(); |
521 | itsLoader->WriteTracks("OVERWRITE"); |
522 | |
523 | itsLoader->UnloadRecPoints(); |
524 | itsLoader->UnloadTracks(); |
525 | tpcLoader->UnloadTracks(); |
526 | return kTRUE; |
527 | } |
528 | |
529 | //_____________________________________________________________________________ |
530 | Bool_t AliMonitorProcess::ReconstructV0s() |
531 | { |
532 | // find V0s |
533 | |
1899848d |
534 | SetStatus(kRecV0s); |
04fa961a |
535 | |
536 | AliITSLoader* itsLoader = (AliITSLoader*) fRunLoader->GetLoader("ITSLoader"); |
537 | if (!itsLoader) { |
538 | Error("ReconstructV0", "no ITS loader found"); |
539 | return kFALSE; |
540 | } |
541 | gSystem->Unlink("ITS.V0s.root"); |
542 | |
543 | // V0 finder |
544 | Info("ReconstructV0s", "reconstructing V0s..."); |
545 | itsLoader->LoadTracks("read"); |
546 | itsLoader->LoadV0s("recreate"); |
547 | AliV0vertexer vertexer; |
548 | TTree* tracks = itsLoader->TreeT(); |
549 | if (!tracks) { |
550 | Error("ReconstructV0s", "no ITS tracks tree found"); |
551 | return kFALSE; |
552 | } |
553 | if (!itsLoader->TreeV0()) itsLoader->MakeTree("V0"); |
554 | TTree* v0s = itsLoader->TreeV0(); |
555 | vertexer.Tracks2V0vertices(tracks, v0s); |
556 | itsLoader->WriteV0s("OVERWRITE"); |
557 | |
558 | itsLoader->UnloadTracks(); |
559 | itsLoader->UnloadV0s(); |
560 | return kTRUE; |
561 | } |
562 | |
97d6eb66 |
563 | //_____________________________________________________________________________ |
564 | #ifdef ALI_HLT |
565 | void AliMonitorProcess::CreateHLT(const char* fileName) |
566 | { |
1899848d |
567 | |
97d6eb66 |
568 | // create the HLT (Level3) object |
569 | |
570 | if (fHLT) delete fHLT; |
571 | |
572 | char name[256]; |
573 | strcpy(name, fileName); |
574 | fHLT = new AliLevel3(name); |
575 | fHLT->Init("./", AliLevel3::kRaw, 1); |
576 | |
577 | fHLT->SetClusterFinderParam(0, 0, kTRUE); |
578 | |
c4bd737c |
579 | Int_t phiSegments = 50; |
580 | Int_t etaSegments = 100; |
97d6eb66 |
581 | Int_t trackletlength = 3; |
1899848d |
582 | Int_t tracklength = 40;//40 or 5 |
97d6eb66 |
583 | Int_t rowscopetracklet = 2; |
584 | Int_t rowscopetrack = 2; |
c4bd737c |
585 | Double_t minPtFit = 0; |
97d6eb66 |
586 | Double_t maxangle = 1.31; |
587 | Double_t goodDist = 5; |
588 | Double_t maxphi = 100; |
589 | Double_t maxeta = 100; |
1899848d |
590 | Double_t hitChi2Cut = 15;//100 or 15 |
591 | Double_t goodHitChi2 = 5;//20 or 5 |
592 | Double_t trackChi2Cut = 10;//50 or 10 |
c4bd737c |
593 | fHLT->SetTrackerParam(phiSegments, etaSegments, |
97d6eb66 |
594 | trackletlength, tracklength, |
595 | rowscopetracklet, rowscopetrack, |
c4bd737c |
596 | minPtFit, maxangle, goodDist, hitChi2Cut, |
97d6eb66 |
597 | goodHitChi2, trackChi2Cut, 50, maxphi, maxeta, kTRUE); |
598 | |
599 | fHLT->WriteFiles("./hlt/"); |
600 | } |
601 | #endif |
602 | |
603 | //_____________________________________________________________________________ |
54a41755 |
604 | Bool_t AliMonitorProcess::ReconstructHLT( |
605 | #ifdef ALI_HLT |
606 | Int_t iEvent |
607 | #else |
9abb4463 |
608 | Int_t /* iEvent */ |
54a41755 |
609 | #endif |
610 | ) |
97d6eb66 |
611 | { |
612 | // run the HLT cluster and track finder |
613 | |
1899848d |
614 | SetStatus(kRecHLT); |
97d6eb66 |
615 | |
616 | #ifndef ALI_HLT |
617 | Warning("ReconstructHLT", "the code was compiled without HLT support"); |
618 | return kTRUE; |
619 | |
620 | #else |
621 | gSystem->Exec("rm -rf hlt"); |
622 | gSystem->MakeDirectory("hlt"); |
623 | if (!fHLT) return kFALSE; |
624 | |
625 | fHLT->ProcessEvent(0, 35, iEvent); |
626 | |
627 | // remove the event number from the file names |
628 | char command[256]; |
629 | sprintf(command, "rename points_%d points hlt/*.raw", iEvent); |
630 | gSystem->Exec(command); |
631 | sprintf(command, "rename tracks_tr_%d tracks_tr hlt/*.raw", iEvent); |
632 | gSystem->Exec(command); |
633 | sprintf(command, "rename tracks_gl_%d tracks_gl hlt/*.raw", iEvent); |
634 | gSystem->Exec(command); |
635 | sprintf(command, "rename tracks_%d tracks hlt/*.raw", iEvent); |
636 | gSystem->Exec(command); |
637 | return kTRUE; |
638 | #endif |
639 | } |
640 | |
04fa961a |
641 | |
642 | //_____________________________________________________________________________ |
643 | Bool_t AliMonitorProcess::WriteHistos() |
644 | { |
645 | // write the monitor tree and the monitor histograms to the file |
646 | // "monitor_<run number>[_<sub_run_number>].root" |
647 | // if at least fNEventsMin events were monitored |
648 | |
1899848d |
649 | SetStatus(kWriting); |
04fa961a |
650 | |
651 | // rename tree file and create a new one |
652 | fFile->cd(); |
653 | fTree->Write(); |
654 | fFile->Close(); |
655 | delete fFile; |
656 | |
657 | char fileName[256]; |
658 | sprintf(fileName, "monitor_tree_%d.root", fRunNumber); |
659 | if (fSubRunNumber > 0) { |
660 | sprintf(fileName, "monitor_tree_%d_%d.root", fRunNumber, fSubRunNumber); |
661 | } |
662 | if (fNEvents < fNEventsMin) { |
663 | gSystem->Unlink("monitor_tree.root"); |
664 | } else { |
665 | gSystem->Rename("monitor_tree.root", fileName); |
666 | } |
667 | |
668 | fFile = TFile::Open("monitor_tree.root", "RECREATE"); |
669 | if (!fFile || !fFile->IsOpen()) { |
670 | Fatal("WriteHistos", "could not open file for tree"); |
671 | } |
672 | fTree = new TTree("MonitorTree", "tree for monitoring"); |
673 | for (Int_t iMonitor = 0; iMonitor < fMonitors.GetEntriesFast(); iMonitor++) { |
674 | ((AliMonitor*) fMonitors[iMonitor])->CreateBranches(fTree); |
675 | } |
676 | gROOT->cd(); |
677 | |
678 | // write the histograms |
679 | if (fNEvents < fNEventsMin) return kTRUE; |
680 | |
681 | if (!fWriteHistoList) { |
682 | TIterator* iFolder = fTopFolder->GetListOfFolders()->MakeIterator(); |
683 | while (TFolder* folder = (TFolder*) iFolder->Next()) { |
684 | TIterator* iHisto = folder->GetListOfFolders()->MakeIterator(); |
685 | while (AliMonitorPlot* histo = (AliMonitorPlot*) iHisto->Next()) { |
686 | histo->ResetList(); |
687 | } |
688 | delete iHisto; |
689 | } |
690 | delete iFolder; |
691 | } |
692 | |
693 | Bool_t result = kTRUE; |
694 | sprintf(fileName, "monitor_%d.root", fRunNumber); |
695 | if (fSubRunNumber > 0) { |
696 | sprintf(fileName, "monitor_%d_%d.root", fRunNumber, fSubRunNumber); |
697 | } |
698 | TFile* file = TFile::Open(fileName, "recreate"); |
699 | if (!file || !file->IsOpen()) { |
700 | Error("WriteHistos", "could not open file %s", fileName); |
701 | result = kFALSE; |
702 | } else { |
703 | fTopFolder->Write(); |
704 | file->Close(); |
705 | } |
706 | if (file) delete file; |
707 | |
708 | return result; |
709 | } |
710 | |
711 | //_____________________________________________________________________________ |
712 | void AliMonitorProcess::StartNewRun() |
713 | { |
714 | // reset the histograms for a new run |
715 | |
1899848d |
716 | SetStatus(kResetting); |
04fa961a |
717 | TIterator* iFolder = fTopFolder->GetListOfFolders()->MakeIterator(); |
718 | while (TFolder* folder = (TFolder*) iFolder->Next()) { |
719 | TIterator* iHisto = folder->GetListOfFolders()->MakeIterator(); |
720 | while (AliMonitorPlot* histo = (AliMonitorPlot*) iHisto->Next()) { |
721 | histo->Reset(); |
722 | } |
723 | delete iHisto; |
724 | } |
725 | delete iFolder; |
726 | |
727 | fNEvents = 0; |
728 | } |
729 | |
730 | |
731 | //_____________________________________________________________________________ |
732 | void AliMonitorProcess::CheckForConnections() |
733 | { |
734 | // check if new clients want to connect and add them to the list of sockets |
735 | |
736 | TMessage message(kMESS_OBJECT); |
737 | message.WriteObject(fTopFolder); |
1899848d |
738 | SetStatus(kConnecting); |
04fa961a |
739 | |
740 | TSocket* socket; |
741 | while ((socket = fServerSocket->Accept()) != (TSocket*)-1) { |
54a41755 |
742 | socket->SetOption(kNoBlock, 1); |
04fa961a |
743 | char socketType[256]; |
744 | if (!socket->Recv(socketType, 255)) continue; |
745 | if (strcmp(socketType, "client") == 0) { |
746 | if ((fNEvents == 0) || (socket->Send(message) >= 0)) { |
747 | fSockets.Add(socket); |
748 | TInetAddress adr = socket->GetInetAddress(); |
749 | Info("CheckForConnections", "new client:\n %s (%s), port %d\n", |
750 | adr.GetHostName(), adr.GetHostAddress(), adr.GetPort()); |
751 | } |
752 | } else if (strcmp(socketType, "display") == 0) { |
753 | if (fDisplaySocket) { |
754 | fDisplaySocket->Close(); |
755 | delete fDisplaySocket; |
756 | } |
757 | fDisplaySocket = socket; |
758 | fDisplaySocket->SetOption(kNoBlock, 1); |
759 | TInetAddress adr = socket->GetInetAddress(); |
760 | Info("CheckForConnections", "new display:\n %s (%s), port %d\n", |
761 | adr.GetHostName(), adr.GetHostAddress(), adr.GetPort()); |
762 | } |
763 | } |
764 | |
765 | for (Int_t iSocket = 0; iSocket < fSockets.GetEntriesFast(); iSocket++) { |
766 | socket = (TSocket*) fSockets[iSocket]; |
767 | if (!socket) continue; |
1899848d |
768 | // remove finished client |
769 | char str[256]; |
770 | if (socket->Recv(str, 255)) { |
771 | TString socketMessage(str); |
772 | if(socketMessage.CompareTo("Finished") == 0) { |
773 | TInetAddress adr = socket->GetInetAddress(); |
774 | Info("CheckForConnections", |
775 | "disconnect finished client:\n %s (%s), port %d\n", |
776 | adr.GetHostName(), adr.GetHostAddress(), adr.GetPort()); |
777 | delete fSockets.RemoveAt(iSocket); |
778 | continue; |
779 | } |
780 | } |
04fa961a |
781 | if (!socket->IsValid()) { |
782 | // remove invalid sockets from the list |
783 | TInetAddress adr = socket->GetInetAddress(); |
784 | Info("BroadcastHistos", "disconnect client:\n %s (%s), port %d\n", |
785 | adr.GetHostName(), adr.GetHostAddress(), adr.GetPort()); |
786 | delete fSockets.RemoveAt(iSocket); |
787 | } |
788 | } |
789 | fSockets.Compress(); |
790 | } |
791 | |
792 | //_____________________________________________________________________________ |
793 | void AliMonitorProcess::BroadcastHistos() |
794 | { |
795 | // send the monitor histograms to the clients |
796 | |
1899848d |
797 | SetStatus(kBroadcasting); |
04fa961a |
798 | TMessage message(kMESS_OBJECT); |
799 | message.WriteObject(fTopFolder); |
800 | |
801 | for (Int_t iSocket = 0; iSocket < fSockets.GetEntriesFast(); iSocket++) { |
802 | TSocket* socket = (TSocket*) fSockets[iSocket]; |
803 | if (!socket) continue; |
54a41755 |
804 | socket->SetOption(kNoBlock, 0); |
04fa961a |
805 | if (!socket->IsValid() || (socket->Send(message) < 0)) { |
806 | // remove the socket from the list if there was an error |
807 | TInetAddress adr = socket->GetInetAddress(); |
808 | Info("BroadcastHistos", "disconnect client:\n %s (%s), port %d\n", |
809 | adr.GetHostName(), adr.GetHostAddress(), adr.GetPort()); |
810 | delete fSockets.RemoveAt(iSocket); |
54a41755 |
811 | } else { |
812 | socket->SetOption(kNoBlock, 1); |
04fa961a |
813 | } |
814 | } |
815 | fSockets.Compress(); |
816 | } |