SdFootprintAnalyzer.h

Go to the documentation of this file.

#ifndef _SdFootprintAnalyzer_h_
#define _SdFootprintAnalyzer_h_

// Offline includes
#include <fwk/VModule.h>

#include <utl/TimeStamp.h>
#include <utl/Point.h>
#include <utl/Vector.h>
#include <utl/CoordinateSystemPtr.h>
#include <utl/Trace.h>

#include <evt/Event.h>

#include <sevt/Station.h>
#include <sevt/SEvent.h>
#include <sevt/SdFootprintData.h>

using namespace std;

namespace SdFootprintAnalyzerNS{   // you should assign your module to a namespace
  
  class Station {
  public:
    Station(int Id, utl::TimeStamp Time, utl::Point Pos, double Signal, double AoP, bool TOT, double Lp) :
      fId(Id),
      fTime(Time),
      fPos(Pos),
      fSignal(Signal),
      fAoP(AoP),
      fTOT(TOT),
      fLp(Lp)
    { }   
    int fId;
    utl::TimeStamp fTime;
    utl::Point fPos;
    double fSignal;
    double fAoP;
    bool fTOT;
    double fLp;
  };
  
  typedef std::vector<SdFootprintAnalyzerNS::Station> StationCollection;
  
  class ByIncreasingDistanceToAxis {
  public:
    ByIncreasingDistanceToAxis(const utl::CoordinateSystemPtr& cs);
    bool operator()(const sevt::Station* const s1, const sevt::Station* const s2) const;
  private:
    utl::CoordinateSystemPtr fShowerCS;
  };
  
  class ByIncreasingTime {
  public:
    bool operator()(SdFootprintAnalyzerNS::Station i, SdFootprintAnalyzerNS::Station j) const;
  };
  
  class SdFootprintAnalyzer : public fwk::VModule{
    
  public:
    
    SdFootprintAnalyzer();
    virtual ~SdFootprintAnalyzer();
    
    fwk::VModule::ResultFlag Init();
    fwk::VModule::ResultFlag Run(evt::Event & event);
    fwk::VModule::ResultFlag Finish();

  private:

    utl::Point ComputeBarycenter(const sevt::SEvent& theSEvent) const;
    StationCollection GetFootprintStations(sevt::SEvent& theSEvent) const;
    void ComputeGroundVariables(evt::Event& event, StationCollection station, sevt::SdFootprintData& theFootprint);
    double ComputeAoPAsymmetry(StationCollection station);
    
    utl::TraceD GetVEMChargeTrace(const sevt::Station& station) const;
 
    void SortStationsAroundBarycenter(evt::Event& event, utl::Point barycenter) const;
    
    int fVerbose;
    double fMinProjDistForSpeed;
    bool fRecomputeBarycenter;
    double fMaxPhiDiffLine;

    // Register the module
    REGISTER_MODULE("SdFootprintAnalyzerOG",SdFootprintAnalyzer);
  };

} // SdFootprintAnalyzer namespace

#endif

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make  -k"
// End: