In muon-related macros in EVE:
[u/mrichter/AliRoot.git] / EVE / alice-macros / muon_raw.C
1 // $Id$
2
3 /**************************************************************************
4  * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
5  * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for          *
6  * full copyright notice.                                                 *
7  **************************************************************************/
8
9 /// \ingroup evemacros
10 /// \file muon_raw.C
11 /// \brief Macro to visualise rootified raw-data from MUON spectrometer 
12 /// (both tracker and trigger).
13 ///
14 /// Use muon_raw() in order to run it
15 ///
16 /// Needs that alieve_init() is already called
17 ///
18 /// \author P. Pillot, L. Aphecetche; Subatech
19
20 #if !defined(__CINT__) || defined(__MAKECINT__)
21
22 #include "AliMUONDigitMaker.h"
23 #include "AliMUONDigitStoreV2R.h"
24
25 #include "AliMpCDB.h"
26
27 #include "AliRawReader.h"
28
29 #include "EveBase/AliEveEventManager.h"
30
31 #include <TEveManager.h>
32 #include <TEveUtil.h>
33 #include <TEveQuadSet.h>
34
35 #include <TStyle.h>
36 #include <Riostream.h>
37 #include <TROOT.h>
38
39 #endif
40
41 void muon_raw()
42 {
43   // load mapping
44   AliMpCDB::LoadAll(kFALSE);
45
46   // load raw data
47   AliRawReader* reader = AliEveEventManager::AssertRawReader();
48   if ( reader->GetEventHeader() ) 
49     cout << "RUN " << reader->GetRunNumber() << " EVENT " << reader->GetEventIndex() << endl;
50   else
51     cout << "NO EVENT HEADER ?" << endl;
52   
53   // convert raw to digits
54   AliMUONDigitMaker digitMaker;
55   digitMaker.SetMakeTriggerDigits(kTRUE);
56   AliMUONDigitStoreV2R digitStore;
57   digitMaker.Raw2Digits(reader,&digitStore);
58   if (digitStore.GetSize() == 0 && !gEve->GetKeepEmptyCont()) return;
59   
60   // container for graphic representation of digits
61   TEveElementList* cont = new TEveElementList("MUON Raw digits");
62   cont->SetTitle(Form("N=%d",digitStore.GetSize()));
63   
64   TEveQuadSet* bending = new TEveQuadSet(TEveQuadSet::kQT_RectangleXY, kFALSE, 32);
65   bending->SetName("Bending");
66   bending->SetRenderMode(TEveDigitSet::kRM_Fill);
67   bending->SetPickable(kFALSE);
68   cont->AddElement(bending);
69   
70   TEveQuadSet* nonBending = new TEveQuadSet(TEveQuadSet::kQT_RectangleXY, kFALSE, 32);
71   nonBending->SetName("Non bending");
72   nonBending->SetRenderMode(TEveDigitSet::kRM_Line);
73   nonBending->SetPickable(kFALSE);
74   cont->AddElement(nonBending);
75   
76   // add digits to the containers
77   TEveUtil::LoadMacro("muon_digits.C+");
78   TIter next(digitStore.CreateIterator());
79   gROOT->ProcessLine(Form("add_muon_digits((TIter*)%p, (TEveQuadSet*)%p, (TEveQuadSet*)%p, kTRUE);",
80                           &next, bending, nonBending));
81   
82   // set containers' title
83   bending->SetTitle(Form("N=%d",bending->GetPlex()->Size()));
84   nonBending->SetTitle(Form("N=%d",nonBending->GetPlex()->Size()));
85   
86   // automatic scaling
87   gStyle->SetPalette(1);
88   bending->AssertPalette();
89   nonBending->AssertPalette();
90   
91   // add graphic containers
92   gEve->DisableRedraw();
93   gEve->AddElement(cont);
94   gEve->EnableRedraw();
95   gEve->Redraw3D();
96 }