df/d58/KalFitAlg_8h_source.html

//------------------------------------------------------------------------

// Description :

// Main file of the module KalFit in charge of :

// 1/ Refit of the Mdc tracks using Kalman filter

// 2/ Backward filter (smoothing)

// 3/ and also several mass hypothesis, Non unif mag field treatment...

//

//------------------------------------------------------------------------


#ifndef _DEFINE_KALFITALG_H_

#define _DEFINE_KALFITALG_H_

#ifndef DBL_MAX

#define DBL_MAX 9999

#endif


class KalFitTrack;

class KalFitHitMdc;

class KalFitHelixSeg;

class Bfield;

//class Helix;


#include "GaudiKernel/Algorithm.h"

#include "MdcCalibFunSvc/IMdcCalibFunSvc.h"

#include "MdcCalibFunSvc/MdcCalibFunSvc.h"

#include "GaudiKernel/IMagneticFieldSvc.h"

#include "MagneticField/MagneticFieldSvc.h"

#include "KalFitAlg/KalFitWire.h"

#include "KalFitAlg/KalFitHelixSeg.h"

#include "KalFitAlg/KalFitCylinder.h"

#include "KalFitAlg/KalFitLayer_Mdc.h"

#include "KalFitAlg/KalFitSuper_Mdc.h"

#include "KalFitAlg/KalFitMaterial.h"

#include "KalFitAlg/KalFitGemHit.h"

#include "MdcTables/MdcTables.h"

#include "GaudiKernel/NTuple.h"

#include "AIDA/IHistogram1D.h"

#include "AIDA/IHistogramFactory.h"

#include "GaudiKernel/IHistogramSvc.h"

#include "MdcRecEvent/RecMdcKalTrack.h"

#include "MdcRecEvent/RecMdcKalHelixSeg.h"

#include "HepPDT/ParticleDataTable.hh"

#include <vector>

//added header files for 6.0.0

#include "CLHEP/Matrix/Vector.h"

#include "CLHEP/Matrix/Matrix.h"

#include "CLHEP/Matrix/SymMatrix.h"

#include "CLHEP/Vector/ThreeVector.h"

#include "G4Material.hh"

#include "G4Tubs.hh"

#include "CgemGeomSvc/ICgemGeomSvc.h"

#include "CgemGeomSvc/CgemGeomSvc.h"


#ifndef ENABLE_BACKWARDS_COMPATIBILITY

typedef HepGeom::Point3D <double > HepPoint3D;

#endif


using namespace std;

using namespace CLHEP;


class KalFitAlg : public Algorithm {


    public:

        /**@name Main member functions*/

        //@{

        /// constructor

        KalFitAlg(const std::string& name, ISvcLocator* pSvcLocator);

        /// destructor

        ~KalFitAlg(void);

        /// initialize

        StatusCode initialize();

        /// event function

        StatusCode execute();

        StatusCode finalize();

        StatusCode beginRun();

        /// hist definition

        void hist_def ( void );

        /// delete C++ objects, necessary to clean before begin_run or inside

        /// destructor

        void clean(void);

        //@}

        /**@name Set up the description of the Mdc */

        //@{

        /// Set up the wires, layers and superlayers...

        void set_Mdc(void);

        /// Initialize the material for Mdc

        void setMaterial_Mdc(void);

        /// Initialize the cylinders (walls and cathodes) for Mdc

        void setCylinder_Mdc(void);

        ///

        void setDchisqCut(void);

        /// initialize for the services

        void setCalibSvc_init(void);

        void setGeomSvc_init(void);

        void setMagneticFieldSvc_init(void);

        ///

        void getEventStarTime(void);


        /**@name Kalman filter method related member functions*/

        //@{

        /// Kalman filter (forward) in Mdc

        void filter_fwd_anal(KalFitTrack& trk, int l_mass, int way, HepSymMatrix& Eakal);

        void filter_fwd_calib(KalFitTrack& trk, int l_mass, int way, HepSymMatrix& Eakal);


        void init_matrix(MdcRec_trk& trk, HepSymMatrix& Ea);

        void init_matrix(int k, MdcRec_trk& trk, HepSymMatrix& Ea);


        void start_seed(KalFitTrack& track, int lead_, int way, MdcRec_trk& trk);


        /// Kalman filter (smoothing or backward part)

        void smoother_anal(KalFitTrack& trk, int way);

        void smoother_calib(KalFitTrack& trk, int way);


        /// Take the inner walls (eloss and mult scat) into account

        void innerwall(KalFitTrack& trk, int l_mass, int way);

        void innerwall(KalFitTrack& trk, int l_mass, int way,int begin,int end);

        void innerwall(KalFitTrack& track, int l_mass, int way, double r1, double r2, int& index);//add by Huang Zhen

        //@{

        /// with results got at the inner Mdc hit

        void fillTds(MdcRec_trk& TrasanTRK, KalFitTrack& track,

                RecMdcKalTrack* trk,int l_mass);

        void fillTds_lead(MdcRec_trk& TrasanTRK, KalFitTrack& track,

                RecMdcKalTrack* trk, int l_mass);

        ///  with results got at the outer Mdc hit

        void fillTds_back(KalFitTrack& track, RecMdcKalTrack* trk,

                MdcRec_trk& TrasanTRK,int l_mass );


        void fillTds_back(KalFitTrack& track, RecMdcKalTrack* trk,

                MdcRec_trk& TrasanTRK,int l_mass,RecMdcKalHelixSegCol*  segcol);


        ///for smoother process

