[u/mrichter/AliRoot.git] / MUON / AliMUONTriggerGUIboard.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.                  *
 **************************************************************************/

// $Id$

//-----------------------------------------------------------------------------
/// \class AliMUONTriggerGUIboard
///
/// Single trigger board object with geometry information, strips and digits
///
/// \author Bogdan Vulpescu, LPC Clermont-Ferrand
//-----------------------------------------------------------------------------

#include <TClonesArray.h>
#include <TBox.h>
#include <TMath.h>

#include "AliMUONGeometryTransformer.h"

#include "AliMUONTriggerGUIboard.h"

/// \cond CLASSIMP
ClassImp(AliMUONTriggerGUIboard)
/// \endcond

//__________________________________________________________________________
AliMUONTriggerGUIboard::AliMUONTriggerGUIboard() 
  : TObject(),
    fName(0),
    fCrateName(0),
    fID(-1),
    fStatus(0),
    fPosition(0),
    fYOver(0),
    fXSix(0),
    fXSiy1(0),
    fXSiy2(0),
    fYSix1(0),
    fYSix2(0),
    fYSiy(0),
    fDetElemId(0),
    fIdCircuit(-1),
    fIsOpen(0),
    fNPadsX(),
    fNPadsY(),
    fPadsX(),
    fPadsY()
{
  /// board main constructor

  fName = new TString("");
  fCrateName = new TString("");

  for (Int_t i = 0; i < kNMT; i++) {
    fXCenter[i] = 0.;
    fYCenter[i] = 0.;
    fZCenter[i] = 0.;
    fXWidth[i]  = 0.;
    fYWidth[i]  = 0.;
    for (Int_t is = 0; is < kNS; is++) {
      fXDig[i][is] = 0;
      fYDig[i][is] = 0;
      fXDigBox[i][is] = new TBox(0,0,0,0);
      fYDigBox[i][is] = new TBox(0,0,0,0);
      fXDigBox[i][is]->SetBit(kCannotPick);
      fYDigBox[i][is]->SetBit(kCannotPick);
      fXDigBox[i][is]->SetFillStyle(1001);
      fYDigBox[i][is]->SetFillStyle(1001);
      fXDigBox[i][is]->SetFillColor(4);
      fYDigBox[i][is]->SetFillColor(4);
    }
  }

  fXSix  = -1;
  fXSiy1 = -1;
  fXSiy2 = -1;

  fYSix1 = -1;
  fYSix2 = -1;
  fYSiy  = -1;

  fDetElemId = -1;
  fIdCircuit = -1;

  fIsOpen = kFALSE;

  fYOver    = 0;
  fPosition = 0;

  for (Int_t i = 0; i < kNMT; i++) {
    fPadsX[i] = new TClonesArray("AliMpPad",16); fNPadsX[i] = 0;
    fPadsY[i] = new TClonesArray("AliMpPad",16); fNPadsY[i] = 0;
  }

}

//__________________________________________________________________________
AliMUONTriggerGUIboard::~AliMUONTriggerGUIboard() 
{
  /// board destructor

  delete fName;

  for (Int_t imt = 0; imt < kNMT; imt++) {
    for (Int_t is = 0; is < kNS; is++) {
      delete fXDigBox[imt][is];
      delete fYDigBox[imt][is];
    }
  }

}

//__________________________________________________________________________
void AliMUONTriggerGUIboard::SetXDigBox(Int_t imt, Int_t is, Double_t x1, Double_t y1, Double_t x2, Double_t y2)
{
  /// set coordinates of "is" x-strip box in chamber "imt"

  fXDigBox[imt][is]->SetX1(x1);
  fXDigBox[imt][is]->SetY1(y1);
  fXDigBox[imt][is]->SetX2(x2);
  fXDigBox[imt][is]->SetY2(y2);

}

//__________________________________________________________________________
void AliMUONTriggerGUIboard::SetYDigBox(Int_t imt, Int_t is, Double_t x1, Double_t y1, Double_t x2, Double_t y2)
{
  /// set coordinates of "is" y-strip box in chamber "imt"

  fYDigBox[imt][is]->SetX1(x1);
  fYDigBox[imt][is]->SetY1(y1);
  fYDigBox[imt][is]->SetX2(x2);
  fYDigBox[imt][is]->SetY2(y2);

}

