by default only barrel tracks for resolution task

[u/mrichter/AliRoot.git] / JETAN / AliJetDummyGeo.cxx

diff --git a/JETAN/AliJetDummyGeo.cxx b/JETAN/AliJetDummyGeo.cxx
index 02994e474f0486679eaff3363a117fb588ad1c01..052cb2240064cd04a7b690273767f61d8d24a0af 100644 (file)
--- a/JETAN/AliJetDummyGeo.cxx
+++ b/JETAN/AliJetDummyGeo.cxx
@@ -18,6 +18,7 @@
 // necessary for the jet finder
 //
 // M. Estienne
+// Magali.Estienne@cern.ch
 //
 
 #include <Riostream.h>
@@ -169,6 +170,7 @@ AliJetDummyGeo::AliJetDummyGeo():
 }
 
 AliJetDummyGeo::AliJetDummyGeo(const AliJetDummyGeo& geom):
+  TObject(geom),  
   fArm1EtaMin(geom.fArm1EtaMin), 
   fArm1EtaMax(geom.fArm1EtaMax), 
   fArm1PhiMin(geom.fArm1PhiMin), 
@@ -274,7 +276,7 @@ AliJetDummyGeo::AliJetDummyGeo(const AliJetDummyGeo& geom):
 AliJetDummyGeo::~AliJetDummyGeo()
 {
   // Destructor
-  delete [] fMatrixOfSM;
+  // delete [] fMatrixOfSM;
 }
 
 //------------------------------------------------------------------------------------
@@ -603,7 +605,8 @@ void AliJetDummyGeo::CreateListOfTrd1Modules()
   // Generate the list of Trd1 modules
   // which will make up the EMCAL
   // geometry
-
+    printf("CreateListOfTrd1Modules() \n");
+    
   AliJetDummyShishKebabTrd1Module *mod=0, *mTmp=0; // current module
   if(fShishKebabTrd1Modules == 0) {
     fShishKebabTrd1Modules = new TList;
@@ -641,6 +644,9 @@ void AliJetDummyGeo::CreateListOfTrd1Modules()
   Int_t ind=0; // this is phi index
   Int_t ieta=0, nModule=0, iphiTemp;
   Double_t xr, zr, theta, phi, eta, r, x,y;
+  xr = 0.;
+  zr = 0.;
+  
   TVector3 vglob;
   Double_t ytCenterModule=0.0, ytCenterCell=0.0;
 
@@ -682,10 +688,10 @@ void AliJetDummyGeo::CreateListOfTrd1Modules()
         trd1->GetCenterOfCellInLocalCoordinateofSM(ic, xr, zr);      // case of 2X2
         GetCellPhiEtaIndexInSModule(0, nModule, 0, ic, iphiTemp, ieta); 
       } if(fNPHIdiv==3) {
-        trd1->GetCenterOfCellInLocalCoordinateofSM_3X3(ic, xr, zr);  // case of 3X3
+        trd1->GetCenterOfCellInLocalCoordinateofSM3X3(ic, xr, zr);  // case of 3X3
         GetCellPhiEtaIndexInSModule(0, nModule, 0, ic, iphiTemp, ieta); 
       } if(fNPHIdiv==1) {
-        trd1->GetCenterOfCellInLocalCoordinateofSM_1X1(xr, zr);      // case of 1X1
+        trd1->GetCenterOfCellInLocalCoordinateofSM1X1(xr, zr);      // case of 1X1
         GetCellPhiEtaIndexInSModule(0, nModule, 0, ic, iphiTemp, ieta); 
       }
       fCentersOfCellsXDir.AddAt(Float_t(xr) - fParSM[0],ieta);