FdProfileConstrainedGeometryFit.cc

Go to the documentation of this file.

#include <evt/Event.h>
#include <evt/ShowerFRecData.h>

#include <fevt/Eye.h>
#include <fevt/EyeHeader.h>
#include <fevt/FEvent.h>
#include <fevt/EyeRecData.h>

#include <det/Detector.h>

#include <fdet/Eye.h>
#include <fdet/FDetector.h>

#include <fwk/RunController.h>

#include "FdProfileConstrainedGeometryFit.h"

#include "ChiZeroRegression.h"
#include "ProfileChi2.h"

using namespace std;
using namespace fevt;
using namespace utl;


namespace FdProfileConstrainedGeometryFit {

fwk::VModule::ResultFlag
FdProfileConstrainedGeometryFit::Init()
{
  fPCGFitter.Init();
  return eSuccess;
}


fwk::VModule::ResultFlag
FdProfileConstrainedGeometryFit::Finish()
{
  return eSuccess;
}


fwk::VModule::ResultFlag
FdProfileConstrainedGeometryFit::Run(evt::Event& event)
{
  if (!event.HasFEvent())
    return eSuccess;
  FEvent& fevent = event.GetFEvent();

  bool success = false;

  for (FEvent::EyeIterator eye = fevent.EyesBegin(ComponentSelector::eHasData),
       end = fevent.EyesEnd(ComponentSelector::eHasData); eye != end; ++eye) {
    if (!eye->HasRecData() ||
        !eye->GetRecData().HasFRecShower() ||
        det::Detector::GetInstance().GetFDetector().GetEye(*eye).IsVirtual() ||
        eye->GetRecData().GetTimeFitNDof() <=0)
      continue;

    cout << "PCGF for eye " << eye->GetId() << " ...\n";
    if (fPCGFitter.Run(*eye)) {
      success = true;
    } else {
      eye->SetStatus(ComponentSelector::eDeSelected);
    }

    cout << "PCGF for eye " << eye->GetId() << " DONE\n";
  }

  if (success) {
    ++fwk::RunController::GetInstance().GetRunData().GetNamedCounters()["FdProfileConstrainedGeometryFit/Good"];
    return eSuccess;
  }
  
  return eContinueLoop;
}


void
FdProfileConstrainedGeometryFit::FitTimingOfGivenGeometry(const Eye& eye)
  const
{
  const EyeRecData& eyeRec = eye.GetRecData();
  ChiZeroRegression chi0reg(eyeRec);

  cout << "Parameters from Axis Finder:\n"
          "Chi0: " << eyeRec.GetChiZero() << "\n"
          "Chi^2: " << eyeRec.GetTimeFitChiSquare() << "\t"
          "Rp: " << eyeRec.GetRp()/m << "\t"
          "T0: " << eyeRec.GetTZero()/ns << "\n";

   double chi2 = 0;
   double rp = 0;
   double t0 = 0;
   chi0reg(eyeRec.GetChiZero(), chi2, rp, t0);

   cout << "scanning Chi_0\n";

   for (int i = 0; i <= 180; ++i) {
     cout << i << " degrees\n";

     chi0reg(i*degree, chi2, rp, t0);
     cout << "Rp: " << rp/m << "\t"
             "T0: " << t0/ns << "\t"
             "Chi^2_t: " << chi2 << "\t";
  }
}

}