Reconstruction/MdcPatRec/MdcTrkRecon/include/MdcTrkRecon/MdcHistItem.h

Go to the documentation of this file.

//
// Description: class for ntuple items
//
 
#ifndef MDCHISTITEM_H
#define MDCHISTITEM_H
 
#include "AIDA/IHistogram1D.h"
#include "AIDA/IHistogram2D.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
 
// histograms
// extern AIDA::IHistogram1D*  g_actHit;
// extern AIDA::IHistogram1D*  g_hitEff;
extern AIDA::IHistogram2D* g_residVsLayer;
extern AIDA::IHistogram2D* g_residVsDoca;
extern AIDA::IHistogram1D* h_mapHit;
extern AIDA::IHistogram1D* h_sfHit;
 
//-----------cut--------------
// extern AIDA::IHistogram1D*  g_segChi2;
// extern AIDA::IHistogram1D*  g_segChi2SlayPat3[43];
// extern AIDA::IHistogram1D*  g_segChi2SlayPat4[43];
extern AIDA::IHistogram1D* g_maxSegChisqAx;
extern AIDA::IHistogram1D* g_maxSegChisqSt;
extern AIDA::IHistogram1D* g_nSigAdd;
extern AIDA::IHistogram1D* g_cirTkChi2;
extern AIDA::IHistogram1D* g_z0Cut;
extern AIDA::IHistogram1D* g_ctCut;
extern AIDA::IHistogram1D* g_delCt;
extern AIDA::IHistogram1D* g_delZ0;
extern AIDA::IHistogram1D* g_delPhi;
extern AIDA::IHistogram1D* g_delSlope;
extern AIDA::IHistogram1D* g_combChi2;
extern AIDA::IHistogram1D* g_phiDiff;
extern AIDA::IHistogram1D* g_slopeDiff;
extern AIDA::IHistogram1D* g_fitNAct;
extern AIDA::IHistogram1D* g_3dTkChi2;
extern AIDA::IHistogram1D* g_pickHitWire;
 
// tuple of MdcTrkRecon
extern NTuple::Tuple* m_tupleMc;
extern NTuple::Tuple* m_tupleMcHit;
extern NTuple::Tuple* m_tuple1;
extern NTuple::Tuple* m_tupleSeg;
extern NTuple::Tuple* m_tupleEvt;
// tuple item of mc truth
extern NTuple::Item<double>  m_tMcEvtNo;
extern NTuple::Item<long>    m_tMcNTk;
extern NTuple::Item<long>    m_tMcTkId;
extern NTuple::Item<long>    m_tMcPid;
extern NTuple::Item<double>  m_tMcPx;
extern NTuple::Item<double>  m_tMcPy;
extern NTuple::Item<double>  m_tMcPz;
extern NTuple::Item<double>  m_tMcD0;
extern NTuple::Item<double>  m_tMcZ0;
extern NTuple::Item<double>  m_tMcTheta;
extern NTuple::Item<double>  m_tMcFiTerm;
extern NTuple::Item<long>    m_tMcNHit;
extern NTuple::Array<double> m_tMcLayer;
extern NTuple::Array<double> m_tMcWire;
extern NTuple::Array<double> m_tMcDrift;
extern NTuple::Array<double> m_tMcX;
extern NTuple::Array<double> m_tMcY;
extern NTuple::Array<double> m_tMcZ;
extern NTuple::Array<double> m_tMcLR;
 
// tuple item of reconstruction results
extern NTuple::Item<double> m_t0;
extern NTuple::Item<double> m_t0Stat;
extern NTuple::Item<double> m_t0Truth;
extern NTuple::Item<double> m_nTdsTk;
extern NTuple::Item<double> m_nMcTk;
extern NTuple::Item<double> m_p;
extern NTuple::Item<double> m_pt;
extern NTuple::Item<double> m_nSlay;
extern NTuple::Item<double> m_pz;
extern NTuple::Item<double> m_d0;
extern NTuple::Item<double> m_phi0;
extern NTuple::Item<double> m_cpa;
extern NTuple::Item<double> m_z0;
extern NTuple::Item<double> m_tanl;
extern NTuple::Item<double> m_q;
extern NTuple::Item<double> m_pocax;
extern NTuple::Item<double> m_pocay;
extern NTuple::Item<double> m_pocaz;
extern NTuple::Item<double> m_evtNo;
extern NTuple::Item<double> m_nSt;
extern NTuple::Item<double> m_nDof;
extern NTuple::Item<double> m_chi2;
extern NTuple::Item<double> m_tkId;
extern NTuple::Item<long>   m_nHit;
extern NTuple::Item<double> m_nAct;
extern NTuple::Item<double> m_layerCount;
extern NTuple::Item<long>   m_nDigi;
 
