[u/mrichter/AliRoot.git] / ACORDE / AliACORDEAlign.cxx

/**************************************************************************
 * 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.                  *
 **************************************************************************/

#include "AliACORDEAlign.h"
#include "TROOT.h"
#include "Riostream.h"
#include "TFile.h"
#include "TMath.h"
#include "TSystem.h"
#include "AliSurveyObj.h"
#include "AliAlignObjParams.h"
#include "AliCDBStorage.h"
#include <TClonesArray.h>
#include <TFile.h>
#include "AliLog.h"
#include "AliCDBManager.h"
#include "AliSurveyPoint.h" 
#include "AliACORDEConstants.h" 

// Class creating the ACORDE aligmnent objects 
// from the surveys done by surveyers at Point2.
// Survey results are fetched from 
// Survey Depot, based on survey results 
// position of ACORDE alignment objects is computed.


ClassImp(AliACORDEAlign)


//________________________________________________________________________
AliACORDEAlign::AliACORDEAlign() :
  TObject(),
  fFileGlob(0x0),
  fRepLoc(0),
  fRepGlob(0),
  fUser(0x0),  
  fX(),
  fAlignACORDEObjArray(0x0),
  fDebug(0)
{
  //
  //  default constructor
  //
}

AliACORDEAlign::AliACORDEAlign(/*Int_t reportloc,*/Int_t reportglob):
  TObject(),
  fFileGlob(0x0),
  fRepLoc(0),
  fRepGlob(0),
  fUser(0x0),
  fX(120,4),
  fAlignACORDEObjArray(0x0),
  fDebug(0)
{
  
//
  // constructor
  //fRepLoc = new reportloc[80];
  //fRepGlob = new reportglob[80];
  Char_t path[50];
   fFileGlob = new Char_t[80];
   fUser = new Char_t[10];
  snprintf(path,50,gSystem->Getenv("ALICE_ROOT")); 
  // 
  snprintf(fFileGlob,80,"%s/ACORDE/Survey_%d_ACORDE.txt",path,reportglob);
  //
 snprintf(fUser,10,gSystem->Getenv("alien_API_USER"));

}


AliACORDEAlign::AliACORDEAlign(const AliACORDEAlign &align):
  TObject(),
  fFileGlob(0x0),
  fRepLoc(0),
  fRepGlob(0),
  fUser(0x0),
  fX(),
  fAlignACORDEObjArray(0x0),
  fDebug(0)
{
  //
  //  default copy constructor
fDebug = align.fDebug;
}

//__________________________________________________________________________
AliACORDEAlign & AliACORDEAlign::operator =(const AliACORDEAlign &align)
 

{
  //
  // assignment operator - dummy

 //
fDebug = align.fDebug;
  return (*this);
}

//__________________________________________________________________________
AliACORDEAlign::~AliACORDEAlign(){
  //
  // destructor
  //
if(fAlignACORDEObjArray) delete fAlignACORDEObjArray;
  if(fFileGlob) delete[] fFileGlob;
  if(fUser) delete[] fUser;
}


