#include <TObject.h>
#include <TRandom.h>
#include <TObjString.h>
-
-#include <AliAlignObj.h>
+#include <AliGeomManager.h>
+class AliSurveyObj;
class AliTRDalignment : public TObject {
void SetSmZero(); // reset to zero supermodule data
void SetChZero(); // reset to zero chamber data
- void SetZero() {SetSmZero(); SetChZero();} // reset to zero both
void SetSm(int sm, const double a[6]) {for (int i = 0; i < 6; i++) fSm[sm][i] = a[i];}
void SetCh(int ch, const double a[6]) {for (int i = 0; i < 6; i++) fCh[ch][i] = a[i];}
void SetSmRandom(double a[6]); // generate random gaussians with sigmas a
void SetChFull(); // set chamber data to initial aka full
void SetSmResidual(); // set supermodule data to final aka residual
void SetChResidual(); // set chamber data to final aka residual
+ void SetZero() {SetSmZero(); SetChZero();}
+ void SetIdeal() {SetZero();}
void SetFull() {SetSmFull(); SetChFull();}
void SetResidual() {SetSmResidual(); SetChResidual();}
void SetComment(char *s) {fComment.SetString(s);}
// reading-in from file
void ReadAscii(char *filename); // read from ascii file
+ void ReadCurrentGeo(); // read from currently loaded geometry
void ReadRoot(char *filename); // read from root file
void ReadDB(char *filename); // read from DB file
void ReadDB(char *db, char *path, int run, int version=-1, int subversion=-1);
- void ReadGeo(char *misaligned); // read from misaligned_geometry.root
+ Bool_t DecodeSurveyPointName(TString pna, Int_t &sm, Int_t &iz,Int_t &ir, Int_t &iphi);
void ReadSurveyReport(char *filename); // read from survey report
+ void ReadSurveyReport(AliSurveyObj *so); // read from survey object
void ReadAny(char *filename); // read from any kind of file
// writing on file
void WriteAscii(char *filename) const; // store data on ascii file
void WriteRoot(char *filename); // store data on root file
void WriteDB(char *filename, int run0, int run1); // store data on a local DB-like file
- void WriteDB(char *db, char *path, int run0, int run1);
- // store data on DB file
+ void WriteDB(char *db, char *pa, int r0, int r1); // store data on DB file
void WriteGeo(char *filename); // apply misalignment and store geometry
// geometry and symbolic names getters
// phi-sector number of chamber ch, 0-17
- int GetSec(int ch) const {return ch/30;}
+ int GetSec(int ch) const {return ch/30;}
// stack number, 0-4
- int GetSta(int ch) const {return ch%30/6;}
+ int GetSta(int ch) const {return ch%30/6;}
// plane number, 0-5
- int GetPla(int ch) const {return ch%30%6;}
+ int GetPla(int ch) const {return ch%30%6;}
// module number, 0-89
- int GetMod(int ch) const {return 5*GetSec(ch)+GetSta(ch); }
+ int GetMod(int ch) const {return 5*GetSec(ch)+GetSta(ch);}
// layer number, 9-14
- int GetLay(int ch) const {return AliGeomManager::kTRD1+GetPla(ch); }
+ int GetLay(int ch) const {return AliGeomManager::kTRD1+GetPla(ch);}
// volume id
- UShort_t GetVoi(int ch) const {return AliGeomManager::LayerToVolUID(GetLay(ch),GetMod(ch)); }
+ UShort_t GetVoi(int ch) const {return AliGeomManager::LayerToVolUID(GetLay(ch),GetMod(ch));}
// symbolic name of a supermodule
- char *GetSmName(int sm) const {return Form("TRD/sm%02d",sm); }
+ char *GetSmName(int sm) const {return Form("TRD/sm%02d",sm);}
// symbolic name of a chamber
- char *GetChName(int ch) const {return Form("TRD/sm%02d/st%d/pl%d",GetSec(ch),GetSta(ch),GetPla(ch)); }
+ char *GetChName(int ch) const {return Form("TRD/sm%02d/st%d/pl%d",GetSec(ch),GetSta(ch),GetPla(ch));}
+ // index of a supermodule
+ int GetSmIndex(const char *name) {for (int i=0; i<18; i++) if (strcmp(name,GetSmName(i))==0) return i; return -1;}
+ // index of a chamber
+ int GetChIndex(const char *name) {for (int i=0; i<540; i++) if (strcmp(name,GetChName(i))==0) return i; return -1;}
// data analysis
void ArToNumbers(TClonesArray *ar); // read ar and fill fSm and fCh
void NumbersToAr(TClonesArray *ar); // build ar using fSm and fCh data
- void LoadIdealGeometry(char *filename); // load ideal geometry from file
- void LoadIdealGeometry() {LoadIdealGeometry("geometry.root");}
+ int IsGeoLoaded(); // check if geometry is loaded
protected:
double fSurveyX[18][2][2][2]; // supermodule survey point X
double fSurveyY[18][2][2][2]; // supermodule survey point Y
double fSurveyZ[18][2][2][2]; // supermodule survey point Z
- double fSurveyE[18][2][2][2]; // supermodule survey point error
+ double fSurveyEX[18][2][2][2]; // supermodule survey point X error
+ double fSurveyEY[18][2][2][2]; // supermodule survey point Y error
+ double fSurveyEZ[18][2][2][2]; // supermodule survey point Z error
double fSurveyX0[2][2][2]; // ideal X position of the survey marks
double fSurveyY0[2][2][2]; // ideal Y position of the survey marks
double fSurveyZ0[2][2][2]; // ideal Z position of the survey marks