VertexFitterSimple.h

Go to the documentation of this file.

#ifndef  VERTEXFITTERSIMPLE_ANALYZERS_H
#define  VERTEXFITTERSIMPLE_ANALYZERS_H
 
#include <cmath>
#include <vector>
 
#include "ROOT/RVec.hxx"
#include "edm4hep/ReconstructedParticleData.h"
#include "edm4hep/TrackState.h"
 
#include "TVectorD.h"
#include "TVector3.h"
#include "TMatrixDSym.h"
#include "TMath.h"
#include "TDecompChol.h"
#include "TMatrixD.h"
 
#include "ReconstructedParticle2Track.h"
#include "ReconstructedParticle2MC.h"
#include "VertexingUtils.h"
 
#include "edm4hep/VertexData.h"
#include "edm4hep/Vertex.h"
 
#include "TrackCovariance/VertexFit.h" // from Delphes - updates Franco, Jul 2022
#include "TrackCovariance/VertexMore.h"
 
namespace FCCAnalyses{

namespace VertexFitterSimple{

VertexingUtils::FCCAnalysesVertex VertexFitter(
    int Primary,
    ROOT::VecOps::RVec<edm4hep::ReconstructedParticleData> recoparticles,
    ROOT::VecOps::RVec<edm4hep::TrackState> alltracks,
    bool BeamSpotConstraint = false, double sigmax = 0., double sigmay = 0.,
    double sigmaz = 0., double bsc_x = 0., double bsc_y = 0., double bsc_z = 0.,
    bool ComputeMomentaAtVertex = true);

VertexingUtils::FCCAnalysesVertex
VertexFitter_Tk(int Primary, ROOT::VecOps::RVec<edm4hep::TrackState> tracks,
                bool BeamSpotConstraint = false, double sigmax = 0.,
                double sigmay = 0., double sigmaz = 0., double bsc_x = 0.,
                double bsc_y = 0., double bsc_z = 0.,
                bool ComputeMomentaAtVertex = true);
 
VertexingUtils::FCCAnalysesVertex
VertexFitter_Tk(int Primary, ROOT::VecOps::RVec<edm4hep::TrackState> tracks,
                const ROOT::VecOps::RVec<edm4hep::TrackState> &alltracks,
                bool BeamSpotConstraint = false, double sigmax = 0.,
                double sigmay = 0., double sigmaz = 0., double bsc_x = 0.,
                double bsc_y = 0., double bsc_z = 0.,
                bool ComputeMomentaAtVertex = true);

   ROOT::VecOps::RVec<edm4hep::TrackState> get_PrimaryTracks(           ROOT::VecOps::RVec<edm4hep::TrackState> tracks,
                                                                        bool BeamSpotConstraint,
                                                                        double bsc_sigmax, double bsc_sigmay, double bsc_sigmaz,
                                                                        double bsc_x, double bsc_y, double bsc_z ) ;
 

  ROOT::VecOps::RVec<edm4hep::TrackState>  get_NonPrimaryTracks( ROOT::VecOps::RVec<edm4hep::TrackState> allTracks,
                                                                 ROOT::VecOps::RVec<edm4hep::TrackState> primaryTracks ) ;

   ROOT::VecOps::RVec<bool> IsPrimary_forTracks( ROOT::VecOps::RVec<edm4hep::TrackState> allTracks,
                                                                 ROOT::VecOps::RVec<edm4hep::TrackState> primaryTracks ) ;
 
 
/*
  Double_t FastRv(TVectorD p1, TVectorD p2) ;
  TMatrixDSym RegInv3(TMatrixDSym &Smat0) ;
  TMatrixD Fill_A(TVectorD par, Double_t phi) ;
  TVectorD Fill_a(TVectorD par, Double_t phi) ;
  TVectorD Fill_x0(TVectorD par) ;
  TVectorD Fill_x(TVectorD par, Double_t phi) ;
 
  TVectorD XPtoPar(TVector3 x, TVector3 p, Double_t Q);
  TVector3 ParToP(TVectorD Par);
 
  TVectorD XPtoPar(TVector3 x, TVector3 p, Double_t Q);
  TVector3 ParToP(TVectorD Par);
*/
 
 
 
}//end NS VertexFitterSimple
 
}//end NS FCCAnalyses
#endif