[u/mrichter/AliRoot.git] / PMD / AliPMDUtility.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.                  *
 **************************************************************************/
//-----------------------------------------------------//
//                                                     //
//                                                     //
//  Date   : August 05 2003                            //
//                                                     //
//  Utility code for ALICE-PMD                         //
//                                                     //
//-----------------------------------------------------//

#include "Riostream.h"
#include "AliPMDUtility.h"
#include "TMath.h"
#include <stdio.h>
#include <math.h>


ClassImp(AliPMDUtility)

AliPMDUtility::AliPMDUtility()
{
  // Default constructor
  fPx    = 0.;
  fPy    = 0.;
  fPz    = 0.;
  fTheta = 0.;
  fEta   = 0.;
  fPhi   = 0.;
}

AliPMDUtility::AliPMDUtility(Float_t px, Float_t py, Float_t pz)
{
  // Constructor
  fPx = px;
  fPy = py;
  fPz = pz;
  fTheta = 0.;
  fEta   = 0.;
  fPhi   = 0.;
}

AliPMDUtility::~AliPMDUtility()
{
  // Default destructor
}
void AliPMDUtility::HexGeomCellPos(Int_t ism, Int_t xpad, Int_t ypad, Float_t &xpos, Float_t &ypos)
{
  // This converts PMD cluster or CELL coordinates
  // to Global coordinates.
  // Written by Prof. S.C. Phatak

  const Float_t  kCellDia = 0.5;
  const Float_t  kPi      = TMath::Pi();    //3.14159;
  const Double_t kSqroot3by2   = 0.8660254;  // sqrth = sqrt(3.)/2.

  Int_t i;
  Int_t j = xpad;
  Int_t k = ypad;

  /*
    Supermodeule number starting from 0
    ism --> supermodule no ( 0 - 26 )
    idet --> detector ( pmd or cpv : not required now )
    j --> xpad ( goes from 1 to 72 )
    k --> ypad ( goes from 1 to 72 )
    xp --> global x coordinate
    yp --> global y coordinate
    
    (xp0,yp0) corner positions of all supermodules in global
    coordinate system. That is the origin
    of the local ( supermodule ) coordinate system.
  */ 
  
  Float_t xp0[27] = 
  {
    -17.9084, 18.2166, 54.3416, -35.9709, 0.154144, 
    36.2791, -54.0334, -17.9084, 18.2166, 36.7791, 
    18.7166, 0.654194, 72.9041, 54.8416, 36.7792, 
    109.029, 90.9666, 72.9042, -18.8708, -36.9334, 
    -54.996, -36.9332, -54.9958, -73.0584, -54.9956, 
    -73.0582, -91.1208
  };

  Float_t yp0[27] = 
  {
    -32.1395, -32.1395, -32.1395, -63.4247, -63.4247, 
    -63.4247, -94.7098, -94.7098, -94.7098, 0.545689, 
    31.8309, 63.1161, 0.545632, 31.8308, 63.116, 
    0.545573, 31.8308, 63.116, 31.5737, 0.288616, 
    -30.9965, 62.859, 31.5738, 0.288733, 94.1442, 
    62.8591, 31.574
  };

  /* 
     angles of rotation for three sets of supermodules
     The angle is same for first nine, next nine and last nine 
     supermodules 
  */
  
  Float_t th[3] = {0., -2.*kPi/3., 2.*kPi/3.};
  Float_t xr, yr, xinit, yinit, cs, sn;
  
  /* 
     xinit and yinit are coordinates of the cell in local coordinate system
  */
  
  xinit = (j)*kCellDia+(k)/2.*kCellDia;
  yinit = kSqroot3by2*(k)/2.;
  i=ism/9;
  cs=cos(th[i]);
  sn=sin(th[i]);
  //
  // rotate first
  //
  xr=cs*xinit+sn*yinit;
  yr=-sn*xinit+cs*yinit;
  //
  // then translate
  //
  xpos=xr+xp0[ism];
  ypos=yr+yp0[ism];

}

