EdbPlateTracking.h

Go to the documentation of this file.

#ifndef ROOT_EdbPlateTracking
#define ROOT_EdbPlateTracking
 
 
#include "EdbRun.h"
#include "EdbPattern.h"
#include "EdbLayer.h"
#include "EdbRunAccess.h"
#include "EdbBrick.h"
 
//class EdbPlateP;
 
//______________________________________________________________________________
class EdbPlateTracking : public TObject {
 
 public:
 
 
  EdbSegP    ePred;              
  EdbPattern eSide1;             
  EdbPattern eSide2;             
  EdbPlateP  ePlate;             
 
 
  EdbSegP eS, eS1, eS2;          
  EdbSegP eNext;                 
  Int_t   eStatus;               
 
 
  EdbScanCond eCondMT;   // conditions for microtracks
  EdbScanCond eCondBT;   // conditions for basetracks
 
  /*   Float_t eDeltaRview;           // (400) */
  /*   Float_t eDeltaTheta;           // (0.15) slope acceptance */
  /*   Float_t ePreliminaryPulsMinMT; // (6)   \_ preliminary cuts to microtracks candidates for */
  Float_t ePreliminaryChi2MaxMT; 
 
  EdbPattern eS1pre;             
  EdbPattern eS2pre;
 
  Float_t eDeltaR;               
 
  EdbPattern eSpre;              
 
  Float_t ePulsMinBT;            
  Float_t ePulsMinDegradMT;      
  Float_t eChi2MaxBT;            
 
  EdbPattern eScnd;              
 
  Float_t ePulsMinMT;            
  Float_t ePulsMinDegradBT;      
  Float_t eChi2MaxMT;            
 
  EdbPattern eS1cnd;             
  EdbPattern eS2cnd;
 
  Float_t eDegradPos;            
  Float_t eDegradSlope;          
 
  Int_t   eIdp[4];               
  Int_t   eIdf[4];
 
  bool    ePredictionScan;       
 
 public:
 
  EdbPlateTracking()
    { 
      Set0(); 
    };
 
  EdbPlateTracking(EdbPattern &S1,  EdbPattern &S2, EdbSegP &prediction, EdbPlateP &plate)
    { 
      
      for(int i=0;i<S1.N();i++)
    eSide1.AddSegment(*(S1.GetSegment(i)));
      for(int i=0;i<S2.N();i++)
    eSide2.AddSegment(*(S2.GetSegment(i)));
 
      eSide1.SetZ(S1.Z());
      eSide2.SetZ(S2.Z());
      
      ePred.Copy(prediction);
      ePlate.Copy(plate);
      Set0(); 
    };
 
  EdbPlateTracking(EdbPattern &S1,  EdbPattern &S2, EdbSegP &prediction)
    { 
      
      for(int i=0;i<S1.N();i++)
    eSide1.AddSegment(*(S1.GetSegment(i)));
      for(int i=0;i<S2.N();i++)
    eSide2.AddSegment(*(S2.GetSegment(i)));
      
      eSide1.SetZ(S1.Z());
      eSide2.SetZ(S2.Z());
 
      ePred.Copy(prediction);
      Set0(); 
    };
 
  virtual ~EdbPlateTracking();
 
  void Set0();
  void SetCondMT(EdbScanCond &cond) { eCondMT = cond; }
  void SetCondBT(EdbScanCond &cond) { eCondBT = cond; }
 
  void SetPred(const EdbSegP &pred);
 
  int UpdateFlag(int flag, int status);
  static int GetBTHoles(int flag) { return(flag/10000);     }
  static int GetMTHoles(int flag) { return((flag/100)%100); }
  int ExtrapolateCond(EdbScanCond &inputcond, int flag, EdbScanCond &outputcond );
 
  int FindCompliments( EdbSegP &s, EdbPattern &pat, TObjArray &found, float chi2max, TArrayF &chiarr );
 
  int FindCandidates( EdbSegP &spred, EdbPattern &fndbt, EdbPattern &fnds1, EdbPattern &fnds2 );
  int FindCandidateMT( EdbPattern &fnds1, EdbPattern &fnds2, EdbSegP &fnd );
  //  int FindCandidateBT(EdbPattern &fndbt, EdbSegP &fnd );
  int FindBestCandidate(EdbPattern &fndbt, EdbSegP &fnd,EdbPattern &cnd, float wmin, float wmindegrad, float chi2max);
  int FindPrediction( EdbSegP &spred, EdbSegP &fndbt, EdbSegP &fnds1, EdbSegP &fnds2, EdbSegP &snewpred );
 
  int FindBestCandidateDS(EdbPattern &fndbt, EdbSegP &fnd,EdbPattern &cnd, float wmin, float wmindegrad, float chi2max, EdbSegP &spred, float maxdmin);
  int FindPredictionDS( EdbSegP &spred, EdbSegP &fndbt, EdbSegP &fnds1, EdbSegP &fnds2, EdbSegP &snewpred, float maxdmin);
  int FindCandidateMTDS( EdbPattern &fnds1, EdbPattern &fnds2, EdbSegP &fnd, EdbSegP &spred, float maxdmin);
 
  //void TransformToPlateRS( EdbPlateP &plate);
  //void TransformFromPlateRS( EdbPlateP &plate);
  void TransformFromPlateRS();
 
  int FindTrack(EdbTrackP &pred,EdbTrackP &found, EdbPlateP &plate);
 
  void Print();
 
  //void Transform(const EdbAffine2D &aff);
  
  static void   CloseSBtree(TTree *tree);
  static TTree *InitSBtree(const char *file_name="sbt.root", const char *mode="RECREATE");
  bool   UpdateSBtree( TTree &tsbt, int idp[4], int idf[4]);
  bool   GetSBtreeEntry( int entry, TTree &tsbt);
    // int stat, 
    //         EdbSegP &ps, EdbSegP &fndbt, EdbSegP &fnds1, EdbSegP &fnds2, EdbSegP &nextpred, 
    //         EdbPattern *scnd=0, EdbPattern *s1cnd=0, EdbPattern *s2cnd=0 );
 
  ClassDef(EdbPlateTracking,1)  // track following in one emulsion plate
};
 
#endif /* ROOT_EdbPlateTracking */