        void fillTds_back(KalFitTrack& track, RecMdcKalTrack* trk,

                MdcRec_trk& TrasanTRK,int l_mass,RecMdcKalHelixSegCol*  segcol, int smoothflag);

        ///  with results got at (0,0,0)

        void fillTds_ip(MdcRec_trk& TrasanTRK, KalFitTrack& track,

                RecMdcKalTrack* trk, int l_mass);


        /// complete the RecMdcKalTrackCol

        void sameas(RecMdcKalTrack* trk, int l_mass, int imain);


        void complete_track(MdcRec_trk& TrasanTRK,

                MdcRec_trk_add& TrasanTRK_add,

                KalFitTrack& track_lead,

                RecMdcKalTrack*  kaltrk,

                RecMdcKalTrackCol* kalcol,RecMdcKalHelixSegCol* segcol,int flagsmooth);


        void complete_track(MdcRec_trk& TrasanTRK,

                MdcRec_trk_add& TrasanTRK_add,

                KalFitTrack& track_lead,

                RecMdcKalTrack*  kaltrk,

                RecMdcKalTrackCol* kalcol,RecMdcKalHelixSegCol* segcol);


        // Careful refit

        void kalman_fitting_anal(void);

        void kalman_fitting_calib(void);

        void kalman_fitting_csmalign(void);

        void kalman_fitting_MdcxReco_Csmc_Sew(void);


        // clear tables by wangdy

        void clearTables( );


        ///

        int getWallMdcNumber(const HepPoint3D& point);


        ///

        void extToAnyPoint(KalFitTrack& trk, const HepPoint3D& point);


        ///

        void setBesFromGdml(void);


        // --- about GEM

        void makeGemHitsCol();

        void fromFHitToInnerWall(KalFitTrack& track, int way);

        void fitGemHits(KalFitTrack& track, int hypo, int way);

        KalFitTrack Cgem_filter_anal(KalFitTrack& track, int hypo, int way);//add by Huang Zhen

        void residual(KalFitTrack& track);//add by Huang Zhen


        /// this usage is used to control the usage of this algorithm ,to be

        /// analysis  or calibration.

        int usage_;

        ///

        int choice_;

        //// flag to calculate path length in each Mdc laye :

        int pathl_;

        /// flag to take account the wire sag into account

        int wsag_;

        /// flag to perform smoothing

        int back_;


        int eventno;


        int Tds_back_no;

        /// value of the pT cut for backward filter

        double pT_;

        /// leading mass assumption

        int lead_;

        ///

        int mhyp_;

        /// value of the momentum cut to decide refit

        double pe_cut_, pmu_cut_, ppi_cut_, pk_cut_, pp_cut_;

        double pt_cut_,theta_cut_;


        /// Flag account to multiple scattering and energy loss,

        /// where lr flag from and whether use active hits only

        int muls_, loss_, lr_, activeonly_, tofflag_, tof_hyp_;

        /// flag to enhance the error matrix at the inner hit of Mdc (cosmic)

        int enhance_;

        double fac_h1_, fac_h2_, fac_h3_, fac_h4_, fac_h5_, matrixg_;

        ///

        double  gain1_,gain2_,gain3_,gain4_,gain5_;

        int steplev_;

        int numfcor_;

        int numf_;

        int inner_steps_;

        int outer_steps_;

        int numf_in_;

        int numf_out_;

        int fitnocut_;

        ///

        int drifttime_choice_;

        /// mass assumption for backward filter (if <0 means use leading mass)

        int i_back_;

        int i_front_;

        /// Debug flag for the track finder part

        int debug_ , debug_kft_;

        /// Fill ntuples of KalFit

        int ntuple_;

        // dir of files

        string matfile_, cylfile_;

        //cut for delta_chi2

        double dchi2cutf_, dchi2cuts_;


        double  dchi2cut_mid1_,dchi2cut_mid2_,dchi2cut_inner_, dchi2cut_outer_,dchi2cut_layid2_,dchi2cut_layid3_;


        /// factor of energy loss straggling for electron

        double fstrag_;

        //wire resoltion flag

        int resolution_;

        int iqual_front_[5],iqual_back_;


        /// propagation correction

        int tprop_;

        int myEventNo;

        int m_usevtxdb;

        int m_csmflag;    //cosmic events flag; for cosmic events tof(y>0) should be minus

        double m_dangcut, m_dphicut;  //for cosmic events cut

        double m_alpha;

    private:

        // --- KalFitCylinders

        std::vector<KalFitCylinder> _BesKalmanFitWalls;

        std::vector<KalFitMaterial> _BesKalmanFitMaterials;

        std::vector<KalFitCylinder> m_innerWall;

        std::vector<KalFitCylinder> m_CGEM;

        std::vector<KalFitMaterial> m_CGEMMaterials;

        //std::vector<G4Tubs>         _BesKalmanFitTubs;


        // --- statistics

        int nTotalTrks;

        int nFailedTrks[5];


        // --- GEM resolution

        double myRphiResGem; // in cm

        double myZResGem;    // in cm


        // --- GEM hit coordinate choice

        bool myUseGemPos3D;


        KalFitWire * _wire;

        KalFitLayer_Mdc * _layer;

        KalFitSuper_Mdc * _superLayer;

        HepPDT::ParticleDataTable* m_particleTable;

        static const double RIW;


        const MdcCalibFunSvc* m_mdcCalibFunSvc_;

        const IMagneticFieldSvc* m_MFSvc_;

        static IMdcGeomSvc* imdcGeomSvc_;

