Dummy default constructor
[u/mrichter/AliRoot.git] / MONITOR / AliMonitorProcess.cxx
CommitLineData
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
56ClassImp(AliMonitorProcess)
57
58
c4bd737c 59const Int_t AliMonitorProcess::fgkPort = 9327;
04fa961a 60
61
62//_____________________________________________________________________________
63AliMonitorProcess::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//_____________________________________________________________________________
150AliMonitorProcess::AliMonitorProcess(const AliMonitorProcess& process) :
151 TObject(process)
152{
153 Fatal("AliMonitorProcess", "copy constructor not implemented");
154}
155
156//_____________________________________________________________________________
157AliMonitorProcess& AliMonitorProcess::operator = (const AliMonitorProcess&
158 /*process*/)
159{
160 Fatal("operator =", "assignment operator not implemented");
161 return *this;
162}
163
04fa961a 164//_____________________________________________________________________________
165AliMonitorProcess::~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//_____________________________________________________________________________
191const char* AliMonitorProcess::GetRevision()
192{
193 return "$Revision$";
194}
195
196
1899848d 197//_____________________________________________________________________________
198void 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//_____________________________________________________________________________
208void 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//_____________________________________________________________________________
236void 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//_____________________________________________________________________________
246void 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//_____________________________________________________________________________
267Bool_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//_____________________________________________________________________________
306Bool_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//_____________________________________________________________________________
387void 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 399UInt_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 407UInt_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 415UInt_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 423Int_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//_____________________________________________________________________________
449Bool_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//_____________________________________________________________________________
486Bool_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//_____________________________________________________________________________
530Bool_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
565void 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 604Bool_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//_____________________________________________________________________________
643Bool_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//_____________________________________________________________________________
712void 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//_____________________________________________________________________________
732void 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//_____________________________________________________________________________
793void 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}