Transition to NewIO
[u/mrichter/AliRoot.git] / PMD / AliPMDv0.cxx
index 40f3968..e5925e2 100644 (file)
  * about the suitability of this software for any purpose. It is          *
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
-/*
-$Log$
-Revision 1.15  2002/10/23 07:36:35  alibrary
-Introducing Riostream.h
-
-Revision 1.14  2001/05/21 17:44:04  hristov
-Backslash to continue strings
-
-Revision 1.13  2001/05/21 10:59:09  morsch
-Printouts in debug mode only.
-
-Revision 1.12  2001/05/21 09:39:28  morsch
-Minor modifications on the geometry. (Tapan Nayak)
 
-Revision 1.11  2001/05/14 14:01:04  morsch
-AliPMDv0 coarse geometry and AliPMDv1 detailed simulation, completely revised versions by Tapan Nayak.
-*/
+/* $Id$ */
 
 //
 ///////////////////////////////////////////////////////////////////////////////
@@ -44,11 +29,14 @@ AliPMDv0 coarse geometry and AliPMDv1 detailed simulation, completely revised ve
 ///////////////////////////////////////////////////////////////////////////////
 ////
 
-#include "AliPMDv0.h"
-#include "AliRun.h"
+#include "Riostream.h"
+
+#include <TVirtualMC.h>
+
 #include "AliConst.h" 
 #include "AliMagF.h" 
-#include "Riostream.h"
+#include "AliPMDv0.h"
+#include "AliRun.h"
  
 static Int_t   kdet, ncell_sm, ncell_hole;
 static Float_t zdist, zdist1;
@@ -218,7 +206,7 @@ void AliPMDv0::CreateSupermodule()
 
   // --- place  inner hex inside outer hex 
 
-  gMC->Gsposp("ECAR", 1, "ECCU", 0., 0., 0., 0, "ONLY", hexd2, 10);
+  gMC->Gspos("ECAR", 1, "ECCU", 0., 0., 0., 0, "ONLY");
 
 // Rhombus shaped supermodules (defined by PARA) 
 
@@ -272,7 +260,7 @@ void AliPMDv0::CreateSupermodule()
     }
     for (i = 1; i <= ncell_sm; ++i) {
       number = i+(j-1)*ncell_sm;
-      gMC->Gsposp("ECCU", number, "EHC1", xb,yb,zb, ihrotm, "ONLY", hexd1,10);
+      gMC->Gspos("ECCU", number, "EHC1", xb,yb,zb, ihrotm, "ONLY");
       xb += (hexd1[6]*2.);
     }
     xrow = xrow+1;
@@ -351,7 +339,7 @@ void AliPMDv0::CreateSupermodule()
     }
     for (i = 1; i <= ncell_sm; ++i) {
       number = i+(j-1)*ncell_sm;
-         gMC->Gsposp("ECCU", number, "EHC2", xb,yb,zb, ihrotm, "ONLY", hexd1,10);
+           gMC->Gspos("ECCU", number, "EHC2", xb,yb,zb, ihrotm, "ONLY");
       xb += (hexd1[6]*2.);
     }
     xrow = xrow+1;
@@ -430,7 +418,7 @@ void AliPMDv0::CreateSupermodule()
     }
     for (i = 1; i <= (ncell_sm - ncell_hole); ++i) {
       number = i+(j-1)*(ncell_sm - ncell_hole);
-         gMC->Gsposp("ECCU", number, "EHC3", xb,yb,zb, ihrotm, "ONLY", hexd1,10);
+      gMC->Gspos("ECCU", number, "EHC3", xb,yb,zb, ihrotm, "ONLY");
       xb += (hexd1[6]*2.);
     }
     xrow = xrow+1;
@@ -722,35 +710,33 @@ void AliPMDv0::CreatePMD()
   Int_t irotate[3] = {0, jhrot12, jhrot13};
 
   num_mod=0;
-  for (j=0; j<3; ++j)
-    {
-      gMC->Gsposp("EALM", j+1, "EPMD", xalm[j],yalm[j], 0., irotate[j], "ONLY", Al_rod, 3);
-      x2=xemm2*TMath::Cos(theta[j]) - yemm2*TMath::Sin(theta[j]);
-      y2=xemm2*TMath::Sin(theta[j]) + yemm2*TMath::Cos(theta[j]);
+  for (j=0; j<3; ++j) {
+     gMC->Gspos("EALM", j+1, "EPMD", xalm[j],yalm[j], 0., irotate[j], "ONLY");
+     x2=xemm2*TMath::Cos(theta[j]) - yemm2*TMath::Sin(theta[j]);
+     y2=xemm2*TMath::Sin(theta[j]) + yemm2*TMath::Cos(theta[j]);
 
-      gMC->Gsposp("EMM2", j+1, "EPMD", x2,y2, 0., irotate[j], "ONLY", dpara_emm2, 6);
+     gMC->Gspos("EMM2", j+1, "EPMD", x2,y2, 0., irotate[j], "ONLY");
 
-      x3=xemm3*TMath::Cos(theta[j]) - yemm3*TMath::Sin(theta[j]);
-      y3=xemm3*TMath::Sin(theta[j]) + yemm3*TMath::Cos(theta[j]);
+     x3=xemm3*TMath::Cos(theta[j]) - yemm3*TMath::Sin(theta[j]);
+     y3=xemm3*TMath::Sin(theta[j]) + yemm3*TMath::Cos(theta[j]);
 
-      gMC->Gsposp("EMM3", j+4, "EPMD", x3,y3, 0., irotate[j], "ONLY", dpara_emm3, 6);
+     gMC->Gspos("EMM3", j+4, "EPMD", x3,y3, 0., irotate[j], "ONLY");
 
-      for (i=1; i<9; ++i)
-       {
-         xpos[i]=xmod[j] + xsup[i]*TMath::Cos(theta[j]) - ysup[i]*TMath::Sin(theta[j]);
-         ypos[i]=ymod[j] + xsup[i]*TMath::Sin(theta[j]) + ysup[i]*TMath::Cos(theta[j]);
-         if(fDebug) 
-             printf("%s: %f %f \n", ClassName(), xpos[i], ypos[i]);
+     for (i=1; i<9; ++i) {
+             xpos[i]=xmod[j] + xsup[i]*TMath::Cos(theta[j]) - ysup[i]*TMath::Sin(theta[j]);
+             ypos[i]=ymod[j] + xsup[i]*TMath::Sin(theta[j]) + ysup[i]*TMath::Cos(theta[j]);
+             if(fDebug) 
+                printf("%s: %f %f \n", ClassName(), xpos[i], ypos[i]);
 
-         num_mod = num_mod+1;
+             num_mod = num_mod+1;
 
-         if(fDebug) 
-             printf("\n%s: Num_mod %d\n",ClassName(),num_mod);
+             if(fDebug) 
+                printf("\n%s: Num_mod %d\n",ClassName(),num_mod);
 
-         gMC->Gsposp("EMM1", num_mod + 6, "EPMD", xpos[i],ypos[i], 0., irotate[j], "ONLY", dpara_emm1, 6);
+             gMC->Gspos("EMM1", num_mod + 6, "EPMD", xpos[i],ypos[i], 0., irotate[j], "ONLY");
 
-       }
-    }
+          }
+  }
 
        
   // place EHOL in the centre of EPMD
@@ -1058,3 +1044,4 @@ void AliPMDv0::GetParameters()
 
 
 
+