void CbmRichID ( Int_t NStart = 0, Int_t NStop = 0){

  TStopwatch timer;
  timer.Start();

  // Special libs...
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();

  gSystem->Load("libGeoBase");
  gSystem->Load("libBase");
  gSystem->Load("libPassive");
  gSystem->Load("libSts");
  gSystem->Load("libTof");
  gSystem->Load("libRich");
  gSystem->Load("libTrd");
  gSystem->Load("libEcal");
  gSystem->Load("libGen");



  // Here we go!!!

  CbmRunAna *fRun = new CbmRunAna();
  CbmRuntimeDb *rtdb = fRun->GetRuntimeDb();
  CbmParRootFileIo *io1 = new CbmParRootFileIo();
  io1->open("parfiles/testparams.root");
  rtdb->setFirstInput(io1);


  fRun->SetInputFile("Electron25GeV.root");
  fRun->AddFriend ( "Electron_hit.root" );
  fRun->AddFriend ( "Electron_ring.root" );
  fRun->AddFriend ( "Electron_cp.root" );

  fRun->SetOutputFile("test.root");


 
 CbmRichID *p= new CbmRichID("rich","rich task");

	fRun->AddTask(p);

//	fRun->LoadGeometry();
	fRun->Init();
 	if(NStart == 0 && NStop == 0)
   	 fRun->Run ();
  	else
    	fRun->Run ( NStart, NStop );

	//fRun->Run();
//        p->Draw();
 	timer.Stop();
	Double_t rtime = timer.RealTime();
	Double_t ctime = timer.CpuTime();
 	delete p;
	printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
 
	
}