]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PHOS/AliPHOSOnlineMonitor.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Class for online monitoring of BeamTests
[u/mrichter/AliRoot.git] / PHOS / AliPHOSOnlineMonitor.cxx
CommitLineData
f4374ddc 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// Class intended to perform online monitoring of PHOS beamtests
20// Being constructed, produces menu with list of available histograms to fill
21// Once histograms are selected, button "Go" should be pressed to start scan of data.
22// Prepared histograms will be periodically updated during scan of the data.
23// Note:
24// 1. To plot most of histograms, a "Connection Table", relating ADC signal index and AbsId
25// of PHOS crystal, should be created beforehand. To do this, call macro
26// $ALICE_ROOT/PHOS/macros/BeamTest/MakeConTableDB.C
27// with apropriate number of raws and columns of prototype.
28// 2. To perform reconstruction and e.g. invariant mass analysis, a "Calibration Database"
29// should be created beforehand. To do this, call macro
30// $ALICE_ROOT/PHOS/macros/BeamTest/MakeConTableDB.C
31// to read calibration parameters from file or use AliPHOSCalibrator to calculate
32// pedestal and gains.
33// 3. Once histograms are filled with "Go" method, they can be written to file
34// with WriteHisto("Filename.root") method.
35//
36//*-- Author : D.Peressounko (RRC KI) after A.V. Kuryakin, (Sarov)
37//////////////////////////////////////////////////////////////////////////////
38
39// --- ROOT system ---
40#include "TROOT.h"
41#include "TClonesArray.h"
42#include "TCanvas.h"
43#include "TH1.h"
44#include "TH2.h"
45#include "TGroupButton.h"
46#include "TFile.h"
47#include "TSystem.h"
48
49
50// --- Standard library ---
51#include "TBenchmark.h"
52#include "Riostream.h"
53
54// --- AliRoot header files ---
55#include "AliPHOSOnlineMonitor.h"
56#include "AliPHOSConTableDB.h"
57#include "AliPHOSGeometry.h"
58#include "AliPHOSRawReaderDate.h"
59#include "AliPHOSRawStream.h"
60#include "AliPHOSDigit.h"
61#include "AliPHOSGetterLight.h"
62#include "AliPHOSClusterizerv1.h"
63#include "AliPHOSTrackSegmentMakerv1.h"
64#include "AliPHOSPIDv1.h"
65#include "AliPHOSCalibrManager.h"
66#include "AliPHOSCalibrationDB.h"
67
68ClassImp(AliPHOSOnlineMonitor)
69
70
71//____________________________________________________________________________
72AliPHOSOnlineMonitor::AliPHOSOnlineMonitor(): TDialogCanvas("PHOS","PHOS",150,300)
73{
74
75 MakeButtons() ;
76 Modified(kTRUE);
77 Update();
78 SetEditable(kFALSE);
79
80 //add this TFitPanel to the list of cleanups such that in case
81 //the referenced object is deleted, its pointer be reset
82 gROOT->GetListOfCleanups()->Add(this);
83
84 fRefObject = this;
85 fRefPad = (TPad*)gROOT->GetSelectedPad();
86
87 fScanPed = 0 ;
88 fScanSig = 0;
89 fReconstruct = 0 ;
90
91 //Prepare list of canvas with histograms
92 fCanvasList = new TList() ;
93 fHistosList = new TList() ;
94
95 fcdb = 0 ;
96 fGeom = AliPHOSGeometry::GetInstance("IHEP","") ;
97 fNUpdate = 1000 ;
98
99}
100//____________________________________________________________________________
101AliPHOSOnlineMonitor::AliPHOSOnlineMonitor(const char * inputfile): TDialogCanvas("PHOS","PHOS",150,300)
102{
103
104 MakeButtons() ;
105 Modified(kTRUE);
106 Update();
107 SetEditable(kFALSE);
108
109 //add this TFitPanel to the list of cleanups such that in case
110 //the referenced object is deleted, its pointer be reset
111 gROOT->GetListOfCleanups()->Add(this);
112
113 fRefObject = this;
114 fRefPad = (TPad*)gROOT->GetSelectedPad();
115
116 fScanPed = 0 ;
117 fScanSig = 0;
118 fReconstruct = 0 ;
119
120 //Prepare list of canvas with histograms
121 fCanvasList = new TList() ;
122 fHistosList = new TList() ;
123
124 fInputFile = inputfile ;
125 fGeom = AliPHOSGeometry::GetInstance("IHEP","") ;
126 fcdb = 0 ;
127 fNUpdate = 1000 ;
128
129}
130//____________________________________________________________________________
131 AliPHOSOnlineMonitor::~AliPHOSOnlineMonitor()
132{
133 //Obvious, but unevoidable comment for destructor: cleans up everething.
134 TIter nextCanvas(fCanvasList);
135 TCanvas * c ;
136 while((c=(TCanvas*)nextCanvas()))
137 delete c ;
138 delete fCanvasList ;
139
140 TIter nextHisto(fHistosList);
141 TH1D * h ;
142 while((h=(TH1D*)nextHisto()))
143 delete h ;
144 delete fHistosList ;
145
146 if(fcdb)
147 delete fcdb ;
148
149}
150//____________________________________________________________________________
151void AliPHOSOnlineMonitor::MakeButtons(void){
152 //Make buttons on graphical menu
153 Int_t nButtons = 16;
154 TGroupButton * b ;
155 Float_t xmin = 0.0;
156 Float_t ymin = 0.01;
157 Float_t xmax = 0.99;
158 Float_t ymax = 0.99;
159 Float_t dy = (ymax-ymin)/nButtons ;
160
161 Float_t y2=ymax ;
162 Float_t y1=y2-dy ;
163 b = new TGroupButton("APPLY","Triggers","",xmin,y1,xmax,y2);
164 b->Draw();
165 y2=y1 ;
166 y1=y1-dy ;
167 b = new TGroupButton("APPLY","Pedestals","",xmin,y1,xmax,y2);
168 b->Draw();
169 y2=y1 ;
170 y1=y1-dy ;
171 b = new TGroupButton("APPLY","Spectrum all","",xmin,y1,xmax,y2);
172 b->Draw();
173 y2=y1 ;
174 y1=y1-dy ;
175 b = new TGroupButton("APPLY","Spectrum g","",xmin,y1,xmax,y2);
176 b->Draw();
177 y2=y1 ;
178 y1=y1-dy ;
179 b = new TGroupButton("APPLY","Inv Mass","",xmin,y1,xmax,y2);
180 b->Draw();
181 for(Int_t i=1; i<=5; i++){
182 y2=y1 ;
183 y1=y1-dy ;
184 char name[10] ;
185 sprintf(name,"Edep(ADC) %d",i) ;
186 b = new TGroupButton("APPLY",name,"",xmin,y1,xmax,y2);
187 b->Draw();
188 }
189 for(Int_t i=1; i<=5; i++){
190 y2=y1 ;
191 y1=y1-dy ;
192 char name[10] ;
193 sprintf(name,"Edep(Cal) %d",i) ;
194 b = new TGroupButton("APPLY",name,"",xmin,y1,xmax,y2);
195 b->Draw();
196 }
197 y2=y1 ;
198 y1=y1-dy ;
199 b = new TGroupButton("APPLY","Go","",xmin,y1,xmax,y2);
200 b->SetTextColor(2);
201 b->Draw();
202}
203//____________________________________________________________________________
204void AliPHOSOnlineMonitor::Apply(const char *action){
205 //Function to handle button actions
206
207 TDialogCanvas::Apply() ;
208
209 TObject *obj;
210 TGroupButton *button;
211 TIter next(fPrimitives);
212
213 if (!strcmp(action,"Triggers")) {
214 DrawTriggers() ;
215 }
216
217 if (!strcmp(action,"Pedestals")) {
218 DrawPedestals() ;
219 }
220 if (!strcmp(action,"Spectrum all")) {
221 DrawSpectrum("all") ;
222 }
223 if (!strcmp(action,"Spectrum g")) {
224 DrawSpectrum("gamma") ;
225 }
226 if (!strcmp(action,"Inv Mass")) {
227 DrawMinv() ;
228 }
229 if(strstr(action,"Edep")){
230 Int_t n ;
231 char tmp[10] ;
232 sscanf(action,"%s %d",tmp,&n) ;
233 char opt[5]="" ;
234 if(strstr(action,"Cal"))
235 sprintf(opt,"Cal") ;
236 DrawEdep(n,opt) ;
237 }
238 if (!strcmp(action,"Go")) {
239 Go() ;
240 }
241
242 //Mark button as pressed
243 if(strcmp(action,"Go")){ //Do not mark "Go" button
244 while ((obj = next())) {
245 if (obj->InheritsFrom(TGroupButton::Class())) {
246 button = (TGroupButton*)obj;
247 if(!strcmp(button->GetTitle(),action)){
248 if (button->GetBorderMode() > 0){
249 button->SetBorderMode(-1) ;
250 button->Modified(kTRUE);
251 }
252 }
253 }
254 }
255 }
256}
257//____________________________________________________________________________
258void AliPHOSOnlineMonitor::SetInputFile(const char * filename){
259 //close previously opened
260
261 fInputFile = filename ;
262}
263//____________________________________________________________________________
264void AliPHOSOnlineMonitor::DrawPedestals(){
265 //Prepare canvas and histograms for drawing pedestals
266
267 TIter nextCanvas(fCanvasList);
268 TCanvas * c ;
269 Bool_t exists = kFALSE ;
270 while((c=(TCanvas*)nextCanvas())){
271 if(!strcmp(c->GetName(),"Pedestals")){
272 exists = kTRUE ;
273 break;
274 }
275 }
276 if(!exists){
277 c = new TDialogCanvas("Pedestals","Pedestals",300,200) ;
278 fCanvasList->AddLast(c) ;
279 }
280
281 TIter nextHisto(fHistosList);
282 TH1D * h ;
283 exists = kFALSE ;
284 while((h=(TH1D*)nextHisto())){
285 if(!strcmp(h->GetName(),"hPedestals")){
286 exists = kTRUE ;
287 break;
288 }
289 }
290 if(!exists){
291 h = new TH1D("hPedestals","Pedestals per event",fGeom->GetNModules()*fGeom->GetNCristalsInModule(),0.,
292 1.*fGeom->GetNModules()*fGeom->GetNCristalsInModule()) ;
293 fHistosList->AddLast(h) ;
294 }
295
296 c->cd() ;
297 h->SetStats(0) ;
298 h->Draw() ;
299 fScanPed = kTRUE ; //We will scan pedestals
300
301}
302//____________________________________________________________________________
303void AliPHOSOnlineMonitor::DrawTriggers(){
304 //Prepare canvas and histogram for drawing triggers
305
306 TIter nextCanvas(fCanvasList);
307 TCanvas * c ;
308 Bool_t exists = kFALSE ;
309 while((c=(TCanvas*)nextCanvas())){
310 if(!strcmp(c->GetName(),"Triggers")){
311 exists = kTRUE ;
312 break;
313 }
314 }
315 if(!exists){
316 c = new TDialogCanvas("Triggers","Triggers",200,200) ;
317 fCanvasList->AddLast(c) ;
318 }
319
320 TIter nextHisto(fHistosList);
321 TH1D * h ;
322 exists = kFALSE ;
323 while((h=(TH1D*)nextHisto())){
324 if(!strcmp(h->GetName(),"hTriggers")){
325 exists = kTRUE ;
326 break;
327 }
328 }
329 if(!exists){
330 h = new TH1D("hTriggers","Triggers",2,0.,2.) ;
331 fHistosList->AddLast(h) ;
332 }
333 //Make Labels
334 h->SetBit(TH1::kCanRebin);
335 h->Fill("LED",0.0000001) ;
336 h->Fill("PUL",0.0000001) ;
337 h->Fill("PED",0.0000001) ;
338 h->Fill("NEL",0.0000001) ;
339 h->Fill("WEL",0.0000001) ;
340 h->Fill("SOB",0.0000001) ;
341 h->Fill("EOB",0.0000001) ;
342 h->Fill("wrong",0.0000001) ;
343 h->LabelsOption("h");
344 h->LabelsDeflate();
345 h->SetStats(0) ;
346 c->cd() ;
347 h->Draw() ;
348}
349//____________________________________________________________________________
350void AliPHOSOnlineMonitor::DrawSpectrum(const char * opt){
351 //Prepare canvas and histograms for drawing spectra of all reconstructed particles or photons
352
353 TString name("Spectrum") ;
354 name+=opt ;
355
356 TIter nextCanvas(fCanvasList);
357 TCanvas * c ;
358 Bool_t exists = kFALSE ;
359 while((c=(TCanvas*)nextCanvas())){
360 if(!strcmp(c->GetName(),name.Data())){
361 exists = kTRUE ;
362 break;
363 }
364 }
365 if(!exists){
366 c = new TDialogCanvas(name,name,250,300) ;
367 fCanvasList->AddLast(c) ;
368 }
369
370 TIter nextHisto(fHistosList);
371 TH1D * h ;
372 exists = kFALSE ;
373 name.Prepend("h") ;
374 while((h=(TH1D*)nextHisto())){
375 if(!strcmp(h->GetName(),name.Data())){
376 exists = kTRUE ;
377 break;
378 }
379 }
380 if(!exists){
381 h = new TH1D(name,name,100,0.,100.) ;
382 fHistosList->AddLast(h) ;
383 }
384
385 h->SetStats(0) ;
386 c->cd() ;
387 h->Draw() ;
388 fReconstruct = kTRUE ;
389 fScanSig = kTRUE ; //We will scan pedestals
390
391}
392//____________________________________________________________________________
393void AliPHOSOnlineMonitor::DrawMinv(){
394 TIter nextCanvas(fCanvasList);
395 TCanvas * c ;
396 Bool_t exists = kFALSE ;
397 while((c=(TCanvas*)nextCanvas())){
398 if(!strcmp(c->GetName(),"InvMass")){
399 exists = kTRUE ;
400 break;
401 }
402 }
403 if(!exists){
404 c = new TDialogCanvas("InvMass","Invariant mass",300,200) ;
405 fCanvasList->AddLast(c) ;
406 }
407
408 TIter nextHisto(fHistosList);
409 TH1D * h ;
410 exists = kFALSE ;
411 while((h=(TH1D*)nextHisto())){
412 if(!strcmp(h->GetName(),"hInvMass")){
413 exists = kTRUE ;
414 break;
415 }
416 }
417 if(!exists){
418 h = new TH1D("hInvMass","hInvMass",1000,0.,1.0) ;
419 fHistosList->AddLast(h) ;
420 }
421
422 c->cd() ;
423 h->Draw() ;
424 h->SetStats(0) ;
425 fReconstruct = kTRUE ;
426 fScanSig = kTRUE ; //We will scan pedestals
427}
428//____________________________________________________________________________
429void AliPHOSOnlineMonitor::DrawEdep(Int_t mod,const char * opt){
430 char name[15] ;
431 sprintf(name,"Edep%s %d",opt,mod) ;
432
433 TIter nextCanvas(fCanvasList);
434 TCanvas * c ;
435 Bool_t exists = kFALSE ;
436 while((c=(TCanvas*)nextCanvas())){
437 if(!strcmp(c->GetName(),name)){
438 exists = kTRUE ;
439 break;
440 }
441 }
442 if(!exists){
443 c = new TDialogCanvas(name,name,300,200) ;
444 fCanvasList->AddLast(c) ;
445 }
446
447 TIter nextHisto(fHistosList);
448 TH2D * h ;
449 exists = kFALSE ;
450 sprintf(name,"hEdep%s%d",opt,mod) ;
451 while((h=(TH2D*)nextHisto())){
452 if(!strcmp(h->GetName(),name)){
453 exists = kTRUE ;
454 break;
455 }
456 }
457 if(!exists){
458 h = new TH2D(name,name,fGeom->GetNPhi(),0.,1.*fGeom->GetNPhi(),fGeom->GetNZ(),0.,1.*fGeom->GetNZ()) ;
459 fHistosList->AddLast(h) ;
460 }
461
462 c->cd() ;
463 h->Draw("col") ;
464 h->SetStats(0) ;
465 fScanSig = kTRUE ; //We will scan signal events
466 if(strstr(opt,"Cal"))
467 fReconstruct = kTRUE ;
468
469}
470//____________________________________________________________________________
471void AliPHOSOnlineMonitor::ScanPedestals(TClonesArray * digits){
472 //This method is called for events with PED trigger
473 //We fill bins with ADC values
474
475 TH1D * h = (TH1D*)gROOT->FindObjectAny("hPedestals");
476 if(!h){
477 Error("ScanPedestals","Can not fild histogram hPedestals") ;
478 return ;
479 }
480 for(Int_t i=0; i<digits->GetEntriesFast(); i++){
481 AliPHOSDigit * dig = static_cast<AliPHOSDigit*>(digits->At(i)) ;
482 h->AddBinContent(dig->GetId(),dig->GetAmp()) ;
483 }
484}
485//____________________________________________________________________________
486void AliPHOSOnlineMonitor::ScanEdep(TClonesArray * digits){
487 //Fill 2D distribution of ADC values in NEL and WEL events
488 AliPHOSGetter * gime = AliPHOSGetter::Instance() ;
489 AliPHOSCalibrationDB *cdb = 0 ;
490 if(gime)
491 cdb = gime->CalibrationDB() ;
492 Int_t mod = 0 ;
493 char name[15] ;
494 TH2D * h = 0 ;
495 TH2D * hCal = 0 ;
496 for(Int_t i=0; i<digits->GetEntriesFast(); i++){
497 AliPHOSDigit * dig = static_cast<AliPHOSDigit*>(digits->At(i)) ;
498 Int_t relId[4] ;
499 fGeom->AbsToRelNumbering(dig->GetId(),relId) ;
500 if(mod != relId[0]){ //new module, look for histograms
501 mod = relId[0] ;
502 sprintf(name,"hEdep%d",mod) ;
503 h = (TH2D*)gROOT->FindObjectAny(name);
504 sprintf(name,"hEdepCal%d",mod) ;
505 hCal = (TH2D*)gROOT->FindObjectAny(name);
506 }
507 if(h)
508 h->Fill(relId[2]-0.1,relId[3]-0.1,1.*dig->GetAmp()) ;
509 if(hCal)
510 hCal->Fill(relId[2]-0.1,relId[3]-0.1,cdb->Calibrate(dig->GetAmp(),dig->GetId())) ;
511 }
512}
513//____________________________________________________________________________
514void AliPHOSOnlineMonitor::ScanRecon(TClonesArray * recParticles){
515 if(!recParticles || recParticles->GetEntries()==0) return ;
516
517 TH1D* hSpectr = (TH1D*)gROOT->FindObjectAny("hSpectrumall");
518 TH1D* hSpectrGam= (TH1D*)gROOT->FindObjectAny("hSpectrumgamma");
519 TH1D* hInvMass = (TH1D*)gROOT->FindObjectAny("hInvMass");
520 for(Int_t i=0; i<recParticles->GetEntriesFast() ; i++){
521 AliPHOSRecParticle * p = (AliPHOSRecParticle *)recParticles->At(i) ;
522 if(hSpectr)hSpectr->Fill(p->Energy()) ;
523 if(hSpectrGam && p->IsPhoton())hSpectrGam->Fill(p->Energy()) ;
524 if(hInvMass){
525 for(Int_t j=i+1; j<recParticles->GetEntriesFast() ; j++){
526 AliPHOSRecParticle * p2 = (AliPHOSRecParticle *)recParticles->At(j) ;
527 Double_t e = p->Energy() + p2->Energy() ;
528 Double_t x = p->Px() + p2->Px() ;
529 Double_t y = p->Py() + p2->Py() ;
530 Double_t z = p->Pz() + p2->Pz() ;
531 Double_t m = e*e-x*x-y*y-z*z ;
532 hInvMass->Fill(m>0?TMath::Sqrt(m): 0. ) ;
533 }
534 }
535 }
536}
537//____________________________________________________________________________
538void AliPHOSOnlineMonitor::ScanTrigger(Int_t trig){
539 //Fills trigger distribution
540
541 TH1D * h = (TH1D*)gROOT->FindObjectAny("hTriggers");
542 if(!h) return ;
543 switch(trig){
544 case AliPHOSRawStream::kLED : h->Fill("LED",1.) ; break ;
545 case AliPHOSRawStream::kPUL : h->Fill("PUL",1.) ; break ;
546 case AliPHOSRawStream::kPED : h->Fill("PED",1.) ; break ;
547 case AliPHOSRawStream::kNEL : h->Fill("NEL",1.) ; break ;
548 case AliPHOSRawStream::kWEL : h->Fill("WEL",1.) ; break ;
549 case AliPHOSRawStream::kSOB : h->Fill("SOB",1.) ; break ;
550 case AliPHOSRawStream::kEOB : h->Fill("EOB",1.) ; break ;
551 default : h->Fill("wrong",1.) ;
552 }
553}
554//____________________________________________________________________________
555void AliPHOSOnlineMonitor::SetConTableDB(const char * filename){
556 //Read ConnectionTableDB from file
557 TFile * file = new TFile(filename) ;
558 AliPHOSConTableDB * tmp = (AliPHOSConTableDB*)file->Get("AliPHOSConTableDB") ;
559 fcdb = new AliPHOSConTableDB(*tmp) ;
560 file->Close() ;
561}
562//____________________________________________________________________________
563void AliPHOSOnlineMonitor::Go(){
564 //Perform scan of curent event
565 gBenchmark->Start("PHOSOnlineMon");
566
567 //First test if we need "Connection table" then open it
568 if(!fcdb){
569 SetConTableDB() ;
570 if(fcdb){
571 Info("Go","Read Connection table from file \"ConTableDB.root\"") ;
572 }else{
573 Error("Go","Please, set connection table with SetConTableDB() method") ;
574 return ;
575 }
576 }
577
578 AliPHOSGetterLight * gime = AliPHOSGetterLight::Instance("PHOS","On Flight") ;
579
580 //Configure CalibrManager to read data from file
581 //Create calibration database and read it
582 AliPHOSCalibrationDB * calibDB = 0 ;
583 if(fScanSig || fReconstruct){ //We will ned calibration parameters
584 AliPHOSCalibrManager::GetInstance("CalibrDB.root","root") ;
585 //If we configured manager to read from ASCII file,
586 //give him connection table. OK, it will not harm in any case.
587 AliPHOSCalibrManager::GetInstance()->SetConTable(fcdb) ;
588
589 calibDB = new AliPHOSCalibrationDB("OnLine") ;
590 calibDB->GetParameters() ; //Read parameters using Manager
591 gime->SetCalibrationDB(calibDB) ;
592 }
593
594 //Now open data file
595 AliPHOSRawReaderDate *rawReader = new AliPHOSRawReaderDate(fInputFile) ;
596 AliPHOSRawStream *rawStream = new AliPHOSRawStream(rawReader) ;
597 rawStream->SetConTableDB(fcdb) ;
598
599 TClonesArray * digits = gime->Digits() ;
600 TClonesArray * recParticles = gime->RecParticles() ;
601 AliPHOSClusterizerv1* clu = 0 ;
602 AliPHOSTrackSegmentMakerv1 * tsm = 0 ;
603 AliPHOSPIDv1 * pid = 0 ;
604 if(fReconstruct){ //We will need calibation parameters
605 clu = new AliPHOSClusterizerv1("PHOS","On Flight") ;
606 clu->SetWriting(0) ; //Do not write to file
607 clu->SetEmcMinE(0.05) ; //Minimal energy of the digit
608 clu->SetEmcLocalMaxCut(0.05) ; //Height of local maximum over environment
609 clu->SetEmcClusteringThreshold(0.2) ; //Minimal energy to start cluster
610// clu->SetUnfolding(kFALSE) ; //Do not unfold
611 tsm = new AliPHOSTrackSegmentMakerv1("PHOS","On Flight") ;
612 tsm->SetWriting(0) ; //Do not write to file
613 pid = new AliPHOSPIDv1("PHOS","On Flight") ;
614 pid->SetWriting(0) ; //Do not write to file
615 }
616
617 fNevents=0 ;
618 //Scan all event in file
619 printf(" ") ;
620 while(rawReader->NextEvent()){
621 //Is it PHYSICAL event
622 if(rawReader->GetType() == PHYSICS_EVENT){
623 fNevents++ ;
624 if(fNevents%100 ==0){
625 printf("\b\b\b\b\b\b%6d",fNevents) ;
626 }
627 if(rawStream->ReadDigits(digits)){
628
629 //Test trigger
630 //Pedestal Event
631 ScanTrigger(rawStream->GetTrigger()) ;
632 if(rawStream->IsPEDevent() && fScanPed){
633 ScanPedestals(digits) ;
634 }
635 if((rawStream->IsNELevent() || rawStream->IsWELevent()) && fScanSig){
636 ScanEdep(digits) ;
637 if(fReconstruct){
638 gime->Clusterizer()->Exec("") ;
639 gime->TrackSegmentMaker()->Exec("") ;
640 gime->PID()->Exec("") ;
641 ScanRecon(recParticles) ;
642 }
643 }
644 }
645
646 if(fNevents%fNUpdate == 0 ){ //upate all histograms
647 TIter nextCanvas(fCanvasList);
648 TCanvas * c ;
649 while((c=(TCanvas*)nextCanvas())){
650 c->Modified() ;
651 c->Update() ;
652 }
653 }
654 gSystem->ProcessEvents();
655 }
656 // if(fNevents>=200)break ;
657 }
658 printf("\n") ;
659 gBenchmark->Stop("PHOSOnlineMon");
660 Float_t time = gBenchmark->GetCpuTime("PHOSOnlineMon") ;
661 printf("took %f seconds for scanning, i.e. %f seconds per event %d \n",
662 time,time/fNevents,fNevents) ;
663
664 //Update canvas with histograms at the end
665 TIter nextCanvas(fCanvasList);
666 TCanvas * c ;
667 while((c=(TCanvas*)nextCanvas())){
668 c->Modified(kTRUE) ;
669 }
670
671 if(clu)delete clu ;
672 if(tsm)delete tsm ;
673 if(pid)delete pid ;
674 printf("delete 1 \n") ;
675 if(calibDB) delete calibDB ;
676 delete rawStream ;
677 delete rawReader ;
678}
679//____________________________________________________________________________
680void AliPHOSOnlineMonitor::Clean(){
681 //Cleans content of all histograms
682
683 TIter nextHisto(fHistosList);
684 TH1D * h ;
685 while((h=(TH1D*)nextHisto())){
686 h->Reset("ISE") ;
687 }
688 TIter nextCanvas(fCanvasList);
689 TCanvas * c ;
690 while((c=(TCanvas*)nextCanvas())){
691 c->Modified() ;
692 }
693}
694//____________________________________________________________________________
695void AliPHOSOnlineMonitor::Reset(){
696 //delets all canvas and histograms,
697 //marks buttons as unpressed
698
699 TIter nextHisto(fHistosList);
700 TH1D * h ;
701 while((h=(TH1D*)nextHisto())){
702 fHistosList->Remove(h) ;
703 delete h ;
704 }
705 TIter nextCanvas(fCanvasList);
706 TCanvas * c ;
707 while((c=(TCanvas*)nextCanvas())){
708 fCanvasList->Remove(c) ;
709 delete c ;
710 }
711 TObject *obj;
712 TGroupButton *button;
713 TIter next(fPrimitives);
714
715 //Mark buttons as anpressed
716 while ((obj = next())) {
717 if (obj->InheritsFrom(TGroupButton::Class())) {
718 button = (TGroupButton*)obj;
719 if (button->GetBorderMode() < 0){
720 button->SetBorderMode(1) ;
721 button->Modified(kTRUE);
722 }
723 }
724 }
725
726}
727//____________________________________________________________________________
728void AliPHOSOnlineMonitor::WriteHistograms(const char * filename){
729 //Write filled histograms to file
730 TFile * file = new TFile(filename,"Update") ;
731 file->cd() ;
732 TIter nextHisto(fHistosList);
733 TH1 * h ;
734 while((h=(TH1*)nextHisto())){
735 h->Write(0,TObject::kOverwrite) ;
736 }
737 file->Close() ;
738}