void muon_clusters()
{
  AliRunLoader* rl =  Alieve::Event::AssertRunLoader();

  rl->LoadRecPoints("MUON");

  TTree* ct = rl->GetTreeR("MUON", false);

  Alieve::MUONDigitsInfo* di = new Alieve::MUONDigitsInfo();
  di->SetRTree(ct);

  gStyle->SetPalette(1, 0);

  gReve->DisableRedraw();

  rl->LoadgAlice();

  char name[128];
  char title[128];

  /* CLUSTERS */

  for (Int_t iSta = 1; iSta <= 5; iSta++) {

    for (Int_t iCha = 1; iCha <= 2; iCha++) {

      sprintf(name,"M-ST%1dCH%1d/Clusters",iSta,iCha);

      Reve::RenderElementList* l = new Reve::RenderElementList(name);
      
      if (iSta <= 5) {
	sprintf(title,"Station %1d chamber %1d (tracking)",iSta,iCha);
      } else {
	sprintf(title,"Station %1d chamber %1d (trigger)",iSta,iCha);
      }
      
      l->SetTitle(title);
      l->SetMainColor((Color_t)4);
      gReve->AddRenderElement(l);
    
      Int_t iChamber = (iSta-1) * 2 + iCha; 

      AliMpDEIterator it;
      for ( it.First(iChamber-1); ! it.IsDone(); it.Next() ) {
	
	Int_t detElemId = it.CurrentDE();
	
	Alieve::MUONModule* m = new Alieve::MUONModule(detElemId, 0, di, 0, 1, (Color_t)2);
	gReve->AddRenderElement(l, m);      
	
      }

    }
    
  }
  
  gReve->EnableRedraw();

}