        const CgemGeomSvc* m_CgemGeomSvc_;


        int useNCGem_;

        int nInnerCFiber_;

        KalFitGemHitCol myGemHitCol[4];


        int testCDC;


        //sort the rec hits by layer

        static bool  order_rechits(const SmartRef<RecMdcHit>& m1, const SmartRef<RecMdcHit>& m2);


        //ntuples

        NTuple::Tuple* m_nt1;  // KalFit track params

        NTuple::Tuple* m_nt2;  // KalFit 2-prong comparison

        NTuple::Tuple* m_nt3;  // PatRec track params

        NTuple::Tuple* m_nt4;  // PatRec 2-prong comparison

        NTuple::Tuple* m_nt5;  // for hit checking and cut

        NTuple::Tuple* m_nt6;  // for helix seg of calib

        NTuple::Tuple* m_nt7;  // for kalmanfit failure  --Long LI


        //for nt1

        NTuple::Item<long>  m_trackid,m_evtid;

        NTuple::Item<double>  m_chi2direct,m_prob;

        NTuple::Matrix<double> m_ndf,m_chisq,m_stat;

        NTuple::Array<double> m_length,m_tof,m_nhits;

        NTuple::Item<double> m_zptot,m_zptote,m_zptotmu,m_zptotk,m_zptotp;

        NTuple::Item<double> m_zpt,m_zpte,m_zptmu,m_zptk,m_zptp;

        NTuple::Item<double> m_fptot,m_fptote,m_fptotmu,m_fptotk,m_fptotp;

        NTuple::Item<double> m_fpt,m_fpte,m_fptmu,m_fptk,m_fptp;

        NTuple::Item<double> m_lptot,m_lptote,m_lptotmu,m_lptotk,m_lptotp;

        NTuple::Item<double> m_lpt,m_lpte,m_lptmu,m_lptk,m_lptp;

        NTuple::Item<double> m_zsigp,m_zsigpe,m_zsigpmu,m_zsigpk,m_zsigpp;

        NTuple::Array<double> m_zhelix,m_zhelixe,m_zhelixmu,m_zhelixk,m_zhelixp;

        NTuple::Array<double> m_fhelix,m_fhelixe,m_fhelixmu,m_fhelixk,m_fhelixp;

        NTuple::Array<double> m_lhelix,m_lhelixe,m_lhelixmu,m_lhelixk,m_lhelixp;

        NTuple::Array<double> m_zerror,m_zerrore,m_zerrormu,m_zerrork,m_zerrorp;

        NTuple::Array<double> m_ferror,m_ferrore,m_ferrormu,m_ferrork,m_ferrorp;

        NTuple::Array<double> m_lerror,m_lerrore,m_lerrormu,m_lerrork,m_lerrorp;

        NTuple::Array<double> m_rGem, m_chi2Gem, m_phiGemExp, m_phiGemHit, m_zGemExp, m_zGemHit;

        //for nt1 single track MCTruth

        NTuple::Array<double> m_mchelix;

        NTuple::Item<double> m_mcptot;

        NTuple::Item<long> m_mcpid;

        //for nt3

        NTuple::Array<double> m_trkhelix, m_trkerror;

        NTuple::Item<double> m_trkndf, m_trkchisq, m_trkptot, m_trksigp;

        //for nt2

        NTuple::Item<double> m_delx,m_dely,m_delz,m_delthe,m_delphi,m_delp;

        NTuple::Item<double> m_delpx,m_delpy,m_delpz;


        //for nt4

        NTuple::Item<double> m_trkdelx,m_trkdely,m_trkdelz;

        NTuple::Item<double> m_trkdelthe,m_trkdelphi,m_trkdelp;

        //for nt5

        NTuple::Item<double>  m_dchi2,m_orichi2,m_fitchi2,m_residest, m_residnew,m_anal_dr, m_anal_phi0, m_anal_kappa, m_anal_dz, m_anal_tanl, m_anal_ea_dr, m_anal_ea_phi0, m_anal_ea_kappa, m_anal_ea_dz, m_anal_ea_tanl;

        NTuple::Item<long> m_masshyp, m_layer;

        //for nt6

        NTuple::Item<double> m_docaInc,m_docaExc, m_tdrift;

        NTuple::Item<long>  m_layerid,m_eventNo;

        NTuple::Item<double> m_residualInc, m_residualExc, m_lr, m_yposition, m_dd;


        //for nt7 --Long LI

        NTuple::Item<int> m_run_kal;

        NTuple::Item<int> m_event_kal;

        NTuple::Item<int> m_trkid_kal;

        NTuple::Item<int> m_dropedHits_kal_e;

        NTuple::Item<double> m_kappa2_kal_e;

        NTuple::Item<int> m_trackNhits_kal_e;

        NTuple::Item<int> m_trackNster_kal_e;

        NTuple::Item<int> m_trackNaxis_kal_e;

        NTuple::Item<double> m_chi2_kal_e;

        NTuple::Item<double> m_Ea00_kal_e;

        NTuple::Item<double> m_Ea11_kal_e;

        NTuple::Item<double> m_Ea22_kal_e;

        NTuple::Item<double> m_Ea33_kal_e;

        NTuple::Item<double> m_Ea44_kal_e;

        NTuple::Item<double> m_iqual_front_kal_e;

        NTuple::Item<int> m_dropedHits_kal_mu;

        NTuple::Item<double> m_kappa2_kal_mu;

        NTuple::Item<int> m_trackNhits_kal_mu;

