AliAnalysisTask*
AddTaskMCParticleFilter() 
{
  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  if (!mgr) {
    Error("AddTaskFMD", "No analysis manager to connect to.");
    return NULL;
  }   

  AliAnalysisTaskMCParticleFilter* mctask = new AliAnalysisTaskMCParticleFilter("mcfilter");
  mgr->AddTask(mctask);
  AliAnalysisDataContainer* histOut = 
    mgr->CreateContainer("mcfilter", TList::Class(), 
			 AliAnalysisManager::kOutputContainer, "mcParticles.root");
  mgr->ConnectInput(mctask, 0, mgr->GetCommonInputContainer());
  mgr->ConnectOutput(mctask, 0, mgr->GetCommonOutputContainer());
  mgr->ConnectOutput(mctask, 1, histOut);

  return mctask;
}