1 // $Id: MUON_display.C 24485 2008-03-13 15:27:38Z mtadel $
2 // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
4 /**************************************************************************
5 * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
6 * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for *
7 * full copyright notice. *
8 **************************************************************************/
11 class AliEveEventManager;
13 AliEveMUONData *g_muon_data = 0;
15 Int_t g_currentEvent = -1;
16 Bool_t g_fromRaw = kFALSE;
18 AliMagFMaps *g_field = 0;
20 void MUON_displaySimu(Bool_t fromRaw = kFALSE, Bool_t showTracks = kTRUE, Bool_t clustersFromESD = kTRUE)
23 // display from simulated digits (or produced raw data)
24 // tracks: ESD, Refs, MC
27 if (!AliMpSegmentation::Instance()) AliMpCDB::LoadMpSegmentation();
28 if (!AliMpDDLStore::Instance()) AliMpCDB::LoadDDLStore();
31 printf("Loading field map...\n");
32 g_field = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG);
33 AliTracker::SetFieldMap(g_field, kFALSE);
34 AliMUONTrackExtrap::SetField(AliTracker::GetFieldMap());
41 if (AliEveEventManager::GetMaster() == 0) {
42 printf("No alieve event: use alieve_init(...) \n");
46 if (g_currentEvent == AliEveEventManager::GetMaster()->GetEventId()) {
47 if (g_fromRaw == fromRaw) {
48 printf("Same event... \n");
52 printf("Same event with digits.\n");
53 AliEveEventManager::GetMaster()->GotoEvent(g_currentEvent);
55 printf("Same event with raw.\n");
56 AliEveEventManager::GetMaster()->GotoEvent(g_currentEvent);
63 TString dataPath = TString(AliEveEventManager::GetMaster()->GetTitle());
64 dataPath.Append("/rawmuon.root");
66 AliRunLoader* rl = AliEveEventManager::AssertRunLoader();
67 g_muon_data = new AliEveMUONData;
70 rl->LoadDigits("MUON");
71 dt = rl->GetTreeD("MUON", false);
73 cout << "No digits produced!" << endl;
75 cout << "With aliroot digits!" << endl;
76 g_muon_data->LoadDigits(dt);
79 if (gSystem->AccessPathName(dataPath.Data(),kFileExists)) {
80 cout << "No raw data produced!" << endl;
82 cout << "With raw digits!" << endl;
83 g_muon_data->LoadRaw(dataPath.Data());
87 TString esdDataPath = TString(AliEveEventManager::GetMaster()->GetTitle());
88 esdDataPath.Append("/AliESDs.root");
89 if (clustersFromESD) {
90 g_muon_data->LoadRecPointsFromESD(esdDataPath.Data());
92 rl->LoadRecPoints("MUON");
93 ct = rl->GetTreeR("MUON", false);
94 g_muon_data->LoadRecPoints(ct);
98 ht = rl->GetTreeH("MUON", false);
99 g_muon_data->LoadHits(ht);
101 g_currentEvent = AliEveEventManager::GetMaster()->GetEventId();
103 gStyle->SetPalette(1, 0);
105 gEve->DisableRedraw();
107 TEveElementList* l = new TEveElementList("MUONChambers");
108 l->SetTitle("MUON chambers");
112 for (Int_t ic = 0; ic < 14; ic++)
114 AliEveMUONChamber* mucha = new AliEveMUONChamber(ic);
116 mucha->SetFrameColor(2);
117 mucha->SetChamberID(ic);
119 mucha->SetDataSource(g_muon_data);
121 gEve->AddElement(mucha, l);
130 gEve->Redraw3D(kTRUE);
131 gEve->EnableRedraw();
134 //______________________________________________________________________________
135 void MUON_ESD_tracks()
137 AliESDEvent* esd = AliEveEventManager::AssertESD();
139 TEveTrackList* lt = new TEveTrackList("ESD-Tracks");
143 gEve->AddElement(lt);
147 Int_t nMuonTracks = esd->GetNumberOfMuonTracks();
149 for (Int_t n = 0; n < nMuonTracks; n++)
151 mt = esd->GetMuonTrack(n);
153 if (mt->GetNHit() == 0) continue;
158 AliEveMUONTrack* track = new AliEveMUONTrack(&rt, lt->GetPropagator());
160 track->MakeESDTrack(mt);
162 gEve->AddElement(track, lt);
167 //______________________________________________________________________________
168 void MUON_Ref_tracks()
170 TString dataPathESD = TString(AliEveEventManager::GetMaster()->GetTitle());
171 dataPathESD.Append("/AliESDs.root");
172 TString dataPathSIM = TString(AliEveEventManager::GetMaster()->GetTitle());
173 dataPathSIM.Append("/");
175 AliMUONRecoCheck recoCheck(dataPathESD.Data(),dataPathSIM.Data());
176 AliMUONVTrackStore* trackRefStore = recoCheck.ReconstructibleTracks(AliEveEventManager::GetMaster()->GetEventId());
177 TIter next(trackRefStore->CreateIterator());
178 AliMUONTrack* trackRef;
180 TEveTrackList* lt = new TEveTrackList("Ref-Tracks");
183 gEve->AddElement(lt);
187 while ( ( trackRef = static_cast<AliMUONTrack*>(next()) ) )
191 AliEveMUONTrack* track = new AliEveMUONTrack(&rt, lt->GetPropagator());
193 track->MakeRefTrack(trackRef);
195 gEve->AddElement(track, lt);
200 //______________________________________________________________________________
201 void MUON_MC_tracks()
203 Double_t RADDEG = 180.0/TMath::Pi();
205 AliRunLoader* rl = AliEveEventManager::AssertRunLoader();
206 rl->LoadKinematics();
207 AliStack* stack = rl->Stack();
209 Int_t nPrimary = stack->GetNprimary();
210 Int_t nTracks = stack->GetNtrack();
212 TEveTrackList* lt = new TEveTrackList("MC-Tracks");
216 gEve->AddElement(lt);
222 Int_t nHitTracks = g_muon_data->GetNTrackList();
224 for (Int_t i = 0; i < nHitTracks; i++)
226 index = g_muon_data->GetTrack(i);
227 if (index >= nTracks) {
228 cout << "TEveHit track index larger than number in stack!" << endl;
232 part = stack->Particle(index);
233 if (part->P() < 0.001) continue; // skip momenta < 1.0 MeV/c
236 AliEveMUONTrack* track = new AliEveMUONTrack(&rt, lt->GetPropagator());
238 track->MakeMCTrack(part);
240 gEve->AddElement(track, lt);
243 rl->UnloadKinematics();