void AliACORDEAlign::LoadSurveyData()
{

//
// Create a new survey object and fill it.
 
AliSurveyObj * s1 = new AliSurveyObj(); 

if(fRepLoc != 0) 
 { 
 // Filling from DCDB (via GRID)
 s1->SetGridUser(fUser);
 s1->Fill("ACORDE",1014872,1,fUser); 
 }
 else
 {
   s1->FillFromLocalFile(fFileGlob);
 }


 //s1->GetEntries();
 //s1->GetUnits();
 //TObjArray* arr = s1->GetData();
 //cout<< "number of entries " << arr->GetEntries() <<endl;
 //arr->UncheckedAt(0)->ClassName();
 //AliSurveyPoint *sp0 = (AliSurveyPoint*) arr->UncheckedAt(0);   
 //cout << "point name " << sp0->GetPointName() << endl  ;
  

//
TString ML= "M" ;
//TString PL= "P";
TString underscore =  "_";
TString  endInner =  "_I";
TString  endOuter =  "_O";
TString  endCenter = "_P";
//
TString surveyname;
TString surveynameInner;
TString surveynameOuter;
TString surveynameCenter;

//TString surveynameAngles;
// 
TString pointNamesInner[60];
TString pointNamesOuter[60];
TString pointNamesCenter[60];

//
Int_t  nid=0;
//
 //for regular modules 
 for (Int_t ncolum=0; ncolum<6; ncolum++)
   {
     for (Int_t nrow=0; nrow<10; nrow++)
       {	
         
         surveyname=ML;
         surveyname+=ncolum;
         surveyname+=underscore;
         surveyname+=nrow;
  
         surveynameInner=surveyname; 
         surveynameInner+=endInner;

         surveynameOuter=surveyname;
         surveynameOuter+=endOuter; 
                          
         surveynameCenter=surveyname;
         surveynameCenter+=endCenter;
  
        pointNamesInner[nid] =  surveynameInner;
        pointNamesOuter[nid] = surveynameOuter; 
        pointNamesCenter[nid] = surveynameCenter;
  	 ++nid; 
       }
   }


//Read  two points 
AliSurveyPoint  *InnerPoint;
AliSurveyPoint *OuterPoint
; 
AliSurveyPoint  *CenterPoint;


 for(Int_t i=0;i<60;i++)
 {

   InnerPoint=0;
   OuterPoint=0; 
   CenterPoint=0;
 
   InnerPoint = (AliSurveyPoint *) s1->GetData()->FindObject(pointNamesInner[i]);
   OuterPoint = (AliSurveyPoint *) s1->GetData()->FindObject(pointNamesOuter[i]);
   CenterPoint = (AliSurveyPoint *) s1->GetData()->FindObject(pointNamesCenter[i]);


  if(InnerPoint && OuterPoint)
   {     
     //Use center if it is available
     if(CenterPoint)
       { 
         fX(i+60,0) =  CenterPoint->GetX()*100;  
         fX(i+60,1) =  CenterPoint->GetY()*100; 
         fX(i+60,2) =  CenterPoint->GetZ()*100;
        }
      else
        {
      //calculate center point 
         fX(i+60,0) = 100*(InnerPoint->GetX() + OuterPoint->GetX())/2.0;  
         fX(i+60,1) = 100*(InnerPoint->GetY() + OuterPoint->GetY())/2.0; 
         fX(i+60,2) = 100*(InnerPoint->GetZ() + OuterPoint->GetZ())/2.0;
        } 
   
      fX(i,0) =  OuterPoint->GetX()*100;  
      fX(i,1) =  OuterPoint->GetY()*100;
      fX(i,2) =  OuterPoint->GetZ()*100;
   }
   else 
   {
       if(InnerPoint && CenterPoint) 
         {

          fX(i+60,0) =  CenterPoint->GetX()*100;  
          fX(i+60,1) =  CenterPoint->GetY()*100; 
          fX(i+60,2) =  CenterPoint->GetZ()*100;

          fX(i,0) =  InnerPoint->GetX()*100;  
          fX(i,1) =  InnerPoint->GetY()*100;
          fX(i,2) =  InnerPoint->GetZ()*100;   
         } 
        else
        {
          if(OuterPoint && CenterPoint)
            { 
        
             fX(i+60,0) =  CenterPoint->GetX()*100;  
             fX(i+60,1) =  CenterPoint->GetY()*100; 
             fX(i+60,2) =  CenterPoint->GetZ()*100;

             fX(i,0) =  OuterPoint->GetX()*100;  
             fX(i,1) =  OuterPoint->GetY()*100;
             fX(i,2) =  OuterPoint->GetZ()*100;   
            }
          else
            { 

             fX(i+60,0) = -99.0;  
             fX(i+60,1) = -99.0; 
             fX(i+60,2) = -99.0;

             fX(i,0) =  -99.0;  
             fX(i,1) =  -99.0;
             fX(i,2) =  -99.0;   
            
            }
        }  
   } 


 }//ends  for

 delete s1;

}