//__________________________________________________________________________
void AliMUONTriggerGUIboard::ClearXDigits()
{
  /// delete the set x-digits

  for (Int_t imt = 0; imt < kNMT; imt++) {
    for (Int_t is = 0; is < kNS; is++) {
      fXDig[imt][is] = 0;
    }
  }

}

//__________________________________________________________________________
void AliMUONTriggerGUIboard::ClearYDigits()
{
  /// delete the set y-digits

  for (Int_t imt = 0; imt < kNMT; imt++) {
    for (Int_t is = 0; is < kNS; is++) {
      fYDig[imt][is] = 0;
    }
  }

}

//__________________________________________________________________________
void AliMUONTriggerGUIboard::MakeGeometry()
{
  /// create the display geometry from the mapping pads

  AliMpPad *pad;

  // circuit number and manu channel (from x-strips)
  for (Int_t ich = 0; ich < kNMT; ich++) {
    if (fNPadsX[ich]) {
      pad = (AliMpPad*)fPadsX[ich]->At(0);
      fIdCircuit = pad->GetLocalBoardId(0);
      break;
    }
  }

  // position index
  if (fName->Length()) {
    if (fName->Contains("12")) fPosition = 1;
    if (fName->Contains("34")) fPosition = 2;
    if (fName->Contains("56")) fPosition = 3;
    if (fName->Contains("78")) fPosition = 4;
  }

  // position index for common y-strip boards
  for (Int_t ich = 0; ich < kNMT; ich++) {
    if (fNPadsY[ich]) {
      pad = (AliMpPad*)fPadsY[ich]->At(0);
      fYOver = pad->GetNofLocations();
      break;
    }
  }

  // pad indices
  Int_t padxIx = -1, padxIy1 = +999, padxIy2 = -999;
  Int_t padyIy = -1, padyIx1 = +999, padyIx2 = -999;
  for (Int_t ip = 0; ip < fNPadsX[0]; ip++) {
    pad = (AliMpPad*)fPadsX[0]->At(ip);
    padxIx = pad->GetIx();
    padxIy1 = TMath::Min(padxIy1,pad->GetIy());
    padxIy2 = TMath::Max(padxIy2,pad->GetIy());
  }
  for (Int_t ip = 0; ip < fNPadsY[0]; ip++) {
    pad = (AliMpPad*)fPadsY[0]->At(ip);
    padyIy = pad->GetIy();
    padyIx1 = TMath::Min(padyIx1,pad->GetIx());
    padyIx2 = TMath::Max(padyIx2,pad->GetIx());
  }
  fXSix  = padxIx;
  fXSiy1 = padxIy1;
  fXSiy2 = padxIy2;
  fYSiy  = padyIy;
  fYSix1 = padyIx1;
  fYSix2 = padyIx2;

  // position and dimension

  AliMUONGeometryTransformer transformer;
  transformer.LoadGeometryData("transform.dat");

  Float_t minX, maxX, minY, maxY;
  Float_t dx, dy;
  Float_t xloc, yloc, xglo, yglo, zglo;
  for (Int_t ich = 0; ich < kNMT; ich++) {
    minX = +9999; maxX = -9999;
    minY = +9999; maxY = -9999;
    for (Int_t ix = 0; ix < fNPadsX[ich]; ix++) {
      pad = (AliMpPad*)fPadsX[ich]->At(ix);
      xloc = pad->GetPositionX();
      yloc = pad->GetPositionY();
      dx = pad->GetDimensionX();
      dy = pad->GetDimensionY();
      transformer.Local2Global((11+ich)*100+GetDetElemId(), xloc, yloc, 0, xglo, yglo, zglo);
      minX = TMath::Min(minX,(xglo-dx));
      maxX = TMath::Max(maxX,(xglo+dx));
      minY = TMath::Min(minY,(yglo-dy));
      maxY = TMath::Max(maxY,(yglo+dy));
    }
    fXCenter[ich] = 0.5*(minX+maxX);
    fYCenter[ich] = 0.5*(minY+maxY);
    fZCenter[ich] = zglo;
    fXWidth[ich]  = maxX-minX;
    fYWidth[ich]  = maxY-minY;
    // truncate to same precision as in the old guimap files
    fXCenter[ich] = 0.01*TMath::Nint(fXCenter[ich]*100.0);
    fYCenter[ich] = 0.01*TMath::Nint(fYCenter[ich]*100.0);
    fXWidth[ich] = 0.01*TMath::Nint(fXWidth[ich]*100.0);
    fYWidth[ich] = 0.01*TMath::Nint(fYWidth[ich]*100.0);

  }

  // delete the pads arrays
  for (Int_t ich = 0; ich < kNMT; ich++) {
    delete fPadsX[ich]; fNPadsX[ich] = 0;
    delete fPadsY[ich]; fNPadsY[ich] = 0;
  }
  
}

