CscSegmentMaker.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRIGL2MUONSA_CSCSEGMENTMAKER_H
#define TRIGL2MUONSA_CSCSEGMENTMAKER_H
 
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ServiceHandle.h"
 
#include "GeoPrimitives/GeoPrimitives.h"
#include "MuonIdHelpers/IMuonIdHelperSvc.h"
#include "StoreGate/ReadCondHandleKey.h"
#include "MuonReadoutGeometry/MuonDetectorManager.h"
 
#include "TrackData.h"
#include "TgcFitResult.h"
#include "MuonRoad.h"
#include "CscData.h"
#include "CscRegUtils.h"
 
#include <string>
#include <vector>
 
namespace TrigL2MuonSA{
 
  class CscSegment;
 
  class CscSegmentMaker: public AthAlgTool{
 
 
  public:
 
 
    typedef struct{
      double loc_x;
      double loc_y;
      double loc_z;
      double error;
      double residual;
      int measphi;
      bool enabled;
      int stationname;
      bool isIP;
    } localCscHit;
 
    typedef struct{
      double a;
      double b;
      double zshift;
      double chi2;
      int nhit;
      int stationname;
      double residual;
      std::vector<localCscHit> localHits;
    }local2dSegment;
 
  public:
    CscSegmentMaker(const std::string& type, const std::string &name, const IInterface* parent);
 
 
    virtual StatusCode initialize() override;
 
 
    ReturnCode  FindSuperPointCsc( const TrigL2MuonSA::CscHits &cscHits, 
                                   std::vector<TrigL2MuonSA::TrackPattern> &v_trackPatterns, 
                                   const TrigL2MuonSA::TgcFitResult &tgcFitResult, 
                                   const TrigL2MuonSA::MuonRoad &muroad) const;
 
 
    ReturnCode make_segment(int mod_hash, TrigL2MuonSA::CscHits clusters[8], CscSegment &cscsegment, CscSegment &cscsegment_noip, const MuonGM::MuonDetectorManager* muDetMgr) const;
 
    ReturnCode make_2dsegment(int measphi, const localCscHit &ip_loc,const std::vector<localCscHit> hits_loc[4], local2dSegment &seg2d_eta,local2dSegment &local2d_noip, int &nhite) const;
 
    ReturnCode make_2dseg4hit(int measphi, const localCscHit &ip_loc,std::vector<localCscHit> hits_loc[4], std::vector<local2dSegment> &seg2d_4hitCollection, int &nhite) const;
 
    ReturnCode make_2dseg3hit(int measphi, const localCscHit &ip_loc, const std::vector<localCscHit> hits_loc[4], std::vector<local2dSegment> &seg2d_4hitCollection, int &nhit) const;
    ReturnCode fit_clusters(int measphi, const std::vector<localCscHit> &hits_fit, local2dSegment &seg2d) const;
 
    ReturnCode make_4dsegment(const local2dSegment &seg2d_eta, const local2dSegment &seg2d_phi,
                  Amg::Vector3D &seg_pos, Amg::Vector3D &seg_dir) const;
 
    ReturnCode getModuleSP(int mod_hash[2], const TrigL2MuonSA::TgcFitResult &tgcFitResult, int phibin, const TrigL2MuonSA::MuonRoad &muroad, const int exist_clusters[32]) const;
 
    ReturnCode display_hits(const std::vector<localCscHit> localHits[4]) const;
 
    CscSegment segmentAtFirstLayer(int mod_hash, TrigL2MuonSA::CscSegment *mu_seg) const;
 
  private:
    UtilTools m_util;
    ToolHandle<CscRegDict> m_cscregdict {
      this, "CscRegDict", "TrigL2MuonSA::CscRegDict", ""};
    ServiceHandle<Muon::IMuonIdHelperSvc> m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};
    SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> m_muDetMgrKey {this, "DetectorManagerKey", "MuonDetectorManager", "Key of input MuonDetectorManager condition data"}; 
 
    //properties
    Gaudi::Property< bool > m_use_geometry {
      this, "UseGeometry", false, ""};
    Gaudi::Property< double > m_max_chisquare {
      this, "max_chisquare", 25., ""};
    Gaudi::Property< double > m_max_residual_eta {
      this, "max_residual_eta", 100., ""};
    Gaudi::Property< double > m_max_residual_phi {
      this, "max_residual_phi", 250., ""};
    Gaudi::Property< double > m_err_eta {
      this, "err_eta", 3., ""};
    Gaudi::Property< double > m_err_phi {
      this, "err_phi", 6., ""};
    Gaudi::Property< double > m_err_ip {
      this, "err_ip", 250., ""};
 
 
  };
 
  class CscSegment{
 
  public:
 
    CscSegment();
 
    double x(){ return m_x; }
    double y(){ return m_y; }
    double z(){ return m_z; }
    double px(){ return m_px; }
    double py(){ return m_py; }
    double pz(){ return m_pz; }
    double slopeRZ(){ return m_slopeRZ; }
    double interceptRZ(){ return m_interceptRZ; }
    double chiSquare(){ return m_chisquare; }
    double chiSquarePhi(){ return m_chisquare_phi; }
 
    int nHitEta(){ return m_nhit_eta; }
    void setNHitEta( int nhite){ m_nhit_eta = nhite; }
    int nHitPhi(){ return m_nhit_phi; }
    void setNHitPhi( int nhitp){ m_nhit_phi = nhitp; }
    bool isClean();
 
    ReturnCode set(double x, double y, double z, double px, double py, double pz, double chisquare, double chisquare_phi);
    ReturnCode set(Amg::Vector3D &seg_pos, Amg::Vector3D &seg_dir, double chisquare, double chisquare_phi);
 
  private:
    double m_x, m_y, m_z, m_px, m_py, m_pz;
    double m_slopeRZ, m_interceptRZ;
    int m_nhit_eta, m_nhit_phi;
    double m_chisquare;    
    double m_chisquare_phi;
 
  };
 
 
}//namespace TrigL2MuonSA
 
 
#endif /* TRIGL2MUONSA_CSCSEGMENTMAKER_H */