FullGridChi2TrackTimeExtractor.cc

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#include <tracking/eventTimeExtraction/findlets/FullGridChi2TrackTimeExtractor.h>
#include <tracking/eventTimeExtraction/findlets/BaseEventTimeExtractor.icc.h>
#include <tracking/eventTimeExtraction/findlets/GridEventTimeExtractor.icc.h>
#include <tracking/eventTimeExtraction/findlets/IterativeEventTimeExtractor.icc.h>
#include <tracking/trackFindingCDC/utilities/StringManipulation.h>
#include <tracking/dataobjects/RecoTrack.h>
#include <framework/core/ModuleParamList.templateDetails.h>
 
using namespace Belle2;
using namespace TrackFindingCDC;
 
template class Belle2::GridEventTimeExtractor<Chi2BasedEventTimeExtractor>;
 
FullGridChi2TrackTimeExtractor::FullGridChi2TrackTimeExtractor()
{
  addProcessingSignalListener(&m_finalExtractor);
  addProcessingSignalListener(&m_gridExtractor);
}

void FullGridChi2TrackTimeExtractor::exposeParameters(ModuleParamList* moduleParamList, const std::string& prefix)
{
  m_gridExtractor.exposeParameters(moduleParamList, prefixed("Grid", prefix));
  m_finalExtractor.exposeParameters(moduleParamList, prefixed("Refiner", prefix));
 
  moduleParamList->getParameter<unsigned int>("GridIterations").setDefaultValue(1);
  moduleParamList->getParameter<bool>("RefinerUseLastEventT0").setDefaultValue(true);
 
  moduleParamList->addParameter(prefixed(prefix, "skipIfSVDEventT0Present"),
                                m_skipIfSVDEventT0Present,
                                "Skip execution of the time consuming FullGridChi2 algorithm if a valid SVD based EventT0 estimate is present. " \
                                "If set to true, this module is only used to get an CDC based EventT0 in the few events where no sufficient SVD information is available to estimate EventT0.",
                                m_skipIfSVDEventT0Present);
}

void FullGridChi2TrackTimeExtractor::apply(std::vector<RecoTrack*>& recoTracks)
{
  // Don't execute this time consuming algorithm if a valid SVD based EventT0 estimate is present
  if (m_skipIfSVDEventT0Present and m_eventT0->hasTemporaryEventT0(Const::EDetector::SVD)) {
    return;
  }
 
  const auto& temporaryT0Extractions = m_eventT0->getTemporaryEventT0s(Const::CDC);
  B2ASSERT("There should only be a single or none extraction at this stage!", temporaryT0Extractions.size() <= 1);
 
  if (temporaryT0Extractions.empty()) {
    m_gridExtractor.apply(recoTracks);
    // if neither the hit based nor the grid based extraction finds something, we can not do anything
    if (not m_gridExtractor.wasSuccessful()) {
      m_wasSuccessful = false;
      return;
    }
  } else {
    // Make sure to actually use the temporary hit based extraction result
    m_eventT0->setEventT0(temporaryT0Extractions[0]);
  }
 
  m_finalExtractor.apply(recoTracks);
  m_wasSuccessful = m_finalExtractor.wasSuccessful();
 
  if (not wasSuccessful()) {
    B2DEBUG(25, "Resetting the event t0 as the final extraction was not successful.");
    resetEventT0();
  }
}