New class AliITSSurveyToAlignSSD to convert the SSD survey data to alignment objects...

diff --git a/ITS/AliITSSurveyToAlignSSD.cxx b/ITS/AliITSSurveyToAlignSSD.cxx
new file mode 100644 (file)
index 0000000..5ac2ddf
--- /dev/null
+++ b/ITS/AliITSSurveyToAlignSSD.cxx
@@ -0,0 +1,223 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+/*
+$Log$
+Revision 1.1  2007/10/01 14:12:45  pchrist
+Class creating the aligmnent object fro the surveyor measurements.
+
+*/ 
+
+//
+//  Creates the SSD align object
+//
+//
+#include "Riostream.h"
+#include "TFile.h"
+#include "TSystem.h"
+#include "TClonesArray.h"
+#include "TGeoManager.h"
+#include "TGeoMatrix.h"
+#include "TGeoPhysicalNode.h"
+
+#include "AliITSSurveyToAlignSSD.h"
+#include "AliSurveyObj.h"
+#include "AliSurveyPoint.h"
+#include "AliAlignObjParams.h"
+
+#include "AliLog.h"
+
+#include "AliCDBManager.h"
+#include "AliCDBEntry.h"
+#include "AliCDBStorage.h"
+
+ClassImp(AliITSSurveyToAlignSSD)
+
+//________________________________________________________________________
+AliITSSurveyToAlignSSD::AliITSSurveyToAlignSSD() :
+  TObject(),
+  fRun(0),
+  fFileLoc(0x0),
+  fFileGlob(0x0),
+  fSurveyPoints(0),
+  fSSDAlignObj(new TClonesArray("AliAlignObjParams",100)),
+  fSSDAlignObjParam(0),
+  fDebug(0) {
+  //
+  //  default constructor
+  //
+}   
+
+//________________________________________________________________________
+AliITSSurveyToAlignSSD::AliITSSurveyToAlignSSD(Int_t run, Int_t reportloc, Int_t reportglob) :
+  TObject(),
+  fRun(0),
+  fFileLoc(0x0),
+  fFileGlob(0x0),
+  fSurveyPoints(0),
+  fSSDAlignObj(new TClonesArray("AliAlignObjParams",100)),
+  fSSDAlignObjParam(0),
+  fDebug(0) {
+  //
+  // constructor - defines data files
+  //
+  fFileLoc = new Char_t[80];
+  fFileGlob = new Char_t[80];
+  Char_t path[50];
+  sprintf(path,gSystem->Getenv("ALICE_ROOT")); 
+  //
+  sprintf(fFileLoc,"%s/ITS/Survey_SSD_%d.txt",path,reportloc);  
+  sprintf(fFileGlob,"%s/ITS/Survey_SSD_%d.txt",path,reportglob);
+  //
+
+}
+//_________________________________________________________________________
+AliITSSurveyToAlignSSD::AliITSSurveyToAlignSSD(const AliITSSurveyToAlignSSD &align) :
+  TObject(),
+  fFileLoc(0x0),
+  fFileGlob(0x0),
+  fSurveyPoints(0),
+  fSSDAlignObj(new TClonesArray("AliAlignObjParams",100)),
+  fSSDAlignObjParam(0),
+  fDebug(0) {
+  //
+  //  copy constructor - dummy
+  
+  fDebug = align.fDebug;
+}
+
+//__________________________________________________________________________
+AliITSSurveyToAlignSSD & AliITSSurveyToAlignSSD::operator =(const AliITSSurveyToAlignSSD & align) {
+  //
+  // assignment operator - dummy
+  //
+  fDebug = align.fDebug;
+  return (*this);
+}
+
+//__________________________________________________________________________
+AliITSSurveyToAlignSSD::~AliITSSurveyToAlignSSD() {
+  //
+  // destructor
+  //
+  if(fSurveyPoints) delete fSurveyPoints;
+  if(fSSDAlignObj) delete fSSDAlignObj;
+  if(fSSDAlignObjParam) delete fSSDAlignObjParam;
+}
+
+//______________________________________________________________________
+void AliITSSurveyToAlignSSD::Run() { 
+  //
+  // runs the full chain
+  //
+  SetDebug(0);
+  Bool_t flag = LoadSurveyData();
+  if(!flag) {
+    cout<<"Missing points"<<endl;
+    return;
+  }
+  CreateAlignObj();
+  StoreAlignObj();
+  //
+}
+
+//__________________________________________________________________________
+Bool_t AliITSSurveyToAlignSSD::LoadSurveyData() {
+  //
+  // for a time being it loads from the local file the surveyed point
+  // and has the ideal points hardwired. I am waiting until Ricardo
+  // completes his job
+  // 
+  
+  //Load survey data from the local file
+  AliSurveyObj * s1 = new AliSurveyObj();
+  if(s1->FillFromLocalFile(fFileGlob))
+    fSurveyPoints = s1->GetData();
+  else 
+    return kFALSE;
+  
+  cout<<"Number of SSD survey points: "<<fSurveyPoints->GetEntries()<<endl;  
+  
+  return kTRUE;
+}
+
+//_______________________________________________________________________
+void AliITSSurveyToAlignSSD::CreateAlignObj(){
+  //
+  // This method creates AliAlignObj and fills it with Euler angles
+  // and shifts. The units are degrees and cm.
+  // 
+  TClonesArray &alobj = *fSSDAlignObj;
+
+  //load ideal geometry from the OCDB
+  AliCDBManager *cdb = AliCDBManager::Instance();
+  cdb->SetDefaultStorage("local://$ALICE_ROOT");
+  cdb->SetRun(fRun);
+  AliCDBEntry* entry = cdb->Get("GRP/Geometry/Data");
+  AliGeomManager::SetGeometry(gGeoManager);
+
+  AliSurveyPoint* pt = 0;
+  TGeoHMatrix* hm = 0;
+  TGeoPNEntry* pne = 0;
+  Double_t* tr;
+  Int_t uid;
+  const char* symname;
+  Double_t sx, sy, sz;
+  Int_t ilayer, imodule;
+
+  for(Int_t imod = 0; imod < fSurveyPoints->GetEntries(); imod++) {
+    pt = (AliSurveyPoint*) fSurveyPoints->At(imod);
+    if(!pt) continue;
+    sx = pt->GetX()*100.;
+    sy = pt->GetY()*100.;
+    sz = pt->GetZ()*100.;
+
+    ilayer = (imod < 748) ? AliGeomManager::kSSD1 : AliGeomManager::kSSD2;
+    imodule = (imod < 748) ? imod : imod - 748;
+
+    uid = AliGeomManager::LayerToVolUID(ilayer,imodule);
+    symname = AliGeomManager::SymName(uid);
+    pne = gGeoManager->GetAlignableEntryByUID(uid);
+    hm = pne->GetGlobalOrig();
+    tr = hm->GetTranslation();
+
+    //Printf("symname %s", symname);
+    //Printf("x,y,z from survey: %f %f %f", sx, sy, sz);
+    //Printf("x,y,z from ideal : %f %f %f", tr[0], tr[1], tr[2]);
+    
+   new(alobj[imod]) AliAlignObjParams(symname, uid, sx-tr[0], sy-tr[1], sz-tr[2], 0., 0., 0., kTRUE);
+  }//module loop
+
+  delete entry;
+}
+
+//_________________________________________________________________________
+void AliITSSurveyToAlignSSD::StoreAlignObj(){
+  // Stores the TClonesArray of AliAlignObj in
+  // $ALICE_ROOT/ITS/Align/Data/SSDfromSurvey.root
+  // In a later version these objects will be merged 
+  // with the SDD ones and will be put in the OCDB
+  const char* filename = "$ALICE_ROOT/ITS/Align/Data/SSDfromSurvey.root";
+  TFile *f = TFile::Open(filename,"RECREATE");
+  if(!f){
+    Error(filename,"cannot open file for output\n");
+    return;
+  }
+  cout<<"Saving alignment objects to the file "<<filename<<endl;
+  f->cd();
+  f->WriteObject(fSSDAlignObj,"SSDAlignObjs","kSingleKey");
+  f->Close();
+
+}