        NTuple::Item<int> m_trackNster_kal_mu;

        NTuple::Item<int> m_trackNaxis_kal_mu;

        NTuple::Item<double> m_chi2_kal_mu;

        NTuple::Item<double> m_Ea00_kal_mu;

        NTuple::Item<double> m_Ea11_kal_mu;

        NTuple::Item<double> m_Ea22_kal_mu;

        NTuple::Item<double> m_Ea33_kal_mu;

        NTuple::Item<double> m_Ea44_kal_mu;

        NTuple::Item<double> m_iqual_front_kal_mu;

        NTuple::Item<int> m_dropedHits_kal_pi;

        NTuple::Item<double> m_kappa2_kal_pi;

        NTuple::Item<int> m_trackNhits_kal_pi;

        NTuple::Item<int> m_trackNster_kal_pi;

        NTuple::Item<int> m_trackNaxis_kal_pi;

        NTuple::Item<double> m_chi2_kal_pi;

        NTuple::Item<double> m_Ea00_kal_pi;

        NTuple::Item<double> m_Ea11_kal_pi;

        NTuple::Item<double> m_Ea22_kal_pi;

        NTuple::Item<double> m_Ea33_kal_pi;

        NTuple::Item<double> m_Ea44_kal_pi;

        NTuple::Item<double> m_iqual_front_kal_pi;

        NTuple::Item<int> m_dropedHits_kal_k;

        NTuple::Item<double> m_kappa2_kal_k;

        NTuple::Item<int> m_trackNhits_kal_k;

        NTuple::Item<int> m_trackNster_kal_k;

        NTuple::Item<int> m_trackNaxis_kal_k;

        NTuple::Item<double> m_chi2_kal_k;

        NTuple::Item<double> m_Ea00_kal_k;

        NTuple::Item<double> m_Ea11_kal_k;

        NTuple::Item<double> m_Ea22_kal_k;

        NTuple::Item<double> m_Ea33_kal_k;

        NTuple::Item<double> m_Ea44_kal_k;

        NTuple::Item<double> m_iqual_front_kal_k;

        NTuple::Item<int> m_dropedHits_kal_p;

        NTuple::Item<double> m_kappa2_kal_p;

        NTuple::Item<int> m_trackNhits_kal_p;

        NTuple::Item<int> m_trackNster_kal_p;

        NTuple::Item<int> m_trackNaxis_kal_p;

        NTuple::Item<double> m_chi2_kal_p;

        NTuple::Item<double> m_Ea00_kal_p;

        NTuple::Item<double> m_Ea11_kal_p;

        NTuple::Item<double> m_Ea22_kal_p;

        NTuple::Item<double> m_Ea33_kal_p;

        NTuple::Item<double> m_Ea44_kal_p;

        NTuple::Item<double> m_iqual_front_kal_p;


        NTuple::Item<int> m_hit_no;

        NTuple::Item<int> m_nCluster;

        NTuple::Array<int> m_layer_hit_e;

        NTuple::Array<double> m_dchi2_hit_e;

        NTuple::Array<double> m_residest_hit_e;

        NTuple::Array<double> m_anal_dr_hit_e;

        NTuple::Array<double> m_anal_phi0_hit_e;

        NTuple::Array<double> m_anal_kappa_hit_e;

        NTuple::Array<double> m_anal_dz_hit_e;

        NTuple::Array<double> m_anal_tanl_hit_e;

        NTuple::Array<double> m_anal_ea_dr_hit_e;

        NTuple::Array<double> m_anal_ea_phi0_hit_e;

        NTuple::Array<double> m_anal_ea_kappa_hit_e;

        NTuple::Array<double> m_anal_ea_dz_hit_e;

        NTuple::Array<double> m_anal_ea_tanl_hit_e;

        NTuple::Array<double> m_residnew_hit_e;

        NTuple::Array<int> m_layer_hit_mu;

        NTuple::Array<double> m_dchi2_hit_mu;

        NTuple::Array<double> m_residest_hit_mu;

        NTuple::Array<double> m_anal_dr_hit_mu;

        NTuple::Array<double> m_anal_phi0_hit_mu;

        NTuple::Array<double> m_anal_kappa_hit_mu;

        NTuple::Array<double> m_anal_dz_hit_mu;

        NTuple::Array<double> m_anal_tanl_hit_mu;

        NTuple::Array<double> m_anal_ea_dr_hit_mu;

        NTuple::Array<double> m_anal_ea_phi0_hit_mu;

        NTuple::Array<double> m_anal_ea_kappa_hit_mu;

        NTuple::Array<double> m_anal_ea_dz_hit_mu;

        NTuple::Array<double> m_anal_ea_tanl_hit_mu;

        NTuple::Array<double> m_residnew_hit_mu;

        NTuple::Array<int> m_layer_hit_pi;

        NTuple::Array<double> m_dchi2_hit_pi;

        NTuple::Array<double> m_residest_hit_pi;

        NTuple::Array<double> m_anal_dr_hit_pi;

        NTuple::Array<double> m_anal_phi0_hit_pi;

        NTuple::Array<double> m_anal_kappa_hit_pi;

        NTuple::Array<double> m_anal_dz_hit_pi;

        NTuple::Array<double> m_anal_tanl_hit_pi;

        NTuple::Array<double> m_anal_ea_dr_hit_pi;

        NTuple::Array<double> m_anal_ea_phi0_hit_pi;

        NTuple::Array<double> m_anal_ea_kappa_hit_pi;

