]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ZDC/AliZDCv3.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
PDG code of incoming particles stored in ZDC hits
[u/mrichter/AliRoot.git] / ZDC / AliZDCv3.cxx
diff --git a/ZDC/AliZDCv3.cxx b/ZDC/AliZDCv3.cxx
index d3f527d7d0c0155775cc4bca6dc40f232f2bfb02..aaca85f05cc3943b3d55c17f48d800ecb3c110da 100644 (file)
--- a/ZDC/AliZDCv3.cxx
+++ b/ZDC/AliZDCv3.cxx
@@ -33,6 +33,7 @@
 #include <TTree.h>
 #include <TVirtualMC.h>
 #include <TGeoManager.h>
+#include <TParticle.h>
 
 // --- AliRoot classes
 #include "AliConst.h"
@@ -58,7 +59,7 @@ AliZDCv3::AliZDCv3() :
   fMedSensZEM(0),
   fMedSensGR(0),
   fMedSensPI(0),
-  fMedSensTDI(0), 
+  fMedSensTDI(0),
   fNalfan(0),
   fNalfap(0),
   fNben(0),  
@@ -68,11 +69,14 @@ AliZDCv3::AliZDCv3() :
   fpLostD1C(0), 
   fpDetectedC(0),
   fnDetectedC(0),
+  fnLumiC(0),
   fpLostITA(0), 
   fpLostD1A(0), 
   fpLostTDI(0), 
   fpDetectedA(0),
-  fnDetectedA(0)
+  fnDetectedA(0),
+  fnLumiA(0),
+  fnTrou(0)
 {
   //
   // Default constructor for Zero Degree Calorimeter
@@ -81,7 +85,7 @@ AliZDCv3::AliZDCv3() :
 }
  
 //_____________________________________________________________________________
-AliZDCv3::AliZDCv3(const char *name, const char *title): 
+AliZDCv3::AliZDCv3(const char *name, const char *title) : 
   AliZDC(name,title),
   fMedSensF1(0),
   fMedSensF2(0),
@@ -90,21 +94,24 @@ AliZDCv3::AliZDCv3(const char *name, const char *title):
   fMedSensZEM(0),
   fMedSensGR(0),
   fMedSensPI(0),
-  fMedSensTDI(0), 
+  fMedSensTDI(0),
   fNalfan(90),
   fNalfap(90),
   fNben(18),  
   fNbep(28), 
+  fZEMLength(0),
   fpLostITC(0), 
   fpLostD1C(0), 
   fpDetectedC(0),
   fnDetectedC(0),
+  fnLumiC(0),
   fpLostITA(0), 
   fpLostD1A(0), 
   fpLostTDI(0), 
   fpDetectedA(0),
