#ifndef PPC_H
#define PPC_H

#include <icetray/I3Frame.h>
#include <icetray/I3Module.h>
#include <simclasses/I3MMCTrack.h>
#include <dataclasses/physics/I3MCTree.h>
#include <dataclasses/physics/I3MCTreeUtils.h>
#include <dataclasses/geometry/I3Geometry.h>

class i3ppc : public I3Module{

 public:
  SET_LOGGER("i3ppc");
  i3ppc(const I3Context& ctx);
  void Physics(I3FramePtr frame);

 private:

  bool verbose;
  map<pair<int,unsigned long long>,const I3MMCTrack *> i3mmctracks;
  void pparticle(const I3Particle&, int, I3MCTreeConstPtr);
  double res(double, double);
};

#endif //PPC_H