#include <TMath.h>
#include "AliTRDpadPlane.h"
-#include "AliTRDgeometry.h"
ClassImp(AliTRDpadPlane)
//_____________________________________________________________________________
AliTRDpadPlane::AliTRDpadPlane()
:TObject()
- ,fGeo(0)
- ,fPla(0)
- ,fCha(0)
+ ,fLayer(0)
+ ,fStack(0)
,fLength(0)
,fWidth(0)
,fLengthRim(0)
,fPadRow(0)
,fPadCol(0)
,fPadRowSMOffset(0)
+ ,fAnodeWireOffset(0)
{
//
// Default constructor
}
//_____________________________________________________________________________
-AliTRDpadPlane::AliTRDpadPlane(Int_t p, Int_t c)
+AliTRDpadPlane::AliTRDpadPlane(Int_t layer, Int_t stack)
:TObject()
- ,fGeo(0)
- ,fPla(0)
- ,fCha(0)
+ ,fLayer(layer)
+ ,fStack(stack)
,fLength(0)
,fWidth(0)
,fLengthRim(0)
,fPadRow(0)
,fPadCol(0)
,fPadRowSMOffset(0)
+ ,fAnodeWireOffset(0)
{
//
- // Constructor that initializes a given pad plane type
+ // Constructor
//
- fGeo = new AliTRDgeometry();
-
- fPla = p;
- fCha = c;
-
- fRowSpacing = 0.0;
- fColSpacing = 0.0;
-
- fLengthRim = 1.0;
- fWidthRim = 0.5;
-
- fNcols = 144;
-
- //
- // The pad plane parameter
- //
- switch (p) {
- case 0:
- if (c == 2) {
- // L0C0 type
- fNrows = 12;
- fLength = 108.0;
- fWidth = 92.2;
- fLengthOPad = 8.0;
- fWidthOPad = 0.515;
- fLengthIPad = 9.0;
- fWidthIPad = 0.635;
- fTiltingAngle = -2.0;
- }
- else {
- // L0C1 type
- fNrows = 16;
- fLength = 122.0;
- fWidth = 92.2;
- fLengthOPad = 7.5;
- fWidthOPad = 0.515;
- fLengthIPad = 7.5;
- fWidthIPad = 0.635;
- fTiltingAngle = -2.0;
- }
- break;
- case 1:
- if (c == 2) {
- // L1C0 type
- fNrows = 12;
- fLength = 108.0;
- fWidth = 96.6;
- fLengthOPad = 8.0;
- fWidthOPad = 0.585;
- fLengthIPad = 9.0;
- fWidthIPad = 0.665;
- fTiltingAngle = 2.0;
- }
- else {
- // L1C1 type
- fNrows = 16;
- fLength = 122.0;
- fWidth = 96.6;
- fLengthOPad = 7.5;
- fWidthOPad = 0.585;
- fLengthIPad = 7.5;
- fWidthIPad = 0.665;
- fTiltingAngle = 2.0;
- }
- break;
- case 2:
- if (c == 2) {
- // L2C0 type
- fNrows = 12;
- fLength = 108.0;
- fWidth = 101.1;
- fLengthOPad = 8.0;
- fWidthOPad = 0.705;
- fLengthIPad = 9.0;
- fWidthIPad = 0.695;
- fTiltingAngle = -2.0;
- }
- else {
- // L2C1 type
- fNrows = 16;
- fLength = 129.0;
- fWidth = 101.1;
- fLengthOPad = 7.5;
- fWidthOPad = 0.705;
- fLengthIPad = 8.0;
- fWidthIPad = 0.695;
- fTiltingAngle = -2.0;
- }
- break;
- case 3:
- if (c == 2) {
- // L3C0 type
- fNrows = 12;
- fLength = 108.0;
- fWidth = 105.5;
- fLengthOPad = 8.0;
- fWidthOPad = 0.775;
- fLengthIPad = 9.0;
- fWidthIPad = 0.725;
- fTiltingAngle = 2.0;
- }
- else {
- // L3C1 type
- fNrows = 16;
- fLength = 136.0;
- fWidth = 105.5;
- fLengthOPad = 7.5;
- fWidthOPad = 0.775;
- fLengthIPad = 8.5;
- fWidthIPad = 0.725;
- fTiltingAngle = 2.0;
- }
- break;
- case 4:
- if (c == 2) {
- // L4C0 type
- fNrows = 12;
- fLength = 108.0;
- fWidth = 109.9;
- fLengthOPad = 8.0;
- fWidthOPad = 0.845;
- fLengthIPad = 9.0;
- fWidthIPad = 0.755;
- fTiltingAngle = -2.0;
- }
- else {
- // L4C1 type
- fNrows = 16;
- fLength = 143.0;
- fWidth = 109.9;
- fLengthOPad = 7.5;
- fWidthOPad = 0.845;
- fLengthIPad = 9.0;
- fWidthIPad = 0.755;
- fTiltingAngle = -2.0;
- }
- break;
- case 5:
- if (c == 2) {
- // L5C0 type
- fNrows = 12;
- fLength = 108.0;
- fWidth = 114.4;
- fLengthOPad = 8.0;
- fWidthOPad = 0.965;
- fLengthIPad = 9.0;
- fWidthIPad = 0.785;
- fTiltingAngle = 2.0;
- }
- else {
- // L5C1 type
- fNrows = 16;
- fLength = 145.0;
- fWidth = 114.4;
- fLengthOPad = 8.5;
- fWidthOPad = 0.965;
- fLengthIPad = 9.0;
- fWidthIPad = 0.785;
- fTiltingAngle = 2.0;
- }
- break;
- };
-
- //
- // Store tilting angle as tangens
- //
- fTiltingTan = TMath::Tan(TMath::Pi()/180.0 * fTiltingAngle);
-
- //
- // The positions of the borders of the pads
- //
- // Row direction
- //
- fPadRow = new Double_t[fNrows];
- Double_t row = fGeo->GetChamberLength(p,c) / 2.0
- - fGeo->RpadW()
- - fLengthRim;
- for (Int_t ir = 0; ir < fNrows; ir++) {
- fPadRow[ir] = row;
- row -= fRowSpacing;
- if (ir == 0) {
- row -= fLengthOPad;
- }
- else {
- row -= fLengthIPad;
- }
- }
- //
- // Column direction
- //
- fPadCol = new Double_t[fNcols];
- Double_t col = fGeo->GetChamberWidth(p) / 2.0
- + fGeo->CroWid()
- - fWidthRim;
- for (Int_t ic = 0; ic < fNcols; ic++) {
- fPadCol[ic] = col;
- col -= fColSpacing;
- if (ic == 0) {
- col -= fWidthOPad;
- }
- else {
- col -= fWidthIPad;
- }
- }
- // Calculate the offset to translate from the local ROC system into
- // the local supermodule system, which is used for clusters
- Double_t rowTmp = fGeo->GetChamberLength(p,0)
- + fGeo->GetChamberLength(p,1)
- + fGeo->GetChamberLength(p,2) / 2.0;
- for (Int_t ic = 0; ic < c; ic++) {
- rowTmp -= fGeo->GetChamberLength(p,ic);
- }
- fPadRowSMOffset = rowTmp - fGeo->GetChamberLength(p,c)/2.0;
-
-}
-
-//_____________________________________________________________________________
-AliTRDpadPlane::AliTRDpadPlane(const AliTRDpadPlane &p)
- :TObject(p)
- ,fGeo(0)
- ,fPla(p.fPla)
- ,fCha(p.fCha)
- ,fLength(p.fLength)
- ,fWidth(p.fWidth)
- ,fLengthRim(p.fLengthRim)
- ,fWidthRim(p.fLengthRim)
- ,fLengthOPad(p.fLengthOPad)
- ,fWidthOPad(p.fWidthOPad)
- ,fLengthIPad(p.fLengthIPad)
- ,fWidthIPad(p.fWidthIPad)
- ,fRowSpacing(p.fRowSpacing)
- ,fColSpacing(p.fColSpacing)
- ,fNrows(p.fNrows)
- ,fNcols(p.fNcols)
- ,fTiltingAngle(p.fTiltingAngle)
- ,fTiltingTan(p.fTiltingTan)
- ,fPadRow(0)
- ,fPadCol(0)
- ,fPadRowSMOffset(p.fPadRowSMOffset)
-{
- //
- // AliTRDpadPlane copy constructor
- //
-
- Int_t iBin = 0;
-
- if (((AliTRDpadPlane &) p).fPadRow) {
- delete [] ((AliTRDpadPlane &) p).fPadRow;
- }
- ((AliTRDpadPlane &) p).fPadRow = new Double_t[fNrows];
- for (iBin = 0; iBin < fNrows; iBin++) {
- ((AliTRDpadPlane &) p).fPadRow[iBin] = fPadRow[iBin];
- }
-
- if (((AliTRDpadPlane &) p).fPadCol) {
- delete [] ((AliTRDpadPlane &) p).fPadCol;
- }
- ((AliTRDpadPlane &) p).fPadCol = new Double_t[fNrows];
- for (iBin = 0; iBin < fNrows; iBin++) {
- ((AliTRDpadPlane &) p).fPadCol[iBin] = fPadCol[iBin];
- }
-
}
//_____________________________________________________________________________
// AliTRDpadPlane destructor
//
- if (fGeo) {
- delete fGeo;
- fGeo = 0;
- }
-
if (fPadRow) {
delete [] fPadRow;
fPadRow = 0;
}
-//_____________________________________________________________________________
-AliTRDpadPlane &AliTRDpadPlane::operator=(const AliTRDpadPlane &p)
-{
- //
- // Assignment operator
- //
-
- if (this != &p) {
- ((AliTRDpadPlane &) p).Copy(*this);
- }
-
- return *this;
-
-}
-
//_____________________________________________________________________________
void AliTRDpadPlane::Copy(TObject &p) const
{
Int_t iBin = 0;
- ((AliTRDpadPlane &) p).fGeo = 0;
+ ((AliTRDpadPlane &) p).fLayer = fLayer;
+ ((AliTRDpadPlane &) p).fStack = fStack;
- ((AliTRDpadPlane &) p).fPla = fPla;
- ((AliTRDpadPlane &) p).fCha = fCha;
+ ((AliTRDpadPlane &) p).fLength = fLength;
+ ((AliTRDpadPlane &) p).fWidth = fWidth;
+ ((AliTRDpadPlane &) p).fLengthRim = fLengthRim;
+ ((AliTRDpadPlane &) p).fWidthRim = fWidthRim;
+ ((AliTRDpadPlane &) p).fLengthOPad = fLengthOPad;
+ ((AliTRDpadPlane &) p).fWidthOPad = fWidthOPad;
+ ((AliTRDpadPlane &) p).fLengthIPad = fLengthIPad;
+ ((AliTRDpadPlane &) p).fWidthIPad = fWidthIPad;
- ((AliTRDpadPlane &) p).fLength = fLength;
- ((AliTRDpadPlane &) p).fWidth = fWidth;
- ((AliTRDpadPlane &) p).fLengthRim = fLengthRim;
- ((AliTRDpadPlane &) p).fWidthRim = fWidthRim;
- ((AliTRDpadPlane &) p).fLengthOPad = fLengthOPad;
- ((AliTRDpadPlane &) p).fWidthOPad = fWidthOPad;
- ((AliTRDpadPlane &) p).fLengthIPad = fLengthIPad;
- ((AliTRDpadPlane &) p).fWidthIPad = fWidthIPad;
+ ((AliTRDpadPlane &) p).fRowSpacing = fRowSpacing;
+ ((AliTRDpadPlane &) p).fColSpacing = fColSpacing;
- ((AliTRDpadPlane &) p).fRowSpacing = fRowSpacing;
- ((AliTRDpadPlane &) p).fColSpacing = fColSpacing;
+ ((AliTRDpadPlane &) p).fNrows = fNrows;
+ ((AliTRDpadPlane &) p).fNcols = fNcols;
- ((AliTRDpadPlane &) p).fNrows = fNrows;
- ((AliTRDpadPlane &) p).fNcols = fNcols;
+ ((AliTRDpadPlane &) p).fTiltingAngle = fTiltingAngle;
+ ((AliTRDpadPlane &) p).fTiltingTan = fTiltingTan;
- ((AliTRDpadPlane &) p).fTiltingAngle = fTiltingAngle;
- ((AliTRDpadPlane &) p).fTiltingTan = fTiltingTan;
-
- ((AliTRDpadPlane &) p).fPadRowSMOffset = fPadRowSMOffset;
+ ((AliTRDpadPlane &) p).fPadRowSMOffset = fPadRowSMOffset;
+ ((AliTRDpadPlane &) p).fAnodeWireOffset = fAnodeWireOffset;
if (((AliTRDpadPlane &) p).fPadRow) {
delete [] ((AliTRDpadPlane &) p).fPadRow;
}
+//_____________________________________________________________________________
+void AliTRDpadPlane::SetTiltingAngle(Double_t t)
+{
+ //
+ // Set the tilting angle of the pads
+ //
+
+ fTiltingAngle = t;
+ fTiltingTan = TMath::Tan(TMath::Pi()/180.0 * fTiltingAngle);
+
+}
+
//_____________________________________________________________________________
Int_t AliTRDpadPlane::GetPadRowNumber(Double_t z) const
{
Int_t AliTRDpadPlane::GetPadColNumber(Double_t rphi) const
{
//
- // Finds the pad column number for a given global rphi-position
+ // Finds the pad column number for a given rphi-position
//
Int_t col = 0;
Int_t nbelow = 0;
Int_t middle = 0;
- if ((rphi > GetCol0() ) ||
- (rphi < GetColEnd())) {
+ if ((rphi < GetCol0() ) ||
+ (rphi > GetColEnd())) {
col = -1;
}
else {
- nabove = fNcols + 1;
+ nabove = fNcols;
nbelow = 0;
while (nabove - nbelow > 1) {
middle = (nabove + nbelow) / 2;
- if (rphi == fPadCol[middle-1]) {
+ if (rphi == fPadCol[middle]) {
col = middle;
}
- if (rphi > fPadCol[middle-1]) {
- nabove = middle;
+ if (rphi > fPadCol[middle]) {
+ nbelow = middle;
}
else {
- nbelow = middle;
+ nabove = middle;
}
}
- col = nbelow - 1;
+ col = nbelow;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff