-// Read rec. hits info
-// Loop on all detectors
-
- if (chamber > 6) return;
-
- printf("Entering LoadRecHits\n");
-
-
- AliRICH *pRICH = (AliRICH*)gAlice->GetModule("RICH");
- AliRICHChamber* iChamber;
-
- TClonesArray *pRICHrechits = pRICH->RecHitsAddress(chamber);
- printf ("Chamber:%d\n", chamber);
- if (pRICHrechits == 0) return;
-
- //RICH->ResetRecHits();
-
-
- Int_t nent=(Int_t)gAlice->TreeR()->GetEntries();
- gAlice->TreeR()->GetEvent(nent-1+cathode-1);
- Int_t nrechits = pRICHrechits->GetEntriesFast();
- printf ("nrechits:%d\n",nrechits);
- if (nrechits == 0) return;
- if (fRecpoints == 0) fRecpoints = new TObjArray(50);
-
- iChamber = &(pRICH->Chamber(chamber));
- AliRICHRecHit *mRec;
- AliRICHPoints *points = 0;
- //AliRICHEllipse *ellipse = 0;
- //
- //loop over all rechits and store their position
-
- points = new AliRICHPoints(nrechits);
- for (Int_t irec=0;irec<nrechits;irec++) {
- mRec = (AliRICHRecHit*)pRICHrechits->UncheckedAt(irec);
- fRecpoints->AddAt(points,irec);
- points->SetMarkerColor(38);
- points->SetMarkerStyle(8);
- points->SetMarkerSize(1.);
- points->SetParticle(-1);
- points->SetHitIndex(-1);
- points->SetTrackIndex(-1);
- points->SetDigitIndex(-1);
- Float_t vectorLoc[3]={mRec->fX,6.276,mRec->fY};
- Float_t vectorGlob[3];
- iChamber->LocaltoGlobal(vectorLoc,vectorGlob);
- points->SetPoint(irec,vectorGlob[0],vectorGlob[1],vectorGlob[2]);
- //Float_t theta = iChamber->GetRotMatrix()->GetTheta();
- //Float_t phi = iChamber->GetRotMatrix()->GetPhi();
- //ellipse=new TEllipse(vectorGlob[0],vectorGlob[2],10,10,0,360,phi);
- printf("Generating ellipse %d\n",irec);
- AliRICHEllipse *ellipse=new AliRICHEllipse(mRec->fX,mRec->fY,mRec->fOmega,mRec->fTheta,mRec->fPhi,mRec->fEmissPoint);
- ellipse->CerenkovRingDrawing(chamber,irec);
- //ellipse->SetFillStyle(1001);
- ellipse->SetMarkerColor(38);
- ellipse->Draw();
- //marker->SetRefObject((TObject*)points);
- //points->Set3DMarker(0, marker);
- }