void AliPMDUtility::RectGeomCellPos(Int_t ism, Int_t ium, Int_t xpad, Int_t ypad, Float_t &xpos, Float_t &ypos)
{
  // This routine finds the cell eta,phi for the new PMD rectangular 
  // geometry in ALICE
  // Authors : Bedanga Mohanty and Dipak Mishra - 29.4.2003
  // modified by B. K. Nnadi for change of coordinate sys
  //
  // SMA  ---> Supermodule Type A           ( SM - 0)
  // SMAR ---> Supermodule Type A ROTATED   ( SM - 1)
  // SMB  ---> Supermodule Type B           ( SM - 2)
  // SMBR ---> Supermodule Type B ROTATED   ( SM - 3)
  //
  // ism   : number of supermodules in one plane = 4
  // ium   : number of unitmodules  in one SM    = 6
  // gbum  : (global) unit module numbering in a supermodule
  //

  Int_t gbum = ism*6 + ium;
  Int_t irow  = xpad;
  Int_t icol  = ypad;

  // Corner positions (x,y) of the 24 unit moudles in ALICE PMD
  
  Double_t xcorner[24] =
  {
    85.15,  60.85,  36.55,  85.15,  60.85,  36.55, //SMA 
    -85.15, -60.85, -36.55, -85.15, -60.85, -36.55, //SMAR
    84.90,  36.60,  84.90,  36.60,  84.90,  36.60, //SMB
    -84.90, -36.60, -84.90, -36.60, -84.90, -36.60  //SMBR
  };
  
  Double_t ycorner[24] =
  { 
    32.45708755,  32.45708755,  32.45708755,        //SMA
    -9.30645245,  -9.30645245,  -9.30645245,        //SMA
    -32.45708755, -32.45708755, -32.45708755,        //SMAR
    9.30645245,   9.30645245,   9.30645245,        //SMAR
    -31.63540818, -31.63540818, -52.61435544,        //SMB
    -52.61435544, -73.59330270, -73.59330270,        //SMB
    31.63540818,  31.63540818,  52.61435544,        //SMBR
    52.61435544,  73.59330270,  73.59330270         //SMBR
  };
  
  const Float_t kSqroot3      = 1.73205;  // sqrt(3.);
  const Float_t kCellRadius   = 0.25;
  
  //
  //Every even row of cells is shifted and placed
  //in geant so this condition
  //
  Float_t shift = 0.0;
  if(irow%2 == 0)
    {
      shift = 0.25;
    }
  else
    {
      shift = 0.0;
    }
  if(ism == 0 || ism == 2)
    {
      ypos = ycorner[gbum] + 
	irow*kCellRadius*kSqroot3;

      xpos = xcorner[gbum] - 
	icol*2.0*kCellRadius - shift;
    }
  else if(ism == 1 || ism == 3)
    {
      ypos = ycorner[gbum] -
	irow*kCellRadius*kSqroot3;

      xpos = xcorner[gbum] +
	icol*2.0*kCellRadius + shift;
    }
}

