Reconstruction/CgemSegmentFitAlg/CgemSegmentFitAlg-00-00-04/CgemSegmentFitAlg/CgemSegmentFitAlg.h

Go to the documentation of this file.

#ifndef CGEMSEGMENTFITALG_H 
#define CGEMSEGMENTFITALG_H 
#include "GaudiKernel/Algorithm.h"
#include "CgemRecEvent/RecCgemCluster.h"
#include "CgemRecEvent/RecCgemSegment.h"
#include "McTruth/McParticle.h"
#include "CgemCalibFunSvc/CgemCalibFunSvc.h"
 
#include "VertexFit/IVertexDbSvc.h"
#include "GaudiKernel/NTuple.h"
#include "GaudiKernel/INTupleSvc.h"
//#include "TrackUtil/Helix.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 <vector>
 
#include <string>
 
using namespace std;
 
class CgemSegmentFitAlg:public Algorithm {
public:   
     CgemSegmentFitAlg(const std::string& name, ISvcLocator* pSvcLocator);
     StatusCode execute();
     StatusCode initialize();
     StatusCode finalize(); 
 
     StatusCode exe_v1();
     StatusCode exe_v2();
 
 
private:     
 
 
     bool CheckClusterID(int layerID, int clusterID, int vec_clusterID[3]);
     int GetChargeOfTrack();
     int GetChargeOfTrack(double *vec_phi);
     int estiChargeOfSeg();
     double IntersectCylinder(double r, KalmanFit::Helix helix);
     double GetMin(double *a);
     double GetMax(double *a);
     int mnTrkFit(Double_t trkpar[], Double_t trkparErr[], Double_t* emat, double& chisq);
 
     int m_debug;
     int m_checkIdx;
     int m_ver;
 
     int m_run;
     int m_evt;
 
     double wide_delta;
     double pitch;
     int m_failed;
     int check_segNum;
     int check_segOver;
     TMinuit* m_pMnTrk;
 
     map<int , RecCgemSegment*> m_seg_map;               
     map<int , RecCgemSegment*>::iterator m_seg_map_it;  
 
     double m_helix_mc[5];
     void getMcPar();
 
     ICgemCalibFunSvc* myCgemCalibSvc;
 
     NTuple::Tuple*  helix_ntuple;
     NTuple::Item<long> cgem_run;
     NTuple::Item<long> cgem_evt;
     NTuple::Item<long> cgem_nsegment;
     NTuple::Item<long> truth_ntrack;
 
     NTuple::Array<double> truth_dr;
     NTuple::Array<double> truth_dz;
     NTuple::Array<double> truth_phi0;
     NTuple::Array<double> truth_kappa;
     NTuple::Array<double> truth_tanl;
     NTuple::Array<double> truth_CosTheta;
     NTuple::Array<int> truth_charge;
 
     NTuple::Array<int> cgem_segmentid;
     NTuple::Array<double> cgem_dr_ini;
     NTuple::Array<double> cgem_dz_ini;
     NTuple::Array<double> cgem_phi0_ini;
     NTuple::Array<double> cgem_kappa_ini;
     NTuple::Array<double> cgem_tanl_ini;
     NTuple::Array<double> cgem_dr;
     NTuple::Array<double> cgem_dz;
     NTuple::Array<double> cgem_phi0;
     NTuple::Array<double> cgem_kappa;
     NTuple::Array<double> cgem_tanl;
     NTuple::Array<double> cgem_chisq;
     NTuple::Array<int> cgem_charge;
     NTuple::Array<int> cgem_fitFlag;
 
     NTuple::Item<int> n_layer;
     NTuple::Array<int> n_stripX;
     NTuple::Array<int> n_stripV;
     NTuple::Array<double> mc_phi;
     NTuple::Array<double> mc_v;
     NTuple::Array<double> rec_phi;
     NTuple::Array<double> rec_v;
 
 
};
 
 
#endif /* CGEMSEGMENTFITALG_H */