[u/mrichter/AliRoot.git] / TRD / AliTRDtransform.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$ */

////////////////////////////////////////////////////////////////////////////
//                                                                        //
//  Transforms clusters into space points with calibrated positions       //
//  defined in the local tracking system                                  //
//                                                                        //
////////////////////////////////////////////////////////////////////////////

#include <TGeoMatrix.h>

#include "AliLog.h"

#include "AliTRDtransform.h"
#include "AliTRDcluster.h"
#include "AliTRDgeometry.h"
#include "AliTRDpadPlane.h"
#include "AliTRDCommonParam.h"
#include "AliTRDcalibDB.h"
#include "Cal/AliTRDCalDet.h"
#include "Cal/AliTRDCalROC.h"

ClassImp(AliTRDtransform)

AliTRDgeometry* AliTRDtransform::fgGeo = NULL;

//_____________________________________________________________________________
AliTRDtransform::AliTRDtransform()
  :TObject()
  ,fDetector(0)
  ,fParam(0x0)
  ,fCalibration(0x0)
  ,fCalVdriftROC(0x0)
  ,fCalT0ROC(0x0)
  ,fCalPRFROC(0x0)
  ,fkCalVdriftDet(0x0)
  ,fkCalT0Det(0x0)
  ,fCalVdriftDetValue(0)
  ,fCalT0DetValue(0)
  ,fSamplingFrequency(0)
  ,fPadPlane(0x0)
  ,fZShiftIdeal(0)
  ,fMatrix(0x0)
{
  //
  // AliTRDtransform default constructor
  //

  if(!fgGeo){
    fgGeo = new AliTRDgeometry();
    if (!fgGeo->CreateClusterMatrixArray()) {
      AliError("Could not get transformation matrices\n");
    }
  }

  fParam             = AliTRDCommonParam::Instance();
  if (!fParam) {
    AliError("Could not get common parameters\n");
  }
  fSamplingFrequency = fParam->GetSamplingFrequency();

  fCalibration       = AliTRDcalibDB::Instance();
  if (!fCalibration) {
    AliError("Cannot find calibration object");
  }

  // Get the calibration objects for the global calibration
  fkCalVdriftDet     = fCalibration->GetVdriftDet();
  fkCalT0Det         = fCalibration->GetT0Det();

}

//_____________________________________________________________________________
AliTRDtransform::AliTRDtransform(Int_t det)
  :TObject()
  ,fDetector(0)
  ,fParam(0x0)
  ,fCalibration(0x0)
  ,fCalVdriftROC(0x0)
  ,fCalT0ROC(0x0)
  ,fCalPRFROC(0x0)
  ,fkCalVdriftDet(0x0)
  ,fkCalT0Det(0x0)
  ,fCalVdriftDetValue(0)
  ,fCalT0DetValue(0)
  ,fSamplingFrequency(0)
  ,fPadPlane(0x0)
  ,fZShiftIdeal(0)
  ,fMatrix(0x0)
{
  //
  // AliTRDtransform constructor for a given detector
  //

  if(!fgGeo){
    fgGeo = new AliTRDgeometry();
    if (!fgGeo->CreateClusterMatrixArray()) {
      AliError("Could not get transformation matrices\n");
    }
  }

  fParam             = AliTRDCommonParam::Instance();
  if (!fParam) {
    AliError("Could not get common parameters\n");
  }
  fSamplingFrequency = fParam->GetSamplingFrequency();

  fCalibration       = AliTRDcalibDB::Instance();
  if (!fCalibration) {
    AliError("Cannot find calibration object");
  }

  // Get the calibration objects for the global calibration
  fkCalVdriftDet     = fCalibration->GetVdriftDet();
  fkCalT0Det         = fCalibration->GetT0Det();

  SetDetector(det);

}