        NTuple::Array<double> m_anal_ea_dz_hit_pi;

        NTuple::Array<double> m_anal_ea_tanl_hit_pi;

        NTuple::Array<double> m_residnew_hit_pi;

        NTuple::Array<int> m_layer_hit_k;

        NTuple::Array<double> m_dchi2_hit_k;

        NTuple::Array<double> m_residest_hit_k;

        NTuple::Array<double> m_anal_dr_hit_k;

        NTuple::Array<double> m_anal_phi0_hit_k;

        NTuple::Array<double> m_anal_kappa_hit_k;

        NTuple::Array<double> m_anal_dz_hit_k;

        NTuple::Array<double> m_anal_tanl_hit_k;

        NTuple::Array<double> m_anal_ea_dr_hit_k;

        NTuple::Array<double> m_anal_ea_phi0_hit_k;

        NTuple::Array<double> m_anal_ea_kappa_hit_k;

        NTuple::Array<double> m_anal_ea_dz_hit_k;

        NTuple::Array<double> m_anal_ea_tanl_hit_k;

        NTuple::Array<double> m_residnew_hit_k;

        NTuple::Array<int> m_layer_hit_p;

        NTuple::Array<double> m_dchi2_hit_p;

        NTuple::Array<double> m_residest_hit_p;

        NTuple::Array<double> m_anal_dr_hit_p;

        NTuple::Array<double> m_anal_phi0_hit_p;

        NTuple::Array<double> m_anal_kappa_hit_p;

        NTuple::Array<double> m_anal_dz_hit_p;

        NTuple::Array<double> m_anal_tanl_hit_p;

        NTuple::Array<double> m_anal_ea_dr_hit_p;

        NTuple::Array<double> m_anal_ea_phi0_hit_p;

        NTuple::Array<double> m_anal_ea_kappa_hit_p;

        NTuple::Array<double> m_anal_ea_dz_hit_p;

        NTuple::Array<double> m_anal_ea_tanl_hit_p;

        NTuple::Array<double> m_residnew_hit_p;


        NTuple::Item<double>  m_dchisq0,m_dchisq1,m_dchisq2,m_dchisq3,m_dchisq4,m_dchisq5,m_dchisq6,m_dchisq7,m_dchisq8,m_dchisq9,m_dchisq10,m_dchisq11,m_dchisq12,m_dchisq13,m_dchisq14,m_dchisq15,m_dchisq16,m_dchisq17,m_dchisq18,m_dchisq19,m_dchisq20,m_dchisq21,m_dchisq22,m_dchisq23,m_dchisq24,m_dchisq25,m_dchisq26,m_dchisq27,m_dchisq28,m_dchisq29,m_dchisq30,m_dchisq31,m_dchisq32,m_dchisq33,m_dchisq34,m_dchisq35,m_dchisq36,m_dchisq37,m_dchisq38,m_dchisq39,m_dchisq40,m_dchisq41,m_dchisq42;

        NTuple::Item<double>  m_dtrack0,m_dtrack1,m_dtrack2,m_dtrack3,m_dtrack4,m_dtrack5,m_dtrack6,m_dtrack7,m_dtrack8,m_dtrack9,m_dtrack10,m_dtrack11,m_dtrack12,m_dtrack13,m_dtrack14,m_dtrack15,m_dtrack16,m_dtrack17,m_dtrack18,m_dtrack19,m_dtrack20,m_dtrack21,m_dtrack22,m_dtrack23,m_dtrack24,m_dtrack25,m_dtrack26,m_dtrack27,m_dtrack28,m_dtrack29,m_dtrack30,m_dtrack31,m_dtrack32,m_dtrack33,m_dtrack34,m_dtrack35,m_dtrack36,m_dtrack37,m_dtrack38,m_dtrack39,m_dtrack40,m_dtrack41,m_dtrack42;

        NTuple::Item<double>  m_dtdc0,m_dtdc1,m_dtdc2,m_dtdc3,m_dtdc4,m_dtdc5,m_dtdc6,m_dtdc7,m_dtdc8,m_dtdc9,m_dtdc10,m_dtdc11,m_dtdc12,m_dtdc13,m_dtdc14,m_dtdc15,m_dtdc16,m_dtdc17,m_dtdc18,m_dtdc19,m_dtdc20,m_dtdc21,m_dtdc22,m_dtdc23,m_dtdc24,m_dtdc25,m_dtdc26,m_dtdc27,m_dtdc28,m_dtdc29,m_dtdc30,m_dtdc31,m_dtdc32,m_dtdc33,m_dtdc34,m_dtdc35,m_dtdc36,m_dtdc37,m_dtdc38,m_dtdc39,m_dtdc40,m_dtdc41,m_dtdc42;


        bool ifProdNt10;

        NTuple::Tuple* m_nt10;

        NTuple::Item<long> m_evt3;

        NTuple::Item<long> m_nGemHits,m_qua;

        NTuple::Array<double> m_meas_r,m_meas_phi,m_meas_z,m_esti1_r,m_esti1_phi,m_esti1_z,m_esti2_r,m_esti2_phi    ,m_esti2_z,m_diff1_phi,m_diff2_z,m_diff1_z,m_diff2_phi,m_Gchi2;

        NTuple::Array<int>    m_GemLayer,m_mass;

        NTuple::Array<double>    m_meas_phierr,m_meas_zerr,m_esti_phierr,m_esti_zerr;

        NTuple::Item<int> ctrk;

        NTuple::Array<double> m_phi1,m_z1,m_phi2,m_z2,m_resphi1,m_resphi2,m_resz1,m_resz2,m_achi2;