extern NTuple::Array<double> m_resid;
extern NTuple::Array<double> m_sigma;
extern NTuple::Array<double> m_driftD;
extern NTuple::Array<double> m_driftT;
extern NTuple::Array<double> m_doca;
extern NTuple::Array<double> m_entra;
extern NTuple::Array<double> m_fltLen;
extern NTuple::Array<double> m_tof;
extern NTuple::Array<double> m_ambig;
extern NTuple::Array<double> m_act;
extern NTuple::Array<double> m_tdc;
extern NTuple::Array<double> m_adc;
extern NTuple::Array<double> m_layer;
extern NTuple::Array<double> m_wire;
extern NTuple::Array<double> m_x;
extern NTuple::Array<double> m_y;
extern NTuple::Array<double> m_z;
extern NTuple::Array<double> m_dx;
extern NTuple::Array<double> m_dy;
extern NTuple::Array<double> m_dz;
extern NTuple::Array<double> m_dDriftD;
extern NTuple::Array<double> m_dlr;
extern NTuple::Array<double> m_cellWidth; // for MdcTrackList::pickHits() goodDriftCut
 
// tuple item of event
extern NTuple::Item<long>    m_t4EvtNo;
extern NTuple::Item<double>  m_t4nMcTk;
extern NTuple::Item<double>  m_t4nRecTk;
extern NTuple::Item<double>  m_t4t0;
extern NTuple::Item<long>    m_t4t0Stat;
extern NTuple::Item<double>  m_t4t0Truth;
extern NTuple::Item<long>    m_t4nDigi;
extern NTuple::Item<long>    m_t4nDigiUnmatch;
extern NTuple::Array<long>   m_t4Layer;
extern NTuple::Array<long>   m_t4Wire;
extern NTuple::Array<double> m_t4Time;
extern NTuple::Array<double> m_t4Charge;
extern NTuple::Array<double> m_t4recHit;
extern NTuple::Array<double> m_t4PhiMid;
extern NTuple::Array<double> m_t4Hot;
 
// tuple item of segment
extern NTuple::Item<long>  m_tsNSeg;
extern NTuple::Item<long>  m_tsEvtNo;
extern NTuple::Item<long>  m_tsNDigi;
extern NTuple::Array<long> m_tsLayer;
extern NTuple::Array<long> m_tsWire;
extern NTuple::Array<long> m_tsInSeg;
extern NTuple::Array<long> m_tsMcTkId;
 
// tuple item of combine ax segs
extern NTuple::Tuple*       g_tupleCombAx;
extern NTuple::Item<double> g_combAxdPhi0;
extern NTuple::Item<double> g_combAxdCurv;
extern NTuple::Item<double> g_combAxSigPhi0;
extern NTuple::Item<double> g_combAxSigCurv;
extern NTuple::Item<double> g_combAxSlSeed;
extern NTuple::Item<double> g_combAxSlTest;
extern NTuple::Item<double> g_combAxQualitySeed;
extern NTuple::Item<double> g_combAxQualityTest;
extern NTuple::Item<double> g_combAxPatSeed;
extern NTuple::Item<double> g_combAxPatTest;
extern NTuple::Item<double> g_combAxNhitSeed;
extern NTuple::Item<double> g_combAxNhitTest;
extern NTuple::Item<double> g_combAxMc;
extern NTuple::Item<double> g_combAxMcPt;
extern NTuple::Item<double> g_combAxMcTheta;
extern NTuple::Item<double> g_combAxMcPhi;
extern NTuple::Item<double> g_combAxMcAmbigSeed;
extern NTuple::Item<double> g_combAxMcAmbigTest;
 
extern NTuple::Tuple*       m_tuplet;
extern NTuple::Item<long>   m_tl;
extern NTuple::Item<long>   m_tw;
extern NTuple::Item<double> m_tphi;
extern NTuple::Item<double> m_tdphi;
extern NTuple::Item<double> m_tdncell;
 