//_____________________________________________________________________________
AliTRDtransform::AliTRDtransform(const AliTRDtransform &t)
  :TObject(t)
  ,fDetector(t.fDetector)
  ,fParam(0x0)
  ,fCalibration(0x0)
  ,fCalVdriftROC(0x0)
  ,fCalT0ROC(0x0)
  ,fCalPRFROC(0x0)
  ,fkCalVdriftDet(0x0)
  ,fkCalT0Det(0x0)
  ,fCalVdriftDetValue(0)
  ,fCalT0DetValue(0)
  ,fSamplingFrequency(0)
  ,fPadPlane(0x0)
  ,fZShiftIdeal(0)
  ,fMatrix(0x0)
{
  //
  // AliTRDtransform copy constructor
  //

  fParam             = AliTRDCommonParam::Instance();
  if (!fParam) {
    AliError("Could not get common parameters\n");
  }
  fSamplingFrequency = fParam->GetSamplingFrequency();

  fCalibration = AliTRDcalibDB::Instance();
  if (!fCalibration) {
    AliError("Cannot find calibration object");
  }
  fkCalVdriftDet     = fCalibration->GetVdriftDet();
  fkCalT0Det         = fCalibration->GetT0Det();

}

//_____________________________________________________________________________
AliTRDtransform::~AliTRDtransform()
{
  //
  // AliTRDtransform destructor
  //

}

//_____________________________________________________________________________
AliTRDtransform &AliTRDtransform::operator=(const AliTRDtransform &t)
{
  //
  // Assignment operator
  //

  if (this != &t) {
    ((AliTRDtransform &) t).Copy(*this);
  }

  return *this;

}

//_____________________________________________________________________________
void AliTRDtransform::Copy(TObject &t) const
{
  //
  // Copy function
  //

  ((AliTRDtransform &) t).fDetector          = fDetector;
  ((AliTRDtransform &) t).fParam             = AliTRDCommonParam::Instance();
  ((AliTRDtransform &) t).fCalibration       = AliTRDcalibDB::Instance();
  if (fCalibration) {
    ((AliTRDtransform &) t).fkCalVdriftDet   = fCalibration->GetVdriftDet();
    ((AliTRDtransform &) t).fkCalT0Det       = fCalibration->GetT0Det();
  }
  else {
    ((AliTRDtransform &) t).fkCalVdriftDet   = 0;
    ((AliTRDtransform &) t).fkCalT0Det       = 0;
  }
  ((AliTRDtransform &) t).fCalVdriftROC      = 0x0;
  ((AliTRDtransform &) t).fCalT0ROC          = 0x0;
  ((AliTRDtransform &) t).fCalPRFROC         = 0x0;
  ((AliTRDtransform &) t).fCalVdriftDetValue = 0; 
  ((AliTRDtransform &) t).fCalT0DetValue     = 0;
  ((AliTRDtransform &) t).fSamplingFrequency = 0;
  ((AliTRDtransform &) t).fPadPlane          = 0x0;
  ((AliTRDtransform &) t).fZShiftIdeal       = 0;
  ((AliTRDtransform &) t).fMatrix            = 0x0;

}

//_____________________________________________________________________________
void AliTRDtransform::SetDetector(Int_t det)
{
  //
  // Set to a new detector number and update the calibration objects
  // and values accordingly
  //

  fDetector          = det;

  // Get the calibration objects for the pad-by-pad calibration
  fCalVdriftROC      = fCalibration->GetVdriftROC(det);
  fCalT0ROC          = fCalibration->GetT0ROC(det);
  fCalPRFROC         = fCalibration->GetPRFROC(det);

  // Get the detector wise defined calibration values
  fCalVdriftDetValue = fkCalVdriftDet->GetValue(det);
  fCalT0DetValue     = fkCalT0Det->GetValue(det);

  // Shift needed to define Z-position relative to middle of chamber
  Int_t layer        = fgGeo->GetLayer(det);
  Int_t stack        = fgGeo->GetStack(det);
  fPadPlane          = fgGeo->GetPadPlane(layer,stack);
  fZShiftIdeal       = 0.5 * (fPadPlane->GetRow0() + fPadPlane->GetRowEnd());

  // Get the current transformation matrix
  fMatrix            = fgGeo->GetClusterMatrix(det);

}