diff --git a/ITS/AliITSSurveyToAlignSSD.h b/ITS/AliITSSurveyToAlignSSD.h
new file mode 100644 (file)
index 0000000..b1ea1c8
--- /dev/null
+++ b/ITS/AliITSSurveyToAlignSSD.h
@@ -0,0 +1,47 @@
+#ifndef ALIITSSURVEYTOALIGNSSD_H
+#define ALIITSSURVEYTOALIGNSSD_H
+
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/* $Id$ */
+
+/////////////////////////////////////////////////////////////////////////
+//     This class creates the alignment object from the survey data    //      
+//     for the SSD layers.                                             //
+//                                                                     //
+//     Author: Panos.Christakoglou@cern.ch - NIKHEF/UU                 //
+/////////////////////////////////////////////////////////////////////////
+#include <TObjArray.h>
+#include <TClonesArray.h>
+
+#include "AliAlignObjParams.h"
+
+
+class AliITSSurveyToAlignSSD : public TObject {
+
+ public:
+  AliITSSurveyToAlignSSD();
+  AliITSSurveyToAlignSSD(Int_t run, Int_t reportloc, Int_t reportglob);
+  AliITSSurveyToAlignSSD(const AliITSSurveyToAlignSSD &align); // copy constructor
+  AliITSSurveyToAlignSSD &operator = (const AliITSSurveyToAlignSSD &align); //assignment operator
+  Bool_t LoadSurveyData();
+  void CreateAlignObj();
+  void Run();
+  void SetDebug(Int_t debug){fDebug=debug;}
+  void StoreAlignObj();
+  virtual   ~AliITSSurveyToAlignSSD();
+  //
+ private:
+  Int_t   fRun;                         // the run number for the OCDB
+  Char_t *fFileLoc;                     // file with ideal points
+  Char_t *fFileGlob;                    // file with surveyed points
+  TObjArray *fSurveyPoints;             // array of survey points
+  TClonesArray *fSSDAlignObj;           // TClonesArray of AliAlignObj
+  AliAlignObjParams *fSSDAlignObjParam; // SSD alignment object param
+  Int_t fDebug;                         // debug flag
+  
+
+  ClassDef(AliITSSurveyToAlignSSD,0);
+};
+#endif