-  fnDetectedA(0)
-  
+  fnDetectedA(0),  
+  fnLumiA(0),
+  fnTrou(0)
 {
   //
   // Standard constructor for Zero Degree Calorimeter 
@@ -147,16 +154,16 @@ AliZDCv3::AliZDCv3(const char *name, const char *title):
   fDimZP[2] = 75.;    
   fPosZNC[0] = 0.;
   fPosZNC[1] = 1.2;
-  fPosZNC[2] = -11650.; 
+  fPosZNC[2] = -11600.; 
   fPosZPC[0] = 23.9;
-  fPosZPA[1] = 0.;
-  fPosZPA[2] = -11600.; 
+  fPosZPC[1] = 0.;
+  fPosZPC[2] = -11600.; 
   fPosZNA[0] = 0.;
   fPosZNA[1] = 1.2;
   fPosZNA[2] = 11620.; 
-  fPosZPC[0] = 24.;
-  fPosZPC[1] = 0.;
-  fPosZPC[2] = 11620.; 
+  fPosZPA[0] = 24.;
+  fPosZPA[1] = 0.;
+  fPosZPA[2] = 11620.; 
   fFibZN[0] = 0.;
   fFibZN[1] = 0.01825;
   fFibZN[2] = 50.;
@@ -220,7 +227,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[6] = -zd1;
   conpar[7] = 0.;
   conpar[8] = 55.;
-  gMC->Gsvolu("ZDCC", "PCON", idtmed[11], conpar, 9);
+  gMC->Gsvolu("ZDCC", "PCON", idtmed[10], conpar, 9);
   gMC->Gspos("ZDCC", 1, "ALIC", 0., 0., 0., 0, "ONLY");
   
 
@@ -467,11 +474,20 @@ void AliZDCv3::CreateBeamLine()
   gMC->Matrix(im1, 90.+0.143, 0., 90., 90., 0.143, 0.); // x<0  
   gMC->Gspos("QT17", 1, "ZDCC", TMath::Sin(angle) * 680.8/ 2. - 9.4, 
              0., -tubpar[2]-zd1, im1, "ONLY"); 
+  //printf("\n QT17-18 pipe from z = %f to z= %f\n",-zd1,-2*tubpar[2]-zd1);
             
   //AliMatrix(im2, 90.-0.143, 0., 90., 90., 0.143, 180.); // x>0 (ZP)
   gMC->Matrix(im2, 90.-0.143, 0., 90., 90., 0.143, 180.); // x>0 (ZP)  
   gMC->Gspos("QT18", 1, "ZDCC", 9.7 - TMath::Sin(angle) * 680.8 / 2., 
              0., -tubpar[2]-zd1, im2, "ONLY"); 
+          
+  // -- Luminometer (Cu box) in front of ZN - side C
+  boxpar[0] = 8.0/2.;
+  boxpar[1] = 8.0/2.;
+  boxpar[2] = 15./2.;
+  gMC->Gsvolu("QLUC", "BOX ", idtmed[6], boxpar, 3);
+  gMC->Gspos("QLUC", 1, "ZDCC", 0., 0.,  fPosZNC[2]+66.+boxpar[2], 0, "ONLY");
+  //printf("\n QLUC LUMINOMETER from z = %f to z= %f\n",  fPosZNC[2]+66., fPosZNC[2]+66.+2*boxpar[2]);
                 
   // --  END OF BEAM PIPE VOLUME DEFINITION FOR SIDE C (RB26 SIDE) 
   // ----------------------------------------------------------------
@@ -513,7 +529,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 6.0/2.;
   tubpar[1] = 6.4/2.;
   tubpar[2] = 386.5/2.;
-  gMC->Gsvolu("QA01", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA01", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA01", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("\n QA01 TUBE from z = %f to z= %f (Inner triplet beg.)\n",zd2,2*tubpar[2]+zd2);
@@ -525,7 +541,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 6.3/2.;
   tubpar[1] = 6.7/2.;
   tubpar[2] = 3541.8/2.;
-  gMC->Gsvolu("QA02", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA02", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA02", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("\n QA02 TUBE from z = %f to z= %f (D1 beg.)\n",zd2,2*tubpar[2]+zd2);
@@ -547,7 +563,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 6.75/2.;
   tubpar[1] = 7.15/2.;
   tubpar[2] = (945.0+80.1)/2.;
-  gMC->Gsvolu("QA03", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA03", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA03", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("\n QA03 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -560,7 +576,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 7.15/2.;
   conpar[3] = 8.0/2.;
   conpar[4] = 8.4/2.;
-  gMC->Gsvolu("QA04", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QA04", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QA04", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   //printf("   QA04 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -569,7 +585,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 8.0/2.;
   tubpar[1] = 8.4/2.;
   tubpar[2] = 43.9/2.;
-  gMC->Gsvolu("QA05", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA05", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA05", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("\n QA05 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -580,7 +596,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 8.0/2.;
   tubpar[1] = 8.4/2.;
   tubpar[2] = 20./2.;
-  gMC->Gsvolu("QA06", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA06", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA06", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA06 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -591,7 +607,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 8.0/2.;
   tubpar[1] = 8.4/2.;
   tubpar[2] = 28.5/2.;
-  gMC->Gsvolu("QA07", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA07", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA07", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA07 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -602,7 +618,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 8.0/2.;
   tubpar[1] = 8.4/2.;
   tubpar[2] = 28.5/2.;
-  gMC->Gsvolu("QAV1", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QAV1", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QAV1", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QAV1 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -615,7 +631,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 8.4/2.;
   conpar[3] = 9.8/2.;
   conpar[4] = 10.2/2.;
-  gMC->Gsvolu("QAV2", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QAV2", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QAV2", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   //printf("   QAV2 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -627,7 +643,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 10.2/2.;
   conpar[3] = 9.0/2.;
   conpar[4] = 9.4/2.;
-  gMC->Gsvolu("QAV3", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QAV3", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QAV3", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   //printf("   QAV3 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -637,7 +653,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 9.0/2.;
   tubpar[1] = 9.4/2.;
   tubpar[2] = 31.0/2.;
-  gMC->Gsvolu("QAV4", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QAV4", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QAV4", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QAV4 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -652,7 +668,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 13.0/2.;
   conpar[3] = 9.6/2.;
   conpar[4] = 13.0/2.;
-  gMC->Gsvolu("Q01T", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("Q01T", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("Q01T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   //printf("   Q01T CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -662,7 +678,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 9.6/2.;
   tubpar[1] = 10.0/2.;
   tubpar[2] = 1.0/2.;
-  gMC->Gsvolu("Q02T", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("Q02T", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("Q02T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   Q02T TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -675,7 +691,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 10.0/2.;
   conpar[3] = 13.8/2.;
   conpar[4] = 14.2/2.;
-  gMC->Gsvolu("Q03T", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("Q03T", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("Q03T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   //printf("   Q03T CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -685,7 +701,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 13.8/2.;
   tubpar[1] = 14.2/2.;
   tubpar[2] = 38.6/2.;
-  gMC->Gsvolu("Q04T", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("Q04T", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("Q04T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   Q04T TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -696,7 +712,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.0/2.;
   tubpar[1] = 21.4/2.;
   tubpar[2] = 100.12/2.;
-  gMC->Gsvolu("Q05T", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("Q05T", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("Q05T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   Q05T TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -707,7 +723,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 13.8/2.;
   tubpar[1] = 14.2/2.;
   tubpar[2] = 38.6/2.;
-  gMC->Gsvolu("Q06T", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("Q06T", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("Q06T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   Q06T TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -720,7 +736,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 14.2/2.;
   conpar[3] = 18.0/2.;
   conpar[4] = 18.4/2.;
-  gMC->Gsvolu("Q07T", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("Q07T", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("Q07T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   //printf("   Q07T CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -735,7 +751,7 @@ void AliZDCv3::CreateBeamLine()
   tubspar[2] = 100.0/2.;
   tubspar[3] = 0.;
   tubspar[4] = 180.;  
-  gMC->Gsvolu("Q08T", "TUBS", idtmed[6], tubspar, 5);
+  gMC->Gsvolu("Q08T", "TUBS", idtmed[7], tubspar, 5);
   // Ch.debug
   //printf("\n upper part : one single phi segment of a tube (Q08T)\n");  
   
@@ -753,7 +769,7 @@ void AliZDCv3::CreateBeamLine()
   tubspar[2] = 100.0/2.;
   tubspar[3] = 180.;
   tubspar[4] = 360.;  
-  gMC->Gsvolu("Q10T", "TUBS", idtmed[6], tubspar, 5);
+  gMC->Gsvolu("Q10T", "TUBS", idtmed[7], tubspar, 5);
   // rectangular beam pipe inside TCDD lower section (Vacuum)  
   boxpar[0] = 7.0/2.;
   boxpar[1] = 2.5/2.;
@@ -770,7 +786,7 @@ void AliZDCv3::CreateBeamLine()
   boxpar[0] = 0.2/2.;
   boxpar[1] = 5.0/2.;
   boxpar[2] = 100./2.;
-  gMC->Gsvolu("Q12T", "BOX ", idtmed[6], boxpar, 3);  
+  gMC->Gsvolu("Q12T", "BOX ", idtmed[7], boxpar, 3);  
   // positioning RF screen at both sides of TCDD
   gMC->Gspos("Q12T", 1, "ZDCA", tubspar[1]+boxpar[0], 0., -100+zd2, 0, "ONLY");  
   gMC->Gspos("Q12T", 2, "ZDCA", -tubspar[1]-boxpar[0], 0., -100+zd2, 0, "ONLY");      
@@ -780,7 +796,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 8.0/2.;
   tubpar[1] = 8.4/2.;
   tubpar[2] = 30.0/2.;
-  gMC->Gsvolu("QA08", "TUBE", idtmed[8], tubpar, 3);
+  gMC->Gsvolu("QA08", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA08", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA08 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -791,7 +807,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 8.0/2.;
   tubpar[1] = 8.4/2.;
   tubpar[2] = 2.0/2.;
-  gMC->Gsvolu("QA09", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA09", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA09", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA09 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -804,7 +820,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 8.4/2.;
   conpar[3] = 21.2/2.;
   conpar[4] = 21.8/2.;
-  gMC->Gsvolu("QA10", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QA10", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QA10", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   //printf("   QA10 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -814,7 +830,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 403.54/2.;
-  gMC->Gsvolu("QA11", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA11", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA11", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA11 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -825,7 +841,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 40.0/2.;
-  gMC->Gsvolu("QA12", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA12", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA12", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA12 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -836,7 +852,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 30.0/2.;
-  gMC->Gsvolu("QA13", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA13", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA13", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA13 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -847,7 +863,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 40.0/2.;
-  gMC->Gsvolu("QA14", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA14", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA14", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA14 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -858,7 +874,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 20.0/2.;
-  gMC->Gsvolu("QA15", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA15", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA15", 1, "ZDCA", 1.10446, 0., tubpar[2]+zd2, irotpipe2, "ONLY");
   // Ch.debug
   //printf("   QA15 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -869,7 +885,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 54.0/2.;
   tubpar[1] = 54.6/2.;
   tubpar[2] = 540.0/2.;
-  gMC->Gsvolu("Q13T", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("Q13T", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("Q13T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   Q13T TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -880,17 +896,17 @@ void AliZDCv3::CreateBeamLine()
   boxpar[0] = 11.0/2.;
   boxpar[1] = 9.0/2.;
   boxpar[2] = 540.0/2.;
-  gMC->Gsvolu("QTD1", "BOX ", idtmed[6], boxpar, 3);
+  gMC->Gsvolu("QTD1", "BOX ", idtmed[7], boxpar, 3);
   gMC->Gspos("QTD1", 1, "Q13T", -3.8, 10.5,  0., 0, "ONLY");
   boxpar[0] = 11.0/2.;
   boxpar[1] = 9.0/2.;
   boxpar[2] = 540.0/2.;
-  gMC->Gsvolu("QTD2", "BOX ", idtmed[6], boxpar, 3);
+  gMC->Gsvolu("QTD2", "BOX ", idtmed[7], boxpar, 3);
   gMC->Gspos("QTD2", 1, "Q13T", -3.8, -10.5,  0., 0, "ONLY");  
   boxpar[0] = 5.1/2.;
   boxpar[1] = 0.2/2.;
   boxpar[2] = 540.0/2.;
-  gMC->Gsvolu("QTD3", "BOX ", idtmed[6], boxpar, 3);
+  gMC->Gsvolu("QTD3", "BOX ", idtmed[7], boxpar, 3);
   gMC->Gspos("QTD3", 1, "Q13T", -3.8+5.5+boxpar[0], 6.1,  0., 0, "ONLY");  
   gMC->Gspos("QTD3", 2, "Q13T", -3.8+5.5+boxpar[0], -6.1,  0., 0, "ONLY"); 
   gMC->Gspos("QTD3", 3, "Q13T", -3.8-5.5-boxpar[0], 6.1,  0., 0, "ONLY");  
@@ -901,14 +917,14 @@ void AliZDCv3::CreateBeamLine()
   tubspar[2] = 540.0/2.;
   tubspar[3] = 90.;
   tubspar[4] = 270.;  
-  gMC->Gsvolu("QTD4", "TUBS", idtmed[6], tubspar, 5);
+  gMC->Gsvolu("QTD4", "TUBS", idtmed[7], tubspar, 5);
   gMC->Gspos("QTD4", 1, "Q13T", -3.8-10.6, 0.,  0., 0, "ONLY");
   tubspar[0] = 12.0/2.;
   tubspar[1] = 12.4/2.;
   tubspar[2] = 540.0/2.;
   tubspar[3] = -90.;
   tubspar[4] = 90.;  
-  gMC->Gsvolu("QTD5", "TUBS", idtmed[6], tubspar, 5);
+  gMC->Gsvolu("QTD5", "TUBS", idtmed[7], tubspar, 5);
   gMC->Gspos("QTD5", 1, "Q13T", -3.8+10.6, 0.,  0., 0, "ONLY"); 
   //---------------- END DEFINING TDI INSIDE Q13T -------------------------------
   
@@ -916,7 +932,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 20.0/2.;
-  gMC->Gsvolu("QA16", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA16", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA16", 1, "ZDCA", 1.10446+2.2, 0., tubpar[2]+zd2, irotpipe2, "ONLY");
   // Ch.debug
   //printf("   QA16 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -927,7 +943,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 40.0/2.;
-  gMC->Gsvolu("QA17", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA17", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA17", 1, "ZDCA", 4.4, 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA17 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -938,7 +954,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 30.0/2.;
-  gMC->Gsvolu("QA18", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA18", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA18", 1, "ZDCA", 4.4, 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA18 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -949,7 +965,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 40.0/2.;
-  gMC->Gsvolu("QA19", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA19", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA19", 1, "ZDCA", 4.4, 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA19 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -960,7 +976,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 50.0/2.;
-  gMC->Gsvolu("QA20", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA20", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA20", 1, "ZDCA", 4.4, 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA20 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -971,7 +987,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 21.2/2.;
   tubpar[1] = 21.8/2.;
   tubpar[2] = 120.0/2.;
-  gMC->Gsvolu("QA21", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA21", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA21", 1, "ZDCA", 4.4, 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA21 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -984,7 +1000,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 21.8/2.;
   conpar[3] = 79.7/2.;
   conpar[4] = 81.3/2.;
-  gMC->Gsvolu("QA22", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QA22", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QA22", 1, "ZDCA", 4.4-2.201, 0., conpar[0]+zd2, irotpipe1, "ONLY");
   //printf("   QA22 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -994,7 +1010,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 79.7/2.;
   tubpar[1] = 81.3/2.;
   tubpar[2] = 2393.05/2.;
-  gMC->Gsvolu("QA23", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA23", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA23", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA23 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1015,7 +1031,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 81.3/2.; // thickness 8 mm  
   conpar[3] = 74.82868/2.;
   conpar[4] = 76.42868/2.; // thickness 8 mm 
-  gMC->Gsvolu("Q24A", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("Q24A", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("Q24A", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   //printf("   Q24A CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -1027,7 +1043,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 75.42868/2.; // thickness 3 mm  
   conpar[3] = 23.19588/2.;
   conpar[4] = 23.79588/2.; // thickness 3 mm 
-  gMC->Gsvolu("QA25", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QA25", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QA25", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");  
   //printf("   QA25 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -1039,7 +1055,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 24.79588/2.;// thickness 8 mm 
   conpar[3] = 19.6/2.;
   conpar[4] = 21.2/2.;// thickness 8 mm 
-  gMC->Gsvolu("QA26", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QA26", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QA26", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");  
   //printf("   QA26 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
 
@@ -1049,7 +1065,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 19.6/2.;
   tubpar[1] = 21.2/2.;
   tubpar[2] = 9.55/2.;
-  gMC->Gsvolu("QA27", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA27", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA27", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA27 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1060,7 +1076,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 19.6/2.;
   tubpar[1] = 25.3/2.;
   tubpar[2] = 4.9/2.;
-  gMC->Gsvolu("QF01", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QF01", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QF01", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QF01  TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1071,7 +1087,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 20.2/2.;
   tubpar[1] = 20.6/2.;
   tubpar[2] = 2.15/2.;
-  gMC->Gsvolu("QA28", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA28", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA28", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA28  TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1083,7 +1099,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 20.6/2.;
   conpar[3] = 23.9/2.;
   conpar[4] = 24.3/2.;
-  gMC->Gsvolu("QA29", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QA29", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QA29", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   // Ch.debug  
   //printf("   QA29 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
@@ -1093,7 +1109,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 23.9/2.;
   tubpar[1] = 25.5/2.;
   tubpar[2] = 17.0/2.;
-  gMC->Gsvolu("QA30", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA30", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA30", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug  
   //printf("   QA30  TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1105,7 +1121,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 24.3/2.;
   conpar[3] = 20.2/2.;
   conpar[4] = 20.6/2.;
-  gMC->Gsvolu("QA31", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QA31", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QA31", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   // Ch.debug  
   //printf("   QA31 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
@@ -1115,7 +1131,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 20.2/2.;
   tubpar[1] = 20.6/2.;
   tubpar[2] = 2.15/2.;
-  gMC->Gsvolu("QA32", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA32", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA32", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug  
   //printf("   QA32  TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1126,7 +1142,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 19.6/2.;
   tubpar[1] = 25.3/2.;
   tubpar[2] = 4.9/2.;
-  gMC->Gsvolu("QF02", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QF02", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QF02", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QF02 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1139,7 +1155,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 19.6/2.;
   tubpar[1] = 20.0/2.;
   tubpar[2] = 3.9/2.;
-  gMC->Gsvolu("QA33", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA33", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA33", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug
   //printf("   QA33  TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1152,7 +1168,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 20.0/2.;
   conpar[3] = 21.66/2.;
   conpar[4] = 22.06/2.;
-  gMC->Gsvolu("QA34", "CONE", idtmed[6], conpar, 5);
+  gMC->Gsvolu("QA34", "CONE", idtmed[7], conpar, 5);
   gMC->Gspos("QA34", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
   // Ch.debug  
   //printf("   QA34 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
@@ -1163,21 +1179,19 @@ void AliZDCv3::CreateBeamLine()
   boxpar[0] = 25.3/2.;
   boxpar[1] = 25.3/2.;
   boxpar[2] = 2.5/2.;
-  gMC->Gsvolu("QF03", "BOX ", idtmed[6], boxpar, 3);
+  gMC->Gsvolu("QF03", "BOX ", idtmed[7], boxpar, 3);
   tubpar[0] = 0.0/2.;
   tubpar[1] = 22.06/2.;
   tubpar[2] = 2.5/2.;
   gMC->Gsvolu("QFV1", "TUBE", idtmed[10], tubpar, 3);
   gMC->Gspos("QFV1", 1, "QF03", 0., 0., 0., 0, "MANY");
   gMC->Gspos("QF03", 1, "ZDCA", 0., 0., 14.3+zd2, 0, "MANY");  
-  // Ch.debug
-  //printf("\n Flange:  first support for the trousers\n");
   
   // tube  
   tubpar[0] = 21.66/2.;
   tubpar[1] = 22.06/2.;
   tubpar[2] = 28.6/2.;
-  gMC->Gsvolu("QA35", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA35", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA35", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
   // Ch.debug 
   //printf("\n QA35  TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);
@@ -1190,8 +1204,8 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 21.6/2.;
   conpar[3] = 0.0/2.;
   conpar[4] = 5.8/2.;
-  gMC->Gsvolu("QAL1", "CONE", idtmed[6], conpar, 5);
-  gMC->Gsvolu("QAL2", "CONE", idtmed[6], conpar, 5); 
+  gMC->Gsvolu("QAL1", "CONE", idtmed[7], conpar, 5);
+  gMC->Gsvolu("QAL2", "CONE", idtmed[7], conpar, 5); 
   gMC->Gspos("QAL1", 1, "ZDCA", -3.45-0.52, 0., (90.1/2.)+zd2, irotpipe5, "MANY");
   gMC->Gspos("QAL2", 1, "ZDCA", 3.45+0.52, 0., (90.1/2.)+zd2, irotpipe6, "MANY");
   
@@ -1211,7 +1225,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 5.4/2.;
   tubpar[1] = 5.8/2.;
   tubpar[2] = 40.0/2.;
-  gMC->Gsvolu("QA36", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA36", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA36", 1, "ZDCA", -15.8/2., 0., tubpar[2]+zd2, 0, "ONLY");
   gMC->Gspos("QA36", 2, "ZDCA",  15.8/2., 0., tubpar[2]+zd2, 0, "ONLY");  
   // Ch.debug
@@ -1225,7 +1239,7 @@ void AliZDCv3::CreateBeamLine()
   conpar[2] = 5.8/2.;
   conpar[3] = 6.3/2.;
   conpar[4] = 7.0/2.;
-  gMC->Gsvolu("QA37", "CONE", idtmed[6], conpar, 5); 
+  gMC->Gsvolu("QA37", "CONE", idtmed[7], conpar, 5); 
   gMC->Gspos("QA37", 1, "ZDCA", -7.9-0.175, 0., conpar[0]+zd2, irotpipe7, "ONLY");
   gMC->Gspos("QA37", 2, "ZDCA", 7.9+0.175, 0., conpar[0]+zd2, irotpipe8, "ONLY");
   //printf("   QA37 CONE from z = %f to z= %f\n",zd2,2*conpar[0]+zd2);
@@ -1236,7 +1250,7 @@ void AliZDCv3::CreateBeamLine()
   boxpar[0] = 25.9/2.;
   boxpar[1] = 9.4/2.;
   boxpar[2] = 1./2.;
-  gMC->Gsvolu("QF04", "BOX ", idtmed[6], boxpar, 3);
+  gMC->Gsvolu("QF04", "BOX ", idtmed[7], boxpar, 3);
   boxpar[0] = 16.5/2.;
   boxpar[1] = 7./2.;
   boxpar[2] = 1./2.;
@@ -1258,19 +1272,25 @@ void AliZDCv3::CreateBeamLine()
   gMC->Gspos("QFV4", 1, "QF04", 16.5/2., 0.,  0., 0, "MANY");
   gMC->Gspos("QF04", 1, "ZDCA", 0., 0.,  18.5+zd2, 0, "MANY");
   
-
   // 2 tubes (ID = 63 mm OD=70 mm)      
   tubpar[0] = 6.3/2.;
   tubpar[1] = 7.0/2.;
   tubpar[2] = 512.9/2.;
-  gMC->Gsvolu("QA38", "TUBE", idtmed[6], tubpar, 3);
+  gMC->Gsvolu("QA38", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("QA38", 1, "ZDCA", -16.5/2., 0., tubpar[2]+zd2, 0, "ONLY");
   gMC->Gspos("QA38", 2, "ZDCA",  16.5/2., 0., tubpar[2]+zd2, 0, "ONLY");
   //printf("   QA38 TUBE from z = %f to z= %f\n",zd2,2*tubpar[2]+zd2);  
 
   zd2 += 2.*tubpar[2];
   //printf("\n END OF BEAM PIPE VOLUME DEFINITION AT z= %f\n",zd2);
-  //printf("   MAGNET DEFINITION FOLLOWS\n\n"); 
+          
+  // -- Luminometer (Cu box) in front of ZN - side A
+  boxpar[0] = 8.0/2.;
+  boxpar[1] = 8.0/2.;
+  boxpar[2] = 15./2.;
+  gMC->Gsvolu("QLUA", "BOX ", idtmed[7], boxpar, 3);
+  gMC->Gspos("QLUA", 1, "ZDCA", 0., 0.,  fPosZNA[2]-66.-boxpar[2], 0, "ONLY");
+  //printf("\n QLUC LUMINOMETER from z = %f to z= %f\n\n",  fPosZNA[2]-66., fPosZNA[2]-66.-2*boxpar[2]);
   
 
   // ----------------------------------------------------------------
@@ -1280,19 +1300,22 @@ void AliZDCv3::CreateBeamLine()
   //           SIDE C - RB26  (dimuon side) 
   // ***************************************************************   
   // --  COMPENSATOR DIPOLE (MBXW)
-  zc = 1921.6;   
+  //zc = 1921.6;   
+  zc = 1972.5;   
   
   // --  GAP (VACUUM WITH MAGNETIC FIELD)
   tubpar[0] = 0.;
   tubpar[1] = 4.5;
-  tubpar[2] = 170./2.;
+  //tubpar[2] = 170./2.;
+  tubpar[2] = 153./2.;
   gMC->Gsvolu("MBXW", "TUBE", idtmed[11], tubpar, 3);
 
   // --  YOKE 
   tubpar[0] = 4.5;
   tubpar[1] = 55.;
-  tubpar[2] = 170./2.;
-  gMC->Gsvolu("YMBX", "TUBE", idtmed[13], tubpar, 3);
+  //tubpar[2] = 170./2.;
+  tubpar[2] = 153./2.;
+  gMC->Gsvolu("YMBX", "TUBE", idtmed[7], tubpar, 3);
 
   gMC->Gspos("MBXW", 1, "ZDCC", 0., 0., -tubpar[2]-zc, 0, "ONLY");
   gMC->Gspos("YMBX", 1, "ZDCC", 0., 0., -tubpar[2]-zc, 0, "ONLY");
@@ -1318,8 +1341,8 @@ void AliZDCv3::CreateBeamLine()
   gMC->Gspos("MQXL", 1, "ZDCC", 0., 0., -tubpar[2]-zq, 0, "ONLY");
   gMC->Gspos("YMQL", 1, "ZDCC", 0., 0., -tubpar[2]-zq, 0, "ONLY");
   
-  gMC->Gspos("MQXL", 2, "ZDCC", 0., 0., -tubpar[2]-zq-2430., 0, "ONLY");
-  gMC->Gspos("YMQL", 2, "ZDCC", 0., 0., -tubpar[2]-zq-2430., 0, "ONLY");
+  gMC->Gspos("MQXL", 2, "ZDCC", 0., 0., -tubpar[2]-zq-2400., 0, "ONLY");
+  gMC->Gspos("YMQL", 2, "ZDCC", 0., 0., -tubpar[2]-zq-2400., 0, "ONLY");
   
   // --  MQX 
   // --  GAP (VACUUM WITH MAGNETIC FIELD) 
@@ -1366,10 +1389,12 @@ void AliZDCv3::CreateBeamLine()
   
   gMC->Gspos("YD1 ", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
   gMC->Gspos("MD1 ", 1, "YD1 ", 0., 0., 0., 0, "ONLY");
+  // Ch debug
+  //printf("\t **** D1 positioned! It goes from z = %1.2f to z = %1.2f cm\n",-zd1, -zd1-2*tubpar[2]); 
   
   // -- DIPOLE D2 
   // --- LHC optics v6.4
-  zd2 = 12147.6;
+  zd2 = 12167.8;
   
   // --  GAP (VACUUM WITH MAGNETIC FIELD) 
   tubpar[0] = 0.;
@@ -1404,7 +1429,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 4.5;
   tubpar[1] = 55.;
   tubpar[2] = 153./2.;
-  gMC->Gsvolu("YMCB", "TUBE", idtmed[8], tubpar, 3);
+  gMC->Gsvolu("YMCB", "TUBE", idtmed[7], tubpar, 3);
   gMC->Gspos("YMCB", 1, "ZDCA", 0., 0., tubpar[2]+1972.5, 0, "ONLY");  
   
   
@@ -1423,7 +1448,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 3.5;
   tubpar[1] = 22.;
   tubpar[2] = 637./2.;
-  gMC->Gsvolu("YMQ1", "TUBE", idtmed[8], tubpar, 3);
+  gMC->Gsvolu("YMQ1", "TUBE", idtmed[7], tubpar, 3);
 
    // -- BEAM SCREEN FOR Q1
    tubpar[0] = 4.78/2.;
@@ -1473,7 +1498,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 3.5;
   tubpar[1] = 22.;
   tubpar[2] = 550./2.;
-  gMC->Gsvolu("YMQ2", "TUBE", idtmed[8], tubpar, 3);
+  gMC->Gsvolu("YMQ2", "TUBE", idtmed[7], tubpar, 3);
   
 
    // -- BEAM SCREEN FOR Q2
@@ -1525,7 +1550,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 7.34/2.; // to be checked
   tubpar[1] = 110./2;
   tubpar[2] = 945./2.;
-  gMC->Gsvolu("YD1L", "TUBE", idtmed[8], tubpar, 3);
+  gMC->Gsvolu("YD1L", "TUBE", idtmed[7], tubpar, 3);
   
   gMC->Gspos("YD1L", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");  
   gMC->Gspos("MD1L", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "MANY");
@@ -1545,7 +1570,7 @@ void AliZDCv3::CreateBeamLine()
   tubpar[0] = 0.;
   tubpar[1] = 55.;
   tubpar[2] = 945./2.;
-  gMC->Gsvolu("YD2L", "TUBE", idtmed[8], tubpar, 3);
+  gMC->Gsvolu("YD2L", "TUBE", idtmed[7], tubpar, 3);
   
   gMC->Gspos("YD2L", 1, "ZDCA", 0., 0., tubpar[2]+zd2l, 0, "ONLY");
   
@@ -1895,7 +1920,7 @@ void AliZDCv3::CreateMaterials()
   z[1] = 30.;
   wmat[0] = .63;
   wmat[1] = .37;
-  AliMixture(2, "BRASS               ", a, z, dens, 2, wmat);
+  AliMixture(2, "BRASS", a, z, dens, 2, wmat);
   
   // --- SiO2 
   dens = 2.64;
@@ -1905,7 +1930,7 @@ void AliZDCv3::CreateMaterials()
   z[1] = 8.;
   wmat[0] = 1.;
   wmat[1] = 2.;
-  AliMixture(3, "SIO2                ", a, z, dens, -2, wmat);  
+  AliMixture(3, "SIO2", a, z, dens, -2, wmat);  
   
   // --- Lead 
   ubuf[0] = 1.12;
@@ -1921,7 +1946,7 @@ void AliZDCv3::CreateMaterials()
   
   // --- Iron (no energy loss)
   ubuf[0] = 1.1;
-  AliMaterial(8, "IRON1", 55.85, 26., 7.87, 1.76, 0., ubuf, 1);
+  AliMaterial(8,  "IRON1", 55.85, 26., 7.87, 1.76, 0., ubuf, 1);
   AliMaterial(13, "IRON2", 55.85, 26., 7.87, 1.76, 0., ubuf, 1);
     
   // ---------------------------------------------------------  
@@ -2067,6 +2092,7 @@ void AliZDCv3::Init()
   gMC->Gstpar(idtmed[i], "CUTELE", .1);
   gMC->Gstpar(idtmed[i], "CUTNEU", 1.);
   gMC->Gstpar(idtmed[i], "CUTHAD", 1.);
+  
   // Avoid too detailed showering along the beam line 
   i = 13; //iron with energy loss (ZIRONN)
   gMC->Gstpar(idtmed[i], "CUTGAM", 1.);
@@ -2217,37 +2243,34 @@ void AliZDCv3::StepManager()
   //
   // Routine called at every step in the Zero Degree Calorimeters
   //
-    
-  Int_t j, vol[2], ibeta=0, ialfa, ibe, nphe;
-  Float_t x[3], xdet[3], destep, hits[10], m, ekin, um[3], ud[3], be, out;
-  //Float_t radius;
-  Float_t xalic[3], z, guiEff;
+  Int_t   j, vol[2]={0,0}, ibeta=0, ialfa=0, ibe=0, nphe=0;
+  Float_t hits[11], x[3], xdet[3], um[3], ud[3];
+  Float_t m=0., ekin=0., destep=0., be=0., out=0.;
   // Parametrization for light guide uniformity
-  // -> OBSOLETE!!!! For guide tilted @ 46 degrees
-  //Float_t guiPar[4]={0.31,-0.0004,0.0197,0.7958};
   // NEW!!! Light guide tilted @ 51 degrees
   Float_t guiPar[4]={0.31,-0.0006305,0.01337,0.8895};
   Double_t s[3], p[3];
   const char *knamed;
   //
-  for (j=0;j<10;j++) hits[j]=-999.;
+  for(j=0;j<11;j++) hits[j]=-999.;
   //
   // --- This part is for no shower developement in beam pipe and TDI
   // If particle interacts with beam pipe or TDI -> return
   if((gMC->CurrentMedium() == fMedSensPI) || (gMC->CurrentMedium() == fMedSensTDI)){ 
      // If option NoShower is set -> StopTrack
-    if(fNoShower==1) {
+    if(fNoShower==1){
       gMC->TrackPosition(s[0],s[1],s[2]);
-      if(gMC->CurrentMedium() == fMedSensPI) {
+      if(gMC->CurrentMedium() == fMedSensPI){
         knamed = gMC->CurrentVolName();
         if(!strncmp(knamed,"YMQ",3)){
          if(s[2]<0) fpLostITC += 1;
          else fpLostITA += 1;
         }
-       if(!strncmp(knamed,"YD1",3)){
+       else if(!strncmp(knamed,"YD1",3)){
          if(s[2]<0) fpLostD1C += 1;
          else fpLostD1A += 1;
        }
+       else if(!strncmp(knamed,"QAL",3)) fnTrou++;
       }
       else if(gMC->CurrentMedium() == fMedSensTDI){ 
         knamed = gMC->CurrentVolName();
@@ -2255,15 +2278,29 @@ void AliZDCv3::StepManager()
          if(s[2]<0) fpLostD1C += 1;
          else  fpLostD1A += 1;
         }
-       if(!strncmp(knamed,"QTD",3)) fpLostTDI += 1;
+       else if(!strncmp(knamed,"QTD",3)) fpLostTDI += 1;
+       else if(!strncmp(knamed,"QLU",3)){
+         if(s[2]<0) fnLumiC ++;
+         else fnLumiA++;
+       }
       }
-      printf("\n\t ---------- Side C ----------\n");
-      printf("\n      # of spectators lost in IT = %d\n",fpLostITC);
-      printf("\n      # of spectators lost in D1  = %d\n",fpLostD1C);
-      printf("\n\t ---------- Side A ----------\n");
-      printf("\n      # of spectators lost in IT = %d\n",fpLostITA);
-      printf("\n      # of spectators lost in D1  = %d\n",fpLostD1A);
-      printf("\n      # of spectators lost in TDI = %d\n\n",fpLostTDI);
+      //
+      //gMC->TrackMomentum(p[0], p[1], p[2], p[3]);
+      //printf("\t Particle: mass = %1.3f, E = %1.3f GeV, pz = %1.2f GeV -> stopped in volume %s\n", 
+      //     gMC->TrackMass(), p[3], p[2], gMC->CurrentVolName());
+      //
+      printf("\n\t **********************************\n");
+      printf("\t ********** Side C **********\n");
+      printf("\t # of spectators in IT = %d\n",fpLostITC);
+      printf("\t # of spectators in D1 = %d\n",fpLostD1C);
+      printf("\t # of spectators in luminometer = %d\n",fnLumiC);
+      printf("\t ********** Side A **********\n");
+      printf("\t # of spectators in IT = %d\n",fpLostITA);
+      printf("\t # of spectators in D1 = %d\n",fpLostD1A);
+      printf("\t # of spectators in TDI = %d\n",fpLostTDI);
+      printf("\t # of spectators in luminometer = %d\n",fnLumiA);
+      printf("\t # of spectators in trousers = %d\n",fnTrou);
+      printf("\t **********************************\n");
       gMC->StopTrack();
     }
     return;
@@ -2274,7 +2311,7 @@ void AliZDCv3::StepManager()
      (gMC->CurrentMedium() == fMedSensGR) || (gMC->CurrentMedium() == fMedSensF1) ||
      (gMC->CurrentMedium() == fMedSensF2) || (gMC->CurrentMedium() == fMedSensZEM)){
 
-  
+    
   //Particle coordinates 
     gMC->TrackPosition(s[0],s[1],s[2]);
     for(j=0; j<=2; j++) x[j] = s[j];
@@ -2399,28 +2436,37 @@ void AliZDCv3::StepManager()
         hits[7] = 0;
         hits[8] = 0;
         hits[9] = 0;
+       //
+       Int_t curTrackN = gAlice->GetMCApp()->GetCurrentTrackNumber();
+        TParticle *part = (gAlice->GetMCApp())->Particle(curTrackN);
+       hits[10] = part->GetPdgCode();
+       //printf("\t PDGCode = %d\n", part->GetPdgCode());
 
-       AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
+       AddHit(curTrackN, vol, hits);
        
        if(fNoShower==1){
          //printf("\t VolName %s -> det %d quad %d - x = %f, y = %f, z = %f\n", 
            //knamed, vol[0], vol[1], x[0], x[1], x[2]);
          if(vol[0]==1){
            fnDetectedC += 1;
-           printf("\n    # of detected neutrons in ZNC = %d\n\n",fnDetectedC);
+           printf("\n    # of particles in ZNC = %d\n\n",fnDetectedC);
          }
          else if(vol[0]==2){
            fpDetectedC += 1;
-           printf("\n    # of detected protons in ZPC = %d\n\n",fpDetectedC);
+           printf("\n    # of particles in ZPC = %d\n\n",fpDetectedC);
          }
          else if(vol[0]==4){
            fnDetectedA += 1;
-           printf("\n    # of detected neutrons in ZNA = %d\n\n",fnDetectedA);     
+           printf("\n    # of particles in ZNA = %d\n\n",fnDetectedA);     
          }
          else if(vol[0]==5){
            fpDetectedA += 1;
-           printf("\n    # of detected protons in ZPA = %d\n\n",fpDetectedA);      
+           printf("\n    # of particles in ZPA = %d\n\n",fpDetectedA);      
          }
+          //
+         //printf("\t Particle: mass = %1.3f, E = %1.3f GeV, pz = %1.2f GeV -> stopped in volume %s\n", 
+          //   gMC->TrackMass(), p[3], p[2], gMC->CurrentVolName());
+         //
          gMC->StopTrack();
          return;
        }
@@ -2540,15 +2586,16 @@ void AliZDCv3::StepManager()
          if(ibe>fNbep) ibe=fNbep;
          out =  charge*charge*fTablep[ibeta][ialfa][ibe];
         gMC->TrackPosition(s[0],s[1],s[2]);
-         for(j=0; j<=2; j++){
+        Float_t xalic[3];
+         for(j=0; j<3; j++){
             xalic[j] = s[j];
          }
         // z-coordinate from ZEM front face 
         // NB-> fPosZEM[2]+fZEMLength = -1000.+2*10.3 = 979.69 cm
-        z = -xalic[2]+fPosZEM[2]+2*fZEMLength-xalic[1];
+        Float_t z = -xalic[2]+fPosZEM[2]+2*fZEMLength-xalic[1];
 //      z = xalic[2]-fPosZEM[2]-fZEMLength-xalic[1]*(TMath::Tan(45.*kDegrad));
 //         printf("\n  fPosZEM[2]+2*fZEMLength = %f", fPosZEM[2]+2*fZEMLength);
-        guiEff = guiPar[0]*(guiPar[1]*z*z+guiPar[2]*z+guiPar[3]);
+        Float_t guiEff = guiPar[0]*(guiPar[1]*z*z+guiPar[2]*z+guiPar[3]);
         out = out*guiEff;
         nphe = gRandom->Poisson(out);
 //         printf("    out*guiEff = %f nphe = %d", out, nphe);