//_____________________________________________________________________________
Bool_t AliTRDtransform::Transform(AliTRDcluster *c)
{
  //
  // Transforms the local cluster coordinates into calibrated 
  // space point positions defined in the local tracking system.
  //
  // Here the calibration for T0, Vdrift and ExB is applied as well.
  //
  // Input: Cluster in the local chamber coordinates
  // Output: Tracking cluster

  if (!fMatrix) return kFALSE;


  // Parameters to adjust the X position of clusters in the alignable volume
  const Double_t kX0shift = AliTRDgeometry::AnodePos(); //[cm]

 
  // Retrieve calibration values
  Int_t col = c->GetPadCol(), row = c->GetPadRow();
  // drift velocity
  Double_t vd  = fCalVdriftDetValue * fCalVdriftROC->GetValue(col,row);
  // t0
  Double_t t0  = fCalT0DetValue     + fCalT0ROC->GetValue(col,row);
  t0 /= fSamplingFrequency;
  // ExB correction
  Double_t exb = AliTRDCommonParam::Instance()->GetOmegaTau(vd);

  Float_t x = c->GetXloc(t0, vd);

  // Pad dimensions
  Double_t rs = fPadPlane->GetRowSize(row);
  Double_t cs = fPadPlane->GetColSize(col);

  // cluster error with diffusion corrections
  Double_t s2  = cs*fCalPRFROC->GetValue(col, row); 
  s2 *= s2; 
  Float_t dl, dt;
  AliTRDCommonParam::Instance()->GetDiffCoeff(dl, dt, vd);

  Double_t y0 = fPadPlane->GetColPos(col) + .5*cs;
  Double_t loc[] = {
    kX0shift-x,                    // Invert the X-position,
    c->GetYloc(y0, s2, cs) - x*exb,// apply ExB correction
    fPadPlane->GetRowPos(row) - .5*rs - fZShiftIdeal // move the Z-position relative to the middle of the chamber
  };

  // Go to tracking coordinates
  Double_t trk[3];
  fMatrix->LocalToMaster(loc, trk);

  // store tracking values
  c->SetX(trk[0]);c->SetY(trk[1]);c->SetZ(trk[2]);
  c->SetSigmaY2(s2, dt, exb, x);
  c->SetSigmaZ2(fPadPlane->GetRowSize(row)*fPadPlane->GetRowSize(row)/12.);
  
  return kTRUE;

}

//_____________________________________________________________________________
void AliTRDtransform::Recalibrate(AliTRDcluster *c, Bool_t setDet)
{
  //
  // Recalibrates the position of a given cluster
  // If <setDet> is TRUE, the detector number is set for each cluster
  // automatically. Otherwise, AliTRDtransform::SetDetector() has to
  // be used.
  //

  if (setDet) SetDetector(c->GetDetector());
  Transform(c);

}
Commit	Line	Data
af26ce80	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	/* $Id$ */
	17
	18	////////////////////////////////////////////////////////////////////////////
	19	// //
	20	// Transforms clusters into space points with calibrated positions //
	21	// defined in the local tracking system //
	22	// //
	23	////////////////////////////////////////////////////////////////////////////
	24
	25	#include <TGeoMatrix.h>
	26
	27	#include "AliLog.h"
af26ce80	28
	29	#include "AliTRDtransform.h"
	30	#include "AliTRDcluster.h"
	31	#include "AliTRDgeometry.h"
	32	#include "AliTRDpadPlane.h"
	33	#include "AliTRDCommonParam.h"
	34	#include "AliTRDcalibDB.h"
	35	#include "Cal/AliTRDCalDet.h"
	36	#include "Cal/AliTRDCalROC.h"
	37
	38	ClassImp(AliTRDtransform)
	39
8c499dbf	40	AliTRDgeometry* AliTRDtransform::fgGeo = NULL;
8c499dbf	41
af26ce80	42	//_____________________________________________________________________________
af26ce80	43	AliTRDtransform::AliTRDtransform()
af26ce80	44	:TObject()
af26ce80	45	,fDetector(0)
	46	,fParam(0x0)
	47	,fCalibration(0x0)
	48	,fCalVdriftROC(0x0)
	49	,fCalT0ROC(0x0)
