PXDDQMEfficiencyNtupleSelftrackModule.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Bjoern Spruck                                            *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
 
#pragma once
 
#include <framework/core/Module.h>
#include <framework/datastore/StoreArray.h>
 
#include <pxd/dataobjects/PXDCluster.h>
 
#include <vxd/geometry/GeoCache.h>
#include <vxd/geometry/SensorInfoBase.h>
#include <mdst/dataobjects/Track.h>
#include <tracking/dataobjects/RecoTrack.h>
#include <tracking/dataobjects/ROIid.h>
 
 
//root stuff
#include "TVector3.h"
 
#include "TFile.h"
#include "TNtuple.h"
 
 
namespace Belle2 {
  class PXDDQMEfficiencyNtupleSelftrackModule : public Module {
 
  public:
 
    PXDDQMEfficiencyNtupleSelftrackModule();
 
  private:
 
    void event() override final;
 
    void initialize() override final;
 
    void terminate() override final;
 
 
  private:
    /* returns the space point in local coordinates where the track hits the sensor:
      sensorInfo: info of the sensor under investigation
      aTrack: the track to be tested
      isgood: flag which is false if some error occured (do not use the point if false)
      du and dv are the uncertainties in u and v on the sensor plane of the fit (local coordinates)
     */
    TVector3 getTrackInterSec(const VXD::SensorInfoBase& pxdSensorInfo, const RecoTrack& aTrack, bool& isgood, double& du, double& dv);
    int findClosestCluster(const VxdID& vxdid, TVector3 intersection);
    bool isCloseToBorder(int u, int v, int checkDistance);
    bool isDeadPixelClose(int u, int v, int checkDistance, const VxdID& moduleID);
 
    bool m_useAlignment;
 
    VXD::GeoCache& m_vxdGeometry;
 
    std::string m_pxdClustersName; 
    std::string m_tracksName; 
    std::string m_recoTracksName; 
    std::string m_ROIsName; 
 
    StoreArray<PXDCluster> m_pxdclusters; 
    StoreArray<Track> m_tracks; 
    StoreArray<RecoTrack> m_recoTracks; 
    StoreArray<ROIid> m_ROIs; 
 
    double m_pcut; 
    double m_momCut; 
    double m_pTCut; 
    unsigned int m_minSVDHits; 
    int m_maskedDistance; 
 
    TFile* m_file{}; 
    TNtuple* m_tuple{}; 
 
  };
}