void  AliACORDEAlign::ComputePosition()
{


//Residuals for rotations

Double_t theta;
Double_t resphi[60]; 
Double_t resiphi;

for (Int_t imod=0; imod<60; imod++)
  {
   if(TMath::Abs(fX(imod+60,0)-fX(imod,0))>=0.000001)
   {
   theta = (fX(imod+60,1)-fX(imod,1))/(fX(imod+60,0)-fX(imod,0));
   resiphi = TMath::ATan(theta)*(180.0/TMath::Pi());
   // calculate the residuals  special  modules 
   if(imod==0 || imod==9 || imod==50 || imod==59 )
    {    
    resphi[imod] = 0.0-resiphi;
    continue; 
    }
   // for module with no measurements 
   if(imod == 42 )
    {
     resphi[imod]= 0.0;
    continue;
    } 
   //face A
   if(imod>0 && imod <20)
    {
    resphi[imod] = resiphi + 45.0;  
    }
   //face B
   if(imod>=20 && imod <40)
    {
    resphi[imod] = -resiphi;  
    }
   //face C
   if(imod>=40 && imod <60)
    {
    resphi[imod] = resiphi - 45.0;  
    }
 
   }

}


//Get the  residuals for translations 

AliCDBManager* cdb = AliCDBManager::Instance();
if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
cdb->SetRun(0);

//AliCDBStorage* storage;
AliGeomManager::LoadGeometry(); 


TString symname;
TString basename = "ACORDE/Array";


// Get the ideal directly from the geometry 
 Double_t *tr;
 TGeoHMatrix *matrix;  
  for (Int_t imod=0; imod<60; imod++)
  {
    symname = basename;
    symname += imod; 
    cout<< symname << endl;
    matrix = AliGeomManager::GetMatrix(symname);
    tr=matrix->GetTranslation();  

    if(imod == 42)
      {
    fX(imod+60,0) = 0.0;  
    fX(imod+60,1) = 0.0;  
    fX(imod+60,2) = 0.0; 
    fX(imod,0) = 0.0;  
    fX(imod,1) = 0.0;  
    fX(imod,2) = 0.0; 
    continue;
      }

    fX(imod+60,0)=fX(imod+60,0)- tr[0];
    fX(imod+60,1)=fX(imod+60,1)- tr[1]- 4.0;
    fX(imod+60,2)=fX(imod+60,2)- tr[2];  
  
    fX(imod,0) = resphi[imod];  
    fX(imod,1) = 0.0;  
    fX(imod,2) = 0.0; 

   }
  

}

//______________________________________________________________________
void AliACORDEAlign::Run(){
  //
  // runs the full chain
  //
  
  //if(!LoadSurveyFromAlienFile("ACORDE",999999,1))
  //{
   // cout<<"Missing points"<<endl;
    //return;
  //}
  //else 
  //{
    //LoadSurveyfromLocalFile("ACORDE",99999,1);
  //} 


  LoadSurveyData();
  ComputePosition();
  //CreateACORDEAlignObjs();
  StoreAlignObj();

}

//_________________________________________________________________________