        bool ifProdNt11;

        NTuple::Tuple* m_nt11;

        NTuple::Item<int> m_evt4,m_ntruth;

        NTuple::Array<double> m_dtphi,m_dtv,m_dtpostphi,m_dtpostz;

        NTuple::Array<int> m_tlayer;


        //calculate residual of 5 helix parameter and momemtum between kalman filter and MCtruth, add by Huang Zhen

        bool ifProdNt12;

        NTuple::Tuple* m_nt12;

        NTuple::Item<long>  m_track,m_evt;

        NTuple::Item<double> m_diff_dr,m_diff_phi0,m_diff_kappa,m_diff_dz,m_diff_tanl,m_diff_p;

};


#endif


CgemGeomSvc.h

ICgemGeomSvc.h

IMdcCalibFunSvc.h

HepPoint3D
HepGeom::Point3D< double > HepPoint3D
Definition KalFitAlg.h:54

KalFitCylinder.h

KalFitGemHit.h

KalFitGemHitCol
vector< KalFitGemHit > KalFitGemHitCol
Definition KalFitGemHit.h:49

KalFitHelixSeg.h

KalFitLayer_Mdc.h

KalFitMaterial.h

KalFitSuper_Mdc.h

KalFitWire.h

MagneticFieldSvc.h

MdcCalibFunSvc.h

MdcTables.h

RecMdcKalHelixSeg.h

RecMdcKalHelixSegCol
ObjectVector< RecMdcKalHelixSeg > RecMdcKalHelixSegCol
Definition RecMdcKalHelixSeg.h:184

RecMdcKalTrack.h

RecMdcKalTrackCol
ObjectVector< RecMdcKalTrack > RecMdcKalTrackCol
Definition RecMdcKalTrack.h:348

Bfield
Definition KalFitAlg/KalFitAlg-00-15-18/KalFitAlg/coil/Bfield.h:43

CgemGeomSvc
Definition CgemGeomSvc.h:25

HepGeom::Point3D< double >

IMagneticFieldSvc
Definition IMagneticFieldSvc.h:29

IMdcGeomSvc
Definition bak_MdcGeomSvc-00-01-39/MdcGeomSvc/IMdcGeomSvc.h:15

KalFitAlg
Definition KalFitAlg.h:61

KalFitAlg::pe_cut_
double pe_cut_
value of the momentum cut to decide refit
Definition KalFitAlg.h:201

KalFitAlg::debug_kft_
int debug_kft_
Definition KalFitAlg.h:227

KalFitAlg::gain1_
double gain1_
Definition KalFitAlg.h:212

KalFitAlg::residual
void residual(KalFitTrack &track)
Definition KalFitAlg.cxx:8091

KalFitAlg::m_alpha
double m_alpha
Definition KalFitAlg.h:249

KalFitAlg::filter_fwd_anal
void filter_fwd_anal(KalFitTrack &trk, int l_mass, int way, HepSymMatrix &Eakal)
Kalman filter (forward) in Mdc.
Definition KalFitAlg.cxx:3191

KalFitAlg::filter_fwd_calib
void filter_fwd_calib(KalFitTrack &trk, int l_mass, int way, HepSymMatrix &Eakal)
Definition KalFitAlg.cxx:3553

KalFitAlg::dchi2cut_mid1_
double dchi2cut_mid1_
Definition KalFitAlg.h:235

KalFitAlg::fstrag_
double fstrag_
factor of energy loss straggling for electron
Definition KalFitAlg.h:238

KalFitAlg::pmu_cut_
double pmu_cut_
Definition KalFitAlg.h:201

KalFitAlg::beginRun
StatusCode beginRun()
Definition KalFitAlg.cxx:342

KalFitAlg::fac_h4_
double fac_h4_
Definition KalFitAlg.h:210

KalFitAlg::m_dangcut
double m_dangcut
Definition KalFitAlg.h:248

KalFitAlg::loss_
int loss_
Definition KalFitAlg.h:207

KalFitAlg::theta_cut_
double theta_cut_
Definition KalFitAlg.h:202

KalFitAlg::pt_cut_
double pt_cut_
Definition KalFitAlg.h:202

KalFitAlg::back_
int back_
flag to perform smoothing
Definition KalFitAlg.h:189

KalFitAlg::gain3_
double gain3_
Definition KalFitAlg.h:212

KalFitAlg::fillTds_lead
void fillTds_lead(MdcRec_trk &TrasanTRK, KalFitTrack &track, RecMdcKalTrack *trk, int l_mass)
Definition KalFitAlg.cxx:1671

KalFitAlg::iqual_front_
int iqual_front_[5]
Definition KalFitAlg.h:241

KalFitAlg::numfcor_
int numfcor_
Definition KalFitAlg.h:214

KalFitAlg::outer_steps_
int outer_steps_
Definition KalFitAlg.h:217

KalFitAlg::eventno
int eventno
Definition KalFitAlg.h:191

KalFitAlg::~KalFitAlg
~KalFitAlg(void)
destructor
Definition KalFitAlg.cxx:182

KalFitAlg::muls_
int muls_
Definition KalFitAlg.h:207

KalFitAlg::fromFHitToInnerWall
void fromFHitToInnerWall(KalFitTrack &track, int way)
Definition KalFitAlg.cxx:6882

KalFitAlg::gain5_
double gain5_
Definition KalFitAlg.h:212

KalFitAlg::dchi2cutf_
double dchi2cutf_
Definition KalFitAlg.h:233