//__________________________________________________________________________
Int_t AliMUONTriggerGUIboard::GetLine() const
{
  /// get detector side
  if (fName->Length() >= 5) {
    const Char_t *name = fName->Data();
    TString sline = TString(name[4]);
    return sline.Atoi();
  }

  return -1;

}

//__________________________________________________________________________
Int_t AliMUONTriggerGUIboard::GetCol() const
{
  /// get detector side
  if (fName->Length() >= 5) {
    const Char_t *name = fName->Data();
    TString scol = TString(name[2]);
    return scol.Atoi();
  }

  return -1;

}

//__________________________________________________________________________
Int_t AliMUONTriggerGUIboard::GetSide() const
{
  /// get detector side
  if (fName->Length() >= 5) {
    const Char_t *name = fName->Data();
    if (!strncmp(name,"L",1)) return 0;
    if (!strncmp(name,"R",1)) return 1;
  }

  return -1;

}

//__________________________________________________________________________
void AliMUONTriggerGUIboard::PrintBoard() const
{
  /// print information on this board

  printf("Name: %s Id %3d Circ %3d DetElemId %2d Pos %1d YOver %1d\n",GetBoardName(),GetNumber(),GetIdCircuit(),GetDetElemId(),GetPosition(),GetYOver());
  printf("NStrips: X %2d Y %2d \n",GetNStripX(),GetNStripY());
  printf("Pad indices: X: ix %3d iy1 %3d iy2 %3d \n",GetXSix(),GetXSiy1(),GetXSiy2());
  printf("Pad indices: Y: iy %3d ix1 %3d ix2 %3d \n",GetYSiy(),GetYSix1(),GetYSix2());
  printf("Position and dimension:\n");
  for (Int_t imt = 0; imt < 4; imt++) {
    printf("MT=%1d: X %9.4f Y %9.4f Z %10.4f \n",imt,GetXCenter(imt),GetYCenter(imt),GetZCenter(imt));
    printf("      DX %7.4f DY %7.4f \n",GetXWidth(imt),GetYWidth(imt));
  }

}
Commit	Line	Data
905654c2	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
37391c63	16	// $Id$
37391c63	17
3d1463c8	18	//-----------------------------------------------------------------------------
37391c63	19	/// \class AliMUONTriggerGUIboard
	20	///
	21	/// Single trigger board object with geometry information, strips and digits
	22	///
	23	/// \author Bogdan Vulpescu, LPC Clermont-Ferrand
3d1463c8	24	//-----------------------------------------------------------------------------
905654c2	25
8b067dfe	26	#include <TClonesArray.h>
905654c2	27	#include <TBox.h>
8b067dfe	28	#include <TMath.h>
	29
	30	#include "AliMUONGeometryTransformer.h"
905654c2	31
	32	#include "AliMUONTriggerGUIboard.h"
	33
	34	/// \cond CLASSIMP
	35	ClassImp(AliMUONTriggerGUIboard)
	36	/// \endcond
	37
	38	//__________________________________________________________________________
