VXDDedxPIDModule.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/Module.h>
#include <framework/gearbox/Const.h>
#include <framework/datastore/StoreArray.h>
 
#include <reconstruction/dataobjects/VXDDedxTrack.h>
#include <reconstruction/dataobjects/VXDDedxLikelihood.h>
#include <reconstruction/dataobjects/DedxConstants.h>
 
#include <mdst/dataobjects/Track.h>
#include <mdst/dataobjects/MCParticle.h>
#include <tracking/dataobjects/RecoTrack.h>
#include <svd/dataobjects/SVDCluster.h>
#include <pxd/dataobjects/PXDCluster.h>
 
#include <framework/database/DBObjPtr.h>
#include <svd/dbobjects/SVDdEdxPDFs.h>
#include <pxd/dbobjects/PXDdEdxPDFs.h>
 
#include <vector>
 
 
namespace Belle2 {
  class HelixHelper;
 
  class VXDDedxPIDModule : public Module {
 
  public:
 
    VXDDedxPIDModule();
 
    virtual ~VXDDedxPIDModule();
 
    virtual void initialize() override;
 
    virtual void event() override;
 
    virtual void terminate() override;
 
  private:
 
    // required input
    StoreArray<Track> m_tracks; 
    StoreArray<RecoTrack> m_recoTracks; 
    // optional input
    StoreArray<MCParticle> m_mcparticles; 
    StoreArray<SVDCluster> m_svdClusters; 
    StoreArray<PXDCluster> m_pxdClusters; 
    // output
    StoreArray<VXDDedxTrack> m_dedxTracks; 
    StoreArray<VXDDedxLikelihood> m_dedxLikelihoods; 
    int m_eventID;
    int m_trackID;
 
    void calculateMeans(double* mean, double* truncatedMean, double* truncatedMeanErr, const std::vector<double>& dedx) const;
 
    static double getTraversedLength(const PXDCluster* hit, const HelixHelper* helix);
 
    static double getTraversedLength(const SVDCluster* hit, const HelixHelper* helix);
 
    template <class HitClass> void saveSiHits(VXDDedxTrack* track, const HelixHelper& helix, const std::vector<HitClass*>& hits) const;
 
    void savePXDLogLikelihood(double(&logl)[Const::ChargedStable::c_SetSize], double p, float dedx) const;
 
    void saveSVDLogLikelihood(double(&logl)[Const::ChargedStable::c_SetSize], double p, float dedx) const;
 
    bool detectorEnabled(Dedx::Detector d) const
    {
      return (d == Dedx::c_PXD and m_usePXD) or (d == Dedx::c_SVD and m_useSVD);
    }
 
    void checkPDFs();
 
    // pdfs for PID
    DBObjPtr<SVDdEdxPDFs> m_SVDDedxPDFs; 
    DBObjPtr<PXDdEdxPDFs> m_PXDDedxPDFs; 
    // parameters: full likelihood vs. truncated mean
    bool m_useIndividualHits; 
    //no need to define highest and lowest truncated value as we laways remove highest 2 dE/dx values from 8 dE/dx value.
 
    //parameters: technical stuff
    double m_trackDistanceThreshhold; 
    //parameters: which particles and detectors to use
    bool m_onlyPrimaryParticles; 
    bool m_usePXD; 
    bool m_useSVD; 
    //parameters: PDF configuration
    bool m_ignoreMissingParticles; 
  };
} // Belle2 namespace