[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 minAngleWrtITSModulePlanes=0.,
		  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 minAngleWrtITSModulePlanes=0.,
			 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,
d6116b3f	38	Float_t minAngleWrtITSModulePlanes=0.,
9a1304f0	39	Float_t minMom=0.3,Float_t maxMom=1.e9,
babbb915	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,
9a1304f0	43	Float_t maxMatchingAngle=0.17, // 10 deg
9a1304f0	44	Bool_t cuts=kTRUE,
d6116b3f	45	Float_t minAngleWrtITSModulePlanes=0.,
9a1304f0	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.);
98937d93	49
98937d93	50	void BuildIndex();
98937d93	51
46ae650f	52	Bool_t ReadAlignObjs(const char alignObjFileName = "AlignObjs.root", const char arrayName = "Alignment");
98937d93	53
	54	void SetTrackFitter(AliTrackFitter *fitter) { fTrackFitter = fitter; }
	55	void SetMinimizer(AliTrackResiduals *minimizer) { fMinimizer = minimizer; }
	56
25be1e5c	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,
98937d93	64	Int_t iterations = 1);
cc345ce3	65	void AlignVolume(UShort_t volId, UShort_t volIdFit,
	66	Int_t iterations);
	67	void AlignVolumes(const TArrayI volids, const TArrayI volidsfit = 0x0,
25be1e5c	68	AliGeomManager::ELayerID layerRangeMin = AliGeomManager::kFirstLayer,
25be1e5c	69	AliGeomManager::ELayerID layerRangeMax = AliGeomManager::kLastLayer,
46ae650f	70	Int_t iterations = 1);
98937d93	71
cc345ce3	72	AliAlignObj* GetAlignObj(UShort_t volid) const {
cc345ce3	73	Int_t iModule;
25be1e5c	74	AliGeomManager::ELayerID iLayer = AliGeomManager::VolUIDToLayer(volid,iModule);
25be1e5c	75	return fAlignObjs[iLayer-AliGeomManager::kFirstLayer][iModule];
cc345ce3	76	}
91264393	77	void SetUpdate(Bool_t update){fDoUpdate = update;}
91264393	78	Bool_t GetUpdate() const { return fDoUpdate;}
98937d93	79	protected:
	80
	81	void InitIndex();
	82	void ResetIndex();
	83	void DeleteIndex();
	84
	85	void InitAlignObjs();
	86	void ResetAlignObjs();
	87	void DeleteAlignObjs();
	88
cc345ce3	89	Int_t LoadPoints(const TArrayI volids, AliTrackPointArray* &points);
98937d93	90	void UnloadPoints(Int_t n, AliTrackPointArray **points);
	91
	92	AliTrackFitter *CreateFitter();
	93	AliTrackResiduals *CreateMinimizer();
	94
46ae650f	95	Bool_t Misalign(const char misalignObjFileName, const char arrayName);
46ae650f	96
98937d93	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
46ae650f	105	AliAlignObj ***fMisalignObjs; // Array with alignment objects used to introduce misalignment of the space-points
98937d93	106	AliTrackFitter *fTrackFitter; // Pointer to the track fitter
98937d93	107	AliTrackResiduals*fMinimizer; // Pointer to track residuals minimizer
91264393	108	Bool_t fDoUpdate; // Indicator - update Alignment object after minimization
75e3794b	109
	110	private:
	111	AliAlignmentTracks(const AliAlignmentTracks & alignment);
	112	AliAlignmentTracks& operator= (const AliAlignmentTracks& alignment);
	113
91264393	114	ClassDef(AliAlignmentTracks,2)
98937d93	115
	116	};
	117
	118	#endif