From 494744182ede7a2e1359e9522fec8d7d45b9d627 Mon Sep 17 00:00:00 2001 From: mtadel Date: Thu, 15 Mar 2007 14:39:37 +0000 Subject: [PATCH] From Bogdan: removed macro file. --- EVE/alice-macros/MUON_digits.C | 271 --------------------------------- 1 file changed, 271 deletions(-) delete mode 100644 EVE/alice-macros/MUON_digits.C diff --git a/EVE/alice-macros/MUON_digits.C b/EVE/alice-macros/MUON_digits.C deleted file mode 100644 index d141ef45aa9..00000000000 --- a/EVE/alice-macros/MUON_digits.C +++ /dev/null @@ -1,271 +0,0 @@ -#include "TGLViewer.h" - -namespace Alieve { -class MUONData; -class Event; -} - -Alieve::MUONData* g_muon_data = 0; -Alieve::Event* g_muon_last_event = 0; - -Int_t g_currentEvent = -1; -Bool_t g_fromRaw = kFALSE; - -void MUON_digits(Bool_t fromRaw = kFALSE, Bool_t showTracks = kTRUE) -{ - - TTree* dt = 0; - - if (Alieve::gEvent == 0) { - printf("No alieve event: use alieve_init(...) \n"); - return; - } - - if (g_currentEvent == Alieve::gEvent->GetEventId()) { - if (g_fromRaw == fromRaw) { - printf("Same event... \n"); - return; - } else { - if (g_fromRaw) { - printf("Same event with digits.\n"); - Alieve::gEvent->GotoEvent(g_currentEvent); - } else { - printf("Same event with raw.\n"); - Alieve::gEvent->GotoEvent(g_currentEvent); - } - } - } - - g_fromRaw = fromRaw; - - AliRunLoader* rl = Alieve::Event::AssertRunLoader(); - TString fileName = rl->GetFileName(); - Int_t length = fileName.Length(); - fileName.Resize(length-11); - fileName.Append("raw.root"); - - g_muon_data = new Alieve::MUONData; - - if (!fromRaw) { - rl->LoadDigits("MUON"); - dt = rl->GetTreeD("MUON", false); - if (dt == 0) { - cout << "No digits produced!" << endl; - } else { - cout << "Display aliroot digits!" << endl; - g_muon_data->LoadDigits(dt); - } - } else { - if (gSystem->AccessPathName(fileName.Data(),kFileExists)) { - cout << "No raw data produced!" << endl; - } else { - cout << "Display raw digits!" << endl; - g_muon_data->LoadRaw(fileName.Data()); - } - } - - g_muon_last_event = Alieve::gEvent; - - g_currentEvent = g_muon_last_event->GetEventId(); - - gStyle->SetPalette(1, 0); - - gReve->DisableRedraw(); - - Reve::RenderElementList* l = new Reve::RenderElementList("MUONChambers"); - l->SetTitle("MUON chambers"); - l->SetMainColor(Color_t(2)); - gReve->AddRenderElement(l); - - for (Int_t ic = 0; ic < 14; ic++) { - - Alieve::MUONChamber* mucha = new Alieve::MUONChamber(); - - mucha->SetFrameColor(2); - mucha->SetChamberID(ic); - mucha->SetDataSource(g_muon_data); - - gReve->AddRenderElement(l,mucha); - - } - - if (showTracks) MUON_tracks(); - - gReve->EnableRedraw(); - - TGLViewer* view = dynamic_cast(gReve->GetGLCanvas()->GetViewer3D()); - view->ResetCamerasAfterNextUpdate(); - gReve->GetGLCanvas()->Modified(); - gReve->GetGLCanvas()->Update(); - -} - -//_____________________________________________________________________________ -void MUON_tracks() { - - AliRunLoader* rl = Alieve::Event::AssertRunLoader(); - rl->LoadTracks("MUON"); - TTree* tt = rl->GetTreeT("MUON", false); - - TClonesArray *tracks = 0; - tt->SetBranchAddress("MUONTrack",&tracks); - tt->GetEntry(0); - - Int_t ntracks = tracks->GetEntriesFast(); - //printf("Found %d tracks. \n",ntracks); - - Reve::TrackList* lt = new Reve::TrackList("M-Tracks"); - lt->SetMainColor(Color_t(6)); - - gReve->AddRenderElement(lt); - - TMatrixD smatrix(2,2); - TMatrixD sums(2,1); - TMatrixD res(2,1); - - Float_t xRec, xRec0; - Float_t yRec, yRec0; - Float_t zRec, zRec0; - - Float_t zg[4] = { -1603.5, -1620.5, -1703.5, -1720.5 }; - - AliMUONTrack *mt; - Reve::RecTrack rt; - Int_t count; - for (Int_t n = 0; n < ntracks; n++) { - - count = 0; - - mt = (AliMUONTrack*) tracks->At(n); - - //printf("Match trigger %d \n",mt->GetMatchTrigger()); - - rt.label = n; - - Reve::Track* track = new Reve::Track(&rt, lt->GetRnrStyle()); - - //PH The line below is replaced waiting for a fix in Root - //PH which permits to use variable siza arguments in CINT - //PH on some platforms (alphalinuxgcc, solariscc5, etc.) - char form[1000]; - if (mt->GetMatchTrigger()) { - //PH track->SetName(Form("MUONTrack %2d (MT)", rt.label)); - sprintf(form,"MUONTrack %2d (MT)", rt.label); - track->SetName(form); - track->SetLineColor(7); - } else { - //PH track->SetName(Form("MUONTrack %2d ", rt.label)); - sprintf(form,"MUONTrack %2d ", rt.label); - track->SetName(form); - track->SetLineColor(6); - } - - AliMUONTrackParam *trackParam = mt->GetTrackParamAtVertex(); - xRec0 = trackParam->GetNonBendingCoor(); - yRec0 = trackParam->GetBendingCoor(); - zRec0 = trackParam->GetZ(); - - track->SetPoint(count,xRec0,yRec0,zRec0); - count++; - - Float_t xr[20], yr[20], zr[20]; - for (Int_t i = 0; i < 10; i++) xr[i]=yr[i]=zr[i]=0.0; - - Int_t nTrackHits = mt->GetNTrackHits(); - //printf("Nhits = %d \n",nTrackHits); - TClonesArray* trackParamAtHit; - for (Int_t iHit = 0; iHit < nTrackHits; iHit++){ - trackParamAtHit = mt->GetTrackParamAtHit(); - trackParam = (AliMUONTrackParam*) trackParamAtHit->At(iHit); - xRec = trackParam->GetNonBendingCoor(); - yRec = trackParam->GetBendingCoor(); - zRec = trackParam->GetZ(); - - //printf("Hit %d x %f y %f z %f \n",iHit,xRec,yRec,zRec); - - xr[iHit] = xRec; - yr[iHit] = yRec; - zr[iHit] = zRec; - - track->SetPoint(count,xRec,yRec,zRec); - count++; - - } - - Float_t xrc[20], yrc[20], zrc[20]; - Int_t nrc = 0; - if (mt->GetMatchTrigger() && 1) { - - for (Int_t i = 0; i < nTrackHits; i++) { - if (TMath::Abs(zr[i]) > 1000.0) { - //printf("Hit %d x %f y %f z %f \n",iHit,xr[i],yr[i],zr[i]); - xrc[nrc] = xr[i]; - yrc[nrc] = yr[i]; - zrc[nrc] = zr[i]; - nrc++; - } - } - - Double_t xv, yv; - Float_t ax, bx, ay, by; - - if (nrc < 2) continue; - - // fit x-z - smatrix.Zero(); - sums.Zero(); - for (Int_t i = 0; i < nrc; i++) { - xv = (Double_t)zrc[i]; - yv = (Double_t)xrc[i]; - //printf("x-z: xv %f yv %f \n",xv,yv); - smatrix(0,0) += 1.0; - smatrix(1,1) += xv*xv; - smatrix(0,1) += xv; - smatrix(1,0) += xv; - sums(0,0) += yv; - sums(1,0) += xv*yv; - } - res = smatrix.Invert() * sums; - ax = res(0,0); - bx = res(1,0); - - // fit y-z - smatrix.Zero(); - sums.Zero(); - for (Int_t i = 0; i < nrc; i++) { - xv = (Double_t)zrc[i]; - yv = (Double_t)yrc[i]; - //printf("y-z: xv %f yv %f \n",xv,yv); - smatrix(0,0) += 1.0; - smatrix(1,1) += xv*xv; - smatrix(0,1) += xv; - smatrix(1,0) += xv; - sums(0,0) += yv; - sums(1,0) += xv*yv; - } - res = smatrix.Invert() * sums; - ay = res(0,0); - by = res(1,0); - - Float_t xtc, ytc, ztc; - for (Int_t ii = 0; ii < 4; ii++) { - - ztc = zg[ii]; - ytc = ay+by*zg[ii]; - xtc = ax+bx*zg[ii]; - - //printf("tc: x %f y %f z %f \n",xtc,ytc,ztc); - - track->SetPoint(count,xtc,ytc,ztc); - count++; - - } - - } // end match trigger - - gReve->AddRenderElement(lt, track); - - } - -} - -- 2.39.3