8b067dfe	39	AliMUONTriggerGUIboard::AliMUONTriggerGUIboard()
905654c2	40	: TObject(),
905654c2	41	fName(0),
8b067dfe	42	fCrateName(0),
905654c2	43	fID(-1),
	44	fStatus(0),
	45	fPosition(0),
	46	fYOver(0),
	47	fXSix(0),
	48	fXSiy1(0),
	49	fXSiy2(0),
	50	fYSix1(0),
	51	fYSix2(0),
	52	fYSiy(0),
	53	fDetElemId(0),
	54	fIdCircuit(-1),
8b067dfe	55	fIsOpen(0),
	56	fNPadsX(),
	57	fNPadsY(),
	58	fPadsX(),
	59	fPadsY()
905654c2	60	{
	61	/// board main constructor
	62
8b067dfe	63	fName = new TString("");
8b067dfe	64	fCrateName = new TString("");
905654c2	65
	66	for (Int_t i = 0; i < kNMT; i++) {
	67	fXCenter[i] = 0.;
	68	fYCenter[i] = 0.;
	69	fZCenter[i] = 0.;
	70	fXWidth[i] = 0.;
	71	fYWidth[i] = 0.;
	72	for (Int_t is = 0; is < kNS; is++) {
	73	fXDig[i][is] = 0;
	74	fYDig[i][is] = 0;
	75	fXDigBox[i][is] = new TBox(0,0,0,0);
	76	fYDigBox[i][is] = new TBox(0,0,0,0);
	77	fXDigBox[i][is]->SetBit(kCannotPick);
	78	fYDigBox[i][is]->SetBit(kCannotPick);
	79	fXDigBox[i][is]->SetFillStyle(1001);
	80	fYDigBox[i][is]->SetFillStyle(1001);
	81	fXDigBox[i][is]->SetFillColor(4);
	82	fYDigBox[i][is]->SetFillColor(4);
	83	}
	84	}
	85
	86	fXSix = -1;
	87	fXSiy1 = -1;
	88	fXSiy2 = -1;
	89
	90	fYSix1 = -1;
	91	fYSix2 = -1;
	92	fYSiy = -1;
	93
	94	fDetElemId = -1;
	95	fIdCircuit = -1;
	96
	97	fIsOpen = kFALSE;
	98
	99	fYOver = 0;
	100	fPosition = 0;
	101
8b067dfe	102	for (Int_t i = 0; i < kNMT; i++) {
	103	fPadsX[i] = new TClonesArray("AliMpPad",16); fNPadsX[i] = 0;
	104	fPadsY[i] = new TClonesArray("AliMpPad",16); fNPadsY[i] = 0;
	105	}
	106
905654c2	107	}
	108
	109	//__________________________________________________________________________
	110	AliMUONTriggerGUIboard::~AliMUONTriggerGUIboard()
	111	{
	112	/// board destructor
	113
	114	delete fName;
	115
	116	for (Int_t imt = 0; imt < kNMT; imt++) {
	117	for (Int_t is = 0; is < kNS; is++) {
	118	delete fXDigBox[imt][is];
	119	delete fYDigBox[imt][is];
	120	}
	121	}
	122
	123	}
	124
905654c2	125	//__________________________________________________________________________
	126	void AliMUONTriggerGUIboard::SetXDigBox(Int_t imt, Int_t is, Double_t x1, Double_t y1, Double_t x2, Double_t y2)
	127	{
	128	/// set coordinates of "is" x-strip box in chamber "imt"
	129
	130	fXDigBox[imt][is]->SetX1(x1);
	131	fXDigBox[imt][is]->SetY1(y1);
	132	fXDigBox[imt][is]->SetX2(x2);
	133	fXDigBox[imt][is]->SetY2(y2);
	134
	135	}
	136
	137	//__________________________________________________________________________
	138	void AliMUONTriggerGUIboard::SetYDigBox(Int_t imt, Int_t is, Double_t x1, Double_t y1, Double_t x2, Double_t y2)
	139	{
	140	/// set coordinates of "is" y-strip box in chamber "imt"
	141
	142	fYDigBox[imt][is]->SetX1(x1);
	143	fYDigBox[imt][is]->SetY1(y1);
	144	fYDigBox[imt][is]->SetX2(x2);
	145	fYDigBox[imt][is]->SetY2(y2);
	146
	147	}
	148
	149	//__________________________________________________________________________
	150	void AliMUONTriggerGUIboard::ClearXDigits()
	151	{
	152	/// delete the set x-digits
	153
	154	for (Int_t imt = 0; imt < kNMT; imt++) {
	155	for (Int_t is = 0; is < kNS; is++) {
	156	fXDig[imt][is] = 0;
	157	}
	158	}
	159
	160	}
	161
	162	//__________________________________________________________________________
	163	void AliMUONTriggerGUIboard::ClearYDigits()
	164	{
	165	/// delete the set y-digits
	166
	167	for (Int_t imt = 0; imt < kNMT; imt++) {
	168	for (Int_t is = 0; is < kNS; is++) {
	169	fYDig[imt][is] = 0;
	170	}
	171	}
	172
	173	}
