STEER/AliAlignmentTracks.h

   1 #ifndef ALIALIGNMENTTRACKS_H
   2 #define ALIALIGNMENTTRACKS_H
   3
   4 //*************************************************************************
   5 // AliAlignmentTracks: main steering class which deals with the alignment *
   6 // procedures based on reconstructed tracks.                              *
   7 // More comments will come with the development of the interfaces and     *
   8 // functionalities of the class.                                          *
   9 //*************************************************************************
  10
  11 #include <TObject.h>
  12
  13 #include "AliAlignObj.h"
  14
  15 class TChain;
  16 class AliTrackPointArray;
  17 class AliAlignObj;
  18 class AliTrackFitter;
  19 class AliTrackResiduals;
  20
  21 class AliAlignmentTracks : public TObject {
  22
  23  public:
  24
  25   AliAlignmentTracks();
  26   AliAlignmentTracks(TChain *esdchain);
  27   AliAlignmentTracks(const char *esdfilename, const char *esdtreename = "esdTree");
  28   virtual ~AliAlignmentTracks();
  29
  30   void AddESD(TChain *esdchain);
  31   void AddESD(const char *esdfilename, const char *esdtreename = "esdTree");
  32
  33   void SetPointsFilename(const char *pointsfilename = "AliTrackPoints.root") { fPointsFilename = pointsfilename; }
  34
  35   void ProcessESD(TSelector *selector);
  36   void ProcessESD(Bool_t onlyITS=kFALSE,Int_t minITSpts=0,
  37                   Bool_t cuts=kTRUE,
  38                   Float_t minAngleWrtITSModulePlanes=0.,
  39                   Float_t minMom=0.3,Float_t maxMom=1.e9,
  40                   Float_t minAbsSinPhi=0.,Float_t maxAbsSinPhi=1.,
  41                   Float_t minSinTheta=0.,Float_t maxSinTheta=1.);
  42   void ProcessESDCosmics(Bool_t onlyITS=kFALSE,Int_t minITSpts=0,
  43                          Float_t maxMatchingAngle=0.17, // 10 deg
  44                          Bool_t cuts=kTRUE,
  45                          Float_t minAngleWrtITSModulePlanes=0.,
  46                          Float_t minMom=0.3,Float_t maxMom=1.e9,
  47                          Float_t minAbsSinPhi=0.,Float_t maxAbsSinPhi=1.,
  48                          Float_t minSinTheta=0.,Float_t maxSinTheta=1.);
  49
  50   void BuildIndex();
  51
  52   Bool_t ReadAlignObjs(const char *alignObjFileName = "AlignObjs.root", const char* arrayName = "Alignment");
  53
  54   void SetTrackFitter(AliTrackFitter *fitter) { fTrackFitter = fitter; }
  55   void SetMinimizer(AliTrackResiduals *minimizer) { fMinimizer = minimizer; }
  56
  57   void AlignDetector(AliGeomManager::ELayerID firstLayer,
  58                      AliGeomManager::ELayerID lastLayer,
  59                      AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
  60                      AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,Int_t iterations = 1);
  61   void AlignLayer(AliGeomManager::ELayerID layer,
  62                   AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
  63                   AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,
  64                   Int_t iterations = 1);
  65   void AlignVolume(UShort_t volId, UShort_t volIdFit,
  66                    Int_t iterations);
  67   void AlignVolumes(const TArrayI *volids, const TArrayI *volidsfit = 0x0,
  68                    AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
  69                    AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,
  70                    Int_t iterations = 1);
  71
  72   AliAlignObj* GetAlignObj(UShort_t volid) const {
  73     Int_t iModule;
  74     AliGeomManager::ELayerID iLayer = AliGeomManager::VolUIDToLayer(volid,iModule);
  75     return fAlignObjs[iLayer-AliGeomManager::kFirstLayer][iModule];
  76   }
  77   void    SetUpdate(Bool_t update){fDoUpdate = update;}
  78   Bool_t  GetUpdate() const { return fDoUpdate;}
  79  protected:
  80
  81   void InitIndex();
  82   void ResetIndex();
  83   void DeleteIndex();
  84
  85   void InitAlignObjs();
  86   void ResetAlignObjs();
  87   void DeleteAlignObjs();
  88
  89   Int_t LoadPoints(const TArrayI *volids, AliTrackPointArray** &points);
  90   void  UnloadPoints(Int_t n, AliTrackPointArray **points);
  91
  92   AliTrackFitter *CreateFitter();
  93   AliTrackResiduals *CreateMinimizer();
  94
  95   Bool_t Misalign(const char *misalignObjFileName, const char* arrayName);
  96
  97   TChain           *fESDChain;       //! Chain with ESDs
  98   TString           fPointsFilename; //  Name of the file containing the track point arrays
  99   TFile            *fPointsFile;     //  File containing the track point arrays
 100   TTree            *fPointsTree;     //  Tree with the track point arrays
 101   Int_t           **fLastIndex;      //! Last filled index in volume arrays
 102   TArrayI        ***fArrayIndex;     //! Volume arrays which contains the tree index
 103   Bool_t            fIsIndexBuilt;   //  Is points tree index built
 104   AliAlignObj    ***fAlignObjs;      //  Array with alignment objects
 105   AliAlignObj    ***fMisalignObjs;   //  Array with alignment objects used to introduce misalignment of the space-points
 106   AliTrackFitter   *fTrackFitter;    //  Pointer to the track fitter
 107   AliTrackResiduals*fMinimizer;      //  Pointer to track residuals minimizer
 108   Bool_t            fDoUpdate;       //  Indicator - update Alignment object after minimization
 109
 110  private:
 111   AliAlignmentTracks(const AliAlignmentTracks & alignment);
 112   AliAlignmentTracks& operator= (const AliAlignmentTracks& alignment);
 113
 114   ClassDef(AliAlignmentTracks,2)
 115
 116 };
 117
 118 #endif