Cleanup
authorcblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 19 Dec 2005 16:38:01 +0000 (16:38 +0000)
committercblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 19 Dec 2005 16:38:01 +0000 (16:38 +0000)
TRD/AliTRDCalROCVdrift.cxx [deleted file]
TRD/AliTRDCalROCVdrift.h [deleted file]
TRD/AliTRDCalVdrift.cxx [deleted file]
TRD/AliTRDCalVdrift.h [deleted file]
TRD/AliTRDclusterizerV0.cxx [deleted file]
TRD/AliTRDclusterizerV0.h [deleted file]
TRD/AliTRDgeometryHole.cxx [deleted file]
TRD/AliTRDgeometryHole.h [deleted file]
TRD/AliTRDtrackingSector.cxx [deleted file]
TRD/AliTRDtrackingSector.h [deleted file]

diff --git a/TRD/AliTRDCalROCVdrift.cxx b/TRD/AliTRDCalROCVdrift.cxx
deleted file mode 100644 (file)
index ece01bd..0000000
+++ /dev/null
@@ -1,111 +0,0 @@
-/**************************************************************************
- * 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$ */
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-//  Calibration base class for a single ROC                                  //
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
-
-#include "AliTRDCalROCVdrift.h"
-
-ClassImp(AliTRDCalROCVdrift)
-
-//_____________________________________________________________________________
-AliTRDCalROCVdrift::AliTRDCalROCVdrift():AliTRDCalROC()
-{
-  //
-  // Default constructor
-  //
-
-  fNchannels    = 0;
-  fVdrift       = 0;
-
-}
-
-//_____________________________________________________________________________
-AliTRDCalROCVdrift::AliTRDCalROCVdrift(Int_t p, Int_t c)
-                   :AliTRDCalROC(p,c)
-{
-  //
-  // Constructor that initializes a given pad plane type
-  //
-
-  fNchannels = fNrows * fNcols;
-  fVdrift    = new Float_t[fNchannels];
-
-}
-
-//_____________________________________________________________________________
-AliTRDCalROCVdrift::AliTRDCalROCVdrift(const AliTRDCalROCVdrift &c)
-                   :AliTRDCalROC(c)
-{
-  //
-  // AliTRDCalROCVdrift copy constructor
-  //
-
-  ((AliTRDCalROCVdrift &) c).Copy(*this);
-
-}
-
-//_____________________________________________________________________________
-AliTRDCalROCVdrift::~AliTRDCalROCVdrift()
-{
-  //
-  // AliTRDCalROCVdrift destructor
-  //
-
-  if (fVdrift) {
-    delete [] fVdrift;
-    fVdrift = 0;
-  }
-
-}
-
-//_____________________________________________________________________________
-AliTRDCalROCVdrift &AliTRDCalROCVdrift::operator=(const AliTRDCalROCVdrift &c)
-{
-  //
-  // Assignment operator
-  //
-
-  if (this != &c) ((AliTRDCalROCVdrift &) c).Copy(*this);
-  return *this;
-
-}
-
-//_____________________________________________________________________________
-void AliTRDCalROCVdrift::Copy(TObject &c) const
-{
-  //
-  // Copy function
-  //
-
-  Int_t iBin = 0;
-
-  ((AliTRDCalROCVdrift &) c).fNchannels = fNchannels;
-
-  if (((AliTRDCalROCVdrift &) c).fVdrift) delete [] ((AliTRDCalROCVdrift &) c).fVdrift;
-  ((AliTRDCalROCVdrift &) c).fVdrift = new Float_t[fNchannels];
-  for (iBin = 0; iBin < fNchannels; iBin++) {
-    ((AliTRDCalROCVdrift &) c).fVdrift[iBin] = fVdrift[iBin];
-  }                                                                             
-
-  AliTRDCalROC::Copy(c);
-
-}
-
diff --git a/TRD/AliTRDCalROCVdrift.h b/TRD/AliTRDCalROCVdrift.h
deleted file mode 100644 (file)
index 72321f2..0000000
+++ /dev/null
@@ -1,46 +0,0 @@
-#ifndef ALITRDCALROCVDRIFT_H
-#define ALITRDCALROCVDRIFT_H
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * See cxx source for full Copyright notice                               */
-
-/* $Id: AliTRDCalROCVdrift.h,v */
-
-///////////////////////////////////////////////////
-//                                               //
-//  TRD calibration class for Vdrift in one ROC  //
-//                                               //
-///////////////////////////////////////////////////
-
-#include "AliTRDCalROC.h"
-
-//_____________________________________________________________________________
-class AliTRDCalROCVdrift : public AliTRDCalROC {
-
- public:
-
-  AliTRDCalROCVdrift();
-  AliTRDCalROCVdrift(Int_t p, Int_t c);
-  AliTRDCalROCVdrift(const AliTRDCalROCVdrift &c);
-  virtual             ~AliTRDCalROCVdrift();
-  AliTRDCalROCVdrift  &operator=(const AliTRDCalROCVdrift &c);
-  virtual void         Copy(TObject &c) const;
-
-  Int_t        GetChannel(Int_t c, Int_t r)     { return r+c*fNrows; };
-  Int_t        GetNchannels()       const       { return fNchannels;   };
-  Float_t      GetVdrift(Int_t ich) const       { return fVdrift[ich]; };
-  Float_t      GetVdrift(Int_t col, Int_t row)  { return fVdrift[GetChannel(col,row)]; };
-
-  void         SetVdrift(Int_t ich, Float_t vd) { fVdrift[ich] = vd;   };
-  void         SetVdrift(Int_t col, Int_t row, Float_t vd) 
-                                                { fVdrift[GetChannel(col,row)] = vd; };
-
- protected:
-
-  Int_t     fNchannels;             //  Number of channels
-  Float_t  *fVdrift;                //[fNchannels] Drift velocities
-
-  ClassDef(AliTRDCalROCVdrift,1)    //  TRD ROC calibration class for Vdrift
-
-};
-
-#endif
diff --git a/TRD/AliTRDCalVdrift.cxx b/TRD/AliTRDCalVdrift.cxx
deleted file mode 100644 (file)
index d5e2d37..0000000
+++ /dev/null
@@ -1,117 +0,0 @@
-/**************************************************************************
- * 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$ */
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-//  TRD calibration class for Vdrift                                         //
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
-
-#include "AliTRDCalVdrift.h"
-#include "AliTRDCalROCVdrift.h"
-
-ClassImp(AliTRDCalVdrift)
-
-//_____________________________________________________________________________
-AliTRDCalVdrift::AliTRDCalVdrift():TNamed()
-{
-  //
-  // AliTRDCalVdrift default constructor
-  //
-
-  for (Int_t idet = 0; idet < kNdet; idet++) {
-    fROCVdrift[idet] = 0;
-  }
-
-}
-
-//_____________________________________________________________________________
-AliTRDCalVdrift::AliTRDCalVdrift(const Text_t *name, const Text_t *title)
-                :TNamed(name,title)
-{
-  //
-  // AliTRDCalVdrift constructor
-  //
-
-  for (Int_t isec = 0; isec < kNsect; isec++) {
-    for (Int_t ipla = 0; ipla < kNplan; ipla++) {
-      for (Int_t icha = 0; icha < kNcham; icha++) {
-        Int_t idet = GetDet(ipla,icha,isec);
-        fROCVdrift[idet] = new AliTRDCalROCVdrift(ipla,icha);
-      }
-    }
-  }
-
-}
-
-
-//_____________________________________________________________________________
-AliTRDCalVdrift::AliTRDCalVdrift(const AliTRDCalVdrift &c):TNamed(c)
-{
-  //
-  // AliTRDCalVdrift copy constructor
-  //
-
-  ((AliTRDCalVdrift &) c).Copy(*this);
-
-}
-
-///_____________________________________________________________________________
-AliTRDCalVdrift::~AliTRDCalVdrift()
-{
-  //
-  // AliTRDCalVdrift destructor
-  //
-
-  for (Int_t idet = 0; idet < kNdet; idet++) {
-    if (fROCVdrift[idet]) {
-      delete fROCVdrift[idet];
-      fROCVdrift[idet] = 0;
-    }
-  }
-
-}
-
-//_____________________________________________________________________________
-AliTRDCalVdrift &AliTRDCalVdrift::operator=(const AliTRDCalVdrift &c)
-{
-  //
-  // Assignment operator
-  //
-
-  if (this != &c) ((AliTRDCalVdrift &) c).Copy(*this);
-  return *this;
-
-}
-
-//_____________________________________________________________________________
-void AliTRDCalVdrift::Copy(TObject &c) const
-{
-  //
-  // Copy function
-  //
-
-  for (Int_t idet = 0; idet < kNdet; idet++) {
-    if (fROCVdrift[idet]) {
-      fROCVdrift[idet]->Copy(*((AliTRDCalVdrift &) c).fROCVdrift[idet]);
-    }
-  }
-
-  TObject::Copy(c);
-
-}
-
diff --git a/TRD/AliTRDCalVdrift.h b/TRD/AliTRDCalVdrift.h
deleted file mode 100644 (file)
index d37373f..0000000
+++ /dev/null
@@ -1,46 +0,0 @@
-#ifndef ALITRDCALVDRIFT_H
-#define ALITRDCALVDRIFT_H
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * See cxx source for full Copyright notice                               */
-
-/* $Id$ */
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-//  TRD calibration class for Vdrift                                         //
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
-
-#include "TNamed.h"
-
-class AliTRDCalROCVdrift;
-
-class AliTRDCalVdrift : public TNamed {
-
- public:
-  enum { kNplan = 6, kNcham = 5, kNsect = 18, kNdet = 540 };
-
-  AliTRDCalVdrift();
-  AliTRDCalVdrift(const Text_t* name, const Text_t* title);
-  AliTRDCalVdrift(const AliTRDCalVdrift &c);   
-  virtual ~AliTRDCalVdrift();
-  AliTRDCalVdrift &operator=(const AliTRDCalVdrift &c); 
-
-  virtual void     Copy(TObject &c) const;
-
-  Int_t               GetDet(Int_t p, Int_t c, Int_t s) { return p+c*kNplan+s*kNplan*kNcham; };
-
-  AliTRDCalROCVdrift *GetCalROCVdrift(Int_t d) { return fROCVdrift[d]; };
-  AliTRDCalROCVdrift *GetCalROCVdrift(Int_t p, Int_t c, Int_t s) 
-                                               { return fROCVdrift[GetDet(p,c,s)]; };
-
- protected:
-
-  AliTRDCalROCVdrift *fROCVdrift[kNdet];           //  Array of ROC Vdrift objects
-
-  ClassDef(AliTRDCalVdrift,1)                      //  TRD calibration class for Vdrift
-
-};
-
-#endif
diff --git a/TRD/AliTRDclusterizerV0.cxx b/TRD/AliTRDclusterizerV0.cxx
deleted file mode 100644 (file)
index e4a5101..0000000
+++ /dev/null
@@ -1,309 +0,0 @@
-/**************************************************************************
- * 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$ */
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-// TRD cluster finder for the fast simulator. It takes the hits from the     //
-// fast simulator (one hit per plane) and transforms them                    //
-// into cluster, by applying position smearing and merging                   //
-// of nearby cluster. The searing is done uniformly in z-direction           //
-// over the length of a readout pad. In rphi-direction a Gaussian            //
-// smearing is applied with a sigma given by fRphiSigma.                     //
-// Clusters are considered as overlapping when they are closer in            //
-// rphi-direction than the value defined in fRphiDist.                       //
-// Use the macro fastClusterCreate.C to create the cluster.                  //
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
-
-#include <TRandom.h>
-#include <TTree.h>
-#include "AliRun.h"
-
-#include "AliTRD.h"
-#include "AliTRDclusterizerV0.h"
-#include "AliTRDhit.h"
-#include "AliTRDgeometry.h"
-#include "AliTRDrecPoint.h"
-#include "AliTRDparameter.h"
-
-ClassImp(AliTRDclusterizerV0)
-
-//_____________________________________________________________________________
-  AliTRDclusterizerV0::AliTRDclusterizerV0():AliTRDclusterizer(), fTRD(NULL)
-{
-  //
-  // AliTRDclusterizerV0 default constructor
-  //
-
-}
-
-//_____________________________________________________________________________
-AliTRDclusterizerV0::AliTRDclusterizerV0(const Text_t* name, const Text_t* title)
-                    :AliTRDclusterizer(name,title), fTRD(NULL)
-{
-  //
-  // AliTRDclusterizerV0 default constructor
-  //
-
-  Init();
-
-}
-
-//_____________________________________________________________________________
-AliTRDclusterizerV0::~AliTRDclusterizerV0()
-{
-  //
-  // AliTRDclusterizerV0 destructor
-  //
-
-}
-
-//_____________________________________________________________________________
-void AliTRDclusterizerV0::Init()
-{
-  //
-  // Initializes the cluster finder
-  //
-
-  // Position resolution in rphi-direction
-  fRphiSigma  = 0.02;
-  // Minimum distance of non-overlapping cluster
-  fRphiDist   = 1.0;
-
-}
-
-//_____________________________________________________________________________
-Bool_t AliTRDclusterizerV0::MakeClusters()
-{
-  //
-  // Generates the cluster
-  //
-
-  // Get the TRD object
-  fTRD = (AliTRD*) gAlice->GetDetector("TRD"); 
-  if (!fTRD) {
-    printf("AliTRDclusterizerV0::makClusters -- ");
-    printf("No TRD detector object found\n");
-    return kFALSE;
-  }
-
-  if (fTRD->IsVersion() != 0) {
-    printf("AliTRDclusterizerV0::MakeCluster -- ");
-    printf("TRD must be version 0 (fast simulator).\n");
-    return kFALSE; 
-  }
-
-  // Get the geometry
-  AliTRDgeometry *geo = fTRD->GetGeometry();
-
-  // Create a default parameter class if none is defined
-  if (!fPar) {
-    fPar = new AliTRDparameter("TRDparameter","Standard TRD parameter");
-    if (fVerbose > 0) {
-      printf("<AliTRDclusterizerV0::MakeCluster> ");
-      printf("Create the default parameter object.\n");
-    }
-  }
-  
-  printf("AliTRDclusterizerV0::MakeCluster -- ");
-  printf("Start creating cluster.\n");
-
-  Int_t nBytes = 0;
-
-  AliTRDhit *hit;
-  
-  // Get the pointer to the hit tree
-  TTree     *hitTree      = fTRD->TreeH();
-  // Get the pointer to the reconstruction tree
-  TTree     *clusterTree  = gAlice->TreeR();
-
-  TObjArray *chamberArray = new TObjArray();
-
-  // Get the number of entries in the hit tree
-  // (Number of primary particles creating a hit somewhere)
-  Int_t nTrack = (Int_t) hitTree->GetEntries();
-
-  // Loop through all the chambers
-  for (Int_t icham = 0; icham < AliTRDgeometry::Ncham(); icham++) {
-    for (Int_t iplan = 0; iplan < AliTRDgeometry::Nplan(); iplan++) {
-      for (Int_t isect = 0; isect < AliTRDgeometry::Nsect(); isect++) {
-
-        Int_t   nColMax     = fPar->GetColMax(iplan);
-        Float_t row0        = fPar->GetRow0(iplan,icham,isect);
-        Float_t col0        = fPar->GetCol0(iplan);
-        Float_t time0       = fPar->GetTime0(iplan);
-
-        Float_t rowPadSize  = fPar->GetRowPadSize(iplan,icham,isect);
-        Float_t colPadSize  = fPar->GetColPadSize(iplan);
-        Float_t timeBinSize = fPar->GetDriftVelocity()
-                            / fPar->GetSamplingFrequency();
-
-        // Loop through all entries in the tree
-        for (Int_t iTrack = 0; iTrack < nTrack; iTrack++) {
-
-          gAlice->ResetHits();
-          nBytes += hitTree->GetEvent(iTrack);
-
-          // Get the number of hits in the TRD created by this particle
-          Int_t nHit = fTRD->Hits()->GetEntriesFast();
-
-          // Loop through the TRD hits  
-          for (Int_t iHit = 0; iHit < nHit; iHit++) {
-
-            if (!(hit = (AliTRDhit *) fTRD->Hits()->UncheckedAt(iHit))) 
-              continue;
-
-            Float_t pos[3];
-                    pos[0]   = hit->X();
-                    pos[1]   = hit->Y();
-                    pos[2]   = hit->Z();
-            Int_t   track    = hit->Track();
-            Int_t   detector = hit->GetDetector();
-            Int_t   plane    = geo->GetPlane(detector);
-            Int_t   sector   = geo->GetSector(detector);
-            Int_t   chamber  = geo->GetChamber(detector);        
-
-            if ((sector  != isect) ||
-                (plane   != iplan) ||
-                (chamber != icham)) 
-              continue;
-
-            // Rotate the sectors on top of each other
-            Float_t rot[3];
-            geo->Rotate(detector,pos,rot);
-
-            // Add this recPoint to the temporary array for this chamber
-            AliTRDrecPoint *recPoint = new AliTRDrecPoint("");
-            recPoint->SetLocalRow(rot[2]);
-            recPoint->SetLocalCol(rot[1]);
-            recPoint->SetLocalTime(rot[0]);
-            recPoint->SetEnergy(0);
-            recPoint->SetDetector(detector);
-            recPoint->AddDigit(track);
-            chamberArray->Add(recPoint);
-
-         }
-
-       }
-  
-        // Loop through the temporary cluster-array
-        for (Int_t iClus1 = 0; iClus1 < chamberArray->GetEntries(); iClus1++) {
-
-          AliTRDrecPoint *recPoint1 = (AliTRDrecPoint *) 
-                                      chamberArray->UncheckedAt(iClus1);
-          Float_t row1  = recPoint1->GetLocalRow();
-          Float_t col1  = recPoint1->GetLocalCol();
-          Float_t time1 = recPoint1->GetLocalTime();
-
-          if (recPoint1->GetEnergy() < 0) continue;        // Skip marked cluster  
-
-          const Int_t kNsave  = 5;
-          Int_t idxSave[kNsave];
-          Int_t iSave = 0;
-
-          const Int_t kNsaveTrack = 3;
-          Int_t tracks[kNsaveTrack];
-          tracks[0] = recPoint1->GetDigit(0);
-
-          // Check the other cluster to see, whether there are close ones
-          for (Int_t iClus2 = iClus1 + 1; iClus2 < chamberArray->GetEntries(); iClus2++) {
-
-            AliTRDrecPoint *recPoint2 = (AliTRDrecPoint *) 
-                                        chamberArray->UncheckedAt(iClus2);
-            Float_t row2 = recPoint2->GetLocalRow();
-            Float_t col2 = recPoint2->GetLocalCol();
-
-            if ((TMath::Abs(row1 - row2) < rowPadSize) ||
-                (TMath::Abs(col1 - col2) <  fRphiDist)) {
-              if (iSave == kNsave) {
-                printf("AliTRDclusterizerV0::MakeCluster -- ");
-                printf("Boundary error: iSave = %d, kNsave = %d.\n"
-                      ,iSave,kNsave);
-             }
-              else {                              
-                idxSave[iSave]  = iClus2;
-                iSave++;
-                if (iSave < kNsaveTrack) tracks[iSave] = recPoint2->GetDigit(0);
-             }
-           }
-         }
-     
-          // Merge close cluster
-          Float_t rowMerge = row1;
-          Float_t colMerge = col1;
-          if (iSave) {
-            for (Int_t iMerge = 0; iMerge < iSave; iMerge++) {
-              AliTRDrecPoint *recPoint2 =
-                (AliTRDrecPoint *) chamberArray->UncheckedAt(idxSave[iMerge]);
-              rowMerge += recPoint2->GetLocalRow();
-              colMerge += recPoint2->GetLocalCol();
-              recPoint2->SetEnergy(-1);     // Mark merged cluster
-           }
-            rowMerge /= (iSave + 1);
-            colMerge /= (iSave + 1);
-          }
-
-          Float_t smear[3];
-
-          // The position smearing in row-direction (uniform over pad width)            
-          Int_t row = (Int_t) ((rowMerge - row0) / rowPadSize);
-          smear[0]  = (row + gRandom->Rndm()) * rowPadSize + row0;
-
-          // The position smearing in rphi-direction (Gaussian)
-          smear[1] = 0;
-          do
-            smear[1] = gRandom->Gaus(colMerge,fRphiSigma);
-          while ((smear[1] < col0                        ) ||
-                 (smear[1] > col0 + nColMax * colPadSize));
-
-          // Time direction stays unchanged
-          smear[2] = time1;
-         
-         // Transform into local coordinates
-          smear[0] = (Int_t) ((smear[0] -  row0) /  rowPadSize);
-          smear[1] = (Int_t) ((smear[1] -  col0) /  colPadSize);
-          smear[2] = (Int_t) ((time0 - smear[2]) / timeBinSize);
-
-          // Add the smeared cluster to the output array 
-          Int_t   detector  = recPoint1->GetDetector();
-         Int_t   tr[9]     = { -1   };
-          Float_t pos[3];
-          Float_t sigma[2]  = {  0.0 };
-          pos[0] = smear[1];
-          pos[1] = smear[0];
-          pos[2] = (time0 - smear[2]) / timeBinSize;
-          AddCluster(pos,detector,0.0,tr,sigma,0);
-
-       }
-
-        // Clear the temporary cluster-array and delete the cluster
-        chamberArray->Delete();
-
-      }
-    }
-  }
-
-  printf("AliTRDclusterizerV0::MakeCluster -- ");
-  printf("Found %d points.\n",RecPoints()->GetEntries());
-  printf("AliTRDclusterizerV0::MakeCluster -- ");
-  printf("Fill the cluster tree.\n");
-  clusterTree->Fill();
-
-  return kTRUE;
-
-}
diff --git a/TRD/AliTRDclusterizerV0.h b/TRD/AliTRDclusterizerV0.h
deleted file mode 100644 (file)
index beb1c63..0000000
+++ /dev/null
@@ -1,42 +0,0 @@
-#ifndef ALITRDCLUSTERIZERV0_H
-#define ALITRDCLUSTERIZERV0_H
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * See cxx source for full Copyright notice                               */
-
-/* $Id$ */
-
-#include "AliTRDclusterizer.h"
-
-///////////////////////////////////////////////////////
-//  Finds and handles cluster (fast simulation)      //
-///////////////////////////////////////////////////////
-
-class AliTRDclusterizerV0 : public AliTRDclusterizer {
-
- public:
-
-  AliTRDclusterizerV0();
-  AliTRDclusterizerV0(const Text_t* name, const Text_t* title);
-  virtual ~AliTRDclusterizerV0();
-
-  virtual void    Init();
-  virtual Bool_t  MakeClusters();
-  
-  virtual void    SetRphiSigma(Float_t sigma) { fRphiSigma = sigma; };
-  virtual void    SetRphiDist(Float_t dist)   { fRphiDist  = dist;  };
-
-  virtual Float_t GetRphiSigma() const        { return fRphiSigma;  };
-  virtual Float_t GetRphiDist()  const        { return fRphiDist;   };
-
- protected:
-
-  AliTRD       *fTRD;                //! The TRD object
-
-  Float_t      fRphiSigma;           // Gaussian position smearing in rphi-direction
-  Float_t      fRphiDist;            // Maximum distance for non-overlapping cluster
-
-  ClassDef(AliTRDclusterizerV0,1)    // TRD-Cluster manager, fast simulator
-
-};
-
-#endif
diff --git a/TRD/AliTRDgeometryHole.cxx b/TRD/AliTRDgeometryHole.cxx
deleted file mode 100644 (file)
index cf58563..0000000
+++ /dev/null
@@ -1,753 +0,0 @@
-/**************************************************************************
- * 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.                  *
- **************************************************************************/
-
-/*
-$Log$
-Revision 1.12  2004/02/02 13:37:52  hristov
-Implementing of new function to check for holes (M.Ivanov)
-
-Revision 1.11  2003/09/18 09:06:07  cblume
-Geometry update, Removal of compiler warnings
-
-Revision 1.9  2002/11/21 22:38:47  alibrary
-Removing AliMC and AliMCProcess
-
-Revision 1.8  2002/10/31 17:45:35  cblume
-New chamber geometry
-
-Revision 1.7  2002/02/11 14:21:16  cblume
-Update of the geometry. Get rid of MANY
-
-Revision 1.6  2001/05/11 07:56:12  hristov
-Consistent declarations needed on Alpha
-
-Revision 1.5  2001/02/14 18:22:26  cblume
-Change in the geometry of the padplane
-
-Revision 1.4  2000/11/01 14:53:21  cblume
-Merge with TRD-develop
-
-Revision 1.1.4.4  2000/10/15 23:40:01  cblume
-Remove AliTRDconst
-
-Revision 1.1.4.3  2000/10/06 16:49:46  cblume
-Made Getters const
-
-Revision 1.1.4.2  2000/10/04 16:34:58  cblume
-Replace include files by forward declarations
-
-Revision 1.1.4.1  2000/09/22 14:43:41  cblume
-Allow the pad/timebin-dimensions to be changed after initialization
-
-Revision 1.3  2000/10/02 21:28:19  fca
-Removal of useless dependecies via forward declarations
-
-Revision 1.2  2000/06/08 18:32:58  cblume
-Make code compliant to coding conventions
-
-Revision 1.1  2000/02/28 19:01:42  cblume
-Add new TRD classes
-
-*/
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-//  TRD geometry with holes                                                  //
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
-
-#include "TVirtualMC.h"
-
-#include "AliTRDgeometryHole.h"
-
-ClassImp(AliTRDgeometryHole)
-
-//_____________________________________________________________________________
-AliTRDgeometryHole::AliTRDgeometryHole():AliTRDgeometry()
-{
-  //
-  // AliTRDgeometryHole default constructor
-  //
-
-  Init();
-
-}
-
-//_____________________________________________________________________________
-AliTRDgeometryHole::~AliTRDgeometryHole()
-{
-  //
-  // AliTRDgeometryHole destructor
-  //
-
-}
-
-Bool_t   AliTRDgeometryHole::IsHole(Int_t iplan, Int_t icham, Int_t isec) const
-{
-  // Position of Holes for PHOS (P) and RICH (R) starting at 6h
-  //                 P  P  P  -  -  R  R  R  -  -  -  -  -  -  -  -  P  P
-  //Int_t cham[18] = {1, 1, 1, 0, 0, 2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1};
-  Int_t cham[18] = {2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0 };  // did we want this?
-
-  if ((cham[17-isec]==1) && (fClengthPH[iplan][icham]<0.00001)) return kTRUE;
-  if ((cham[17-isec]==2) &&(fClengthRH[iplan][icham]<0.000001))  return kTRUE;
-  return kFALSE;
-  
-}
-
-//_____________________________________________________________________________
-void AliTRDgeometryHole::Init()
-{
-  //
-  // Initializes the geometry parameter
-  //
-
-  Int_t iplan;
-  Int_t icham;
-
-  // The outer lengths of the chambers for the sectors with holes for the PHOS
-  Float_t lengthPH[kNplan][kNcham] = { {   0.0,   0.0,   0.0, 117.0, 124.0 }
-                                    , {   0.0,   0.0,   0.0, 124.0, 131.0 }
-                                    , {   0.0,   0.0,   0.0, 131.0, 138.0 }
-                                    , {   0.0,   0.0,   0.0, 138.0, 145.0 }
-                                    , {   0.0,   0.0,   0.0, 140.0, 147.0 }
-                                     , {   0.0,   0.0,   0.0, 140.0, 147.0 } };
-
-  // The outer lengths of the chambers for the sectors with holes for the RICH
-  Float_t lengthRH[kNplan][kNcham] = { {   0.0,   0.0,   0.0,   0.0,  87.5 }
-                                    , {   0.0,   0.0,   0.0,   0.0, 101.5 }
-                                    , {   0.0,   0.0,   0.0,   0.0, 115.5 }
-                                    , {   0.0,   0.0,   0.0,   0.0, 129.5 }
-                                    , {   0.0,   0.0,   0.0,   0.0, 133.5 }
-                                     , {   0.0,   0.0,   0.0,   0.0, 133.5 } };
-
-  for (icham = 0; icham < kNcham; icham++) {
-    for (iplan = 0; iplan < kNplan; iplan++) {
-      fClengthPH[iplan][icham] = lengthPH[iplan][icham];
-      fClengthRH[iplan][icham] = lengthRH[iplan][icham];
-    }
-  }
-
-}
-
-//_____________________________________________________________________________
-void AliTRDgeometryHole::CreateGeometry(Int_t *idtmed)
-{
-  //
-  // Create the TRD geometry with holes
-  //
-  // Names of the TRD volumina (xx = detector number):
-  //
-  //      Lower part of the readout chambers (gas volume + radiator)
-  //
-  //        UAxx    Aluminum frames             (Al)
-  //        UBxx    G10 frames                  (C)
-  //        UCxx    Inner volumes               (Air)
-  //
-  //      Upper part of the readout chambers (readout plane + fee)
-  //
-  //        UDxx    G10 frames                  (C)
-  //        UExx    Inner volumes of the G10    (Air)
-  //        UFxx    Aluminum frames             (Al)
-  //        UGxx    Inner volumes of the Al     (Air)
-  //
-  //      Inner material layers
-  //
-  //        UHxx    Radiator                    (Rohacell)
-  //        UIxx    Entrance window             (Mylar)
-  //        UJxx    Drift volume                (Xe/CO2)
-  //        UKxx    Amplification volume        (Xe/CO2)
-  //        ULxx    Pad plane                   (Cu)
-  //        UMxx    Support structure           (Rohacell)
-  //        UNxx    FEE + signal lines          (Cu)
-  //        UOxx    Cooling device              (Al)
-  //        UPxx    Cooling device              (Water)
-  //
-
-  const Int_t kNdet    = kNplan * kNcham;
-
-  const Int_t kNparTrd = 4;
-  const Int_t kNparCha = 3;
-
-  Float_t xpos, ypos, zpos;
-
-  Float_t parTrd[kNparTrd];
-  Float_t parCha[kNparCha];
-
-  Char_t  cTagV[5];
-  Char_t  cTagM[5];
-
-  AliTRDgeometry::CreateGeometry(idtmed);
-
-  // The TRD mother volume for one sector (Air), full length in z-direction
-  parTrd[0] = fgkSwidth1/2.;
-  parTrd[1] = fgkSwidth2/2.;
-  parTrd[2] = fgkSlenTR1/2.;
-  parTrd[3] = fgkSheight/2.;
-  gMC->Gsvolu("UTR1","TRD1",idtmed[1302-1],parTrd,kNparTrd);
-
-  // The TRD mother volume for one sector (Air), leaving hole for PHOS
-  parTrd[0] = fgkSwidth1/2.;
-  parTrd[1] = fgkSwidth2/2.;
-  parTrd[2] = fgkSlenTR2/2.;
-  parTrd[3] = fgkSheight/2.;
-  gMC->Gsvolu("UTR2","TRD1",idtmed[1302-1],parTrd,kNparTrd);
-
-  // The TRD mother volume for one sector (Air), leaving hole for RICH
-  parTrd[0] = fgkSwidth1/2.;
-  parTrd[1] = fgkSwidth2/2.;
-  parTrd[2] = fgkSlenTR3/2.;
-  parTrd[3] = fgkSheight/2.;
-  gMC->Gsvolu("UTR3","TRD1",idtmed[1302-1],parTrd,kNparTrd);
-
-  for (Int_t icham = 0; icham < kNcham; icham++) {
-    for (Int_t iplan = 0; iplan < kNplan; iplan++) {  
-
-      Int_t iDet = GetDetectorSec(iplan,icham);
-
-      // The lower part of the readout chambers (gas volume + radiator) 
-      // The aluminum frames 
-      sprintf(cTagV,"UA%02d",iDet);
-      parCha[0] = fCwidth[iplan]/2.;
-      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
-      parCha[2] = fgkCraH/2. + fgkCdrH/2.;
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-      // The G10 frames 
-      sprintf(cTagV,"UB%02d",iDet);
-      parCha[0] = fCwidth[iplan]/2. - fgkCalT; 
-      parCha[1] = -1.;
-      parCha[2] = -1.;
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
-      // The inner part (air)
-      sprintf(cTagV,"UC%02d",iDet);
-      parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT; 
-      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCclfT;
-      parCha[2] = -1.;
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-      if (fClengthPH[iplan][icham] > 0.0) {
-        // The aluminum frames 
-        sprintf(cTagV,"UA%02d",iDet+kNdet);
-        parCha[0] = fCwidth[iplan]/2.;
-        parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.;
-        parCha[2] = fgkCraH/2. + fgkCdrH/2.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-        // The G10 frames 
-        sprintf(cTagV,"UB%02d",iDet+kNdet);
-        parCha[0] = fCwidth[iplan]/2. - fgkCalT; 
-        parCha[1] = -1.;
-        parCha[2] = -1.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
-        // The inner part (air)
-        sprintf(cTagV,"UC%02d",iDet+kNdet);
-        parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT; 
-        parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.- fgkCclfT;
-        parCha[2] = -1.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-      }
-      if (fClengthRH[iplan][icham] > 0.0) {
-        // The aluminum frames 
-        sprintf(cTagV,"UA%02d",iDet+2*kNdet);
-        parCha[0] = fCwidth[iplan]/2.;
-        parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.;
-        parCha[2] = fgkCraH/2. + fgkCdrH/2.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-        // The G10 frames 
-        sprintf(cTagV,"UB%02d",iDet+2*kNdet);
-        parCha[0] = fCwidth[iplan]/2. - fgkCalT; 
-        parCha[1] = -1.;
-        parCha[2] = -1.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
-        // The inner part (air)
-        sprintf(cTagV,"UC%02d",iDet+2*kNdet);
-        parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT; 
-        parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.- fgkCclfT;
-        parCha[2] = -1.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-      }
-
-      // The upper part of the readout chambers (readout plane + fee)
-      // The G10 frames
-      sprintf(cTagV,"UD%02d",iDet);
-      parCha[0] = fCwidth[iplan]/2. + fgkCroW;
-      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
-      parCha[2] = fgkCamH/2.;
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
-      // The inner part of the G10 frame (air)
-      sprintf(cTagV,"UE%02d",iDet);
-      parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT; 
-      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCcuT;
-      parCha[2] = -1.;
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-      // The aluminum frames
-      sprintf(cTagV,"UF%02d",iDet);
-      parCha[0] = fCwidth[iplan]/2. + fgkCroW;
-      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
-      parCha[2] = fgkCroH/2.;
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-      // The inner part of the aluminum frames
-      sprintf(cTagV,"UG%02d",iDet);
-      parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT; 
-      parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCauT;
-      parCha[2] = -1.;
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-      if (fClengthPH[iplan][icham] > 0.0) {
-        sprintf(cTagV,"UD%02d",iDet+kNdet);
-        parCha[0] = fCwidth[iplan]/2. + fgkCroW;
-        parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.;
-        parCha[2] = fgkCamH/2.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
-        // The inner part of the G10 frame (air)
-        sprintf(cTagV,"UE%02d",iDet+kNdet);
-        parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT; 
-        parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.- fgkCcuT;
-        parCha[2] = -1.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-        // The aluminum frames
-        sprintf(cTagV,"UF%02d",iDet+kNdet);
-        parCha[0] = fCwidth[iplan]/2. + fgkCroW;
-        parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.;
-        parCha[2] = fgkCroH/2.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-        // The inner part of the aluminum frames
-        sprintf(cTagV,"UG%02d",iDet+kNdet);
-        parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT; 
-        parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.- fgkCauT;
-        parCha[2] = -1.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-      }
-      if (fClengthRH[iplan][icham] > 0.0) {
-        sprintf(cTagV,"UD%02d",iDet+2*kNdet);
-        parCha[0] = fCwidth[iplan]/2. + fgkCroW;
-        parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.;
-        parCha[2] = fgkCamH/2.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
-        // The inner part of the G10 frame (air)
-        sprintf(cTagV,"UE%02d",iDet+2*kNdet);
-        parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT; 
-        parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.- fgkCcuT;
-        parCha[2] = -1.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-        // The aluminum frames
-        sprintf(cTagV,"UF%02d",iDet+2*kNdet);
-        parCha[0] = fCwidth[iplan]/2. + fgkCroW;
-        parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.;
-        parCha[2] = fgkCroH/2.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-        // The inner part of the aluminum frames
-        sprintf(cTagV,"UG%02d",iDet+2*kNdet);
-        parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT; 
-        parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.- fgkCauT;
-        parCha[2] = -1.;
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
-      }
-
-      // The material layers inside the chambers
-      parCha[0] = -1.;
-      parCha[1] = -1.;
-      // Rohacell layer (radiator)
-      parCha[2] = fgkRaThick/2;
-      sprintf(cTagV,"UH%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
-      // Mylar layer (entrance window + HV cathode) 
-      parCha[2] = fgkMyThick/2;
-      sprintf(cTagV,"UI%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha);
-      // Xe/Isobutane layer (drift volume) 
-      parCha[2] = fgkDrThick/2.;
-      sprintf(cTagV,"UJ%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
-      // Xe/Isobutane layer (amplification volume)
-      parCha[2] = fgkAmThick/2.;
-      sprintf(cTagV,"UK%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);  
-      // Cu layer (pad plane)
-      parCha[2] = fgkCuThick/2;
-      sprintf(cTagV,"UL%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
-      // G10 layer (support structure / honeycomb)
-      parCha[2] = fgkSuThick/2;
-      sprintf(cTagV,"UM%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
-      // Cu layer (FEE + signal lines)
-      parCha[2] = fgkFeThick/2;
-      sprintf(cTagV,"UN%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
-      // Al layer (cooling devices)
-      parCha[2] = fgkCoThick/2;
-      sprintf(cTagV,"UO%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-      // Water layer (cooling)
-      parCha[2] = fgkWaThick/2;
-      sprintf(cTagV,"UP%02d",iDet);
-      gMC->Gsvolu(cTagV,"BOX ",idtmed[1314-1],parCha,kNparCha);
-      if (fClengthPH[iplan][icham] > 0.0) {
-        // Rohacell layer (radiator)
-        parCha[2] = fgkRaThick/2;
-        sprintf(cTagV,"UH%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
-        // Mylar layer (entrance window + HV cathode) 
-        parCha[2] = fgkMyThick/2;
-        sprintf(cTagV,"UI%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha);
-        // Xe/Isobutane layer (drift volume) 
-        parCha[2] = fgkDrThick/2.;
-        sprintf(cTagV,"UJ%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
-        // Xe/Isobutane layer (amplification volume)
-        parCha[2] = fgkAmThick/2.;
-        sprintf(cTagV,"UK%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);  
-        // Cu layer (pad plane)
-        parCha[2] = fgkCuThick/2;
-        sprintf(cTagV,"UL%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
-        // G10 layer (support structure / honeycomb)
-        parCha[2] = fgkSuThick/2;
-        sprintf(cTagV,"UM%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
-        // Cu layer (FEE + signal lines)
-        parCha[2] = fgkFeThick/2;
-        sprintf(cTagV,"UN%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
-        // Al layer (cooling devices)
-        parCha[2] = fgkCoThick/2;
-        sprintf(cTagV,"UO%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-        // Water layer (cooling)
-        parCha[2] = fgkWaThick/2;
-        sprintf(cTagV,"UP%02d",iDet+kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1314-1],parCha,kNparCha);
-      }
-      if (fClengthRH[iplan][icham] > 0.0) {
-        // Rohacell layer (radiator)
-        parCha[2] = fgkRaThick/2;
-        sprintf(cTagV,"UH%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
-        // Mylar layer (entrance window + HV cathode) 
-        parCha[2] = fgkMyThick/2;
-        sprintf(cTagV,"UI%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha);
-        // Xe/Isobutane layer (drift volume) 
-        parCha[2] = fgkDrThick/2.;
-        sprintf(cTagV,"UJ%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
-        // Xe/Isobutane layer (amplification volume)
-        parCha[2] = fgkAmThick/2.;
-        sprintf(cTagV,"UK%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);  
-        // Cu layer (pad plane)
-        parCha[2] = fgkCuThick/2;
-        sprintf(cTagV,"UL%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
-        // G10 layer (support structure / honeycomb)
-        parCha[2] = fgkSuThick/2;
-        sprintf(cTagV,"UM%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
-        // Cu layer (FEE + signal lines)
-        parCha[2] = fgkFeThick/2;
-        sprintf(cTagV,"UN%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
-        // Al layer (cooling devices)
-        parCha[2] = fgkCoThick/2.;
-        sprintf(cTagV,"UO%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
-        // Water layer (cooling)
-        parCha[2] = fgkWaThick/2;
-        sprintf(cTagV,"UP%02d",iDet+2*kNdet);
-        gMC->Gsvolu(cTagV,"BOX ",idtmed[1314-1],parCha,kNparCha);
-      }
-
-      // Position the layers in the chambers
-      xpos = 0;
-      ypos = 0;
-      // Lower part
-      // Rohacell layer (radiator)
-      zpos = fgkRaZpos;
-      sprintf(cTagV,"UH%02d",iDet);
-      sprintf(cTagM,"UC%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // Mylar layer (entrance window + HV cathode)   
-      zpos = fgkMyZpos;
-      sprintf(cTagV,"UI%02d",iDet);
-      sprintf(cTagM,"UC%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // Xe/Isobutane layer (drift volume) 
-      zpos = fgkDrZpos;
-      sprintf(cTagV,"UJ%02d",iDet);
-      sprintf(cTagM,"UC%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // Upper part
-      // Xe/Isobutane layer (amplification volume)
-      zpos = fgkAmZpos;
-      sprintf(cTagV,"UK%02d",iDet);
-      sprintf(cTagM,"UE%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // Readout part
-      // Cu layer (pad plane)
-      zpos = fgkCuZpos; 
-      sprintf(cTagV,"UL%02d",iDet);
-      sprintf(cTagM,"UG%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // G10 layer (support structure)
-      zpos = fgkSuZpos;
-      sprintf(cTagV,"UM%02d",iDet);
-      sprintf(cTagM,"UG%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // Cu layer (FEE + signal lines)
-      zpos = fgkFeZpos; 
-      sprintf(cTagV,"UN%02d",iDet);
-      sprintf(cTagM,"UG%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // Al layer (cooling devices)
-      zpos = fgkCoZpos;
-      sprintf(cTagV,"UO%02d",iDet);
-      sprintf(cTagM,"UG%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // Water layer (cooling)
-      zpos = fgkWaZpos;
-      sprintf(cTagV,"UP%02d",iDet);
-      sprintf(cTagM,"UG%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      if (fClengthPH[iplan][icham] > 0.0) {
-        // Lower part
-        // Rohacell layer (radiator)
-        zpos = fgkRaZpos;
-        sprintf(cTagV,"UH%02d",iDet+kNdet);
-        sprintf(cTagM,"UC%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Mylar layer (entrance window + HV cathode)   
-        zpos = fgkMyZpos;
-        sprintf(cTagV,"UI%02d",iDet+kNdet);
-        sprintf(cTagM,"UC%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Xe/Isobutane layer (drift volume) 
-        zpos = fgkDrZpos;
-        sprintf(cTagV,"UJ%02d",iDet+kNdet);
-        sprintf(cTagM,"UC%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Upper part
-        // Xe/Isobutane layer (amplification volume)
-        zpos = fgkAmZpos;
-        sprintf(cTagV,"UK%02d",iDet+kNdet);
-        sprintf(cTagM,"UE%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Readout part
-        // Cu layer (pad plane)
-        zpos = fgkCuZpos; 
-        sprintf(cTagV,"UL%02d",iDet+kNdet);
-        sprintf(cTagM,"UG%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // G10 layer (support structure)
-        zpos = fgkSuZpos;
-        sprintf(cTagV,"UM%02d",iDet+kNdet);
-        sprintf(cTagM,"UG%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Cu layer (FEE + signal lines)
-        zpos = fgkFeZpos; 
-        sprintf(cTagV,"UN%02d",iDet+kNdet);
-        sprintf(cTagM,"UG%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Al layer (cooling devices)
-        zpos = fgkCoZpos;
-        sprintf(cTagV,"UO%02d",iDet+kNdet);
-        sprintf(cTagM,"UG%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Water layer (cooling)
-        zpos = fgkWaZpos;
-        sprintf(cTagV,"UP%02d",iDet+kNdet);
-        sprintf(cTagM,"UG%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      }
-      if (fClengthRH[iplan][icham] > 0.0) {
-        // Lower part
-        // Rohacell layer (radiator)
-        zpos = fgkRaZpos;
-        sprintf(cTagV,"UH%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UC%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Mylar layer (entrance window + HV cathode)   
-        zpos = fgkMyZpos;
-        sprintf(cTagV,"UI%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UC%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Xe/Isobutane layer (drift volume) 
-        zpos = fgkDrZpos;
-        sprintf(cTagV,"UJ%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UC%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Upper part
-        // Xe/Isobutane layer (amplification volume)
-        zpos = fgkAmZpos;
-        sprintf(cTagV,"UK%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UE%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Readout part
-        // Cu layer (pad plane)
-        zpos = fgkCuZpos; 
-        sprintf(cTagV,"UL%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UG%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // G10 layer (support structure)
-        zpos = fgkSuZpos;
-        sprintf(cTagV,"UM%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UG%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Cu layer (FEE + signal lines)
-        zpos = fgkFeZpos; 
-        sprintf(cTagV,"UN%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UG%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Al layer (cooling devices)
-        zpos = fgkCoZpos;
-        sprintf(cTagV,"UO%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UG%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // Water layer (cooling)
-        zpos = fgkWaZpos;
-        sprintf(cTagV,"UP%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UG%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      }
-
-      // Position the inner volumes of the chambers in the frames
-      xpos      = 0.0;
-      ypos      = 0.0;
-      zpos      = 0.0;
-      // The inside of the lower G10 frame
-      sprintf(cTagV,"UC%02d",iDet);
-      sprintf(cTagM,"UB%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // The lower G10 frame inside the aluminum frame
-      sprintf(cTagV,"UB%02d",iDet);
-      sprintf(cTagM,"UA%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // The inside of the upper G10 frame
-      sprintf(cTagV,"UE%02d",iDet);
-      sprintf(cTagM,"UD%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-      // The inside of the upper aluminum frame
-      sprintf(cTagV,"UG%02d",iDet);
-      sprintf(cTagM,"UF%02d",iDet);
-      gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");      
-      if (fClengthPH[iplan][icham] > 0.0) {
-        // The inside of the lower G10 frame
-        sprintf(cTagV,"UC%02d",iDet+kNdet);
-        sprintf(cTagM,"UB%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // The lower G10 frame inside the aluminum frame
-        sprintf(cTagV,"UB%02d",iDet+kNdet);
-        sprintf(cTagM,"UA%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // The inside of the upper G10 frame
-        sprintf(cTagV,"UE%02d",iDet+kNdet);
-        sprintf(cTagM,"UD%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // The inside of the upper aluminum frame
-        sprintf(cTagV,"UG%02d",iDet+kNdet);
-        sprintf(cTagM,"UF%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");      
-      }
-      if (fClengthRH[iplan][icham] > 0.0) {
-        // The inside of the lower G10 frame
-        sprintf(cTagV,"UC%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UB%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // The lower G10 frame inside the aluminum frame
-        sprintf(cTagV,"UB%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UA%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // The inside of the upper G10 frame
-        sprintf(cTagV,"UE%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UD%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
-        // The inside of the upper aluminum frame
-        sprintf(cTagV,"UG%02d",iDet+2*kNdet);
-        sprintf(cTagM,"UF%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");      
-      }
-
-      // Position the frames of the chambers in the TRD mother volume
-      xpos  = 0.;
-      ypos  = - fClength[iplan][0] - fClength[iplan][1] - fClength[iplan][2]/2.;
-      for (Int_t ic = 0; ic < icham; ic++) {
-        ypos += fClength[iplan][ic];        
-      }
-      ypos += fClength[iplan][icham]/2.;
-      zpos  = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace);
-      // The lower aluminum frame, radiator + drift region
-      sprintf(cTagV,"UA%02d",iDet);
-      gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
-      // The upper G10 frame, amplification region
-      sprintf(cTagV,"UD%02d",iDet);
-      zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.;
-      gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
-      // The upper aluminum frame
-      sprintf(cTagV,"UF%02d",iDet);
-      zpos += fgkCroH/2. + fgkCamH/2.;
-      gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
-      if (fClengthPH[iplan][icham] > 0.0) {
-        xpos  = 0.;
-        ypos  = - fgkSlenTR2/2.;
-        for (Int_t ic = 0; ic < icham; ic++) {
-          ypos += fClengthPH[iplan][ic];        
-        }
-        ypos += fClengthPH[iplan][icham]/2.;
-        zpos  = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace);
-        // The lower aluminum frame, radiator + drift region
-        sprintf(cTagV,"UA%02d",iDet+kNdet);
-        gMC->Gspos(cTagV,1,"UTR2",xpos,ypos,zpos,0,"ONLY");
-        // The upper G10 frame, amplification region
-        sprintf(cTagV,"UD%02d",iDet+kNdet);
-        zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.;
-        gMC->Gspos(cTagV,1,"UTR2",xpos,ypos,zpos,0,"ONLY");
-        // The upper aluminum frame
-        sprintf(cTagV,"UF%02d",iDet+kNdet);
-        zpos += fgkCroH/2. + fgkCamH/2.;
-        gMC->Gspos(cTagV,1,"UTR2",xpos,ypos,zpos,0,"ONLY");
-      }
-      if (fClengthRH[iplan][icham] > 0.0) {
-        xpos  = 0.;
-        ypos  = - fgkSlenTR3/2.;
-        for (Int_t ic = 0; ic < icham; ic++) {
-          ypos += fClengthRH[iplan][ic];        
-        }
-        ypos += fClengthRH[iplan][icham]/2.;
-        zpos  = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace);
-        // The lower aluminum frame, radiator + drift region
-        sprintf(cTagV,"UA%02d",iDet+2*kNdet);
-        gMC->Gspos(cTagV,1,"UTR3",xpos,ypos,zpos,0,"ONLY");
-        // The upper G10 frame, amplification region
-        sprintf(cTagV,"UD%02d",iDet+2*kNdet);
-        zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.;
-        gMC->Gspos(cTagV,1,"UTR3",xpos,ypos,zpos,0,"ONLY");
-        // The upper aluminum frame
-        sprintf(cTagV,"UF%02d",iDet+2*kNdet);
-        zpos += fgkCroH/2. + fgkCamH/2.;
-        gMC->Gspos(cTagV,1,"UTR3",xpos,ypos,zpos,0,"ONLY");
-      }
-
-    }
-  }
-
-  xpos = 0.;
-  ypos = 0.;
-  zpos = 0.;
-  gMC->Gspos("UTR1",1,"BTR1",xpos,ypos,zpos,0,"ONLY");
-  gMC->Gspos("UTR2",2,"BTR2",xpos,ypos,zpos,0,"ONLY");
-  gMC->Gspos("UTR3",3,"BTR3",xpos,ypos,zpos,0,"ONLY");
-
-}
-
-
diff --git a/TRD/AliTRDgeometryHole.h b/TRD/AliTRDgeometryHole.h
deleted file mode 100644 (file)
index 78d3c28..0000000
+++ /dev/null
@@ -1,41 +0,0 @@
-#ifndef ALITRDGEOMETRYHOLE_H
-#define ALITRDGEOMETRYHOLE_H
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * See cxx source for full Copyright notice                               */
-
-/* $Id$ */
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-//  TRD geometry with holes                                                  //
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
-
-#include "AliTRDgeometry.h"
-
-class AliTRDgeometryHole : public AliTRDgeometry {
-
- public:
-
-  AliTRDgeometryHole();
-  virtual ~AliTRDgeometryHole();
-
-          void    CreateGeometry(Int_t *idtmed); 
-          Int_t   IsVersion() const   { return 0; };
-          void    Init();
-
-          void    SetPHOShole()       { };
-          void    SetRICHhole()       { };
-         Bool_t   IsHole(Int_t iplan, Int_t icham, Int_t /*isec*/) const;
-  virtual void    SetOldGeometry();
-
-          Bool_t  GetPHOShole() const { return kTRUE; };
-          Bool_t  GetRICHhole() const { return kTRUE; };
-
- protected:
-
-  ClassDef(AliTRDgeometryHole,2)              // TRD geometry with hole
-
-};
-
-#endif
diff --git a/TRD/AliTRDtrackingSector.cxx b/TRD/AliTRDtrackingSector.cxx
deleted file mode 100644 (file)
index 5852d6f..0000000
+++ /dev/null
@@ -1,202 +0,0 @@
-/**************************************************************************
- * 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$ */
-
-/////////////////////////////////////////////////////////////////////////
-//                                                                     //
-//  Tracking sector                                                    //
-//                                                                     //
-/////////////////////////////////////////////////////////////////////////                       
-                                
-#include <TObject.h>
-
-#include "AliRun.h"
-
-#include "AliTRD.h" 
-#include "AliTRDgeometry.h" 
-#include "AliTRDcluster.h" 
-#include "AliTRDtimeBin.h" 
-#include "AliTRDtrackingSector.h" 
-#include "AliTRDparameter.h"
-
-ClassImp(AliTRDtrackingSector) 
-
-//_______________________________________________________
-AliTRDtrackingSector::~AliTRDtrackingSector()
-{
-  //
-  // Destructor
-  //
-
-  delete fTimeBin;
-
-}
-
-//_______________________________________________________
-AliTRDtimeBin &AliTRDtrackingSector::operator[](Int_t i)
-{
-  //
-  // Index operator 
-  //
-
-  return *(fTimeBin+i);
-
-}
-
-//_______________________________________________________
-void AliTRDtrackingSector::SetUp()
-{ 
-  //
-  // Initialization
-  //
-
-  AliTRD *trd = (AliTRD*) gAlice->GetDetector("TRD");
-  fGeom = trd->GetGeometry();
-
-  fTimeBinSize = fGeom->GetTimeBinSize();
-
-  fN = AliTRDgeometry::Nplan() * (Int_t(AliTRDgeometry::DrThick()
-                                       /fTimeBinSize) + 1);
-
-  fTimeBin = new AliTRDtimeBin[fN]; 
-
-  if (!fPar) {
-    fPar = new AliTRDparameter("TRDparameter","Standard TRD parameter");
-  }
-
-  fTimeBinSize = fPar->GetTimeBinSize();
-
-}
-
-//______________________________________________________
-Double_t AliTRDtrackingSector::GetX(Int_t tb) const
-{
-  //
-  // Get global x coordinate
-  //
-
-  if( (tb<0) || (tb>fN-1)) {
-    fprintf(stderr,"AliTRDtrackingSector::GetX: TimeBin index is out of range !\n");
-    return -99999.;
-  }
-  else { 
-    
-    Int_t tbPerPlane = fN/AliTRDgeometry::Nplan();
-    Int_t localTb = tbPerPlane - tb%tbPerPlane - 1;
-
-    Int_t plane = tb/tbPerPlane;
-    Float_t t0 = fPar->GetTime0(plane);
-    Double_t x = t0 - (localTb + 0.5) * fTimeBinSize;
-
-    return x;
-
-  }
-
-}
-
-//______________________________________________________
-Int_t AliTRDtrackingSector::GetTimeBinNumber(Double_t x) const
-{
-  //
-  // Returns the time bin number
-  //
-
-  Float_t rOut = fPar->GetTime0(AliTRDgeometry::Nplan()-1); 
-  Float_t rIn  = fPar->GetTime0(0) - AliTRDgeometry::DrThick();
-
-
-  if(x >= rOut) return fN-1;
-  if(x <= rIn)  return 0;
-
-  Int_t plane;
-  for (plane = AliTRDgeometry::Nplan()-1; plane >= 0; plane--) {
-    if(x > (fPar->GetTime0(plane) - AliTRDgeometry::DrThick())) break;
-  }  
-  Int_t tbPerPlane = fN/AliTRDgeometry::Nplan();
-  Int_t localTb = Int_t((fPar->GetTime0(plane)-x)/fTimeBinSize);
-
-  if((localTb < 0) || (localTb >= tbPerPlane)) {
-    printf("AliTRDtrackingSector::GetTimeBinNumber: \n");
-    printf("local time bin %d is out of bounds [0..%d]: x = %f \n",
-          localTb,tbPerPlane-1,x);
-    return -1;
-  }
-      
-  Int_t timeBin = (plane + 1) * tbPerPlane - 1 - localTb;
-
-  return timeBin;
-}
-
-//______________________________________________________
-Int_t AliTRDtrackingSector::GetTimeBin(Int_t det, Int_t localTb) const 
-{
-  //
-  // Time bin
-  //
-
-  Int_t plane = fGeom->GetPlane(det);
-
-  Int_t tbPerPlane = fN/AliTRDgeometry::Nplan();
-
-  Int_t timeBin = (plane + 1) * tbPerPlane - 1 - localTb;
-
-  return timeBin;
-
-}
-
-
-//______________________________________________________
-
-Bool_t AliTRDtrackingSector::TECframe(Int_t tb, Double_t y, Double_t z) const
-{
-  //  
-  // Returns <true> if point defined by <x(tb),y,z> is within 
-  // the TEC G10 frame, otherwise returns <false>  
-  //  
-
-  if((tb > (fN-1)) || (tb < 0)) return kFALSE; 
-
-  Int_t tbPerPlane = fN/AliTRDgeometry::Nplan();
-  Int_t plane = tb/tbPerPlane;
-  
-  Double_t x = GetX(tb);
-  y = TMath::Abs(y);
-
-  if((y > fGeom->GetChamberWidth(plane)/2.) &&
-     (y < x*TMath::Tan(0.5*AliTRDgeometry::GetAlpha()))) return kTRUE; 
-
-  Double_t zmin, zmax;
-  Float_t  fRowPadSize, fRow0;
-  Int_t    nPadRows;
-
-  for(Int_t iCha = 1; iCha < AliTRDgeometry::Ncham(); iCha++) {
-
-    fRow0 = fPar->GetRow0(plane,iCha-1,0);
-    fRowPadSize = fPar->GetRowPadSize(plane,iCha-1,0);
-    nPadRows = fPar->GetRowMax(plane,iCha-1,0);
-    zmin = fRow0 - fRowPadSize/2 + fRowPadSize * nPadRows;
-
-    fRow0 = fPar->GetRow0(plane,iCha,0);
-    fRowPadSize = fPar->GetRowPadSize(plane,iCha,0);
-    zmax = fRow0 - fRowPadSize/2;
-
-    if((z > zmin) && (z < zmax)) return kTRUE;     
-  }
-
-  return kFALSE;
-
-}
diff --git a/TRD/AliTRDtrackingSector.h b/TRD/AliTRDtrackingSector.h
deleted file mode 100644 (file)
index 49996d4..0000000
+++ /dev/null
@@ -1,54 +0,0 @@
-#ifndef ALITRDTRACKINGSECTOR_H
-#define ALITRDTRACKINGSECTOR_H
-
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * See cxx source for full Copyright notice                               */
-
-/* $Id: AliTRDtrackingSector.h,v */
-
-/////////////////////////////////////////////////////////////////////////
-//                                                                     //
-//  Tracking sector                                                    //
-//                                                                     //
-/////////////////////////////////////////////////////////////////////////                       
-
-#include <TObject.h>
-
-class AliTRDtimeBin;
-class AliTRDgeometry;
-class AliTRDparameter;
-
-class AliTRDtrackingSector : public TObject {
-
-// Provides tools to address clusters which lay within one sector
-
-public:
-
-  AliTRDtrackingSector() {fN=0; fTimeBin=0; fGeom=0; fTimeBinSize=0;}
-  virtual ~AliTRDtrackingSector();
-  virtual void SetUp();
-  AliTRDtimeBin& operator[](Int_t i);
-  Int_t GetNtimeBins() const { return fN; }
-  Double_t GetX(Int_t tb) const;
-  Int_t   GetTimeBinNumber(Double_t x) const;
-  Int_t   GetTimeBin(Int_t det, Int_t local_tb) const;
-  Bool_t  TECframe(Int_t tb, Double_t y, Double_t z) const;
-
-  virtual void     SetParameter(AliTRDparameter *par)      { fPar           = par; };
-  AliTRDparameter *GetParameter()                    const { return fPar;          };
-
-protected:
-
-  Int_t fN;                                // ???????
-  AliTRDgeometry          *fGeom;          // Pointer to TRD geometry
-  AliTRDtimeBin           *fTimeBin;       // Pointer to array of AliTRDtimeBin
-  Float_t                  fTimeBinSize;   // Time bin size in cm  
-  AliTRDparameter         *fPar;           // TRD parameter
-                                             
-  ClassDef(AliTRDtrackingSector,2)  // Provides tools to address clusters which lay within one sector
-
-}; 
-
-
-#endif