[u/mrichter/AliRoot.git] / PHOS / macros / AlignmentDB / MakeFinalAlignment.C

/*
**********No misalignment*******
i=0 
matrix global_1 - tr=1  rot=1  refl=0  scl=0
  0.766044    0.000000    0.642788    Tx = 310.215729
  0.642788    0.000000   -0.766044    Ty = -369.700684
  0.000000    1.000000    0.000000    Tz =   0.000000
i=1 
matrix global_1 - tr=1  rot=1  refl=0  scl=0
  0.939693    0.000000    0.342020    Tx = 165.062332
  0.342020    0.000000   -0.939693    Ty = -453.505035
  0.000000    1.000000    0.000000    Tz =   0.000000
i=2 
matrix global_1 - tr=1  rot=1  refl=0  scl=0
  1.000000    0.000000    0.000000    Tx =   0.000000
  0.000000    0.000000   -1.000000    Ty = -482.609985
  0.000000    1.000000    0.000000    Tz =   0.000000
i=3 
matrix global_1 - tr=1  rot=1  refl=0  scl=0
  0.939693    0.000000   -0.342020    Tx = -165.062332
 -0.342020    0.000000   -0.939693    Ty = -453.505035
  0.000000    1.000000    0.000000    Tz =   0.000000

*************
=========Should be========

mod 4 (i=0)
matrix  - tr=1  rot=1  refl=0  scl=0
  0.767337   -0.004312    0.641229    Tx = 313.381746
  0.641189   -0.007942   -0.767342    Ty = -372.024764
  0.008402    0.999959   -0.003330    Tz =  -2.354855

mod 3 (i=1)
matrix  - tr=1  rot=1  refl=0  scl=0
  0.940798    0.002136    0.338961    Tx = 165.078518
  0.338959    0.001227   -0.940800    Ty = -454.148853
 -0.002426    0.999997    0.000431    Tz =  -2.801675

mod 3 (i=2)
matrix  - tr=1  rot=1  refl=0  scl=0
  0.999976   -0.000756   -0.006876    Tx =  -0.083626
 -0.006877   -0.001172   -0.999976    Ty = -482.172358
  0.000748    0.999999   -0.001178    Tz =  -1.906440

mod 3 (i=3)
matrix  - tr=1  rot=1  refl=0  scl=0
  0.934365   -0.002092   -0.356311    Tx = -167.233337
 -0.356313   -0.000801   -0.934366    Ty = -460.901723
  0.001670    0.999997   -0.001494    Tz =  -1.821245
  
------------Correct angles, but zero offsets------
  0.767338   -0.004312    0.641229    Tx = 309.463421
  0.641188   -0.007942   -0.767343    Ty = -370.327157
  0.008401    0.999959   -0.003329    Tz =  -1.606798
i=1 
matrix global_1 - tr=1  rot=1  refl=0  scl=0
  0.940798    0.002136    0.338961    Tx = 163.586079
  0.338959    0.001227   -0.940800    Ty = -454.039582
 -0.002425    0.999997    0.000430    Tz =   0.207687
i=2 
matrix global_1 - tr=1  rot=1  refl=0  scl=0
  0.999976   -0.000756   -0.006876    Tx =  -3.318427
 -0.006877   -0.001172   -0.999976    Ty = -482.598242
  0.000748    0.999999   -0.001177    Tz =  -0.568114
i=3 
matrix global_1 - tr=1  rot=1  refl=0  scl=0
  0.934365   -0.002092   -0.356311    Tx = -171.959374
 -0.356314   -0.000801   -0.934366    Ty = -450.934413
  0.001669    0.999997   -0.001494    Tz =  -0.720939

  
  
  
*/

void MakeFinalAlignment(){
  // Create ideal (no misalignment) object for PHOS

  
  const AliPHOSGeometry *phosGeom = AliPHOSGeometry::GetInstance("Run2", "Run2");
  if (!phosGeom) {
    Error(macroName, "Cannot obtain AliPHOSGeometry singleton.\n");
    return;
  }

  //Activate CDB storage and load geometry from CDB
  //[Part of code, taken from ITS version of MakexxxFullMisalignment
  AliCDBManager * cdb = AliCDBManager::Instance();
  cdb->SetDefaultStorage("local://./OCDB");
  cdb->SetRun(0);

  AliPHOSEMCAGeometry * emca = phosGeom->GetEMCAGeometry();
  TClonesArray alobj("AliAlignObjParams", 16);
   
  const Double_t dpsi = 0., dtheta = 0., dphi = 0.;
  Int_t iIndex = 0; //let all modules have index=0 in a layer with no LUT
  const AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer;
  UShort_t volid = AliGeomManager::LayerToVolUID(iLayer, iIndex);
  Int_t i = 0;

  // Alignment for 5 PHOS modules
 Double_t ideal1[9] = {0.766044,    0.000000,   0.642788,  
		       0.642788,    0.000000,  -0.766044,    
		       0.000000,    1.000000,   0.000000} ;    
 Double_t final1[9] = {0.767337,   -0.004312,   0.641229,   
		       0.641189,   -0.007942,  -0.767342, 
		       0.008402,    0.999959,  -0.003330} ;
  TGeoRotation rot1;
  rot1.SetMatrix(ideal1);
  TGeoRotation inv =rot1.Inverse() ;
  TGeoRotation rotF1;
  rotF1.SetMatrix(final1);
  rotF1.MultiplyBy(&inv) ;
  Double_t dX= 313.381746-309.463421; 
  Double_t dY=-372.024764+370.327157; 
  Double_t dZ=-2.354855 +1.606798 ;

 AliAlignObjParams * mod1 = new(alobj[i++]) AliAlignObjParams("PHOS/Module1",volid, dX, dY, dZ, 0., 0., 0., kTRUE);
  mod1->SetRotation(rotF1);

 
  Double_t ideal2[9] = {0.939693,    0.000000,    0.342020,
		        0.342020,    0.000000,   -0.939693,
		        0.000000,    1.000000,   0.000000} ;    
  Double_t final2[9] = {0.940798,    0.002136,    0.338961,
		        0.338959,    0.001227,   -0.940800, 
		       -0.002426,    0.999997,    0.000431} ;
  TGeoRotation rot2;
  rot2.SetMatrix(ideal2);
  TGeoRotation inv2 =rot2.Inverse() ;
  TGeoRotation rotF2;
  rotF2.SetMatrix(final2);
  rotF2.MultiplyBy(&inv2) ;
  dX= 165.078518-163.586079; 
  dY=-454.148853+454.039582; 
  dZ=-2.801675-0.207687 ;

  AliAlignObjParams * mod2 = new(alobj[i++]) AliAlignObjParams("PHOS/Module2",volid, dX, dY, dZ, 0., 0., 0., kTRUE);
  mod2->SetRotation(rotF2);
  
  Double_t ideal3[9] = {1.000000,    0.000000,    0.000000,
		        0.000000,    0.000000,   -1.000000,
		        0.000000,    1.000000,    0.000000} ;    
  Double_t final3[9] = {0.999976,   -0.000756,   -0.006876,
		       -0.006877,   -0.001172,   -0.999976,
		        0.000748,    0.999999,   -0.001178} ;
  TGeoRotation rot3;
  rot3.SetMatrix(ideal3);
  TGeoRotation inv3 =rot3.Inverse() ;
  TGeoRotation rotF3;
  rotF3.SetMatrix(final3);
  rotF3.MultiplyBy(&inv3) ;
  dX= -0.083626+3.318427; 
  dY=-482.172358+482.598242; 
  dZ=-1.906440+0.568114 ;

  AliAlignObjParams * mod3 = new(alobj[i++]) AliAlignObjParams("PHOS/Module3",volid, dX, dY, dZ, 0., 0., 0., kTRUE);
  mod3->SetRotation(rotF3);
  
 Double_t ideal4[9] = {0.939693,    0.000000,  -0.342020,
		      -0.342020,    0.000000,   -0.939693,
		        0.000000,    1.000000,    0.000000} ;    
  Double_t final4[9] = {0.934365,   -0.002092,   -0.356311,
		       -0.356313,   -0.000801,   -0.934366,
		        0.001670,    0.999997,   -0.001494} ;
  TGeoRotation rot4;
  rot4.SetMatrix(ideal4);
  TGeoRotation inv4 =rot4.Inverse() ;
  TGeoRotation rotF4;
  rotF4.SetMatrix(final4);
  rotF4.MultiplyBy(&inv4) ;
  dX=-167.233337+171.959374; 
  dY=-460.901723+450.934413; 
  dZ=-1.821245+0.720939 ;
  
  AliAlignObjParams * mod4 = new(alobj[i++]) AliAlignObjParams("PHOS/Module4",volid, dX, dY, dZ, 0., 0., 0., kTRUE);
  mod4->SetRotation(rotF4);

  new(alobj[i++]) AliAlignObjParams("PHOS/Module5",
	  volid, 0., 0., 0., 0., 0., 0., kTRUE);

  const Double_t dx = 0., dy = 0., dz = 0. ;
  // Alignment of CPV modules
  new(alobj[i++]) AliAlignObjParams("PHOS/Module1/CPV",
        volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  new(alobj[i++]) AliAlignObjParams("PHOS/Module2/CPV",
        volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  new(alobj[i++]) AliAlignObjParams("PHOS/Module3/CPV",
        volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  new(alobj[i++]) AliAlignObjParams("PHOS/Module4/CPV",
        volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
  new(alobj[i++]) AliAlignObjParams("PHOS/Module5/CPV",
        volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
 
  Double_t displacement=0 ;
  // Alignment for PHOS cradle
  new(alobj[i++]) AliAlignObjParams("PHOS/Cradle0",
	  volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
  new(alobj[i++]) AliAlignObjParams("PHOS/Cradle1",
	  volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);

  // Alignment for cradle wheels
  new(alobj[i++]) AliAlignObjParams("PHOS/Wheel0",
	  volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
  new(alobj[i++]) AliAlignObjParams("PHOS/Wheel1",
	  volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
  new(alobj[i++]) AliAlignObjParams("PHOS/Wheel2",
	  volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
  new(alobj[i++]) AliAlignObjParams("PHOS/Wheel3",
	  volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);

  AliCDBMetaData md;
  md.SetResponsible("Dmitri Peressounko");
  md.SetComment("Ideal alignment objects for PHOS");
  md.SetAliRootVersion(gSystem->Getenv("ARVERSION"));
  AliCDBId id("PHOS/Align/Data",0,AliCDBRunRange::Infinity());
  cdb->Put(&alobj, id, &md);

  alobj.Delete();
}
Commit	Line	Data
	1	/*
	2	********No misalignment*****
	3	i=0
	4	matrix global_1 - tr=1 rot=1 refl=0 scl=0
	5	0.766044 0.000000 0.642788 Tx = 310.215729
	6	0.642788 0.000000 -0.766044 Ty = -369.700684
	7	0.000000 1.000000 0.000000 Tz = 0.000000
	8	i=1
	9	matrix global_1 - tr=1 rot=1 refl=0 scl=0
	10	0.939693 0.000000 0.342020 Tx = 165.062332
	11	0.342020 0.000000 -0.939693 Ty = -453.505035
	12	0.000000 1.000000 0.000000 Tz = 0.000000
	13	i=2
	14	matrix global_1 - tr=1 rot=1 refl=0 scl=0
	15	1.000000 0.000000 0.000000 Tx = 0.000000
	16	0.000000 0.000000 -1.000000 Ty = -482.609985
	17	0.000000 1.000000 0.000000 Tz = 0.000000
	18	i=3
	19	matrix global_1 - tr=1 rot=1 refl=0 scl=0
	20	0.939693 0.000000 -0.342020 Tx = -165.062332
	21	-0.342020 0.000000 -0.939693 Ty = -453.505035
	22	0.000000 1.000000 0.000000 Tz = 0.000000
	23
	24	*************
	25	=========Should be========
	26
	27	mod 4 (i=0)
	28	matrix - tr=1 rot=1 refl=0 scl=0
	29	0.767337 -0.004312 0.641229 Tx = 313.381746
	30	0.641189 -0.007942 -0.767342 Ty = -372.024764
	31	0.008402 0.999959 -0.003330 Tz = -2.354855
	32
	33	mod 3 (i=1)
	34	matrix - tr=1 rot=1 refl=0 scl=0
	35	0.940798 0.002136 0.338961 Tx = 165.078518
	36	0.338959 0.001227 -0.940800 Ty = -454.148853
	37	-0.002426 0.999997 0.000431 Tz = -2.801675
	38
	39	mod 3 (i=2)
	40	matrix - tr=1 rot=1 refl=0 scl=0
	41	0.999976 -0.000756 -0.006876 Tx = -0.083626
	42	-0.006877 -0.001172 -0.999976 Ty = -482.172358
	43	0.000748 0.999999 -0.001178 Tz = -1.906440
	44
	45	mod 3 (i=3)
	46	matrix - tr=1 rot=1 refl=0 scl=0
	47	0.934365 -0.002092 -0.356311 Tx = -167.233337
	48	-0.356313 -0.000801 -0.934366 Ty = -460.901723
	49	0.001670 0.999997 -0.001494 Tz = -1.821245
	50
	51	------------Correct angles, but zero offsets------
	52	0.767338 -0.004312 0.641229 Tx = 309.463421
	53	0.641188 -0.007942 -0.767343 Ty = -370.327157
	54	0.008401 0.999959 -0.003329 Tz = -1.606798
	55	i=1
	56	matrix global_1 - tr=1 rot=1 refl=0 scl=0
	57	0.940798 0.002136 0.338961 Tx = 163.586079
	58	0.338959 0.001227 -0.940800 Ty = -454.039582
	59	-0.002425 0.999997 0.000430 Tz = 0.207687
	60	i=2
	61	matrix global_1 - tr=1 rot=1 refl=0 scl=0
	62	0.999976 -0.000756 -0.006876 Tx = -3.318427
	63	-0.006877 -0.001172 -0.999976 Ty = -482.598242
	64	0.000748 0.999999 -0.001177 Tz = -0.568114
	65	i=3
	66	matrix global_1 - tr=1 rot=1 refl=0 scl=0
	67	0.934365 -0.002092 -0.356311 Tx = -171.959374
	68	-0.356314 -0.000801 -0.934366 Ty = -450.934413
	69	0.001669 0.999997 -0.001494 Tz = -0.720939
	70
	71
	72
	73
	74	*/
	75
	76	void MakeFinalAlignment(){
	77	// Create ideal (no misalignment) object for PHOS
	78
	79
	80	const AliPHOSGeometry *phosGeom = AliPHOSGeometry::GetInstance("Run2", "Run2");
	81	if (!phosGeom) {
	82	Error(macroName, "Cannot obtain AliPHOSGeometry singleton.\n");
	83	return;
	84	}
	85
	86	//Activate CDB storage and load geometry from CDB
	87	//[Part of code, taken from ITS version of MakexxxFullMisalignment
	88	AliCDBManager * cdb = AliCDBManager::Instance();
	89	cdb->SetDefaultStorage("local://./OCDB");
	90	cdb->SetRun(0);
	91
	92	AliPHOSEMCAGeometry * emca = phosGeom->GetEMCAGeometry();
	93	TClonesArray alobj("AliAlignObjParams", 16);
	94
	95	const Double_t dpsi = 0., dtheta = 0., dphi = 0.;
	96	Int_t iIndex = 0; //let all modules have index=0 in a layer with no LUT
	97	const AliGeomManager::ELayerID iLayer = AliGeomManager::kInvalidLayer;
	98	UShort_t volid = AliGeomManager::LayerToVolUID(iLayer, iIndex);
	99	Int_t i = 0;
	100
	101	// Alignment for 5 PHOS modules
	102	Double_t ideal1[9] = {0.766044, 0.000000, 0.642788,
	103	0.642788, 0.000000, -0.766044,
	104	0.000000, 1.000000, 0.000000} ;
	105	Double_t final1[9] = {0.767337, -0.004312, 0.641229,
	106	0.641189, -0.007942, -0.767342,
	107	0.008402, 0.999959, -0.003330} ;
	108	TGeoRotation rot1;
	109	rot1.SetMatrix(ideal1);
	110	TGeoRotation inv =rot1.Inverse() ;
	111	TGeoRotation rotF1;
	112	rotF1.SetMatrix(final1);
	113	rotF1.MultiplyBy(&inv) ;
	114	Double_t dX= 313.381746-309.463421;
	115	Double_t dY=-372.024764+370.327157;
	116	Double_t dZ=-2.354855 +1.606798 ;
	117
	118	AliAlignObjParams * mod1 = new(alobj[i++]) AliAlignObjParams("PHOS/Module1",volid, dX, dY, dZ, 0., 0., 0., kTRUE);
	119	mod1->SetRotation(rotF1);
	120
	121
	122	Double_t ideal2[9] = {0.939693, 0.000000, 0.342020,
	123	0.342020, 0.000000, -0.939693,
	124	0.000000, 1.000000, 0.000000} ;
	125	Double_t final2[9] = {0.940798, 0.002136, 0.338961,
	126	0.338959, 0.001227, -0.940800,
	127	-0.002426, 0.999997, 0.000431} ;
	128	TGeoRotation rot2;
	129	rot2.SetMatrix(ideal2);
	130	TGeoRotation inv2 =rot2.Inverse() ;
	131	TGeoRotation rotF2;
	132	rotF2.SetMatrix(final2);
	133	rotF2.MultiplyBy(&inv2) ;
	134	dX= 165.078518-163.586079;
	135	dY=-454.148853+454.039582;
	136	dZ=-2.801675-0.207687 ;
	137
	138	AliAlignObjParams * mod2 = new(alobj[i++]) AliAlignObjParams("PHOS/Module2",volid, dX, dY, dZ, 0., 0., 0., kTRUE);
	139	mod2->SetRotation(rotF2);
	140
	141	Double_t ideal3[9] = {1.000000, 0.000000, 0.000000,
	142	0.000000, 0.000000, -1.000000,
	143	0.000000, 1.000000, 0.000000} ;
	144	Double_t final3[9] = {0.999976, -0.000756, -0.006876,
	145	-0.006877, -0.001172, -0.999976,
	146	0.000748, 0.999999, -0.001178} ;
	147	TGeoRotation rot3;
	148	rot3.SetMatrix(ideal3);
	149	TGeoRotation inv3 =rot3.Inverse() ;
	150	TGeoRotation rotF3;
	151	rotF3.SetMatrix(final3);
	152	rotF3.MultiplyBy(&inv3) ;
	153	dX= -0.083626+3.318427;
	154	dY=-482.172358+482.598242;
	155	dZ=-1.906440+0.568114 ;
	156
	157	AliAlignObjParams * mod3 = new(alobj[i++]) AliAlignObjParams("PHOS/Module3",volid, dX, dY, dZ, 0., 0., 0., kTRUE);
	158	mod3->SetRotation(rotF3);
	159
	160	Double_t ideal4[9] = {0.939693, 0.000000, -0.342020,
	161	-0.342020, 0.000000, -0.939693,
	162	0.000000, 1.000000, 0.000000} ;
	163	Double_t final4[9] = {0.934365, -0.002092, -0.356311,
	164	-0.356313, -0.000801, -0.934366,
	165	0.001670, 0.999997, -0.001494} ;
	166	TGeoRotation rot4;
	167	rot4.SetMatrix(ideal4);
	168	TGeoRotation inv4 =rot4.Inverse() ;
	169	TGeoRotation rotF4;
	170	rotF4.SetMatrix(final4);
	171	rotF4.MultiplyBy(&inv4) ;
	172	dX=-167.233337+171.959374;
	173	dY=-460.901723+450.934413;
	174	dZ=-1.821245+0.720939 ;
	175
	176	AliAlignObjParams * mod4 = new(alobj[i++]) AliAlignObjParams("PHOS/Module4",volid, dX, dY, dZ, 0., 0., 0., kTRUE);
	177	mod4->SetRotation(rotF4);
	178
	179	new(alobj[i++]) AliAlignObjParams("PHOS/Module5",
	180	volid, 0., 0., 0., 0., 0., 0., kTRUE);
	181
	182	const Double_t dx = 0., dy = 0., dz = 0. ;
	183	// Alignment of CPV modules
	184	new(alobj[i++]) AliAlignObjParams("PHOS/Module1/CPV",
	185	volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
	186	new(alobj[i++]) AliAlignObjParams("PHOS/Module2/CPV",
	187	volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
	188	new(alobj[i++]) AliAlignObjParams("PHOS/Module3/CPV",
	189	volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
	190	new(alobj[i++]) AliAlignObjParams("PHOS/Module4/CPV",
	191	volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
	192	new(alobj[i++]) AliAlignObjParams("PHOS/Module5/CPV",
	193	volid, dx, dy, dz, dpsi, dtheta, dphi, kTRUE);
	194
	195	Double_t displacement=0 ;
	196	// Alignment for PHOS cradle
	197	new(alobj[i++]) AliAlignObjParams("PHOS/Cradle0",
	198	volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
	199	new(alobj[i++]) AliAlignObjParams("PHOS/Cradle1",
	200	volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
	201
	202	// Alignment for cradle wheels
	203	new(alobj[i++]) AliAlignObjParams("PHOS/Wheel0",
	204	volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
	205	new(alobj[i++]) AliAlignObjParams("PHOS/Wheel1",
	206	volid, 0., 0., -displacement, dpsi, dtheta, dphi, kTRUE);
	207	new(alobj[i++]) AliAlignObjParams("PHOS/Wheel2",
	208	volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
	209	new(alobj[i++]) AliAlignObjParams("PHOS/Wheel3",
	210	volid, 0., 0., +displacement, dpsi, dtheta, dphi, kTRUE);
	211
	212	AliCDBMetaData md;
	213	md.SetResponsible("Dmitri Peressounko");
	214	md.SetComment("Ideal alignment objects for PHOS");
	215	md.SetAliRootVersion(gSystem->Getenv("ARVERSION"));
	216	AliCDBId id("PHOS/Align/Data",0,AliCDBRunRange::Infinity());
	217	cdb->Put(&alobj, id, &md);
	218
	219	alobj.Delete();
	220	}
	221
	222