AlgFitTrack3.h

Go to the documentation of this file.

// $Id: AlgFitTrack3.h,v 1.4 2003/11/10 20:19:42 brebel Exp $
//
// AlgFitTrack3.h
//
// AlgFitTrack3 is a concrete FitTrack3 Algorithm class.
//
// Author:  R. Nichol 05/2003
//

#ifndef ALGFITTRACK3_H
#define ALGFITTRACK3_H
#include <map>

#include "Algorithm/AlgBase.h"
#include "RecoBase/AlgReco.h"
#include "RecoBase/AlgTrack.h"
#include "Validity/VldContext.h"

#include "TMatrixD.h"
#include "TObjArray.h"
#include "TFile.h"
#include "TTree.h"
#include "TH1F.h"


using namespace std;

typedef map<int,double> IntDoubleMap;
typedef map<int,int> IntIntMap;


class AlgFitTrack3 : public AlgBase, public AlgReco, public AlgTrack
{

public:
   AlgFitTrack3();
   virtual ~AlgFitTrack3();
   virtual void RunAlg(AlgConfig &ac, CandHandle &ch, CandContext &cx);
   virtual void Trace(const char *c) const;
   


private:
   void CleanUp();

  Int_t fParmMaxIterate;
  Double_t fParmMisalignmentError;


   IntDoubleMap fRangeThisPlane; //Need it to calculate momentum from range
   IntDoubleMap fDistanceFromStart; //Relates to D(i,k);
   IntDoubleMap fTheta0i; // mean square scattering angle.
   IntDoubleMap fThicknessi; // thickness of steel plane at i.
   IntDoubleMap fSwamU;
   IntDoubleMap fSwamV;
   IntDoubleMap fSwamdU;
   IntDoubleMap fSwamdV;
   IntDoubleMap fErrorOnMeasure;
   IntDoubleMap fMeasuredU;
   IntDoubleMap fMeasuredV;
   IntDoubleMap fActualZ;
   

  Int_t fFirstPlane;
  Int_t fLastPlane;
  Int_t fLowestPlane;
  Int_t fHighestPlane;
  Int_t fDirection;
  
   Int_t fIsNaturalP0;
   Double_t fMomFromRange;
  Int_t fCharge;
  Double_t fZ0;
  Double_t fU0;
  Double_t fV0;
  Double_t fdUdZ0;
  Double_t fdVdZ0;
  Double_t fP0;

  Int_t fLastCharge;
  Double_t fLastZ0;
  Double_t fLastU0;
  Double_t fLastV0;
  Double_t fLastdUdZ0;
  Double_t fLastdVdZ0;
  Double_t fLastP0;

   Double_t fDP0;
   Double_t fDU0;
   Double_t fDV0;
   Double_t fDdUdZ0;
   Double_t fDdVdZ0;
   
   


  Double_t fChiSqU;
  Double_t fChiSqV;

   Int_t fIterations;
   Int_t fSizeOfUWeightMatrix;
   Int_t fSizeOfVWeightMatrix;
   IntIntMap fReverseMapUIndex;
   IntIntMap fReverseMapVIndex;

  Double_t CalculateD(int i, int k); //Takes plane numbers 
  Double_t CalculateP(int k, int n); //Takes 1-N not plane numbers.

  void FillWeightMatrix();
  void FillMatricesAandC();
  Double_t GetWiju(int row, int col);
  Double_t GetWijv(int row, int col);
  void InvertMatrixAndGetParams();
  void SwimAndFillMaps();
  int InitializeTrkClstList(const CandTrackHandle *track0);
  

   Int_t fTrackNum;

   //TMatrixD *fWeightMatrix;
   TMatrixD *fUWeightMatrix;
   TMatrixD *fVWeightMatrix;
   TMatrixD fMatrixA;
   TMatrixD fMatrixC;
   
   TObjArray fTrkClstList;
   VldContext *fMyVC;

   TFile *fFile;
   TTree *fSillyTree;
   TH1F *fReasonsForNotFitting;
   TTree *fReasonsTree;
   Int_t fR_totalU;
   Int_t fR_totalV;
   Int_t fR_showerU;
   Int_t fR_showerV;
   Int_t fR_validU;
   Int_t fR_validV;


   Double_t fS_z;
   Double_t fS_measuredU;
   Double_t fS_measuredV;
   Double_t fS_swamU;
   Double_t fS_swamV;
   Int_t fS_isU;

ClassDef(AlgFitTrack3,2)                            // FitTrack3 Algorithm Class

};

#endif                                                     // ALGFITTRACK3_H

---