]>
Commit | Line | Data |
---|---|---|
fe88cf1b | 1 | red() |
2 | { | |
3 | gSystem->Load("libMinuit.so"); | |
4 | gSystem->Load("libVMC.so"); | |
5 | gSystem->Load("libESD.so"); | |
6 | gSystem->Load("libSTEER.so"); | |
7 | gSystem->Load("libCDB.so"); | |
8 | ||
9 | gSystem->Load("libGed.so"); | |
10 | gSystem->Load("libRGL.so"); | |
11 | gSystem->Load("libGeom.so"); | |
12 | gSystem->Load("libVMC.so"); | |
13 | ||
14 | gSystem->Load("libReve.so"); | |
15 | ||
16 | gSystem->Load("libHMPIDbase.so"); | |
17 | gSystem->Load("libHMPIDsim.so"); | |
18 | gSystem->Load("libHMPIDrec.so"); | |
19 | ||
20 | TGeoManager::Import("geometry.root"); | |
21 | ||
22 | AliHMPIDParam *pParam=AliHMPIDParam::Instance(); | |
23 | ||
24 | AliRunLoader *pAL=AliRunLoader::Open(); pAL->LoadHits ("HMPID"); TTree *pHitT=pAL->GetTreeH("HMPID", false); | |
25 | pAL->LoadDigits("HMPID"); TTree *pDigT=pAL->GetTreeD("HMPID", false); | |
26 | ||
27 | Reve::PointSet *pHitPnt=new Reve::PointSet("Hits")); | |
28 | TPolyMarker3D *pDigPnt=new TPolyMarker3D; pDigPnt->SetName("Digits"); pDigPnt->SetMarkerColor(kGreen);iPntCnt=0; | |
29 | ||
30 | TPointSelector ps(pHitT,pHitPnt,"fX:fY:fZ",""); ps.Select(); | |
31 | ||
32 | TClonesArray *pDigLst=new TClonesArray("AliHMPIDDigit"); //this is tmp dig list per chamber | |
33 | ||
34 | for(Int_t iCh=0;iCh<7;iCh++){ | |
35 | pDigT->SetBranchAddress(Form("HMPID%i",iCh+1),&pDigLst); | |
36 | pDigT->GetEntry(0); | |
37 | for(Int_t iDig=0;iDig<pDigLst->GetEntries();iDig++){ | |
38 | AliHMPIDDigit *pDig=(AliHMPIDDigit*)pDigLst->At(iDig); | |
39 | TVector2 lors=pParam->Pad2Loc(pDig->PadX(),pDig->PadY()); | |
40 | TVector3 mars=pParam->Lors2Mars(iCh,lors.X(),lors.Y()); | |
41 | pDigPnt->SetPoint(iPntCnt++,mars.X(),mars.Y(),mars.Z()); | |
42 | }//digits loop for chamber | |
43 | }//chambers loop | |
44 | ||
45 | if(!gReve) new Reve::RGTopFrame(0,0,0,2); | |
46 | gReve->AddGlobalRenderElement(new Reve::RenderElementObjPtr(pDigPnt)); | |
47 | gReve->AddRenderElement(pHitPnt); | |
48 | gReve->Redraw3D(); | |
49 | } |