myTrTrackPlus.h

Go to the documentation of this file.

// Authors: M.Duranti - INFN di Perugia
#ifndef myTrTrackPlus_h
#define myTrTrackPlus_h

#include <map>
#include "myTrTrackBase.h"
#include "myPoint.h"

#define DOHOWEVERALLFITS

#ifdef _WITHGBATCH_
//#define PATHONB
#endif //#ifdef _WITHGBATCH_

//---------------------------------------------

struct fitresult {
   Double32_t Rigidity;
   Double32_t ErrRinv;
   Double32_t Chisq;
   Double32_t NormChisqX;
   Double32_t NormChisqY;
};

class myTrTrackPlus: public myTrTrackBase {

 public:
  myTrTrackPlus();
  ~myTrTrackPlus();

  Int_t  QStatusL[9];
  myPoint GlobalCoordinateAL[9];
  myPoint GlobalCoordinateNL[9];
  myPoint NoiseMinDist[5];

 protected:
  std::map<int, fitresult> _fitresultmap;

 public:
  static int ComputeIndex(int algo=0, int pattern=0, int refit=0);
  virtual int iTrTrackPar(int algo=0, int pattern=0, int refit=0);
  virtual Double32_t GetRigidity(int ittp=0);
  virtual Double32_t GetErrRinv(int ittp=0);
  virtual Double32_t GetChisq(int ittp=0);
  virtual Double32_t GetNormChisqX(int ittp=0);
  virtual Double32_t GetNormChisqY(int ittp=0);
  virtual Double32_t GetRigidityCorrectedByMagnField(int ittp=0);

  Int_t BitPatternJ;
  Int_t kDef;

  Short_t span();

 protected:
  Int_t kDefNoMS;
  Int_t kDefSameWeight;
  Int_t k89;
  Int_t k8;
  Int_t k9;
  Int_t kinn;
  Int_t k89NoMS;
  Int_t k8NoMS;
  Int_t k9NoMS;
  Int_t kinnNoMS;
  Int_t k89SameWeight;
  Int_t k8SameWeight;
  Int_t k9SameWeight;
  Int_t kinnSameWeight;
  Int_t kUp;
  Int_t kLow;
  Int_t kExt;
  Int_t kAl;
  Int_t kCIE;
  Int_t kPGCIE;
  Int_t Refit;

 public:
  /* Double32_t B; */
  /* Double32_t SenThetaX; */
  /* Double32_t SenThetaY; */
  /* Double32_t L; */
  /* Double32_t pMDR; */
  /* Double32_t pMDR_known; */

  Double32_t LH;

  myPoint PointTrackCogEcal;
  myDir DirTrackCogEcal;
  myPoint PointTrackEntryEcal;
  myDir DirTrackEntryEcal;
  myPoint PointTrackExitEcal;
  myDir DirTrackExitEcal;
  
  myPoint PointTrackTopTrd;
  myDir DirTrackTopTrd;
  myPoint PointTrackCenterTrd;
  myDir DirTrackCenterTrd;
  myPoint PointTrackBottomTrd;
  myDir DirTrackBottomTrd;

  Double32_t TruncMeanX;
  Double32_t TruncMeanY;
  Double32_t TruncMeanInnX;
  Double32_t TruncMeanInnY;
  
  Int_t Q_NPoints;
  Double32_t Q_RMS;
  Double32_t InnerQ_RMS;

  Double32_t factor;
  Int_t btempcor_status;

  //------------------------

  Double32_t distHMx_inn;
  Double32_t distHMy_inn;
  Double32_t distHMx_L1;
  Double32_t distHMy_L1;
  Double32_t distHMx_L9;
  Double32_t distHMy_L9;
  Double32_t distHMx_L8;
  Double32_t distHMy_L8;

  Double32_t distNHx_inn;
  Double32_t distNHy_inn;
  Double32_t distNHx_L1;
  Double32_t distNHy_L1;
  Double32_t distNHx_L9;
  Double32_t distNHy_L9;
  Double32_t distNHx_L8;
  Double32_t distNHy_L8;
  
 public:
  virtual void Clear(Option_t* option="");
 protected:
  virtual void init();

 public:
  ClassDef(myTrTrackPlus,9);
};

#ifdef _WITHGBATCH_

class myTrTrackPlusFiller: public myTrTrackPlus {  
 public:
  myTrTrackPlusFiller();
  ~myTrTrackPlusFiller();
  
  TrTrackR* track;
  
  virtual void Fill(TrTrackR* _track, short int ipart, EcalShowerR* shower, Double_t beta=999, bool kMC=false, bool heavy=true);
  virtual void FillCharge(TrTrackR* _track, Double_t beta);
  virtual void FillMoreHeavy(TrTrackR* _track);
 private:
  virtual void Fits(Double_t beta=999, bool kMC=false, bool heavy=true);
  
  virtual int AddResults(int algo=0, int pattern=0, int refit=0, float mass=TrTrackR::DefaultMass, float chrg=TrTrackR::DefaultCharge, float beta=999, float fixrig=0.);
  inline virtual void AddResults(fitresult& fr, int algo=0, int pattern=0, int refit=0);
  
 public:
  virtual void Clear(Option_t* option="");
 private:
  virtual void init();
  
 public:
  ClassDef(myTrTrackPlusFiller,1);
};

Double_t getRigidity(TrTrackR* track, int par, bool rigcorr=false);
Double_t getChisq(TrTrackR* track, int par);
Double_t getNormChisqX(TrTrackR* track, int par);
Double_t getNormChisqY(TrTrackR* track, int par);
Double_t getLogExtChisq(TrTrackR *trk, int k89, int k8, int k9);

Double_t HitCloseTrack(int layer, int side, const char* type, AMSEventR* fEvent, TrTrackR* track, bool kMC=false, int TrackId=0);
void BuildHitCloseTrack(AMSEventR* fEvent, Double_t** distHM, Double_t** distNH, TrTrackR* track, int TrackId, std::vector<TrClusterR*> TrClusUsed, bool kMC);
void BuildVectorUsedClusters(AMSEventR* fEvent, std::vector<TrClusterR*>& TrClusUsed);
class myPoint;//"forward declaration"
void GetMinDistanceNoiseHits(TrTrackR *track, myPoint mindistance[]);

void BuildResiduals(bool kMC, Double_t** res, TrTrackR* track, int TrackId);
Double_t LocalResidual(int layer, int side, bool kMC, AMSEventR* fEvent, TrTrackR* track, int TrackId=0);

void PathOnB(TrTrackR* tr, int kDef, MagField* magfield, Double_t& B, Double_t& SenThetaX, Double_t& SenThetaY, Double_t& L, Double_t& pMDR, Double_t& pMDR_known, Double_t step=0.1, bool modulus=false);

#endif //#ifdef _WITHGBATCH_

#endif