a60b4183	50	,fCalPRFROC(0x0)
c893147d	51	,fkCalVdriftDet(0x0)
c893147d	52	,fkCalT0Det(0x0)
af26ce80	53	,fCalVdriftDetValue(0)
	54	,fCalT0DetValue(0)
	55	,fSamplingFrequency(0)
	56	,fPadPlane(0x0)
	57	,fZShiftIdeal(0)
	58	,fMatrix(0x0)
	59	{
	60	//
	61	// AliTRDtransform default constructor
	62	//
	63
8c499dbf	64	if(!fgGeo){
	65	fgGeo = new AliTRDgeometry();
	66	if (!fgGeo->CreateClusterMatrixArray()) {
	67	AliError("Could not get transformation matrices\n");
	68	}
225096c8	69	}
	70
	71	fParam = AliTRDCommonParam::Instance();
	72	if (!fParam) {
	73	AliError("Could not get common parameters\n");
	74	}
	75	fSamplingFrequency = fParam->GetSamplingFrequency();
	76
	77	fCalibration = AliTRDcalibDB::Instance();
	78	if (!fCalibration) {
	79	AliError("Cannot find calibration object");
	80	}
	81
	82	// Get the calibration objects for the global calibration
	83	fkCalVdriftDet = fCalibration->GetVdriftDet();
	84	fkCalT0Det = fCalibration->GetT0Det();
	85
af26ce80	86	}
	87
	88	//_____________________________________________________________________________
af26ce80	89	AliTRDtransform::AliTRDtransform(Int_t det)
af26ce80	90	:TObject()
af26ce80	91	,fDetector(0)
	92	,fParam(0x0)
	93	,fCalibration(0x0)
	94	,fCalVdriftROC(0x0)
	95	,fCalT0ROC(0x0)
a60b4183	96	,fCalPRFROC(0x0)
c893147d	97	,fkCalVdriftDet(0x0)
c893147d	98	,fkCalT0Det(0x0)
af26ce80	99	,fCalVdriftDetValue(0)
	100	,fCalT0DetValue(0)
	101	,fSamplingFrequency(0)
	102	,fPadPlane(0x0)
	103	,fZShiftIdeal(0)
	104	,fMatrix(0x0)
	105	{
	106	//
	107	// AliTRDtransform constructor for a given detector
	108	//
	109
8c499dbf	110	if(!fgGeo){
	111	fgGeo = new AliTRDgeometry();
	112	if (!fgGeo->CreateClusterMatrixArray()) {
	113	AliError("Could not get transformation matrices\n");
	114	}
a938cf15	115	}
af26ce80	116
	117	fParam = AliTRDCommonParam::Instance();
	118	if (!fParam) {
	119	AliError("Could not get common parameters\n");
	120	}
	121	fSamplingFrequency = fParam->GetSamplingFrequency();
	122
	123	fCalibration = AliTRDcalibDB::Instance();
	124	if (!fCalibration) {
	125	AliError("Cannot find calibration object");
	126	}
	127
	128	// Get the calibration objects for the global calibration
c893147d	129	fkCalVdriftDet = fCalibration->GetVdriftDet();
c893147d	130	fkCalT0Det = fCalibration->GetT0Det();
af26ce80	131
	132	SetDetector(det);
	133
	134	}
	135
	136	//_____________________________________________________________________________
af26ce80	137	AliTRDtransform::AliTRDtransform(const AliTRDtransform &t)
af26ce80	138	:TObject(t)
af26ce80	139	,fDetector(t.fDetector)
	140	,fParam(0x0)
	141	,fCalibration(0x0)
	142	,fCalVdriftROC(0x0)
	143	,fCalT0ROC(0x0)