8b067dfe	174
	175	//__________________________________________________________________________
	176	void AliMUONTriggerGUIboard::MakeGeometry()
	177	{
	178	/// create the display geometry from the mapping pads
	179
	180	AliMpPad *pad;
	181
	182	// circuit number and manu channel (from x-strips)
	183	for (Int_t ich = 0; ich < kNMT; ich++) {
	184	if (fNPadsX[ich]) {
	185	pad = (AliMpPad*)fPadsX[ich]->At(0);
168e9c4d	186	fIdCircuit = pad->GetLocalBoardId(0);
8b067dfe	187	break;
	188	}
	189	}
	190
	191	// position index
	192	if (fName->Length()) {
	193	if (fName->Contains("12")) fPosition = 1;
	194	if (fName->Contains("34")) fPosition = 2;
	195	if (fName->Contains("56")) fPosition = 3;
	196	if (fName->Contains("78")) fPosition = 4;
	197	}
	198
	199	// position index for common y-strip boards
	200	for (Int_t ich = 0; ich < kNMT; ich++) {
	201	if (fNPadsY[ich]) {
2129750f	202	pad = (AliMpPad*)fPadsY[ich]->At(0);
8b067dfe	203	fYOver = pad->GetNofLocations();
	204	break;
	205	}
	206	}
	207
	208	// pad indices
	209	Int_t padxIx = -1, padxIy1 = +999, padxIy2 = -999;
	210	Int_t padyIy = -1, padyIx1 = +999, padyIx2 = -999;
	211	for (Int_t ip = 0; ip < fNPadsX[0]; ip++) {
	212	pad = (AliMpPad*)fPadsX[0]->At(ip);
168e9c4d	213	padxIx = pad->GetIx();
	214	padxIy1 = TMath::Min(padxIy1,pad->GetIy());
	215	padxIy2 = TMath::Max(padxIy2,pad->GetIy());
8b067dfe	216	}
	217	for (Int_t ip = 0; ip < fNPadsY[0]; ip++) {
	218	pad = (AliMpPad*)fPadsY[0]->At(ip);
168e9c4d	219	padyIy = pad->GetIy();
	220	padyIx1 = TMath::Min(padyIx1,pad->GetIx());
	221	padyIx2 = TMath::Max(padyIx2,pad->GetIx());
8b067dfe	222	}
	223	fXSix = padxIx;
	224	fXSiy1 = padxIy1;
	225	fXSiy2 = padxIy2;
	226	fYSiy = padyIy;
	227	fYSix1 = padyIx1;
	228	fYSix2 = padyIx2;
	229
	230	// position and dimension
	231
	232	AliMUONGeometryTransformer transformer;
	233	transformer.LoadGeometryData("transform.dat");
	234
	235	Float_t minX, maxX, minY, maxY;
	236	Float_t dx, dy;
	237	Float_t xloc, yloc, xglo, yglo, zglo;
	238	for (Int_t ich = 0; ich < kNMT; ich++) {
	239	minX = +9999; maxX = -9999;
	240	minY = +9999; maxY = -9999;
	241	for (Int_t ix = 0; ix < fNPadsX[ich]; ix++) {
	242	pad = (AliMpPad*)fPadsX[ich]->At(ix);
6e97fbb8	243	xloc = pad->GetPositionX();
	244	yloc = pad->GetPositionY();
	245	dx = pad->GetDimensionX();
	246	dy = pad->GetDimensionY();
8b067dfe	247	transformer.Local2Global((11+ich)*100+GetDetElemId(), xloc, yloc, 0, xglo, yglo, zglo);
	248	minX = TMath::Min(minX,(xglo-dx));
	249	maxX = TMath::Max(maxX,(xglo+dx));
	250	minY = TMath::Min(minY,(yglo-dy));
	251	maxY = TMath::Max(maxY,(yglo+dy));
	252	}
	253	fXCenter[ich] = 0.5*(minX+maxX);
	254	fYCenter[ich] = 0.5*(minY+maxY);
	255	fZCenter[ich] = zglo;
	256	fXWidth[ich] = maxX-minX;
	257	fYWidth[ich] = maxY-minY;
	258	// truncate to same precision as in the old guimap files
	259	fXCenter[ich] = 0.01TMath::Nint(fXCenter[ich]100.0);
	260	fYCenter[ich] = 0.01TMath::Nint(fYCenter[ich]100.0);
	261	fXWidth[ich] = 0.01TMath::Nint(fXWidth[ich]100.0);
	262	fYWidth[ich] = 0.01TMath::Nint(fYWidth[ich]100.0);
	263
	264	}
	265
	266	// delete the pads arrays
	267	for (Int_t ich = 0; ich < kNMT; ich++) {
	268	delete fPadsX[ich]; fNPadsX[ich] = 0;
	269	delete fPadsY[ich]; fNPadsY[ich] = 0;
	270	}
	271
	272	}
	273
	274	//__________________________________________________________________________
	275	Int_t AliMUONTriggerGUIboard::GetLine() const
	276	{
	277	/// get detector side
	278	if (fName->Length() >= 5) {
	279	const Char_t *name = fName->Data();
	280	TString sline = TString(name[4]);
	281	return sline.Atoi();
	282	}
	283
	284	return -1;
	285
	286	}
	287
	288	//__________________________________________________________________________
	289	Int_t AliMUONTriggerGUIboard::GetCol() const
	290	{
	291	/// get detector side
	292	if (fName->Length() >= 5) {
	293	const Char_t *name = fName->Data();
	294	TString scol = TString(name[2]);
	295	return scol.Atoi();
	296	}
	297
	298	return -1;
	299
	300	}
	301
	302	//__________________________________________________________________________
	303	Int_t AliMUONTriggerGUIboard::GetSide() const
	304	{
	305	/// get detector side
	306	if (fName->Length() >= 5) {
	307	const Char_t *name = fName->Data();
	308	if (!strncmp(name,"L",1)) return 0;
	309	if (!strncmp(name,"R",1)) return 1;
	310	}
311
312	return -1;
313
314	}
315
316	//__________________________________________________________________________
317	void AliMUONTriggerGUIboard::PrintBoard() const
318	{
319	/// print information on this board
320
321	printf("Name: %s Id %3d Circ %3d DetElemId %2d Pos %1d YOver %1d\n",GetBoardName(),GetNumber(),GetIdCircuit(),GetDetElemId(),GetPosition(),GetYOver());
322	printf("NStrips: X %2d Y %2d \n",GetNStripX(),GetNStripY());
323	printf("Pad indices: X: ix %3d iy1 %3d iy2 %3d \n",GetXSix(),GetXSiy1(),GetXSiy2());
324	printf("Pad indices: Y: iy %3d ix1 %3d ix2 %3d \n",GetYSiy(),GetYSix1(),GetYSix2());
325	printf("Position and dimension:\n");
326	for (Int_t imt = 0; imt < 4; imt++) {
327	printf("MT=%1d: X %9.4f Y %9.4f Z %10.4f \n",imt,GetXCenter(imt),GetYCenter(imt),GetZCenter(imt));
328	printf(" DX %7.4f DY %7.4f \n",GetXWidth(imt),GetYWidth(imt));
329	}
330
331	}
332