In AliMUONTriggerQADataMakerRec:
[u/mrichter/AliRoot.git] / ZDC / AliZDCv3.cxx
index a986bd2..2565e0e 100644 (file)
@@ -2128,6 +2128,13 @@ void AliZDCv3::Init()
   gMC->Gstpar(idtmed[i], "CUTNEU", 1.);
   gMC->Gstpar(idtmed[i], "CUTHAD", 1.);
   
+  // Thresholds for showering in the luminometer
+  i = 9; //copper
+  gMC->Gstpar(idtmed[i], "CUTGAM", .001);
+  gMC->Gstpar(idtmed[i], "CUTELE", .001);
+  gMC->Gstpar(idtmed[i], "CUTNEU", .01);
+  gMC->Gstpar(idtmed[i], "CUTHAD", .01);
+  
   // Avoid too detailed showering along the beam line 
   i = 7; //iron with energy loss (ZIRON)
   gMC->Gstpar(idtmed[i], "CUTGAM", .1);
@@ -2231,12 +2238,17 @@ void AliZDCv3::InitTables()
      return;
   }
   
+  int read=1;
   for(k=0; k<fNalfan; k++){
      for(j=0; j<fNben; j++){
-       fscanf(fp1,"%f",&fTablen[0][k][j]);
-       fscanf(fp2,"%f",&fTablen[1][k][j]);
-       fscanf(fp3,"%f",&fTablen[2][k][j]);
-       fscanf(fp4,"%f",&fTablen[3][k][j]);
+       read = fscanf(fp1,"%f",&fTablen[0][k][j]);
+       if(read==0) AliDebug(3, " Error in reading light table 1");
+       read = fscanf(fp2,"%f",&fTablen[1][k][j]);
+       if(read==0) AliDebug(3, " Error in reading light table 2");
+       read = fscanf(fp3,"%f",&fTablen[2][k][j]);
+       if(read==0) AliDebug(3, " Error in reading light table 3");
+       read = fscanf(fp4,"%f",&fTablen[3][k][j]);
+       if(read==0) AliDebug(3, " Error in reading light table 4");
      } 
   }
   fclose(fp1);
@@ -2268,10 +2280,14 @@ void AliZDCv3::InitTables()
   
   for(k=0; k<fNalfap; k++){
      for(j=0; j<fNbep; j++){
-       fscanf(fp5,"%f",&fTablep[0][k][j]);
-       fscanf(fp6,"%f",&fTablep[1][k][j]);
-       fscanf(fp7,"%f",&fTablep[2][k][j]);
-       fscanf(fp8,"%f",&fTablep[3][k][j]);
+       read = fscanf(fp5,"%f",&fTablep[0][k][j]);
+       if(read==0) AliDebug(3, " Error in reading light table 5");
+       read = fscanf(fp6,"%f",&fTablep[1][k][j]);
+       if(read==0) AliDebug(3, " Error in reading light table 6");
+       read = fscanf(fp7,"%f",&fTablep[2][k][j]);
+       if(read==0) AliDebug(3, " Error in reading light table 7");
+       read = fscanf(fp8,"%f",&fTablep[3][k][j]);
+       if(read==0) AliDebug(3, " Error in reading light table 8");
      } 
   }
   fclose(fp5);
@@ -2286,12 +2302,12 @@ void AliZDCv3::StepManager()
   // Routine called at every step in the Zero Degree Calorimeters
   //
   Int_t   j, vol[2]={0,0}, ibeta=0, ialfa=0, ibe=0, nphe=0;
-  Float_t hits[12], x[3], xdet[3], um[3], ud[3];
-  Float_t m=0., ekin=0., destep=0., be=0., out=0.;
+  Float_t hits[13], x[3], xdet[3], um[3], ud[3];
+  Float_t destep=0., be=0., out=0.;
   Double_t s[3], p[3];
   const char *knamed;
   //