diff --git a/ITS/ITSbaseLinkDef.h b/ITS/ITSbaseLinkDef.h
index f2e000d29cb5571fae93445a954032d18334562d..9757ac05665135536231b656938ae38227bccd6e 100644 (file)
--- a/ITS/ITSbaseLinkDef.h
+++ b/ITS/ITSbaseLinkDef.h
@@ -96,6 +96,6 @@
 #pragma link C++ class AliITSIntMap+;
 #pragma link C++ class AliITSIntMapNode+;
 #pragma link C++ class AliITSPedestalSSD+;
-
+#pragma link C++ class AliITSSurveyToAlignSSD+;
 
 #endif

diff --git a/ITS/libITSbase.pkg b/ITS/libITSbase.pkg
index 26c44cfc0ed7ee42ba5b083409f36d76e5f04b39..dd1e1373bcbe3aaf109d8578aa1da85b2a34f3a5 100644 (file)
--- a/ITS/libITSbase.pkg
+++ b/ITS/libITSbase.pkg
@@ -63,7 +63,8 @@ SRCS =        AliITSgeom.cxx \
                AliITSDCSDataSDD.cxx \
                 AliITSIntMap.cxx \
                AliITSIntMapNode.cxx \
-               AliITSPedestalSSD.cxx
+               AliITSPedestalSSD.cxx \
+               AliITSSurveyToAlignSSD.cxx
 
 HDRS:=  $(SRCS:.cxx=.h)