a60b4183	144	,fCalPRFROC(0x0)
c893147d	145	,fkCalVdriftDet(0x0)
c893147d	146	,fkCalT0Det(0x0)
af26ce80	147	,fCalVdriftDetValue(0)
	148	,fCalT0DetValue(0)
	149	,fSamplingFrequency(0)
	150	,fPadPlane(0x0)
	151	,fZShiftIdeal(0)
	152	,fMatrix(0x0)
	153	{
	154	//
	155	// AliTRDtransform copy constructor
	156	//
	157
af26ce80	158	fParam = AliTRDCommonParam::Instance();
	159	if (!fParam) {
	160	AliError("Could not get common parameters\n");
	161	}
	162	fSamplingFrequency = fParam->GetSamplingFrequency();
	163
	164	fCalibration = AliTRDcalibDB::Instance();
	165	if (!fCalibration) {
	166	AliError("Cannot find calibration object");
	167	}
c893147d	168	fkCalVdriftDet = fCalibration->GetVdriftDet();
c893147d	169	fkCalT0Det = fCalibration->GetT0Det();
af26ce80	170
	171	}
	172
	173	//_____________________________________________________________________________
	174	AliTRDtransform::~AliTRDtransform()
	175	{
	176	//
	177	// AliTRDtransform destructor
	178	//
	179
af26ce80	180	}
af26ce80	181
fc1c9d79	182	//_____________________________________________________________________________
	183	AliTRDtransform &AliTRDtransform::operator=(const AliTRDtransform &t)
	184	{
	185	//
	186	// Assignment operator
	187	//
	188
	189	if (this != &t) {
	190	((AliTRDtransform &) t).Copy(*this);
	191	}
	192
	193	return *this;
	194
	195	}
	196
	197	//_____________________________________________________________________________
	198	void AliTRDtransform::Copy(TObject &t) const
	199	{
	200	//
	201	// Copy function
	202	//
	203
	204	((AliTRDtransform &) t).fDetector = fDetector;
	205	((AliTRDtransform &) t).fParam = AliTRDCommonParam::Instance();
	206	((AliTRDtransform &) t).fCalibration = AliTRDcalibDB::Instance();
	207	if (fCalibration) {
	208	((AliTRDtransform &) t).fkCalVdriftDet = fCalibration->GetVdriftDet();
	209	((AliTRDtransform &) t).fkCalT0Det = fCalibration->GetT0Det();
	210	}
	211	else {
	212	((AliTRDtransform &) t).fkCalVdriftDet = 0;
	213	((AliTRDtransform &) t).fkCalT0Det = 0;
	214	}
	215	((AliTRDtransform &) t).fCalVdriftROC = 0x0;
	216	((AliTRDtransform &) t).fCalT0ROC = 0x0;
	217	((AliTRDtransform &) t).fCalPRFROC = 0x0;
	218	((AliTRDtransform &) t).fCalVdriftDetValue = 0;
	219	((AliTRDtransform &) t).fCalT0DetValue = 0;
	220	((AliTRDtransform &) t).fSamplingFrequency = 0;
	221	((AliTRDtransform &) t).fPadPlane = 0x0;
	222	((AliTRDtransform &) t).fZShiftIdeal = 0;
	223	((AliTRDtransform &) t).fMatrix = 0x0;
	224
	225	}
	226
af26ce80	227	//_____________________________________________________________________________
	228	void AliTRDtransform::SetDetector(Int_t det)
	229	{
	230	//
	231	// Set to a new detector number and update the calibration objects
	232	// and values accordingly
	233	//
	234
	235	fDetector = det;
	236
	237	// Get the calibration objects for the pad-by-pad calibration
	238	fCalVdriftROC = fCalibration->GetVdriftROC(det);
	239	fCalT0ROC = fCalibration->GetT0ROC(det);
a60b4183	240	fCalPRFROC = fCalibration->GetPRFROC(det);
af26ce80	241
af26ce80	242	// Get the detector wise defined calibration values
c893147d	243	fCalVdriftDetValue = fkCalVdriftDet->GetValue(det);
c893147d	244	fCalT0DetValue = fkCalT0Det->GetValue(det);
af26ce80	245
af26ce80	246	// Shift needed to define Z-position relative to middle of chamber
8c499dbf	247	Int_t layer = fgGeo->GetLayer(det);
	248	Int_t stack = fgGeo->GetStack(det);
	249	fPadPlane = fgGeo->GetPadPlane(layer,stack);
af26ce80	250	fZShiftIdeal = 0.5 * (fPadPlane->GetRow0() + fPadPlane->GetRowEnd());
	251
	252	// Get the current transformation matrix
8c499dbf	253	fMatrix = fgGeo->GetClusterMatrix(det);
af26ce80	254
	255	}
	256
	257	//_____________________________________________________________________________
