/*************************************************************************** * 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 #include 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::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 xcorner[24] = { 74.8833, 53.0045, 31.1255, //Type-A 74.8833, 53.0045, 31.1255, //Type-A -74.8833, -53.0044, -31.1255, //Type-AR -74.8833, -53.0044, -31.1255, //Type-AR 8.9165, -33.7471, //Type-B 8.9165, -33.7471, //Type-B 8.9165, -33.7471, //Type-B -8.9165, 33.7471, //Type-BR -8.9165, 33.7471, //Type-BR -8.9165, 33.7471, //Type-BR }; double ycorner[24] = { 86.225, 86.225, 86.225, //Type-A 37.075, 37.075, 37.075, //Type-A -86.225, -86.225, -86.225, //Type-AR -37.075, -37.075, -37.075, //Type-AR 86.225, 86.225, //Type-B 61.075, 61.075, //Type-B 35.925, 35.925, //Type-B -86.225, -86.225, //Type-BR -61.075, -61.075, //Type-BR -35.925, -35.925 //Type-BR }; 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 cellRadius = 0.25; Float_t shift = 0.0; if(irow%2 == 0) { shift = -cellRadius/2.0; } else { shift = 0.0; } if(ism == 0) { ypos = ycorner[gbum] - irow*kCellRadius*2.0 + shift; xpos = xcorner[gbum] - icol*kSqroot3*kCellRadius; } else if(ism == 1) { ypos = ycorner[gbum] + irow*kCellRadius*2.0 + shift; xpos = xcorner[gbum] + icol*kSqroot3*kCellRadius; } else if(ism == 2) { ypos = ycorner[gbum] - irow*kCellRadius*2.0 + shift; xpos = xcorner[gbum] - icol*kSqroot3*kCellRadius; } else if(ism == 3) { ypos = ycorner[gbum] + irow*kCellRadius*2.0 + shift; xpos = xcorner[gbum] + icol*kSqroot3*kCellRadius; } } 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 xcorner[24] = { 74.8833, 53.0045, 31.1255, //Type-A 74.8833, 53.0045, 31.1255, //Type-A -74.8833, -53.0044, -31.1255, //Type-AR -74.8833, -53.0044, -31.1255, //Type-AR 8.9165, -33.7471, //Type-B 8.9165, -33.7471, //Type-B 8.9165, -33.7471, //Type-B -8.9165, 33.7471, //Type-BR -8.9165, 33.7471, //Type-BR -8.9165, 33.7471, //Type-BR }; double ycorner[24] = { 86.225, 86.225, 86.225, //Type-A 37.075, 37.075, 37.075, //Type-A -86.225, -86.225, -86.225, //Type-AR -37.075, -37.075, -37.075, //Type-AR 86.225, 86.225, //Type-B 61.075, 61.075, //Type-B 35.925, 35.925, //Type-B -86.225, -86.225, //Type-BR -61.075, -61.075, //Type-BR -35.925, -35.925 //Type-BR }; 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 cellRadius = 0.25; Float_t shift = 0.0; Int_t iirow = (Int_t) (irow+0.5); if(iirow%2 == 0) { shift = -cellRadius/2.0; } else { shift = 0.0; } if(ism == 0) { ypos = ycorner[gbum] - irow*kCellRadius*2.0 + shift; xpos = xcorner[gbum] - icol*kSqroot3*kCellRadius; } else if(ism == 1) { ypos = ycorner[gbum] + irow*kCellRadius*2.0 + shift; xpos = xcorner[gbum] + icol*kSqroot3*kCellRadius; } else if(ism == 2) { ypos = ycorner[gbum] - irow*kCellRadius*2.0 + shift; xpos = xcorner[gbum] - icol*kSqroot3*kCellRadius; } else if(ism == 3) { ypos = ycorner[gbum] + irow*kCellRadius*2.0 + shift; xpos = xcorner[gbum] + icol*kSqroot3*kCellRadius; } } 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; }