#include "CalorimeterReconstruction.h"

void reconstruct(Event * event) 
{
  std::vector< CalCell > cells = event -> CalData() ;
  
  // Cree une particule 
  CalorimeterReconstruction reco( cells ) ;
  CalParticle theCandidate = reco.CreateParticle() ;

  event -> setParticle( theCandidate ) ;
  event -> seteReco( theCandidate.energy() ) ;

  event -> setRecoPosition( theCandidate.x() , 
                            theCandidate.y() ) ;

//   std::cout << theCandidate.energy() 
//             << " " << event -> xTrue() 
//             << " " << event -> yTrue() 
//             << " " << event -> xReco() 
//             << " " << event -> yReco() << std::endl ;
}