Risetime1000LLL.h

Go to the documentation of this file.

#ifndef _Risetime1000_Risetime1000LLL_h_
#define _Risetime1000_Risetime1000LLL_h_

// Offline headers
#include <fwk/VModule.h>
#include <evt/Event.h>
#include <sevt/SEvent.h>
#include <sevt/Station.h>
#include <sevt/StationConstants.h>

// ROOT headers
#include <TGraphErrors.h>
#include <TF1.h>

// Standard c++ headers
#include <vector>


namespace Risetime1000 {

  struct StationRisetimeDataLLL {
    unsigned int fStationId;
    double fCorrRisetime;
    double fCorrRisetimeError;
    double fDistance;
    double fDistanceError;
    unsigned short fRejectCode;
  };


  struct RisetimeResultsLLL {
    ~RisetimeResultsLLL() { for (auto p : fStationData) delete p; }

    std::vector<StationRisetimeDataLLL*> fStationData;
    double fFitPar0 = 0;
    double fFitPar1 = 0;
    double fRisetime1000 = -1;
    double fRisetime1000Error = -1;
    double fRisetime1000Chi2 = -1;
    double fRisetime1000NDF = -1;
    double fXmax = -1;
    double fXmaxErrorUp = -1;
    double fXmaxErrorDown = -1;
  };


  class Risetime1000LLL : public fwk::VModule {

  public:
    Risetime1000LLL() = default;
    virtual ~Risetime1000LLL();

    // Init, Run, and Finish functions
    fwk::VModule::ResultFlag Init() override;
    fwk::VModule::ResultFlag Run(evt::Event& event) override;
    fwk::VModule::ResultFlag Finish() override { return eSuccess; }

    // Static public members for access outside this module
    static double fgRisetime;
    static double fgRisetimeError;
    static double fgRisetimeReducedChi2;
    static RisetimeResultsLLL* fgRisetimeResults;

  protected:
    // Variables read from the xml file to configure the
    // risetime recalculation (if used)
    float fRiseTimeStartFraction = 0.1;  // 10%
    float fRiseTimeStopFraction = 0.5;  // 50%

  private:
    // Calculate the Risetime and return that value
    double FitEventRiseTime();
    double RecalculateRiseTime(const sevt::Station& station);
    bool fUseDLECorrectedTrace;

    // Variables used for the riestime fit
    TFormula* fRTWeights = nullptr;
    unsigned int fRejectedStations = 0;

    // Variables read from the xml file to configure the risetime routines
    float fMinimumSignalForRiseTimeFit = 10;  // Signal size in VEM
    float fMinimumDistanceForRiseTimeFit = 0*utl::meter;  // Distance in meters
    float fMaximumDistanceForRiseTimeFit = 1600*utl::meter;  // Distance in meters
    float fDistanceToInterpolateFitResult = 1000*utl::meter;  // Distance in meters to evaluate the fit at
    bool fIncludeSaturatedForRiseTimeFit = false;
    bool fDoRiseTimeFit = true;
    bool fForceRiseTimeRecalculation = false;  // Don't trust others to calculate it
    std::string fRTWeightingFunction;  // Function describing the
                                       // weight as a function of
                                       // distance/signal size

    // Derived variables
    sevt::StationConstants::SignalComponent fComponent = sevt::StationConstants::eTotal;

    REGISTER_MODULE("Risetime1000LLL", Risetime1000LLL);

  };

}


#endif