b72f4eaf	258	Bool_t AliTRDtransform::Transform(AliTRDcluster *c)
af26ce80	259	{
	260	//
	261	// Transforms the local cluster coordinates into calibrated
	262	// space point positions defined in the local tracking system.
	263	//
	264	// Here the calibration for T0, Vdrift and ExB is applied as well.
	265	//
b72f4eaf	266	// Input: Cluster in the local chamber coordinates
b72f4eaf	267	// Output: Tracking cluster
af26ce80	268
b72f4eaf	269	if (!fMatrix) return kFALSE;
af26ce80	270
af26ce80	271
b72f4eaf	272	// Parameters to adjust the X position of clusters in the alignable volume
ec3f0161	273	const Double_t kX0shift = AliTRDgeometry::AnodePos(); //[cm]
1814b374	274
1814b374	275
b72f4eaf	276	// Retrieve calibration values
	277	Int_t col = c->GetPadCol(), row = c->GetPadRow();
	278	// drift velocity
	279	Double_t vd = fCalVdriftDetValue * fCalVdriftROC->GetValue(col,row);
	280	// t0
	281	Double_t t0 = fCalT0DetValue + fCalT0ROC->GetValue(col,row);
566bf887	282	t0 /= fSamplingFrequency;
b72f4eaf	283	// ExB correction
	284	Double_t exb = AliTRDCommonParam::Instance()->GetOmegaTau(vd);
	285
	286	Float_t x = c->GetXloc(t0, vd);
	287
b72f4eaf	288	// Pad dimensions
	289	Double_t rs = fPadPlane->GetRowSize(row);
	290	Double_t cs = fPadPlane->GetColSize(col);
1fd9389f	291
	292	// cluster error with diffusion corrections
	293	Double_t s2 = cs*fCalPRFROC->GetValue(col, row);
	294	s2 *= s2;
	295	Float_t dl, dt;
	296	AliTRDCommonParam::Instance()->GetDiffCoeff(dl, dt, vd);
	297
b72f4eaf	298	Double_t y0 = fPadPlane->GetColPos(col) + .5*cs;
	299	Double_t loc[] = {
	300	kX0shift-x, // Invert the X-position,
	301	c->GetYloc(y0, s2, cs) - x*exb,// apply ExB correction
	302	fPadPlane->GetRowPos(row) - .5*rs - fZShiftIdeal // move the Z-position relative to the middle of the chamber
	303	};
	304
	305	// Go to tracking coordinates
	306	Double_t trk[3];
	307	fMatrix->LocalToMaster(loc, trk);
	308
	309	// store tracking values
	310	c->SetX(trk[0]);c->SetY(trk[1]);c->SetZ(trk[2]);
1fd9389f	311	c->SetSigmaY2(s2, dt, exb, x);
b72f4eaf	312	c->SetSigmaZ2(fPadPlane->GetRowSize(row)*fPadPlane->GetRowSize(row)/12.);
	313
	314	return kTRUE;
c893147d	315
af26ce80	316	}
	317
	318	//_____________________________________________________________________________
	319	void AliTRDtransform::Recalibrate(AliTRDcluster *c, Bool_t setDet)
	320	{
	321	//
	322	// Recalibrates the position of a given cluster
	323	// If <setDet> is TRUE, the detector number is set for each cluster
	324	// automatically. Otherwise, AliTRDtransform::SetDetector() has to
	325	// be used.
	326	//
	327
b72f4eaf	328	if (setDet) SetDetector(c->GetDetector());
b72f4eaf	329	Transform(c);
c893147d	330
af26ce80	331	}