SVDDQMEfficiencyModule.h

/**************************************************************************
 * 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.                  *
 **************************************************************************/
 
 
#pragma once
 
#include <framework/core/HistoModule.h>
#include <framework/datastore/StoreArray.h>
 
#include <svd/dataobjects/SVDCluster.h>
#include <svd/dataobjects/SVDSummaryPlots.h>
 
#include <vxd/geometry/GeoCache.h>
#include <vxd/dataobjects/VxdID.h>
 
#include <tracking/dataobjects/SVDIntercept.h>
#include <tracking/dataobjects/RecoTrack.h>
 
#include "TH2D.h"
 
namespace Belle2 {
  class SVDDQMEfficiencyModule : public HistoModule {
 
  public:
 
    SVDDQMEfficiencyModule();
 
  private:
 
    void event() override final;
 
    void initialize() override final;
 
    void defineHisto() override final;
 
 
    bool isGoodIntercept(SVDIntercept* inter); 
    VXD::GeoCache& m_geoCache; 
    std::string m_svdClustersName; 
    std::string m_interceptsName; 
    StoreArray<SVDCluster> m_svdClusters; 
    StoreArray<SVDIntercept> m_intercepts; 
    StoreArray<RecoTrack> m_recoTracks; 
    float m_fiducialU; 
    float m_fiducialV; 
    float m_maxResidU; 
    float m_maxResidV; 
    //    double m_d0cut; /**<d0-Cut for tracks */
    //    double m_z0cut; /**<z0-Cut for tracks */
    double m_pcut; 
    double m_momCut; 
    double m_ptCut; 
    unsigned int m_minSVDHits; 
    unsigned int m_minCDCHits; 
    /* Histograms to later determine efficiency */
 
    std::string m_histogramDirectoryName; 
    bool m_3Samples = false;  
    bool m_saveExpertHistos; 
    int m_u_bins; 
    int m_v_bins; 
    std::map<VxdID, TH2D*> m_h_track_hits; 
    std::map<VxdID, TH2D*> m_h_matched_clusterU; 
    std::map<VxdID, TH2D*> m_h_matched_clusterV; 
    std::map<VxdID, TH2D*> m_h_matched3_clusterU; 
    std::map<VxdID, TH2D*> m_h_matched3_clusterV; 
    std::map<VxdID, TH2D*> m_h_matched6_clusterU; 
    std::map<VxdID, TH2D*> m_h_matched6_clusterV; 
    SVDSummaryPlots* m_TrackHits = nullptr; 
    SVDSummaryPlots* m_TrackHits3Sample = nullptr; 
    SVDSummaryPlots* m_TrackHits6Sample = nullptr; 
    SVDSummaryPlots* m_MatchedHits = nullptr; 
    SVDSummaryPlots* m_MatchedHits3Sample = nullptr; 
    SVDSummaryPlots* m_MatchedHits6Sample = nullptr; 
  };
}