void AliACORDEAlign::StoreAlignObj()
{
  //
  // Storing ACORDE alignment objects 
  //

  AliCDBManager* cdb = AliCDBManager::Instance();
  if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://$ALICE_ROOT");

  TClonesArray *array = new TClonesArray("AliAlignObjParams",60);
  //
  // storing either in the OCDB or local file
  //

  TString symname;
  TString basename = "ACORDE/Array";
  Int_t iIndex=0; 

  AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer;
  UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iIndex);

  Double_t dx=0., dy=0., dz=0., dpsi=0., dtheta=0., dphi=0.0;
  
   for (Int_t imod=0; imod<60; imod++)
     {
       
       dphi = fX(imod,0);
       dtheta = fX(imod,1);
       dpsi = fX(imod,2);
       dx = fX(imod+60,0);
       dy = fX(imod+60,1);
       dz = fX(imod+60,2);  
       symname = basename;
       symname +=  imod;         
       new((*array)[imod]) AliAlignObjParams(symname,volid,dx,dy,dz,dpsi,dtheta,dphi,kFALSE);     
     }


  if( TString(gSystem->Getenv("TOCDB"))!= TString("kTRUE") )
   {
   
    
 // save on file
    const char* filename = "ACORDESurveyMisalignment.root";
    Char_t fullname[80];

    
    snprintf(fullname,80,filename);
       
   
    TFile *f = new TFile(fullname,"RECREATE");

    
    if(!f)
      {
	AliError("cannot open file for output\n");
	return;
      }
    AliInfo(Form("Saving alignment objects to the file %s", filename));
    f->cd();
    f->WriteObject(array,"ACORDEAlignObjs","kSingleKey");
    f->Close();
  }
  else
    {
      // save in CDB storage
      AliCDBStorage* storage;
      //
       TString Storage = gSystem->Getenv("STORAGE");
       if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://"))
	 {
	   AliError(Form("STORAGE variable set to %s is not valid. Exiting\n",Storage.Data()));
	   return;
	 }
       storage = cdb->GetStorage(Storage.Data());
       if(!storage)
	 {
	   AliError(Form("Unable to open storage %s\n",Storage.Data()));
	   return;
	 }
       //
       AliCDBMetaData* md = new AliCDBMetaData();
       md->SetResponsible("Pedro Podesta");
       md->SetComment("Full misalignment of ACORDE from surveyors");
       AliCDBId id("ACORDE/Align/Data",0,AliCDBRunRange::Infinity());
       storage->Put(fAlignACORDEObjArray,id,md);
    }

}
Commit	Line	Data
091f76f8	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	#include "AliACORDEAlign.h"
	17	#include "TROOT.h"
	18	#include "Riostream.h"
	19	#include "TFile.h"
	20	#include "TMath.h"
	21	#include "TSystem.h"
	22	#include "AliSurveyObj.h"
	23	#include "AliAlignObjParams.h"
	24	#include "AliCDBStorage.h"
	25	#include <TClonesArray.h>
	26	#include <TFile.h>
	27	#include "AliLog.h"
	28	#include "AliCDBManager.h"
	29	#include "AliSurveyPoint.h"
	30	#include "AliACORDEConstants.h"
	31
	32	// Class creating the ACORDE aligmnent objects
	33	// from the surveys done by surveyers at Point2.
	34	// Survey results are fetched from
	35	// Survey Depot, based on survey results
	36	// position of ACORDE alignment objects is computed.
	37
	38
	39	ClassImp(AliACORDEAlign)
	40
	41
	42	//________________________________________________________________________
	43	AliACORDEAlign::AliACORDEAlign() :
	44	TObject(),
	45	fFileGlob(0x0),
	46	fRepLoc(0),
	47	fRepGlob(0),
	48	fUser(0x0),
	49	fX(),
	50	fAlignACORDEObjArray(0x0),
	51	fDebug(0)
	52	{
	53	//
	54	// default constructor
	55	//
	56	}
	57
8f63aa5d	58	AliACORDEAlign::AliACORDEAlign(/Int_t reportloc,/Int_t reportglob):
091f76f8	59	TObject(),
8f63aa5d	60	fFileGlob(0x0),
	61	fRepLoc(0),
	62	fRepGlob(0),
	63	fUser(0x0),
