TCConvertersTestModule.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2014  Belle II Collaboration                              *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Thomas Madlener                                          *
 *                                                                        *
 **************************************************************************/
 
#pragma once
 
#include <framework/core/Module.h>
 
#include <genfit/TrackCand.h>
#include <tracking/spacePointCreation/SpacePointTrackCand.h>
 
#include <string>
#include <array>
#include <vector>
#include <tuple>
#include <algorithm> // for find_if
 
namespace Belle2 {
  class TCConvertersTestModule : public Module {
 
  public:
 
    TCConvertersTestModule(); 
    void initialize() override; /*<< initialize: check if all required StoreArrays are present, initialize counters, etc. */
 
    void event() override; 
    void terminate() override; 
    typedef std::tuple<int, int, int, double> trackCandHit;
 
  protected:
 
    std::string m_PXDClusterName; 
    std::string m_SVDClusterName; 
    std::vector<std::string> m_genfitTCNames; 
    std::string m_SPTCName; 
    std::vector<std::string> m_SpacePointArrayNames; 
    // counter variables
 
    int m_genfitTCCtr; 
    int m_convertedTCCtr; 
    int m_SpacePointTCCtr; 
    int m_failedNoSPTC;
 
    int m_failedNoGFTC;
 
    int m_failedOther; 
    int m_failedNotSameHits;
 
    int m_failedWrongOrder; 
    int m_failedWrongSortingParam; 
    int m_failedNoRelationOrig; 
    int m_failedNoRelationConv; 
    unsigned int m_lessHitsCtr; 
    unsigned int m_moreHitsCtr; 
    unsigned int m_differButOKCtr; 
    void initializeCounters(); 
    bool analyzeMisMatch(const genfit::TrackCand* origTC, const genfit::TrackCand* convTC, const Belle2::SpacePointTrackCand* spTC);
 
    std::vector<trackCandHit> getTrackCandHits(const genfit::TrackCand* trackCand);
 
    std::array<bool, 4> checkHits(const std::vector<trackCandHit>& origHits, const std::vector<trackCandHit>& convHits);
 
    template<int p1, int p2, int p3 = 0>
    std::vector<trackCandHit>::const_iterator checkEntries(const std::vector<trackCandHit>& hits, const trackCandHit& hit)
    {
      return std::find_if(hits.begin(), hits.end(),
      [&hit](const trackCandHit & oHit) {
        return std::get<p1>(hit) == std::get<p1>(oHit) &&
               std::get<p2>(hit) == std::get<p2>(oHit) &&
               std::get<p3>(hit) == std::get<p3>(oHit);
      }
                         );
    }
  };
}