#include "CalorimeterConstants.h"
#include "CalorimeterSimulation.h"
#include "CalorimeterDigi.h"
#include "TRandom.h"

#include <iostream>

void simulate(Event * event) 
{
  double genEnergie = 50. ;
  // simulate the event
  event->seteTrue( genEnergie ); // fixed true energy

  CalorimeterSimulation calsim;
  //  float xstart = gRandom->Uniform(CalConst::XYMin, CalConst::XYMax); 
  //  float ystart = gRandom->Uniform(CalConst::XYMin, CalConst::XYMax); 

  float xstart = gRandom->Uniform( 0. , CalConst::XYSize ) ; 
  float ystart = gRandom->Uniform( 0. , CalConst::XYSize ) ; 

  event->setTruePosition( xstart , ystart ) ;

  calsim.SimulateShower(xstart,ystart,genEnergie);
  std::vector<CalCell> cells;
  calsim.CalorimeterData(cells);

  // Applique digitisation
  CalorimeterDigi caldigi( 0.1 ) ;
  caldigi.ReadCalorimeterTrueCells( cells ) ;
  std::vector< CalCell > digis = caldigi.CreateDigits() ;

  event->setCalData( digis ) ;  
}