KalFitAlg::clean
void clean(void)
Definition KalFitAlg.cxx:190

KalFitAlg::tof_hyp_
int tof_hyp_
Definition KalFitAlg.h:207

KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew
void kalman_fitting_MdcxReco_Csmc_Sew(void)
Definition KalFitAlg.cxx:5074

KalFitAlg::pp_cut_
double pp_cut_
Definition KalFitAlg.h:201

KalFitAlg::setCalibSvc_init
void setCalibSvc_init(void)
initialize for the services
Definition KalFitAlg2.cxx:247

KalFitAlg::numf_in_
int numf_in_
Definition KalFitAlg.h:218

KalFitAlg::smoother_anal
void smoother_anal(KalFitTrack &trk, int way)
Kalman filter (smoothing or backward part)
Definition KalFitAlg.cxx:2838

KalFitAlg::iqual_back_
int iqual_back_
Definition KalFitAlg.h:241

KalFitAlg::i_back_
int i_back_
mass assumption for backward filter (if <0 means use leading mass)
Definition KalFitAlg.h:224

KalFitAlg::complete_track
void complete_track(MdcRec_trk &TrasanTRK, MdcRec_trk_add &TrasanTRK_add, KalFitTrack &track_lead, RecMdcKalTrack *kaltrk, RecMdcKalTrackCol *kalcol, RecMdcKalHelixSegCol *segcol, int flagsmooth)
Definition KalFitAlg.cxx:6125

KalFitAlg::wsag_
int wsag_
flag to take account the wire sag into account
Definition KalFitAlg.h:187

KalFitAlg::fac_h3_
double fac_h3_
Definition KalFitAlg.h:210

KalFitAlg::dchi2cut_mid2_
double dchi2cut_mid2_
Definition KalFitAlg.h:235

KalFitAlg::kalman_fitting_calib
void kalman_fitting_calib(void)
Definition KalFitAlg.cxx:4485

KalFitAlg::lr_
int lr_
Definition KalFitAlg.h:207

KalFitAlg::sameas
void sameas(RecMdcKalTrack *trk, int l_mass, int imain)
complete the RecMdcKalTrackCol
Definition KalFitAlg.cxx:2819

KalFitAlg::fillTds_ip
void fillTds_ip(MdcRec_trk &TrasanTRK, KalFitTrack &track, RecMdcKalTrack *trk, int l_mass)
with results got at (0,0,0)
Definition KalFitAlg.cxx:1766

KalFitAlg::makeGemHitsCol
void makeGemHitsCol()
Definition KalFitAlg.cxx:7264

KalFitAlg::tofflag_
int tofflag_
Definition KalFitAlg.h:207

KalFitAlg::mhyp_
int mhyp_
Definition KalFitAlg.h:199

KalFitAlg::fitGemHits
void fitGemHits(KalFitTrack &track, int hypo, int way)
Definition KalFitAlg.cxx:7369

KalFitAlg::matfile_
string matfile_
Definition KalFitAlg.h:231

KalFitAlg::fillTds
void fillTds(MdcRec_trk &TrasanTRK, KalFitTrack &track, RecMdcKalTrack *trk, int l_mass)
with results got at the inner Mdc hit
Definition KalFitAlg.cxx:1501

KalFitAlg::innerwall
void innerwall(KalFitTrack &trk, int l_mass, int way)
Take the inner walls (eloss and mult scat) into account.
Definition KalFitAlg.cxx:3751

KalFitAlg::fac_h1_
double fac_h1_
Definition KalFitAlg.h:210

KalFitAlg::m_dphicut
double m_dphicut
Definition KalFitAlg.h:248

KalFitAlg::Cgem_filter_anal
KalFitTrack Cgem_filter_anal(KalFitTrack &track, int hypo, int way)
Definition KalFitAlg.cxx:7777

KalFitAlg::fac_h2_
double fac_h2_
Definition KalFitAlg.h:210

KalFitAlg::fitnocut_
int fitnocut_
Definition KalFitAlg.h:220

KalFitAlg::debug_
int debug_
Debug flag for the track finder part.
Definition KalFitAlg.h:227

KalFitAlg::cylfile_
string cylfile_
Definition KalFitAlg.h:231

KalFitAlg::usage_
int usage_
Definition KalFitAlg.h:181

KalFitAlg::pathl_
int pathl_
Definition KalFitAlg.h:185

KalFitAlg::start_seed
void start_seed(KalFitTrack &track, int lead_, int way, MdcRec_trk &trk)
Definition KalFitAlg.cxx:6646

KalFitAlg::dchi2cut_layid2_
double dchi2cut_layid2_
Definition KalFitAlg.h:235

KalFitAlg::getEventStarTime
void getEventStarTime(void)

KalFitAlg::setBesFromGdml
void setBesFromGdml(void)
Definition KalFitReadGdml.cxx:12

KalFitAlg::setGeomSvc_init
void setGeomSvc_init(void)
Definition KalFitAlg2.cxx:259

KalFitAlg::KalFitAlg
KalFitAlg(const std::string &name, ISvcLocator *pSvcLocator)
constructor
Definition KalFitAlg.cxx:99

KalFitAlg::kalman_fitting_anal
void kalman_fitting_anal(void)
Definition KalFitAlg.cxx:3849

KalFitAlg::i_front_
int i_front_
Definition KalFitAlg.h:225

KalFitAlg::setDchisqCut
void setDchisqCut(void)
Definition KalFitAlg.cxx:815