void AliPMDUtility::RectGeomCellPos(Int_t ism, Int_t ium, Float_t xpad, Float_t ypad, Float_t &xpos, Float_t &ypos)
{
  // If the xpad and ypad inputs are float, then 0.5 is added to it
  // to find the layer which is shifted.
  // This routine finds the cell eta,phi for the new PMD rectangular 
  // geometry in ALICE
  // Authors : Bedanga Mohanty and Dipak Mishra - 29.4.2003
  // modified by B. K. Nnadi for change of coordinate sys
  //
  // SMA  ---> Supermodule Type A           ( SM - 0)
  // SMAR ---> Supermodule Type A ROTATED   ( SM - 1)
  // SMB  ---> Supermodule Type B           ( SM - 2)
  // SMBR ---> Supermodule Type B ROTATED   ( SM - 3)
  //
  // ism   : number of supermodules in one plane = 4
  // ium   : number of unitmodules  in one SM    = 6
  // gbum  : (global) unit module numbering in a supermodule
  //

  Int_t gbum    = ism*6 + ium;
  Float_t irow  = xpad;
  Float_t icol  = ypad;

  // Corner positions (x,y) of the 24 unit moudles in ALICE PMD
  
  Double_t xcorner[24] =
  {
    85.15,  60.85,  36.55,  85.15,  60.85,  36.55, //SMA 
    -85.15, -60.85, -36.55, -85.15, -60.85, -36.55, //SMAR
    84.90,  36.60,  84.90,  36.60,  84.90,  36.60, //SMB
    -84.90, -36.60, -84.90, -36.60, -84.90, -36.60  //SMBR
  };
  
  Double_t ycorner[24] =
  { 
    32.45708755,  32.45708755,  32.45708755,        //SMA
    -9.30645245,  -9.30645245,  -9.30645245,        //SMA
    -32.45708755, -32.45708755, -32.45708755,        //SMAR
    9.30645245,   9.30645245,   9.30645245,        //SMAR
    -31.63540818, -31.63540818, -52.61435544,        //SMB
    -52.61435544, -73.59330270, -73.59330270,        //SMB
    31.63540818,  31.63540818,  52.61435544,        //SMBR
    52.61435544,  73.59330270,  73.59330270         //SMBR
  };
  
  const Float_t kSqroot3    = 1.73205;  // sqrt(3.);
  const Float_t kCellRadius = 0.25;
  
  //
  //Every even row of cells is shifted and placed
  //in geant so this condition
  //
  Float_t shift = 0.0;
  Int_t iirow = (Int_t) (irow+0.5);
  if(iirow%2 == 0)
    {
      shift = 0.25;
    }
  else
    {
      shift = 0.0;
    }
  if(ism == 0 || ism == 2)
    {
      ypos = ycorner[gbum] + 
	irow*kCellRadius*kSqroot3;

      xpos = xcorner[gbum] - 
	icol*2.0*kCellRadius - shift;
    }
  else if(ism == 1 || ism == 3)
    {
      ypos = ycorner[gbum] -
	irow*kCellRadius*kSqroot3;

      xpos = xcorner[gbum] +
	icol*2.0*kCellRadius + shift;
    }
}

void AliPMDUtility::SetPxPyPz(Float_t px, Float_t py, Float_t pz)
{
  fPx = px;
  fPy = py;
  fPz = pz;
}

void AliPMDUtility::SetXYZ(Float_t xpos, Float_t ypos, Float_t zpos)
{
  fPx = xpos;
  fPy = ypos;
  fPz = zpos;
}
void AliPMDUtility::CalculateEta()
{
  Float_t rpxpy, theta, eta;

  rpxpy  = TMath::Sqrt(fPx*fPx + fPy*fPy);
  theta  = TMath::ATan2(rpxpy,fPz);
  eta    = -TMath::Log(TMath::Tan(0.5*theta));
  fTheta = theta;
  fEta   = eta;
}
void AliPMDUtility::CalculatePhi()
{
  Float_t pybypx, phi = 0., phi1;

  if(fPx==0)
    {
      if(fPy>0) phi = 90.;
      if(fPy<0) phi = 270.;
    }
  if(fPx != 0)
    {
      pybypx = fPy/fPx;
      if(pybypx < 0) pybypx = - pybypx;
      phi1 = TMath::ATan(pybypx)*180./3.14159;

      if(fPx > 0 && fPy > 0) phi = phi1;        // 1st Quadrant
      if(fPx < 0 && fPy > 0) phi = 180 - phi1;  // 2nd Quadrant
      if(fPx < 0 && fPy < 0) phi = 180 + phi1;  // 3rd Quadrant
      if(fPx > 0 && fPy < 0) phi = 360 - phi1;  // 4th Quadrant

    }
  phi = phi*3.14159/180.;

  fPhi = phi;

}
void AliPMDUtility::CalculateEtaPhi()
{
  Float_t rpxpy, theta, eta;
  Float_t pybypx, phi = 0., phi1;

  rpxpy = TMath::Sqrt(fPx*fPx + fPy*fPy);
  theta = TMath::ATan2(rpxpy,fPz);
  eta   = -TMath::Log(TMath::Tan(0.5*theta));
  
  if(fPx==0)
    {
      if(fPy>0) phi = 90.;
      if(fPy<0) phi = 270.;
    }
  if(fPx != 0)
    {
      pybypx = fPy/fPx;
      if(pybypx < 0) pybypx = - pybypx;
      phi1 = TMath::ATan(pybypx)*180./3.14159;
      if(fPx > 0 && fPy > 0) phi = phi1;        // 1st Quadrant
      if(fPx < 0 && fPy > 0) phi = 180 - phi1;  // 2nd Quadrant
      if(fPx < 0 && fPy < 0) phi = 180 + phi1;  // 3rd Quadrant
      if(fPx > 0 && fPy < 0) phi = 360 - phi1;  // 4th Quadrant

    }
  phi = phi*3.14159/180.;

  fTheta = theta;
  fEta   = eta;
  fPhi   = phi;
}
Float_t AliPMDUtility::GetTheta() const
{
  return fTheta;
}
Float_t AliPMDUtility::GetEta() const
{
  return fEta;
}
Float_t AliPMDUtility::GetPhi() const
{
  return fPhi;
}
Commit	Line	Data
a918d77a	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	**************************************************************************/
638f6e9b	15	//-----------------------------------------------------//
	16	// //
	17	// //
	18	// Date : August 05 2003 //
	19	// //
	20	// Utility code for ALICE-PMD //
	21	// //
	22	//-----------------------------------------------------//
	23
