RdTopDownStationSelector.h

Go to the documentation of this file.

#ifndef _RdTopDownStationSelector_RdTopDownStationSelector_h_
#define _RdTopDownStationSelector_RdTopDownStationSelector_h_

#include <fwk/VModule.h>

namespace RdTopDownStationSelector {

class RdTopDownStationSelector : public fwk::VModule {

 public:
  fwk::VModule::ResultFlag Init() override;
  fwk::VModule::ResultFlag Run(evt::Event& event) override;
  fwk::VModule::ResultFlag Finish() override;

 private:
  int fCurrentRunId = 0;
  int fCurrentEventId = 0;
  std::string fCurrentHeaderId;

  double fIteration = 0;  // iteration of time calibration
  int fLastAddedStationId = 0;

  unsigned int fMaxNumberOfStations = 999;
  unsigned int fInitialNumberOfSignalStations = 3;
  unsigned int fCurrentNumberOfStations = 0;

  // chi2 per ndf of wave fit of previous station selection
  double fChiSquaredNDF = 1;

  // the probability denotes the probability that an observed Chi-squared exceeds
  // the value chi2 by chance, even for a correct model.
  double fChiSquareProbabilityCutValue = 0.1;

  // stop adding station after one station war rejected?
  bool fStopAtFirstStation = true;

  // check if the incoming direction reconstructed from the closest 3 stations
  // roughly matches the SD reconstruction
  bool fInitialSDCheck = true;

  double fMaxIncomingDirectionDifference = 0;
  double fAdditionalSystematicError = 0;

  REGISTER_MODULE("RdTopDownStationSelector", RdTopDownStationSelector);
};

}

#endif