KalFitAlg::m_usevtxdb
int m_usevtxdb
Definition KalFitAlg.h:246

KalFitAlg::tprop_
int tprop_
propagation correction
Definition KalFitAlg.h:244

KalFitAlg::dchi2cut_layid3_
double dchi2cut_layid3_
Definition KalFitAlg.h:235

KalFitAlg::fac_h5_
double fac_h5_
Definition KalFitAlg.h:210

KalFitAlg::choice_
int choice_
Definition KalFitAlg.h:183

KalFitAlg::getWallMdcNumber
int getWallMdcNumber(const HepPoint3D &point)

KalFitAlg::set_Mdc
void set_Mdc(void)
Set up the wires, layers and superlayers...
Definition KalFitAlg2.cxx:114

KalFitAlg::numf_out_
int numf_out_
Definition KalFitAlg.h:219

KalFitAlg::setMagneticFieldSvc_init
void setMagneticFieldSvc_init(void)

KalFitAlg::drifttime_choice_
int drifttime_choice_
Definition KalFitAlg.h:222

KalFitAlg::m_csmflag
int m_csmflag
Definition KalFitAlg.h:247

KalFitAlg::dchi2cut_outer_
double dchi2cut_outer_
Definition KalFitAlg.h:235

KalFitAlg::enhance_
int enhance_
flag to enhance the error matrix at the inner hit of Mdc (cosmic)
Definition KalFitAlg.h:209

KalFitAlg::myEventNo
int myEventNo
Definition KalFitAlg.h:245

KalFitAlg::execute
StatusCode execute()
event function
Definition KalFitAlg.cxx:918

KalFitAlg::gain2_
double gain2_
Definition KalFitAlg.h:212

KalFitAlg::inner_steps_
int inner_steps_
Definition KalFitAlg.h:216

KalFitAlg::lead_
int lead_
leading mass assumption
Definition KalFitAlg.h:197

KalFitAlg::resolution_
int resolution_
Definition KalFitAlg.h:240

KalFitAlg::dchi2cuts_
double dchi2cuts_
Definition KalFitAlg.h:233

KalFitAlg::smoother_calib
void smoother_calib(KalFitTrack &trk, int way)
Definition KalFitAlg.cxx:3053

KalFitAlg::Tds_back_no
int Tds_back_no
Definition KalFitAlg.h:193

KalFitAlg::init_matrix
void init_matrix(MdcRec_trk &trk, HepSymMatrix &Ea)
Definition KalFitAlg.cxx:6548

KalFitAlg::setCylinder_Mdc
void setCylinder_Mdc(void)
Initialize the cylinders (walls and cathodes) for Mdc.

KalFitAlg::numf_
int numf_
Definition KalFitAlg.h:215

KalFitAlg::pk_cut_
double pk_cut_
Definition KalFitAlg.h:201

KalFitAlg::pT_
double pT_
value of the pT cut for backward filter
Definition KalFitAlg.h:195

KalFitAlg::steplev_
int steplev_
Definition KalFitAlg.h:213

KalFitAlg::matrixg_
double matrixg_
Definition KalFitAlg.h:210

KalFitAlg::initialize
StatusCode initialize()
initialize
Definition KalFitAlg.cxx:204

KalFitAlg::finalize
StatusCode finalize()
Definition KalFitAlg.cxx:320

KalFitAlg::kalman_fitting_csmalign
void kalman_fitting_csmalign(void)
Definition KalFitAlg.cxx:5563

KalFitAlg::gain4_
double gain4_
Definition KalFitAlg.h:212

KalFitAlg::ntuple_
int ntuple_
Fill ntuples of KalFit.
Definition KalFitAlg.h:229

KalFitAlg::ppi_cut_
double ppi_cut_
Definition KalFitAlg.h:201

KalFitAlg::hist_def
void hist_def(void)
hist definition
Definition KalFitAlg.cxx:378

KalFitAlg::fillTds_back
void fillTds_back(KalFitTrack &track, RecMdcKalTrack *trk, MdcRec_trk &TrasanTRK, int l_mass)
with results got at the outer Mdc hit
Definition KalFitAlg.cxx:2031

KalFitAlg::setMaterial_Mdc
void setMaterial_Mdc(void)
Initialize the material for Mdc.

KalFitAlg::dchi2cut_inner_
double dchi2cut_inner_
Definition KalFitAlg.h:235

KalFitAlg::activeonly_
int activeonly_
Definition KalFitAlg.h:207

KalFitAlg::extToAnyPoint
void extToAnyPoint(KalFitTrack &trk, const HepPoint3D &point)

KalFitAlg::clearTables
void clearTables()
Definition KalFitAlg.cxx:6755

KalFitHelixSeg
Definition KalFitHelixSeg.h:16

KalFitHitMdc
Description of a Hit in Mdc.
Definition KalFitHitMdc.h:17

KalFitLayer_Mdc
Definition KalFitLayer_Mdc.h:13

KalFitSuper_Mdc
Definition KalFitSuper_Mdc.h:8

KalFitTrack
Description of a track class (<- Helix.cc)
Definition KalFitTrack.h:36

KalFitWire
Description of a Wire class.
Definition KalFitWire.h:46

MdcCalibFunSvc
Definition MdcCalibFunSvc.h:22

MdcRec_trk_add
Definition MdcTables.h:727

MdcRec_trk
Definition MdcTables.h:355

RecMdcKalTrack
Definition RecMdcKalTrack.h:30

std
Definition RootEventData_rootcint.cxx:16