a918d77a	24	#include "Riostream.h"
638f6e9b	25	#include "AliPMDUtility.h"
	26	#include "TMath.h"
	27	#include <stdio.h>
a918d77a	28	#include <math.h>
a918d77a	29
638f6e9b	30
	31	ClassImp(AliPMDUtility)
	32
	33	AliPMDUtility::AliPMDUtility()
	34	{
a918d77a	35	// Default constructor
638f6e9b	36	fPx = 0.;
	37	fPy = 0.;
	38	fPz = 0.;
	39	fTheta = 0.;
	40	fEta = 0.;
	41	fPhi = 0.;
	42	}
	43
a918d77a	44	AliPMDUtility::AliPMDUtility(Float_t px, Float_t py, Float_t pz)
638f6e9b	45	{
a918d77a	46	// Constructor
	47	fPx = px;
	48	fPy = py;
	49	fPz = pz;
638f6e9b	50	fTheta = 0.;
	51	fEta = 0.;
	52	fPhi = 0.;
	53	}
	54
	55	AliPMDUtility::~AliPMDUtility()
	56	{
a918d77a	57	// Default destructor
638f6e9b	58	}
afb8e3a0	59	void AliPMDUtility::HexGeomCellPos(Int_t ism, Int_t xpad, Int_t ypad, Float_t &xpos, Float_t &ypos)
afb8e3a0	60	{
a918d77a	61	// This converts PMD cluster or CELL coordinates
	62	// to Global coordinates.
	63	// Written by Prof. S.C. Phatak
	64
	65	const Float_t kCellDia = 0.5;
	66	const Float_t kPi = TMath::Pi(); //3.14159;
	67	const Double_t kSqroot3by2 = 0.8660254; // sqrth = sqrt(3.)/2.
afb8e3a0	68
a918d77a	69	Int_t i;
afb8e3a0	70	Int_t j = xpad;
	71	Int_t k = ypad;
	72
afb8e3a0	73	/*
a918d77a	74	Supermodeule number starting from 0
afb8e3a0	75	ism --> supermodule no ( 0 - 26 )
	76	idet --> detector ( pmd or cpv : not required now )
	77	j --> xpad ( goes from 1 to 72 )
	78	k --> ypad ( goes from 1 to 72 )
	79	xp --> global x coordinate
	80	yp --> global y coordinate
	81
	82	(xp0,yp0) corner positions of all supermodules in global
	83	coordinate system. That is the origin
	84	of the local ( supermodule ) coordinate system.
a918d77a	85	*/
afb8e3a0	86
	87	Float_t xp0[27] =
	88	{
	89	-17.9084, 18.2166, 54.3416, -35.9709, 0.154144,
	90	36.2791, -54.0334, -17.9084, 18.2166, 36.7791,
	91	18.7166, 0.654194, 72.9041, 54.8416, 36.7792,
	92	109.029, 90.9666, 72.9042, -18.8708, -36.9334,
	93	-54.996, -36.9332, -54.9958, -73.0584, -54.9956,
	94	-73.0582, -91.1208
	95	};
	96
	97	Float_t yp0[27] =
	98	{
	99	-32.1395, -32.1395, -32.1395, -63.4247, -63.4247,
	100	-63.4247, -94.7098, -94.7098, -94.7098, 0.545689,
	101	31.8309, 63.1161, 0.545632, 31.8308, 63.116,
	102	0.545573, 31.8308, 63.116, 31.5737, 0.288616,
	103	-30.9965, 62.859, 31.5738, 0.288733, 94.1442,
	104	62.8591, 31.574
	105	};
	106
	107	/*
	108	angles of rotation for three sets of supermodules
	109	The angle is same for first nine, next nine and last nine
	110	supermodules
	111	*/
	112
a918d77a	113	Float_t th[3] = {0., -2.kPi/3., 2.kPi/3.};
afb8e3a0	114	Float_t xr, yr, xinit, yinit, cs, sn;
	115
	116	/*
	117	xinit and yinit are coordinates of the cell in local coordinate system
	118	*/
	119
a918d77a	120	xinit = (j)kCellDia+(k)/2.kCellDia;
a918d77a	121	yinit = kSqroot3by2*(k)/2.;
afb8e3a0	122	i=ism/9;
	123	cs=cos(th[i]);
	124	sn=sin(th[i]);
	125	//
	126	// rotate first
	127	//
	128	xr=csxinit+snyinit;
	129	yr=-snxinit+csyinit;
	130	//
	131	// then translate
	132	//
	133	xpos=xr+xp0[ism];
	134	ypos=yr+yp0[ism];
	135
	136	}
	137
	138	void AliPMDUtility::RectGeomCellPos(Int_t ism, Int_t ium, Int_t xpad, Int_t ypad, Float_t &xpos, Float_t &ypos)
	139	{
	140	// This routine finds the cell eta,phi for the new PMD rectangular
	141	// geometry in ALICE
	142	// Authors : Bedanga Mohanty and Dipak Mishra - 29.4.2003
	143	// modified by B. K. Nnadi for change of coordinate sys
	144	//
	145	// SMA ---> Supermodule Type A ( SM - 0)
	146	// SMAR ---> Supermodule Type A ROTATED ( SM - 1)
	147	// SMB ---> Supermodule Type B ( SM - 2)
	148	// SMBR ---> Supermodule Type B ROTATED ( SM - 3)
	149	//
	150	// ism : number of supermodules in one plane = 4
	151	// ium : number of unitmodules in one SM = 6
a918d77a	152	// gbum : (global) unit module numbering in a supermodule
afb8e3a0	153	//
afb8e3a0	154
a918d77a	155	Int_t gbum = ism*6 + ium;
afb8e3a0	156	Int_t irow = xpad;
	157	Int_t icol = ypad;
	158
	159	// Corner positions (x,y) of the 24 unit moudles in ALICE PMD
	160
a918d77a	161	Double_t xcorner[24] =
afb8e3a0	162	{
	163	85.15, 60.85, 36.55, 85.15, 60.85, 36.55, //SMA
	164	-85.15, -60.85, -36.55, -85.15, -60.85, -36.55, //SMAR
	165	84.90, 36.60, 84.90, 36.60, 84.90, 36.60, //SMB
	166	-84.90, -36.60, -84.90, -36.60, -84.90, -36.60 //SMBR
	167	};
	168
a918d77a	169	Double_t ycorner[24] =
afb8e3a0	170	{
	171	32.45708755, 32.45708755, 32.45708755, //SMA
	172	-9.30645245, -9.30645245, -9.30645245, //SMA
	173	-32.45708755, -32.45708755, -32.45708755, //SMAR
	174	9.30645245, 9.30645245, 9.30645245, //SMAR
	175	-31.63540818, -31.63540818, -52.61435544, //SMB
	176	-52.61435544, -73.59330270, -73.59330270, //SMB
	177	31.63540818, 31.63540818, 52.61435544, //SMBR
	178	52.61435544, 73.59330270, 73.59330270 //SMBR
	179	};
	180
a918d77a	181	const Float_t kSqroot3 = 1.73205; // sqrt(3.);
a918d77a	182	const Float_t kCellRadius = 0.25;
afb8e3a0	183
	184	//
	185	//Every even row of cells is shifted and placed
	186	//in geant so this condition
	187	//
a2441c6e	188	Float_t shift = 0.0;
afb8e3a0	189	if(irow%2 == 0)
	190	{
	191	shift = 0.25;
	192	}
	193	else
	194	{
	195	shift = 0.0;
	196	}
	197	if(ism == 0 \|\| ism == 2)
	198	{
a918d77a	199	ypos = ycorner[gbum] +
a918d77a	200	irowkCellRadiuskSqroot3;
afb8e3a0	201
a918d77a	202	xpos = xcorner[gbum] -
a918d77a	203	icol2.0kCellRadius - shift;
afb8e3a0	204	}
	205	else if(ism == 1 \|\| ism == 3)
	206	{
a918d77a	207	ypos = ycorner[gbum] -
a918d77a	208	irowkCellRadiuskSqroot3;
afb8e3a0	209
a918d77a	210	xpos = xcorner[gbum] +
a918d77a	211	icol2.0kCellRadius + shift;
a2441c6e	212	}
	213	}
	214
	215	void AliPMDUtility::RectGeomCellPos(Int_t ism, Int_t ium, Float_t xpad, Float_t ypad, Float_t &xpos, Float_t &ypos)
	216	{
	217	// If the xpad and ypad inputs are float, then 0.5 is added to it
	218	// to find the layer which is shifted.
	219	// This routine finds the cell eta,phi for the new PMD rectangular
	220	// geometry in ALICE
	221	// Authors : Bedanga Mohanty and Dipak Mishra - 29.4.2003
	222	// modified by B. K. Nnadi for change of coordinate sys
	223	//
	224	// SMA ---> Supermodule Type A ( SM - 0)
	225	// SMAR ---> Supermodule Type A ROTATED ( SM - 1)
	226	// SMB ---> Supermodule Type B ( SM - 2)
	227	// SMBR ---> Supermodule Type B ROTATED ( SM - 3)
	228	//
	229	// ism : number of supermodules in one plane = 4
	230	// ium : number of unitmodules in one SM = 6
a918d77a	231	// gbum : (global) unit module numbering in a supermodule
a2441c6e	232	//
a2441c6e	233
a918d77a	234	Int_t gbum = ism*6 + ium;
a2441c6e	235	Float_t irow = xpad;
	236	Float_t icol = ypad;
	237
	238	// Corner positions (x,y) of the 24 unit moudles in ALICE PMD
	239
a918d77a	240	Double_t xcorner[24] =
a2441c6e	241	{
	242	85.15, 60.85, 36.55, 85.15, 60.85, 36.55, //SMA
	243	-85.15, -60.85, -36.55, -85.15, -60.85, -36.55, //SMAR
	244	84.90, 36.60, 84.90, 36.60, 84.90, 36.60, //SMB
	245	-84.90, -36.60, -84.90, -36.60, -84.90, -36.60 //SMBR
	246	};
	247
a918d77a	248	Double_t ycorner[24] =
a2441c6e	249	{
	250	32.45708755, 32.45708755, 32.45708755, //SMA
	251	-9.30645245, -9.30645245, -9.30645245, //SMA
	252	-32.45708755, -32.45708755, -32.45708755, //SMAR
	253	9.30645245, 9.30645245, 9.30645245, //SMAR
	254	-31.63540818, -31.63540818, -52.61435544, //SMB
	255	-52.61435544, -73.59330270, -73.59330270, //SMB
	256	31.63540818, 31.63540818, 52.61435544, //SMBR
	257	52.61435544, 73.59330270, 73.59330270 //SMBR
	258	};
	259
a918d77a	260	const Float_t kSqroot3 = 1.73205; // sqrt(3.);
a918d77a	261	const Float_t kCellRadius = 0.25;
a2441c6e	262
	263	//
	264	//Every even row of cells is shifted and placed
	265	//in geant so this condition
	266	//
	267	Float_t shift = 0.0;
	268	Int_t iirow = (Int_t) (irow+0.5);
	269	if(iirow%2 == 0)
	270	{
	271	shift = 0.25;
	272	}
	273	else
	274	{
	275	shift = 0.0;
	276	}
	277	if(ism == 0 \|\| ism == 2)
	278	{
a918d77a	279	ypos = ycorner[gbum] +
a918d77a	280	irowkCellRadiuskSqroot3;
a2441c6e	281
a918d77a	282	xpos = xcorner[gbum] -
a918d77a	283	icol2.0kCellRadius - shift;
a2441c6e	284	}
	285	else if(ism == 1 \|\| ism == 3)
	286	{
a918d77a	287	ypos = ycorner[gbum] -
a918d77a	288	irowkCellRadiuskSqroot3;
a2441c6e	289
a918d77a	290	xpos = xcorner[gbum] +
a918d77a	291	icol2.0kCellRadius + shift;
afb8e3a0	292	}
afb8e3a0	293	}
638f6e9b	294
a918d77a	295	void AliPMDUtility::SetPxPyPz(Float_t px, Float_t py, Float_t pz)
638f6e9b	296	{
a918d77a	297	fPx = px;
	298	fPy = py;
	299	fPz = pz;
638f6e9b	300	}
638f6e9b	301
a918d77a	302	void AliPMDUtility::SetXYZ(Float_t xpos, Float_t ypos, Float_t zpos)
638f6e9b	303	{
a918d77a	304	fPx = xpos;
	305	fPy = ypos;
	306	fPz = zpos;
638f6e9b	307	}
	308	void AliPMDUtility::CalculateEta()
	309	{
	310	Float_t rpxpy, theta, eta;
	311
	312	rpxpy = TMath::Sqrt(fPxfPx + fPyfPy);
	313	theta = TMath::ATan2(rpxpy,fPz);
	314	eta = -TMath::Log(TMath::Tan(0.5*theta));
	315	fTheta = theta;
	316	fEta = eta;
	317	}
	318	void AliPMDUtility::CalculatePhi()
	319	{
	320	Float_t pybypx, phi = 0., phi1;
	321
	322	if(fPx==0)
	323	{
	324	if(fPy>0) phi = 90.;
	325	if(fPy<0) phi = 270.;
	326	}
	327	if(fPx != 0)
	328	{
	329	pybypx = fPy/fPx;
	330	if(pybypx < 0) pybypx = - pybypx;
	331	phi1 = TMath::ATan(pybypx)*180./3.14159;
afb8e3a0	332
	333	if(fPx > 0 && fPy > 0) phi = phi1; // 1st Quadrant
	334	if(fPx < 0 && fPy > 0) phi = 180 - phi1; // 2nd Quadrant
	335	if(fPx < 0 && fPy < 0) phi = 180 + phi1; // 3rd Quadrant
	336	if(fPx > 0 && fPy < 0) phi = 360 - phi1; // 4th Quadrant
	337
638f6e9b	338	}
	339	phi = phi*3.14159/180.;
	340
	341	fPhi = phi;
	342
	343	}
	344	void AliPMDUtility::CalculateEtaPhi()
	345	{
	346	Float_t rpxpy, theta, eta;
	347	Float_t pybypx, phi = 0., phi1;
	348
	349	rpxpy = TMath::Sqrt(fPxfPx + fPyfPy);
	350	theta = TMath::ATan2(rpxpy,fPz);
	351	eta = -TMath::Log(TMath::Tan(0.5*theta));
	352
	353	if(fPx==0)
	354	{
	355	if(fPy>0) phi = 90.;
	356	if(fPy<0) phi = 270.;
	357	}
	358	if(fPx != 0)
	359	{
	360	pybypx = fPy/fPx;
	361	if(pybypx < 0) pybypx = - pybypx;
	362	phi1 = TMath::ATan(pybypx)*180./3.14159;
afb8e3a0	363	if(fPx > 0 && fPy > 0) phi = phi1; // 1st Quadrant
	364	if(fPx < 0 && fPy > 0) phi = 180 - phi1; // 2nd Quadrant
	365	if(fPx < 0 && fPy < 0) phi = 180 + phi1; // 3rd Quadrant
	366	if(fPx > 0 && fPy < 0) phi = 360 - phi1; // 4th Quadrant
	367
638f6e9b	368	}
	369	phi = phi*3.14159/180.;
	370
	371	fTheta = theta;
	372	fEta = eta;
	373	fPhi = phi;
	374	}
	375	Float_t AliPMDUtility::GetTheta() const
	376	{
	377	return fTheta;
	378	}
	379	Float_t AliPMDUtility::GetEta() const
	380	{
	381	return fEta;
	382	}
	383	Float_t AliPMDUtility::GetPhi() const
	384	{
	385	return fPhi;
	386	}
	387