-  for(j=0;j<12;j++) hits[j]=-999.;
+  for(j=0;j<13;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
@@ -2327,7 +2343,7 @@ void AliZDCv3::StepManager()
       //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());
       //
-      if(ipr!=0){
+      /*if(ipr!=0){
         printf("\n\t **********************************\n");
         printf("\t ********** Side C **********\n");
         printf("\t # of spectators in IT = %d\n",fpLostITC);
@@ -2337,7 +2353,7 @@ void AliZDCv3::StepManager()
         printf("\t # of spectators in D1 = %d\n",fpLostD1A);
         printf("\t # of spectators in TDI = %d\n",fpLostTDI);
         printf("\t **********************************\n");
-      }
+      }*/
       gMC->StopTrack();
     }
     return;
@@ -2453,80 +2469,80 @@ void AliZDCv3::StepManager()
     //printf("\t *** det %d vol %d xdet(%f, %f)\n",vol[0], vol[1], xdet[0], xdet[1]);
     
     
-  // Store impact point and kinetic energy of the ENTERING particle
+    // Store impact point and kinetic energy of the ENTERING particle
     
-      if(gMC->IsTrackEntering()){
-        //Particle energy
-        gMC->TrackMomentum(p[0],p[1],p[2],p[3]);
-        hits[3] = p[3];
-        // Impact point on ZDC
-       // X takes into account the LHC x-axis sign
-       // which is opposite to positive x on detcetor front face
-       // for side A detectors (ZNA and ZPA)  
-        if(vol[0]==4 || vol[0]==5) hits[4] = -xdet[0];
-       else hits[4] = xdet[0];
-        hits[5] = xdet[1];
-       hits[6] = 0;
-        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());
-        //
-       hits[11] = 1.0e09*gMC->TrackTime();
-       //printf("\t TrackTime = %f\n", hits[11]);
-
-       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("      # of particles in ZNC = %d\n\n",fnDetectedC);
-         }
-         else if(vol[0]==2){
-           fpDetectedC += 1;
-           printf("      # of particles in ZPC = %d\n\n",fpDetectedC);
-         }
-         else if(vol[0]==4){
-           fnDetectedA += 1;
-           printf("      # of particles in ZNA = %d\n\n",fnDetectedA);     
-         }
-         else if(vol[0]==5){
-           fpDetectedA += 1;
-           printf("      # 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;
-       }
+    if(gMC->IsTrackEntering()){
+      //Particle energy
+      gMC->TrackMomentum(p[0],p[1],p[2],p[3]);
+      hits[3] = p[3];
+      // Impact point on ZDC
+      // X takes into account the LHC x-axis sign
+      // which is opposite to positive x on detcetor front face
+      // for side A detectors (ZNA and ZPA)  
+      if(vol[0]==4 || vol[0]==5){
+        hits[4] = -xdet[0];
       }
-             
-      // Charged particles -> Energy loss
-      if((destep=gMC->Edep())){
-         if(gMC->IsTrackStop()){
-           gMC->TrackMomentum(p[0],p[1],p[2],p[3]);
-          m = gMC->TrackMass();
-          ekin = p[3]-m;
-          hits[9] = ekin;
-          hits[7] = 0.;
-          hits[8] = 0.;
-          AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
-        }
-        else{
-          hits[9] = destep;
-          hits[7] = 0.;
-          hits[8] = 0.;
-          AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
-        }
+      else{
+        hits[4] = xdet[0];
       }
+      hits[5] = xdet[1];
+      hits[6] = 0;
+      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());
+      //
+      Int_t imo = part->GetFirstMother();
+      if(imo>0){
+       TParticle * pmot = gAlice->GetMCApp()->Particle(imo);
+       hits[11] = pmot->GetPdgCode();
+      }
+      else hits[11]=0;
+      //
+      hits[12] = 1.0e09*gMC->TrackTime(); // in ns!
+      //printf("\t TrackTime = %f\n", hits[12]);
+
+      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("     # of particles in ZNC = %d\n\n",fnDetectedC);
+        }
+        else if(vol[0]==2){
+          fpDetectedC += 1;
+          printf("     # of particles in ZPC = %d\n\n",fpDetectedC);
+        }
+        else if(vol[0]==4){
+          fnDetectedA += 1;
+          printf("     # of particles in ZNA = %d\n\n",fnDetectedA);     
+        }
+        else if(vol[0]==5){
+          fpDetectedA += 1;
+          printf("     # 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;
+      }
+    }
+          
+    // Particle energy loss
+    if(gMC->Edep() != 0){
+      hits[9] = gMC->Edep();
+      hits[7] = 0.;
+      hits[8] = 0.;
+      AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
+    }
   }