[u/mrichter/AliRoot.git] / STEER / AliAlignmentTracks.h

#ifndef ALIALIGNMENTTRACKS_H
#define ALIALIGNMENTTRACKS_H

//*************************************************************************
// AliAlignmentTracks: main steering class which deals with the alignment *
// procedures based on reconstructed tracks.                              *
// More comments will come with the development of the interfaces and     *
// functionalities of the class.                                          *
//*************************************************************************

#include <TObject.h>

#include "AliAlignObj.h"

class TChain;
class AliTrackPointArray;
class AliAlignObj;
class AliTrackFitter;
class AliTrackResiduals;

class AliAlignmentTracks : public TObject {

 public:

  AliAlignmentTracks();
  AliAlignmentTracks(TChain *esdchain);
  AliAlignmentTracks(const char *esdfilename, const char *esdtreename = "esdTree");
  virtual ~AliAlignmentTracks();

  void AddESD(TChain *esdchain);
  void AddESD(const char *esdfilename, const char *esdtreename = "esdTree");

  void SetPointsFilename(const char *pointsfilename = "AliTrackPoints.root") { fPointsFilename = pointsfilename; }

  void ProcessESD(TSelector *selector);
  void ProcessESD(Bool_t onlyITS=kFALSE,Int_t minITSpts=0,
		  Bool_t cuts=kTRUE,
		  Float_t minMom=0.3,Float_t maxMom=1.e9,
		  Float_t minAbsSinPhi=0.,Float_t maxAbsSinPhi=1.,
		  Float_t minSinTheta=0.,Float_t maxSinTheta=1.);
  void ProcessESDCosmics(Bool_t onlyITS=kFALSE,Int_t minITSpts=0,
			 Float_t maxMatchingAngle=0.17, // 10 deg
			 Bool_t cuts=kTRUE,
			 Float_t minMom=0.3,Float_t maxMom=1.e9,
			 Float_t minAbsSinPhi=0.,Float_t maxAbsSinPhi=1.,
			 Float_t minSinTheta=0.,Float_t maxSinTheta=1.);

  void BuildIndex();

  Bool_t ReadAlignObjs(const char *alignObjFileName = "AlignObjs.root", const char* arrayName = "Alignment");

  void SetTrackFitter(AliTrackFitter *fitter) { fTrackFitter = fitter; }
  void SetMinimizer(AliTrackResiduals *minimizer) { fMinimizer = minimizer; }

  void AlignDetector(AliGeomManager::ELayerID firstLayer,
		     AliGeomManager::ELayerID lastLayer,
		     AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
		     AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,Int_t iterations = 1);
  void AlignLayer(AliGeomManager::ELayerID layer,
		  AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
		  AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,
		  Int_t iterations = 1);
  void AlignVolume(UShort_t volId, UShort_t volIdFit,
		   Int_t iterations);
  void AlignVolumes(const TArrayI *volids, const TArrayI *volidsfit = 0x0,
		   AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
		   AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,
		   Int_t iterations = 1);

  AliAlignObj* GetAlignObj(UShort_t volid) const {
    Int_t iModule;
    AliGeomManager::ELayerID iLayer = AliGeomManager::VolUIDToLayer(volid,iModule);
    return fAlignObjs[iLayer-AliGeomManager::kFirstLayer][iModule];
  }
  void    SetUpdate(Bool_t update){fDoUpdate = update;}
  Bool_t  GetUpdate() const { return fDoUpdate;}
 protected:

  void InitIndex();
  void ResetIndex();
  void DeleteIndex();

  void InitAlignObjs();
  void ResetAlignObjs();
  void DeleteAlignObjs();

  Int_t LoadPoints(const TArrayI *volids, AliTrackPointArray** &points);
  void  UnloadPoints(Int_t n, AliTrackPointArray **points);

  AliTrackFitter *CreateFitter();
  AliTrackResiduals *CreateMinimizer();

  Bool_t Misalign(const char *misalignObjFileName, const char* arrayName);

  TChain           *fESDChain;       //! Chain with ESDs
  TString           fPointsFilename; //  Name of the file containing the track point arrays
  TFile            *fPointsFile;     //  File containing the track point arrays
  TTree            *fPointsTree;     //  Tree with the track point arrays
  Int_t           **fLastIndex;      //! Last filled index in volume arrays
  TArrayI        ***fArrayIndex;     //! Volume arrays which contains the tree index
  Bool_t            fIsIndexBuilt;   //  Is points tree index built
  AliAlignObj    ***fAlignObjs;      //  Array with alignment objects
  AliAlignObj    ***fMisalignObjs;   //  Array with alignment objects used to introduce misalignment of the space-points
  AliTrackFitter   *fTrackFitter;    //  Pointer to the track fitter
  AliTrackResiduals*fMinimizer;      //  Pointer to track residuals minimizer
  Bool_t            fDoUpdate;       //  Indicator - update Alignment object after minimization

