Removing AliITSgeom dependencies from the old ITS clusterer V2 and the corresponding...
authorcvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 2 Mar 2007 16:00:53 +0000 (16:00 +0000)
committercvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 2 Mar 2007 16:00:53 +0000 (16:00 +0000)
HLT/ITS/AliHLTITSclusterer.cxx
HLT/ITS/AliHLTITSclusterer.h
ITS/AliITSclustererV2.cxx
ITS/AliITSclustererV2.h
MONITOR/AliMonitorProcess.cxx
MONITOR/monitorGDC.cxx
RAW/AliHoughFilter.cxx

index 79ebc6b..54e08d8 100644 (file)
@@ -23,8 +23,8 @@
 //-------------------------------------------------------------------------
 
 #include "AliHLTITSclusterer.h"
+#include "AliITSgeomTGeo.h"
 #include "AliRawReader.h"
-#include "AliITSgeom.h"
 #include "AliITSRawStreamSPD.h"
 #include "AliITSRawStreamSDD.h"
 #include "AliITSRawStreamSSD.h"
@@ -33,9 +33,9 @@
 
 ClassImp(AliHLTITSclusterer)
 
-AliHLTITSclusterer::AliHLTITSclusterer(const AliITSgeom *geom):AliITSclustererV2(geom)
+AliHLTITSclusterer::AliHLTITSclusterer(const Char_t *geom):AliITSclustererV2(geom)
 {
-fNModule = geom->GetIndexMax();
+fNModule = AliITSgeomTGeo::GetNModules();
 }
 
 void AliHLTITSclusterer::Digits2Clusters(AliRawReader* rawReader,TTree *cTree)
index ba8e3fd..d3a56d7 100644 (file)
 #include "AliITSclustererV2.h"
 
 class TTree;
