small fixes for Kalman fitter

[u/mrichter/AliRoot.git] / EVE / alice-macros / trd_friend_tracks.C

void trd_friend_tracks(TEveElement *cont = 0)
{

  // Link data containers
  AliESDfriend *eventESDfriend = 0x0;
  if(!(eventESDfriend = AliEveEventManager::AssertESDfriend())){
    Warning("trd_friend_tracks", "AliESDfriend not found");
    return 0x0;
  }

  AliESDEvent* esd = AliEveEventManager::AssertESD();

  AliEveEventManager::AssertGeometry();

  AliMagFMaps *field = new AliMagFMaps("Maps","Maps", 2, 1., 10., AliMagFMaps::k5kG);
  AliTracker::SetFieldMap(field, kTRUE);

  AliTRDReconstructor *reco = new AliTRDReconstructor();
  reco->SetRecoParam(AliTRDrecoParam::GetLowFluxParam());
  reco->SetOption("!nn");

  AliEveTRDTrackList *tracks = new AliEveTRDTrackList("TRD Tracks");
  for (Int_t n=0; n<esd->GetNumberOfTracks(); n++){
    AliESDtrack* esdTrack = esd->GetTrack(n);
    AliESDfriendTrack *friendTrack = eventESDfriend->GetTrack(n);
  
    TObject *cal = 0x0;
    Int_t ical = 0;
    while(cal = friendTrack->GetCalibObject(ical++)){
      if(strcmp(cal->IsA()->GetName(), "AliTRDtrackV1") != 0) continue;
      AliTRDtrackV1 *trackObj = dynamic_cast<AliTRDtrackV1 *>(cal);
      trackObj->SetReconstructor(reco);
      AliEveTRDTrack *trackEve = new AliEveTRDTrack(trackObj);
      tracks->AddElement(trackEve);
      trackEve->SetESDstatus(esdTrack->GetStatus());
    }
  }
	
  tracks->SetTitle(Form("Tracks %d", tracks->NumChildren()));
  tracks->StampObjProps();
  gEve->AddElement(tracks, cont);
  
  gEve->Redraw3D();

  TGLViewer *v = gEve->GetGLViewer();
  v->SetCurrentCamera(TGLViewer::kCameraOrthoXOY);
  ((TGLOrthoCamera&)v->CurrentCamera()).SetEnableRotate(kTRUE);
  v->UpdateScene();
  
  return;
}