#ifdef MDCPATREC_TIMETEST
extern NTuple::Tuple*       m_tupleTime;
extern NTuple::Item<double> m_eventTime;
extern NTuple::Item<long>   m_recTkNum;
extern NTuple::Item<long>   m_mcTkNum;
extern NTuple::Item<long>   m_t5EvtNo;
extern NTuple::Item<long>   m_t5nHit;
extern NTuple::Item<long>   m_t5nDigi;
extern NTuple::Item<long>   m_t5nSt;
extern NTuple::Item<double> m_t5fitTime;
extern NTuple::Item<double> m_t5p;
extern NTuple::Item<double> m_t5pt;
#endif
 
#ifdef MDCPATREC_RESLAYER
int                         m_resLayer;
extern NTuple::Tuple*       g_tupleres;
extern NTuple::Item<double> g_resLayer;
extern NTuple::Item<double> g_t6Layer;
#endif
 
extern NTuple::Tuple*       g_tupleFindSeg;
extern NTuple::Item<double> g_findSegChi2;
extern NTuple::Item<double> g_findSegIntercept;
extern NTuple::Item<double> g_findSegSlope;
extern NTuple::Item<double> g_findSegChi2Refit;
extern NTuple::Item<double> g_findSegChi2Add;
extern NTuple::Item<int>    g_findSegPat;
extern NTuple::Item<int>    g_findSegNhit;
extern NTuple::Item<int>    g_findSegPat34;
extern NTuple::Item<int>    g_findSegSl;
extern NTuple::Item<double> g_findSegMc;
extern NTuple::Item<double> g_findSegAmbig;
 
extern NTuple::Tuple*        m_tuplePick;
extern NTuple::Item<long>    m_pickLayer;
extern NTuple::Item<long>    m_pickNCell;
extern NTuple::Item<long>    m_pickNCellWithDigi;
extern NTuple::Array<long>   m_pickWire;
extern NTuple::Array<double> m_pickPredDrift;
extern NTuple::Array<double> m_pickDrift;
extern NTuple::Array<double> m_pickDriftTruth;
extern NTuple::Array<int>    m_pickPhizOk;
extern NTuple::Array<double> m_pickPhiMaxDiff;
extern NTuple::Array<double> m_pickPhiMinDiff;
extern NTuple::Array<double> m_pickZ;
extern NTuple::Array<double> m_pickResid;
extern NTuple::Array<double> m_pickSigma;
extern NTuple::Array<double> m_pickPhiLowCut;
extern NTuple::Array<double> m_pickPhiHighCut;
extern NTuple::Array<double> m_pickDriftCut;
extern NTuple::Array<long>   m_pickMcTk;
extern NTuple::Array<long>   m_pickIs2D;
extern NTuple::Array<double> m_pickPt;
extern NTuple::Array<double> m_pickCurv;
 
extern NTuple::Tuple*       m_tupleWireEff;
extern NTuple::Item<long>   m_we_tkId;
extern NTuple::Item<long>   m_we_pdg;
extern NTuple::Item<long>   m_we_primary;
extern NTuple::Item<long>   m_we_layer;
extern NTuple::Item<long>   m_we_wire;
extern NTuple::Item<long>   m_we_gwire;
extern NTuple::Item<double> m_we_pt;
extern NTuple::Item<double> m_we_theta;
extern NTuple::Item<double> m_we_phi;
// extern NTuple::Item<double>             m_we_tdc;
extern NTuple::Item<long> m_we_poisoned;
extern NTuple::Item<long> m_we_seg;
extern NTuple::Item<long> m_we_track;
 
extern NTuple::Tuple*     m_tupleWireEffPoi;
extern NTuple::Item<long> m_we_poi_tkId;
extern NTuple::Item<long> m_we_poi_pdg;
extern NTuple::Item<long> m_we_poi_primary;
extern NTuple::Item<long> m_we_poi_layer;
extern NTuple::Item<long> m_we_poi_wire;
 
extern int    haveDigiTk[43][288];
extern double haveDigiPt[43][288];
extern double haveDigiTheta[43][288];
extern double haveDigiPhi[43][288];
extern double haveDigiDrift[43][288];
extern int    haveDigiAmbig[43][288];
extern int    recHitMap[43][288];
#endif