#include <iostream>
#include <exception>
#include "SealKernel/Exception.h"

class PoolApplication
{
public:
  // Constructor. Initializes all the services.
  PoolApplication();
  // Destructor
  ~PoolApplication();
  // main method
  void run();
};


// The program main.
int main( int, char** )
{
  std::cout << "Reader application starts." << std::endl;
  try {
    PoolApplication app;
    app.run();
  }
  catch ( seal::Exception& e ) {
    std::cerr << e.what() << std::endl;
    return 1;
  }
  catch ( std::exception& e ) {
    std::cerr << e.what() << std::endl;
    return 1;
  }
  catch ( ... ) {
    std::cerr << "Funny error" << std::endl;
    return 1;
  }
  std::cout << "Reader application successfully finished." << std::endl;
  return 0;
}



#include <stdexcept>
#include "SealBase/SharedLibrary.h"
#include "SealBase/SharedLibraryError.h"
#include "POOLCore/POOLContext.h"

// Constructor of PoolApplication
PoolApplication::PoolApplication()
{
  // Loads the seal message stream
  pool::POOLContext::loadComponent( "SEAL/Services/MessageService" );
  // Set the verbosity threshold to warnings
  pool::POOLContext::setMessageVerbosityLevel( seal::Msg::Warning );

  // Loads the dictionary
  const std::string dictlibrary = "EventModelDict";
  try {
    seal::SharedLibrary::load( seal::SharedLibrary::libname( dictlibrary ) );
  }
  catch ( seal::SharedLibraryError *error) {
    throw std::runtime_error( error->explainSelf().c_str() );
  }
  std::cout << "Successfully loaded the services" << std::endl;
}


// Destructor of PoolApplication
PoolApplication::~PoolApplication()
{}


#include <memory>
#include <sstream>

#include "FileCatalog/URIParser.h"
#include "FileCatalog/IFileCatalog.h"

#include "PersistencySvc/ITransaction.h"
#include "PersistencySvc/ISession.h"

#include "DataSvc/IDataSvc.h"
#include "DataSvc/DataSvcFactory.h"

#include "Collection/Collection.h"
#include "AttributeList/AttributeList.h"

#include "Event.h"


// Main method of PoolApplication
void
PoolApplication::run()
{
  // Create a file catalog and connect to it
  pool::URIParser p;
  p.parse();
  std::auto_ptr<pool::IFileCatalog> fileCatalog( new pool::IFileCatalog );
  fileCatalog->addReadCatalog( p.contactstring() );
  fileCatalog->connect();

  // Create a data service and tell it which catalog to use
  std::auto_ptr<pool::IDataSvc> dataSvc( pool::DataSvcFactory::instance( fileCatalog.get() ) );

  // Start a transaction for the file catalog
  fileCatalog->start();

  // Start a read transaction
  dataSvc->transaction().start( pool::ITransaction::READ );

  // Create a new explicit collection
  pool::Collection< pool_tutorial::Event > eventCollection( dataSvc.get(),
							    "RootCollection",
							    "eventCollection", "",
							    pool::ICollection::READ );

  // Select the events with HCAL deposited energy greater than 0.
  pool::Collection<pool_tutorial::Event>::Iterator event = eventCollection.select( "E_hcal > 0" );
  
  // For every event read the number of clusters in ECAL. Calculate the average.
  float nClusters = 0;
  double deEcal = 0;
  int nEvents = 0;
  while ( event.next() ) {
    const pool_tutorial::Sector& ecal = event->ecalData();
    nClusters += ecal.numberOfClusters();
    double de = 0;
    event.attributeList()[ "E_ecal" ].getValue<double>( de );
    deEcal += de;
    ++nEvents;
  }
  std::cout << "Read back " << nEvents << " events " << std::endl;
  std::cout << "Average number of clusters in ECal : " << nClusters / nEvents << std::endl;
  std::cout << "Average energy per cluster in ECal : " << deEcal / nClusters << std::endl;
  
  // Close the collection
  // eventCollection.close();

  // commit the transaction
  dataSvc->transaction().commit();

  // Commit the file catalog transaction
  fileCatalog->commit();
}