 private:
  AliAlignmentTracks(const AliAlignmentTracks & alignment);
  AliAlignmentTracks& operator= (const AliAlignmentTracks& alignment);

  ClassDef(AliAlignmentTracks,2)

};

#endif
Commit	Line	Data
98937d93	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");
98937d93	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);
babbb915	36	void ProcessESD(Bool_t onlyITS=kFALSE,Int_t minITSpts=0,
babbb915	37	Bool_t cuts=kTRUE,
9a1304f0	38	Float_t minMom=0.3,Float_t maxMom=1.e9,
babbb915	39	Float_t minAbsSinPhi=0.,Float_t maxAbsSinPhi=1.,
	40	Float_t minSinTheta=0.,Float_t maxSinTheta=1.);
	41	void ProcessESDCosmics(Bool_t onlyITS=kFALSE,Int_t minITSpts=0,
9a1304f0	42	Float_t maxMatchingAngle=0.17, // 10 deg
	43	Bool_t cuts=kTRUE,
	44	Float_t minMom=0.3,Float_t maxMom=1.e9,
	45	Float_t minAbsSinPhi=0.,Float_t maxAbsSinPhi=1.,
	46	Float_t minSinTheta=0.,Float_t maxSinTheta=1.);
98937d93	47
98937d93	48	void BuildIndex();
98937d93	49
46ae650f	50	Bool_t ReadAlignObjs(const char alignObjFileName = "AlignObjs.root", const char arrayName = "Alignment");
98937d93	51
	52	void SetTrackFitter(AliTrackFitter *fitter) { fTrackFitter = fitter; }
	53	void SetMinimizer(AliTrackResiduals *minimizer) { fMinimizer = minimizer; }
	54
25be1e5c	55	void AlignDetector(AliGeomManager::ELayerID firstLayer,
	56	AliGeomManager::ELayerID lastLayer,
	57	AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
	58	AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,Int_t iterations = 1);
	59	void AlignLayer(AliGeomManager::ELayerID layer,
	60	AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
	61	AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,
98937d93	62	Int_t iterations = 1);
cc345ce3	63	void AlignVolume(UShort_t volId, UShort_t volIdFit,
	64	Int_t iterations);
	65	void AlignVolumes(const TArrayI volids, const TArrayI volidsfit = 0x0,
25be1e5c	66	AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
25be1e5c	67	AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,
46ae650f	68	Int_t iterations = 1);
98937d93	69
cc345ce3	70	AliAlignObj* GetAlignObj(UShort_t volid) const {
cc345ce3	71	Int_t iModule;
25be1e5c	72	AliGeomManager::ELayerID iLayer = AliGeomManager::VolUIDToLayer(volid,iModule);
25be1e5c	73	return fAlignObjs[iLayer-AliGeomManager::kFirstLayer][iModule];
cc345ce3	74	}
91264393	75	void SetUpdate(Bool_t update){fDoUpdate = update;}
91264393	76	Bool_t GetUpdate() const { return fDoUpdate;}
98937d93	77	protected:
	78
	79	void InitIndex();
	80	void ResetIndex();
	81	void DeleteIndex();
	82
	83	void InitAlignObjs();
	84	void ResetAlignObjs();
	85	void DeleteAlignObjs();
	86
cc345ce3	87	Int_t LoadPoints(const TArrayI volids, AliTrackPointArray* &points);
98937d93	88	void UnloadPoints(Int_t n, AliTrackPointArray **points);
	89
	90	AliTrackFitter *CreateFitter();
	91	AliTrackResiduals *CreateMinimizer();
	92
46ae650f	93	Bool_t Misalign(const char misalignObjFileName, const char arrayName);
46ae650f	94
98937d93	95	TChain *fESDChain; //! Chain with ESDs
	96	TString fPointsFilename; // Name of the file containing the track point arrays
	97	TFile *fPointsFile; // File containing the track point arrays
	98	TTree *fPointsTree; // Tree with the track point arrays
	99	Int_t **fLastIndex; //! Last filled index in volume arrays
	100	TArrayI ***fArrayIndex; //! Volume arrays which contains the tree index
	101	Bool_t fIsIndexBuilt; // Is points tree index built
	102	AliAlignObj ***fAlignObjs; // Array with alignment objects
46ae650f	103	AliAlignObj ***fMisalignObjs; // Array with alignment objects used to introduce misalignment of the space-points
98937d93	104	AliTrackFitter *fTrackFitter; // Pointer to the track fitter
98937d93	105	AliTrackResiduals*fMinimizer; // Pointer to track residuals minimizer
91264393	106	Bool_t fDoUpdate; // Indicator - update Alignment object after minimization
75e3794b	107
	108	private:
	109	AliAlignmentTracks(const AliAlignmentTracks & alignment);
	110	AliAlignmentTracks& operator= (const AliAlignmentTracks& alignment);
	111
91264393	112	ClassDef(AliAlignmentTracks,2)
98937d93	113
	114	};
	115
	116	#endif