-class AliITSgeom;
 class AliRawReader;
 
 //-------------------------------------------------------------------------
 class AliHLTITSclusterer : public AliITSclustererV2 {
 public:
   AliHLTITSclusterer():AliITSclustererV2(){fNModule = 0;}
-  AliHLTITSclusterer(const AliITSgeom *geom);
+  AliHLTITSclusterer(const Char_t *geom);
 
   void Digits2Clusters(AliRawReader* rawReader,TTree *cTree);
 
index c572f22..72260b4 100644 (file)
@@ -19,7 +19,7 @@
 #include <TFile.h>
 #include <TTree.h>
 #include <TClonesArray.h>
-#include "AliITSgeom.h"
+#include "AliITSgeomTGeo.h"
 #include "AliITSdigitSPD.h"
 #include "AliITSdigitSDD.h"
 #include "AliITSdigitSSD.h"
@@ -56,11 +56,11 @@ fTanP(0),
 fTanN(0){
    //default constructor
  }
-AliITSclustererV2::AliITSclustererV2(const AliITSgeom *geom):
-fNModules(0),
+AliITSclustererV2::AliITSclustererV2(const Char_t *geom):
+fNModules(AliITSgeomTGeo::GetNModules()),
 fEvent(0),
 fI(0),
-fLastSPD1(0),
+fLastSPD1(AliITSgeomTGeo::GetModuleIndex(2,1,1)-1),
 fNySPD(256),
 fNzSPD(160),
 fYpitchSPD(0.0050),
@@ -75,7 +75,7 @@ fZpitchSDD(0.02940),
 fHwSDD(3.5085),
 fHlSDD(3.7632),
 fYoffSDD(0.0425),
-fLastSSD1(0),
+fLastSSD1(AliITSgeomTGeo::GetModuleIndex(6,1,1)-1),
 fYpitchSSD(0.0095),
 fHwSSD(3.65),
 fHlSSD(2.00),
@@ -84,28 +84,20 @@ fTanN(0.0075) {
   //------------------------------------------------------------
   // Standard constructor
   //------------------------------------------------------------
-  AliITSgeom *g=(AliITSgeom*)geom;
-
-  Int_t mmax=geom->GetIndexMax();
-  if (mmax>2200) {
-     Fatal("AliITSclustererV2","Too many ITS subdetectors !"); 
+  if (geom) {
+    AliWarning("\"geom\" is actually a dummy argument !");
   }
   Int_t m;
-  for (m=0; m<mmax; m++) {
-     Int_t lay,lad,det; g->GetModuleId(m,lay,lad,det);
-     Float_t x,y,z;     g->GetTrans(lay,lad,det,x,y,z); 
-     Double_t rot[9];   g->GetRotMatrix(lay,lad,det,rot);
-     Double_t alpha=TMath::ATan2(rot[1],rot[0])+TMath::Pi();
-     Double_t ca=TMath::Cos(alpha), sa=TMath::Sin(alpha);
-     fYshift[m] = x*ca + y*sa;
-     fZshift[m] = (Double_t)z;
-     fNdet[m] = (lad-1)*g->GetNdetectors(lay) + (det-1);
+  for (m=0; m<fNModules; m++) {
+     Int_t lay,lad,det; AliITSgeomTGeo::GetModuleId(m,lay,lad,det);
+     const TGeoHMatrix *tm=AliITSgeomTGeo::GetTracking2LocalMatrix(m);
+     fYshift[m] = (tm->Inverse()).GetTranslation()[1];
+     fZshift[m] = (tm->Inverse()).GetTranslation()[2];
+     fNdet[m] = (lad-1)*AliITSgeomTGeo::GetNDetectors(lay) + (det-1);
      fNlayer[m] = lay-1;
   }
-  fNModules = g->GetIndexMax();
 
   //SPD geometry  
-  fLastSPD1=g->GetModuleIndex(2,1,1)-1;
   fYSPD[0]=0.5*fYpitchSPD;
   for (m=1; m<fNySPD; m++) fYSPD[m]=fYSPD[m-1]+fYpitchSPD; 
   fZSPD[0]=fZ1pitchSPD;
@@ -122,9 +114,6 @@ fTanN(0.0075) {
     fZSPD[m]-=dz;
   }
 
-  //SSD geometry
-  fLastSSD1=g->GetModuleIndex(6,1,1)-1;
-
 }
 
 
index 028403f..7c54b16 100644 (file)
@@ -15,7 +15,6 @@ class TFile;
 class TTree;
 class TClonesArray;
 
-class AliITSgeom;
 class AliITSclusterV2;
 class AliRawReader;
 class AliITSRawStream;
@@ -23,7 +22,7 @@ class AliITSRawStream;
 class AliITSclustererV2 : public TObject {
 public:
   AliITSclustererV2();
-  AliITSclustererV2(const AliITSgeom *geom);
+  AliITSclustererV2(const Char_t *geom);
 
   void SetEvent(Int_t event) { fEvent=event; }
   Int_t Digits2Clusters(TTree *in, TTree *out);
index 18121b8..779e067 100644 (file)
@@ -647,7 +647,7 @@ Bool_t AliMonitorProcess::ReconstructITS(AliRawReader* rawReader, AliESD* esd)
   // cluster finder
   AliDebug(1, "reconstructing clusters...");
   itsLoader->LoadRecPoints("recreate");
-  AliITSclustererV2 clusterer(fITSgeom);
+  AliITSclustererV2 clusterer(0);
   itsLoader->MakeRecPointsContainer();
   clusterer.Digits2Clusters(rawReader);
 
index c5e0c50..c813b7c 100644 (file)
@@ -170,7 +170,7 @@ int main(int argc, char** argv)
 
        // ITS clusterer and vertexer
        fBenchmark->Start("ITS Clusterer");
-       AliHLTITSclusterer clusterer(geom);
+       AliHLTITSclusterer clusterer(0);
        AliRawReader *itsrawreader=new AliRawReaderDate(ptr);
        TTree* treeClusters = new TTree("TreeL3ITSclusters"," "); //make a tree
        clusterer.Digits2Clusters(itsrawreader,treeClusters);
index 78c640e..5c1df45 100644 (file)
@@ -113,9 +113,7 @@ void AliHoughFilter::RunITSclusterer(AliRawEvent* event, TTree *treeClusters)
   TStopwatch timer;
   timer.Start();
 
-  if(!fITSgeom)
-    AliError("ITS geometry not created!");
-  AliHLTITSclusterer clusterer(fITSgeom);
+  AliHLTITSclusterer clusterer(0);
   AliRawReader *itsrawreader=new AliRawReaderRoot(event);
   clusterer.Digits2Clusters(itsrawreader,treeClusters);
   delete itsrawreader;