091f76f8	64	fX(120,4),
	65	fAlignACORDEObjArray(0x0),
	66	fDebug(0)
	67	{
	68
	69	//
	70	// constructor
	71	//fRepLoc = new reportloc[80];
	72	//fRepGlob = new reportglob[80];
	73	Char_t path[50];
8f63aa5d	74	fFileGlob = new Char_t[80];
8f63aa5d	75	fUser = new Char_t[10];
091f76f8	76	snprintf(path,50,gSystem->Getenv("ALICE_ROOT"));
	77	//
	78	snprintf(fFileGlob,80,"%s/ACORDE/Survey_%d_ACORDE.txt",path,reportglob);
	79	//
10e05d9d	80	snprintf(fUser,10,gSystem->Getenv("alien_API_USER"));
091f76f8	81
	82	}
	83
	84
	85
	86
	87	AliACORDEAlign::AliACORDEAlign(const AliACORDEAlign &align):
	88	TObject(),
	89	fFileGlob(0x0),
	90	fRepLoc(0),
	91	fRepGlob(0),
	92	fUser(0x0),
	93	fX(),
	94	fAlignACORDEObjArray(0x0),
	95	fDebug(0)
	96	{
	97	//
	98	// default copy constructor
	99	fDebug = align.fDebug;
	100	}
	101
	102	//__________________________________________________________________________
	103	AliACORDEAlign & AliACORDEAlign::operator =(const AliACORDEAlign &align)
	104
	105
	106	{
	107	//
	108	// assignment operator - dummy
	109
	110	//
	111	fDebug = align.fDebug;
	112	return (*this);
	113	}
	114
	115	//__________________________________________________________________________
	116	AliACORDEAlign::~AliACORDEAlign(){
	117	//
	118	// destructor
	119	//
	120	if(fAlignACORDEObjArray) delete fAlignACORDEObjArray;
	121	if(fFileGlob) delete[] fFileGlob;
	122	if(fUser) delete[] fUser;
	123	}
	124
	125
	126	void AliACORDEAlign::LoadSurveyData()
	127	{
	128
	129	//
	130	// Create a new survey object and fill it.
	131
	132	AliSurveyObj * s1 = new AliSurveyObj();
	133
	134	if(fRepLoc != 0)
	135	{
	136	// Filling from DCDB (via GRID)
	137	s1->SetGridUser(fUser);
	138	s1->Fill("ACORDE",1014872,1,fUser);
	139	}
	140	else
	141	{
	142	s1->FillFromLocalFile(fFileGlob);
	143	}
	144
145
146	//s1->GetEntries();
147	//s1->GetUnits();
0415bcff	148	//TObjArray* arr = s1->GetData();
091f76f8	149	//cout<< "number of entries " << arr->GetEntries() <<endl;
	150	//arr->UncheckedAt(0)->ClassName();
	151	//AliSurveyPoint sp0 = (AliSurveyPoint) arr->UncheckedAt(0);
	152	//cout << "point name " << sp0->GetPointName() << endl ;
	153
	154
	155	//
	156	TString ML= "M" ;
	157	//TString PL= "P";
	158	TString underscore = "_";
	159	TString endInner = "_I";
	160	TString endOuter = "_O";
	161	TString endCenter = "_P";
	162	//
	163	TString surveyname;
	164	TString surveynameInner;
	165	TString surveynameOuter;
	166	TString surveynameCenter;
	167
	168	//TString surveynameAngles;
	169	//
	170	TString pointNamesInner[60];
	171	TString pointNamesOuter[60];
	172	TString pointNamesCenter[60];
	173
	174	//
	175	Int_t nid=0;
	176	//
	177	//for regular modules
	178	for (Int_t ncolum=0; ncolum<6; ncolum++)
	179	{
	180	for (Int_t nrow=0; nrow<10; nrow++)
	181	{
	182
	183	surveyname=ML;
	184	surveyname+=ncolum;
	185	surveyname+=underscore;
	186	surveyname+=nrow;
	187
	188	surveynameInner=surveyname;
	189	surveynameInner+=endInner;
	190
	191	surveynameOuter=surveyname;
	192	surveynameOuter+=endOuter;
	193
	194	surveynameCenter=surveyname;
	195	surveynameCenter+=endCenter;
	196
	197	pointNamesInner[nid] = surveynameInner;
	198	pointNamesOuter[nid] = surveynameOuter;
	199	pointNamesCenter[nid] = surveynameCenter;
	200	++nid;
	201	}
	202	}
	203
	204
	205	//Read two points
	206	AliSurveyPoint *InnerPoint;
8f63aa5d	207	AliSurveyPoint *OuterPoint
8f63aa5d	208	;
091f76f8	209	AliSurveyPoint *CenterPoint;
	210
	211
	212	for(Int_t i=0;i<60;i++)
	213	{
	214
	215	InnerPoint=0;
	216	OuterPoint=0;
	217	CenterPoint=0;
	218
	219	InnerPoint = (AliSurveyPoint *) s1->GetData()->FindObject(pointNamesInner[i]);
	220	OuterPoint = (AliSurveyPoint *) s1->GetData()->FindObject(pointNamesOuter[i]);
	221	CenterPoint = (AliSurveyPoint *) s1->GetData()->FindObject(pointNamesCenter[i]);
	222
	223
	224
	225	if(InnerPoint && OuterPoint)
	226	{
	227	//Use center if it is available
	228	if(CenterPoint)
	229	{
	230	fX(i+60,0) = CenterPoint->GetX()*100;
	231	fX(i+60,1) = CenterPoint->GetY()*100;
	232	fX(i+60,2) = CenterPoint->GetZ()*100;
	233	}
	234	else
	235	{
	236	//calculate center point
	237	fX(i+60,0) = 100*(InnerPoint->GetX() + OuterPoint->GetX())/2.0;
	238	fX(i+60,1) = 100*(InnerPoint->GetY() + OuterPoint->GetY())/2.0;
	239	fX(i+60,2) = 100*(InnerPoint->GetZ() + OuterPoint->GetZ())/2.0;
	240	}
	241
	242	fX(i,0) = OuterPoint->GetX()*100;
	243	fX(i,1) = OuterPoint->GetY()*100;
	244	fX(i,2) = OuterPoint->GetZ()*100;
	245	}
	246	else
	247	{
	248	if(InnerPoint && CenterPoint)
	249	{
	250
	251	fX(i+60,0) = CenterPoint->GetX()*100;
	252	fX(i+60,1) = CenterPoint->GetY()*100;
	253	fX(i+60,2) = CenterPoint->GetZ()*100;
	254
	255	fX(i,0) = InnerPoint->GetX()*100;
	256	fX(i,1) = InnerPoint->GetY()*100;
	257	fX(i,2) = InnerPoint->GetZ()*100;
	258	}
	259	else
	260	{
	261	if(OuterPoint && CenterPoint)
	262	{
	263
	264	fX(i+60,0) = CenterPoint->GetX()*100;
	265	fX(i+60,1) = CenterPoint->GetY()*100;
	266	fX(i+60,2) = CenterPoint->GetZ()*100;
	267
	268	fX(i,0) = OuterPoint->GetX()*100;
	269	fX(i,1) = OuterPoint->GetY()*100;
	270	fX(i,2) = OuterPoint->GetZ()*100;
	271	}
	272	else
273	{
274
275	fX(i+60,0) = -99.0;
276	fX(i+60,1) = -99.0;
277	fX(i+60,2) = -99.0;
278
279	fX(i,0) = -99.0;
280	fX(i,1) = -99.0;
281	fX(i,2) = -99.0;
282
283	}
284	}
285	}
286
287
288	}//ends for
289
290	delete s1;
291
292	}
293
294	void AliACORDEAlign::ComputePosition()
295	{
296
297
298	//Residuals for rotations
299
300	Double_t theta;
301	Double_t resphi[60];
302	Double_t resiphi;
303
304	for (Int_t imod=0; imod<60; imod++)
305	{
306	if(TMath::Abs(fX(imod+60,0)-fX(imod,0))>=0.000001)
307	{
308	theta = (fX(imod+60,1)-fX(imod,1))/(fX(imod+60,0)-fX(imod,0));
309	resiphi = TMath::ATan(theta)*(180.0/TMath::Pi());
310	// calculate the residuals special modules
311	if(imod==0 \|\| imod==9 \|\| imod==50 \|\| imod==59 )
312	{
313	resphi[imod] = 0.0-resiphi;
314	continue;
315	}
316	// for module with no measurements
317	if(imod == 42 )
318	{
319	resphi[imod]= 0.0;
320	continue;
321	}
322	//face A
323	if(imod>0 && imod <20)
324	{
325	resphi[imod] = resiphi + 45.0;
326	}
327	//face B
328	if(imod>=20 && imod <40)
329	{
330	resphi[imod] = -resiphi;
331	}
332	//face C
333	if(imod>=40 && imod <60)
334	{
335	resphi[imod] = resiphi - 45.0;
336	}
337
338	}
339
340	}
341
342
343	//Get the residuals for translations
344
345	AliCDBManager* cdb = AliCDBManager::Instance();
346	if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
347	cdb->SetRun(0);
348
349	//AliCDBStorage* storage;
350	AliGeomManager::LoadGeometry();
351
352
353	TString symname;
354	TString basename = "ACORDE/Array";
355
356
357	// Get the ideal directly from the geometry
358	Double_t *tr;
359	TGeoHMatrix *matrix;
360	for (Int_t imod=0; imod<60; imod++)
361	{
362	symname = basename;
363	symname += imod;
364	cout<< symname << endl;
365	matrix = AliGeomManager::GetMatrix(symname);
366	tr=matrix->GetTranslation();
367
368	if(imod == 42)
369	{
370	fX(imod+60,0) = 0.0;
371	fX(imod+60,1) = 0.0;
372	fX(imod+60,2) = 0.0;
373	fX(imod,0) = 0.0;
374	fX(imod,1) = 0.0;
375	fX(imod,2) = 0.0;
376	continue;
377	}
378
379	fX(imod+60,0)=fX(imod+60,0)- tr[0];
380	fX(imod+60,1)=fX(imod+60,1)- tr[1]- 4.0;
381	fX(imod+60,2)=fX(imod+60,2)- tr[2];
382
383	fX(imod,0) = resphi[imod];
384	fX(imod,1) = 0.0;
385	fX(imod,2) = 0.0;
386
387	}
388
389
390
391
392	}
393
394	//______________________________________________________________________
395	void AliACORDEAlign::Run(){
396	//
397	// runs the full chain
398	//
399
400	//if(!LoadSurveyFromAlienFile("ACORDE",999999,1))
401	//{
402	// cout<<"Missing points"<<endl;
403	//return;
404	//}
405	//else
406	//{
407	//LoadSurveyfromLocalFile("ACORDE",99999,1);
408	//}
409
410
411	LoadSurveyData();
412	ComputePosition();
413	//CreateACORDEAlignObjs();
414	StoreAlignObj();
415
416	}
417
418	//_________________________________________________________________________
419
420	void AliACORDEAlign::StoreAlignObj()
421	{
422	//
423	// Storing ACORDE alignment objects
424	//
425
426	AliCDBManager* cdb = AliCDBManager::Instance();
427	if(!cdb->IsDefaultStorageSet()) cdb->SetDefaultStorage("local://$ALICE_ROOT");
428
429	TClonesArray *array = new TClonesArray("AliAlignObjParams",60);
430	//
431	// storing either in the OCDB or local file
432	//
433
434	TString symname;
435	TString basename = "ACORDE/Array";
436	Int_t iIndex=0;
437
438	AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer;
439	UShort_t volid = AliGeomManager::LayerToVolUID(iLayer,iIndex);
440
441	Double_t dx=0., dy=0., dz=0., dpsi=0., dtheta=0., dphi=0.0;
442
443	for (Int_t imod=0; imod<60; imod++)
444	{
445
446	dphi = fX(imod,0);
447	dtheta = fX(imod,1);
448	dpsi = fX(imod,2);
449	dx = fX(imod+60,0);
450	dy = fX(imod+60,1);
451	dz = fX(imod+60,2);
452	symname = basename;
453	symname += imod;
454	new((*array)[imod]) AliAlignObjParams(symname,volid,dx,dy,dz,dpsi,dtheta,dphi,kFALSE);
455	}
456
457
458	if( TString(gSystem->Getenv("TOCDB"))!= TString("kTRUE") )
459	{
460
461
462
463	// save on file
464	const char* filename = "ACORDESurveyMisalignment.root";
465	Char_t fullname[80];
466
467
468
0415bcff	469	snprintf(fullname,80,filename);
091f76f8	470
	471
	472	TFile *f = new TFile(fullname,"RECREATE");
	473
	474
	475	if(!f)
	476	{
	477	AliError("cannot open file for output\n");
	478	return;
	479	}
	480	AliInfo(Form("Saving alignment objects to the file %s", filename));
	481	f->cd();
	482	f->WriteObject(array,"ACORDEAlignObjs","kSingleKey");
	483	f->Close();
	484	}
	485	else
	486	{
	487	// save in CDB storage
	488	AliCDBStorage* storage;
	489	//
	490	TString Storage = gSystem->Getenv("STORAGE");
	491	if(!Storage.BeginsWith("local://") && !Storage.BeginsWith("alien://"))
	492	{
	493	AliError(Form("STORAGE variable set to %s is not valid. Exiting\n",Storage.Data()));
	494	return;
	495	}
	496	storage = cdb->GetStorage(Storage.Data());
	497	if(!storage)
	498	{
	499	AliError(Form("Unable to open storage %s\n",Storage.Data()));
	500	return;
	501	}
	502	//
	503	AliCDBMetaData* md = new AliCDBMetaData();
	504	md->SetResponsible("Pedro Podesta");
	505	md->SetComment("Full misalignment of ACORDE from surveyors");
	506	AliCDBId id("ACORDE/Align/Data",0,AliCDBRunRange::Infinity());
	507	storage->Put(fAlignACORDEObjArray,id,md);
	508	}
	509
	510	}