CgemMdcCombAlg.h

Go to the documentation of this file.

#ifndef CGEMMDCCOMBALG_H 
#define CGEMMDCCOMBALG_H 
#include "GaudiKernel/Algorithm.h"
#include "CgemRecEvent/RecCgemCluster.h"
#include "CgemRecEvent/RecCgemSegment.h"
#include "McTruth/McParticle.h"
 
#include "VertexFit/IVertexDbSvc.h"
#include "GaudiKernel/NTuple.h"
#include "GaudiKernel/INTupleSvc.h"
#include "KalFitAlg/helix/Helix.h"
#include "CLHEP/Geometry/Point3D.h"
#include "CLHEP/Matrix/Vector.h"
 
#include "TROOT.h"
#include "TObjArray.h"
#include "TH1F.h"
#include "TProfile.h"
#include "TFile.h"
#include "TFolder.h"
#include "TMinuit.h"
#include "TGraph.h"
 
#include <vector>
 
#include <string>
 
using namespace std;
 
class CgemMdcCombAlg:public Algorithm {
public:
     CgemMdcCombAlg(const std::string& name, ISvcLocator* pSvcLocator);
     StatusCode execute();
     StatusCode initialize();
     StatusCode finalize();
 
private:
      int m_debug, m_combineTracking, m_dropHot, m_keepBadTdc, m_keepUnmatch;
      int m_prodNT;
      int unSegNum;
      int unTrkNum;
      int EvInWindow;
      int EvOutWindow;
 
      int m_evt, m_run;
 
      StatusCode exe_v1();
      StatusCode exe_v2();
 
      void getMcPar();
 
 
      double GetSigChisq(double x, double y, double err);
      double GetTotChisq(double seg[5], double trk[5], double err[15]);
      double GetWindowChisq(double pt, double a0, double a1, double a2, double a3);
      double GetWindowChisq(double pt, vector<double> Pt, vector<double> w_p);
      double GetWindowChisq(double pt, int i_par);
 
      void printMdcTrk();
 
      vector<double> m_dr_par;
      vector<double> m_phi0_par;
      vector<double> m_kappa_par;
      vector<double> m_dz_par;
      vector<double> m_tanl_par;
      
      vector<double> m_Pt;
      vector<double> m_w_dr;
      vector<double> m_w_phi0;
      vector<double> m_w_kappa;
      vector<double> m_w_dz;
      vector<double> m_w_tanl;
 
      TGraph* mg_mean_dr;
      TGraph* mg_mean_phi0;
      TGraph* mg_mean_kappa;
      TGraph* mg_mean_dz;
      TGraph* mg_mean_tanL;
      TGraph* mg_sig_dr;
      TGraph* mg_sig_phi0;
      TGraph* mg_sig_kappa;
      TGraph* mg_sig_dz;
      TGraph* mg_sig_tanL;
      int m_nPt;
 
      double Chi2Match(double delta, double pt, int iPar);
 
      vector<double> mv_Pt;
      vector<double> mv_mean_dr;
      vector<double> mv_mean_phi0;
      vector<double> mv_mean_kappa;
      vector<double> mv_mean_dz;
      vector<double> mv_mean_tanL;
      vector<double> mv_sig_dr;
      vector<double> mv_sig_phi0;
      vector<double> mv_sig_kappa;
      vector<double> mv_sig_dz;
      vector<double> mv_sig_tanL;
      vector<double> mv_matchChi2Cut;
 
 
 
      NTuple::Tuple*  match_ntuple;   
      NTuple::Item<long> evt;
      NTuple::Item<long> run;
      NTuple::Item<double> costhetaMC;
      NTuple::Item<double> phiMC;
      NTuple::Item<int> pdgMC;
      NTuple::Item<double> truth_dr;
      NTuple::Item<double> truth_dz;
      NTuple::Item<double> truth_phi0;
      NTuple::Item<double> truth_kappa;
      NTuple::Item<double> truth_tanl;
      NTuple::Item<long> match_nsegment;
      NTuple::Item<long> track_size;
      NTuple::Array<int> CgemSegMatched;
      NTuple::Array<double> dr_avg;
      NTuple::Array<double> dz_avg;
      NTuple::Array<double> tanL_avg;
      NTuple::Item<long> segment_size;
      NTuple::Array<int> match_failed;
      NTuple::Array<double> match_dr_chisq;
      NTuple::Array<double> match_dz_chisq;
      NTuple::Array<double> match_phi0_chisq;
      NTuple::Array<double> match_kappa_chisq;
      NTuple::Array<double> match_tanl_chisq;
      NTuple::Array<double> match_IsTruth;
      NTuple::Array<double> track_dr;
      NTuple::Array<double> track_dz;
      NTuple::Array<double> track_phi0;
      NTuple::Array<double> track_kappa;
      NTuple::Array<double> track_tanl;
      NTuple::Array<double> track_dr_err;
      NTuple::Array<double> track_dz_err;
      NTuple::Array<double> track_phi0_err;
      NTuple::Array<double> track_kappa_err;
      NTuple::Array<double> track_tanl_err;
      NTuple::Array<double> seg_dr;
      NTuple::Array<double> seg_dz;
      NTuple::Array<double> seg_phi0;
      NTuple::Array<double> seg_kappa;
      NTuple::Array<double> seg_tanl;
      NTuple::Array<double> seg_dr_err;
      NTuple::Array<double> seg_dz_err;
      NTuple::Array<double> seg_phi0_err;
      NTuple::Array<double> seg_kappa_err;
      NTuple::Array<double> seg_tanl_err;
      NTuple::Array<double> track[5];
 
};
 
 
#endif /* CGEMMDCCOMBALG_H */