Some function moved to AliZDC
authorcoppedis <coppedis@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 26 Sep 2001 16:07:55 +0000 (16:07 +0000)
committercoppedis <coppedis@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 26 Sep 2001 16:07:55 +0000 (16:07 +0000)
ZDC/AliZDCv1.cxx
ZDC/AliZDCv1.h
ZDC/AliZDCv2.h

index 4a55b9f1c0e4ef1ebd75f8e48599f7479421f303..ff273ff642348b93fd52b25c9ebbcc98ca44dae1 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.26  2001/06/13 11:10:55  coppedis
+Minor changes
+
 Revision 1.25  2001/06/12 13:45:05  coppedis
 TDI in correct position and minor correction
 
@@ -102,10 +105,9 @@ Introduction of the Copyright and cvs Log
 
 ///////////////////////////////////////////////////////////////////////////////
 //                                                                           //
-//  Zero Degree Calorimeter                                                  //
-//  This class contains the basic functions for the ZDC                      //
-//  Functions specific to one particular geometry are                        //
-//  contained in the derived classes                                         //
+//             AliZDCv1 --- ZDC geometry as designed in TDR (obsolete!)     //
+//                     with the EM ZDC at 116 m from IP                     //
+//  Just one set of ZDC is inserted, on the same side of the dimuon arm      //
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
@@ -124,7 +126,6 @@ Introduction of the Copyright and cvs Log
 // --- AliRoot classes
 #include "AliZDCv1.h"
 #include "AliZDCHit.h"
-#include "AliZDCDigit.h"
 #include "AliRun.h"
 #include "AliDetector.h"
 #include "AliMagF.h"
@@ -137,13 +138,6 @@ Introduction of the Copyright and cvs Log
  
 ClassImp(AliZDCv1)
  
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-//  Zero Degree Calorimeter version 1                                        //
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
-
 //_____________________________________________________________________________
 AliZDCv1::AliZDCv1() : AliZDC()
 {
@@ -171,10 +165,10 @@ AliZDCv1::AliZDCv1(const char *name, const char *title)
   //
   // Check that DIPO, ABSO, DIPO and SHIL is there (otherwise tracking is wrong!!!)
   
-  AliModulePIPE=gAlice->GetModule("PIPE");
-  AliModuleABSO=gAlice->GetModule("ABSO");
-  AliModuleDIPO=gAlice->GetModule("DIPO");
-  AliModuleSHIL=gAlice->GetModule("SHIL");
+  AliModule *PIPE=gAlice->GetModule("PIPE");
+  AliModule *ABSO=gAlice->GetModule("ABSO");
+  AliModule *DIPO=gAlice->GetModule("DIPO");
+  AliModule *SHIL=gAlice->GetModule("SHIL");
   if((!PIPE) || (!ABSO) || (!DIPO) || (!SHIL)) {
     Error("Constructor","ZDC needs PIPE, ABSO, DIPO and SHIL!!!\n");
     exit(1);
@@ -217,7 +211,7 @@ AliZDCv1::AliZDCv1(const char *name, const char *title)
   fDimZP[1] = 6.;
   fDimZP[2] = 75.;    
   fPosZN[0] = 0.;
-  fPosZN[1] = -1.2;
+  fPosZN[1] = 1.2;
   fPosZN[2] = 11650.;
   fPosZP[0] = -24.;
   fPosZP[1] = 0.;
@@ -234,8 +228,6 @@ AliZDCv1::AliZDCv1(const char *name, const char *title)
   fPosZEM[1] = 5.8;
   fPosZEM[2] = 11600.;
   
-
-  fDigits = new TClonesArray("AliZDCDigit",1000);
 }
  
 //_____________________________________________________________________________
@@ -385,8 +377,7 @@ void AliZDCv1::CreateBeamLine()
   // QT09 is 10 cm longer to accomodate TDI
   tubpar[2] = 515.4/2.;
   gMC->Gsvolu("QT09", "TUBE", idtmed[7], tubpar, 3);
-  gMC->Gspos("QT09", 1, "ZDC ", 0., 0., tubpar[2] + zd1, 0, "ONLY");
-  
+  gMC->Gspos("QT09", 1, "ZDC ", 0., 0., tubpar[2] + zd1, 0, "ONLY"); 
   
   // --- Insert TDI (inside ZDC volume)
   
@@ -1245,165 +1236,6 @@ void AliZDCv1::InitTables()
   fclose(fp7);
   fclose(fp8);
 }
-
-//_____________________________________________________________________________
-Int_t AliZDCv1::Digitize(Int_t Det, Int_t Quad, Int_t Light)
-{
-  // Evaluation of the ADC channel corresponding to the light yield Light
-
-  if(fDebug == 1){
-    printf("\n Digitize -> Det = %d, Quad = %d, Light = %d\n", Det, Quad, Light);
-  }   
-  
-  // Parameters for conversion of light yield in ADC channels
-  Float_t fPMGain[3][5];      // PM gain
-  Float_t fADCRes;            // ADC conversion factor
-  
-  Int_t j,i;
-  for(i=0; i<3; i++){
-     for(j=0; j<5; j++){
-        fPMGain[i][j]   = 100000.;
-     }
-  }
-  fADCRes   = 0.00000064; // ADC Resolution: 250 fC/ADCch
-  
-  Int_t ADCch = Int_t(Light*fPMGain[Det-1][Quad]*fADCRes);
-     
-  return ADCch;
-}
-
-
-//_____________________________________________________________________________
-void AliZDCv1::SDigits2Digits()
-{
-   Hits2Digits(gAlice->GetNtrack());
-}
-
-//_____________________________________________________________________________
-void AliZDCv1::Hits2Digits(Int_t ntracks)
-{
-  AliZDCDigit *newdigit;
-  AliZDCHit   *hit;
-
-  Int_t PMCZN = 0, PMCZP = 0, PMQZN[4], PMQZP[4], PMZEM = 0;
-  
-  Int_t i;
-  for(i=0; i<4; i++){
-     PMQZN[i] =0;
-     PMQZP[i] =0;
-  }
-  
-  Int_t itrack = 0;
-  for(itrack=0; itrack<ntracks; itrack++){
-     gAlice->ResetHits();
-     gAlice->TreeH()->GetEvent(itrack);
-     for(i=0; i<fHits->GetEntries(); i++){
-        hit = (AliZDCHit*)fHits->At(i);
-        Int_t det   = hit->GetVolume(0);
-        Int_t quad  = hit->GetVolume(1);
-        Int_t lightQ = Int_t(hit->GetLightPMQ());
-        Int_t lightC = Int_t(hit->GetLightPMC());
-       if(fDebug == 1)
-          printf("        \n itrack = %d, fNhits = %d, det = %d, quad = %d,"
-         "lightC = %d lightQ = %d\n", itrack, fNhits, det, quad, lightC, lightQ);
-           
-        if(det == 1){   //ZN 
-          PMCZN = PMCZN + lightC;
-          PMQZN[quad-1] = PMQZN[quad-1] + lightQ;
-        }
-
-        if(det == 2){   //ZP 
-          PMCZP = PMCZP + lightC;
-          PMQZP[quad-1] = PMQZP[quad-1] + lightQ;
-        }
-
-        if(det == 3){   //ZEM 
-          PMZEM = PMZEM + lightC;
-        }
-     } // Hits loop
-  
-  } // Tracks loop
-  
-     if(fDebug == 1){
-       printf("\n       PMCZN = %d, PMQZN[0] = %d, PMQZN[1] = %d, PMQZN[2] = %d, PMQZN[3] = %d\n"
-           , PMCZN, PMQZN[0], PMQZN[1], PMQZN[2], PMQZN[3]);
-       printf("\n       PMCZP = %d, PMQZP[0] = %d, PMQZP[1] = %d, PMQZP[2] = %d, PMQZP[3] = %d\n"
-           , PMCZP, PMQZP[0], PMQZP[1], PMQZP[2], PMQZP[3]);
-       printf("\n       PMZEM = %d\n", PMZEM);
-     }
-
-  // ------------------------------------    Hits2Digits
-  // Digits for ZN
-     newdigit = new AliZDCDigit(1, 0, Digitize(1, 0, PMCZN));
-     new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
-     fNdigits++;
-     delete newdigit;
-  
-     Int_t j;
-     for(j=0; j<4; j++){
-        newdigit = new AliZDCDigit(1, j+1, Digitize(1, j+1, PMQZN[j]));
-        new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
-        fNdigits++;
-        delete newdigit;
-     }
-  
-     // Digits for ZP
-     newdigit = new AliZDCDigit(2, 0, Digitize(2, 0, PMCZP));
-     new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
-     fNdigits++;
-     delete newdigit;
-  
-     Int_t k;
-     for(k=0; k<4; k++){
-        newdigit = new AliZDCDigit(2, k+1, Digitize(2, k+1, PMQZP[k]));
-        new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
-        fNdigits++;
-        delete newdigit;
-     }
-  
-     // Digits for ZEM
-     newdigit = new AliZDCDigit(3, 0, Digitize(3, 0, PMZEM));
-     new((*fDigits)[fNdigits]) AliZDCDigit(*newdigit);
-     fNdigits++;
-     delete newdigit;
-      
-  
-  gAlice->TreeD()->Fill();
-  gAlice->TreeD()->Write(0,TObject::kOverwrite);
-
-//  if(fDebug == 1){
-//    printf("\n  Event Digits -----------------------------------------------------\n");  
-//    fDigits->Print("");
-//  }
-  
-}
-//_____________________________________________________________________________
- void AliZDCv1::MakeBranch(Option_t *opt, const char *file)
-{
-  //
-  // Create a new branch in the current Root Tree
-  //
-
-  AliDetector::MakeBranch(opt);
-  
-  Char_t branchname[10];
-  sprintf(branchname,"%s",GetName());
-  const char *cD = strstr(opt,"D");
-
-  if (gAlice->TreeD() && cD) {
-
-    // Creation of the digits from hits 
-
-    if(fDigits!=0) fDigits->Clear();
-    else fDigits = new TClonesArray ("AliZDCDigit",1000);
-    char branchname[10];
-    sprintf(branchname,"%s",GetName());
-    MakeBranchInTree(gAlice->TreeD(), 
-                     branchname, &fDigits, fBufferSize, file) ;
-    printf("* AliZDCv1::MakeBranch    * Making Branch %s for digits\n\n",branchname);
-  }
-       
-}
 //_____________________________________________________________________________
 void AliZDCv1::StepManager()
 {
@@ -1413,6 +1245,7 @@ void AliZDCv1::StepManager()
 
   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, radius, out;
+  Float_t xalic[3], z, GuiEff, GuiPar[4]={0.31,-0.0004,0.0197,0.7958};
   TLorentzVector s, p;
   const char *knamed;
 
@@ -1421,6 +1254,8 @@ void AliZDCv1::StepManager()
   if((gMC->GetMedium() == fMedSensZN) || (gMC->GetMedium() == fMedSensZP) ||
      (gMC->GetMedium() == fMedSensGR) || (gMC->GetMedium() == fMedSensF1) ||
      (gMC->GetMedium() == fMedSensF2) || (gMC->GetMedium() == fMedSensZEM)){
+
+//   --- This part is for no shower developement in beam pipe and TDI
 //     (gMC->GetMedium() == fMedSensPI) || (gMC->GetMedium() == fMedSensTDI)){
        
   // If particle interacts with beam pipe -> return
@@ -1452,14 +1287,19 @@ void AliZDCv1::StepManager()
 
   // Determine in which ZDC the particle is
     knamed = gMC->CurrentVolName();
-    if(!strncmp(knamed,"ZN",2))vol[0]=1;
-    if(!strncmp(knamed,"ZP",2))vol[0]=2;
-    if(!strncmp(knamed,"ZE",2))vol[0]=3;
+    if(!strncmp(knamed,"ZN",2)){
+      vol[0]=1;
+    }
+    else if(!strncmp(knamed,"ZP",2)){
+      vol[0]=2;
+    }
+    else if(!strncmp(knamed,"ZE",2)){
+      vol[0]=3;
+    }
   
   // Determine in which quadrant the particle is
-    
-    //Quadrant in ZN
-    if(vol[0]==1){
+       
+    if(vol[0]==1){     //Quadrant in ZN
       xdet[0] = x[0]-fPosZN[0];
       xdet[1] = x[1]-fPosZN[1];
       if((xdet[0]<=0.) && (xdet[1]>=0.))  vol[1]=1;
@@ -1467,9 +1307,7 @@ void AliZDCv1::StepManager()
       if((xdet[0]<0.)  && (xdet[1]<0.))   vol[1]=3;
       if((xdet[0]>0.)  && (xdet[1]<0.))   vol[1]=4;
     }
-    
-    //Quadrant in ZP
-    if(vol[0]==2){
+    else if(vol[0]==2){        //Quadrant in ZP
       xdet[0] = x[0]-fPosZP[0];
       xdet[1] = x[1]-fPosZP[1];
       if(xdet[0]>fDimZP[0])xdet[0]=fDimZP[0]-0.01;
@@ -1478,14 +1316,11 @@ void AliZDCv1::StepManager()
       for(int i=1; i<=4; i++){
          if(xqZP>=(i-3) && xqZP<(i-2)){
           vol[1] = i;
-          if(i==0) printf("\n!!! vol[1] = 0 -> xqZP = %f\n", xqZP);
           break;
         }
       }
     }
-    
-    //ZEM has only 1 quadrant
-    if(vol[0] == 3){
+    else if(vol[0] == 3){      //ZEM has only 1 quadrant
       vol[1] = 1;
       xdet[0] = x[0]-fPosZEM[0];
       xdet[1] = x[1]-fPosZEM[1];
@@ -1551,11 +1386,21 @@ void AliZDCv1::StepManager()
        gMC->TrackMomentum(p);
        Float_t ptot=TMath::Sqrt(p[0]*p[0]+p[1]*p[1]+p[2]*p[2]);
        Float_t beta =  ptot/p[3];
-       if(beta<0.67) return;
-       if((beta>=0.67) && (beta<=0.75)) ibeta = 0;
-       if((beta>0.75)  && (beta<=0.85)) ibeta = 1;
-       if((beta>0.85)  && (beta<=0.95)) ibeta = 2;
-       if(beta>0.95)   ibeta = 3;
+       if(beta<0.67){
+         return;
+       }
+       else if((beta>=0.67) && (beta<=0.75)){
+         ibeta = 0;
+       }
+       if((beta>0.75)  && (beta<=0.85)){
+         ibeta = 1;
+       }
+       if((beta>0.85)  && (beta<=0.95)){
+         ibeta = 2;
+       }
+       if(beta>0.95){
+         ibeta = 3;
+       }
  
        // Angle between particle trajectory and fibre axis
        // 1 -> Momentum directions
@@ -1583,15 +1428,18 @@ void AliZDCv1::StepManager()
          be = TMath::Abs(ud[0]);
        }
  
-       if((vol[0]==1)) radius = fFibZN[1];
-       if((vol[0]==2)) radius = fFibZP[1];
+       if((vol[0]==1)){
+         radius = fFibZN[1];
+       }
+       else if((vol[0]==2)){
+         radius = fFibZP[1];
+       }
        ibe = Int_t(be*1000.+1);
  
        //Looking into the light tables 
        Float_t charge = gMC->TrackCharge();
        
-       // (1)  ZN
-       if((vol[0]==1)) {
+       if((vol[0]==1)) {       // (1)  ZN fibres
          if(ibe>fNben) ibe=fNben;
          out =  charge*charge*fTablen[ibeta][ialfa][ibe];
         nphe = gRandom->Poisson(out);
@@ -1610,9 +1458,7 @@ void AliZDCv1::StepManager()
           AddHit(gAlice->CurrentTrack(), vol, hits);
         }
        } 
-       
-       // (2) ZP
-       if((vol[0]==2)) {
+       else if((vol[0]==2)) {  // (2) ZP fibres
          if(ibe>fNbep) ibe=fNbep;
          out =  charge*charge*fTablep[ibeta][ialfa][ibe];
         nphe = gRandom->Poisson(out);
@@ -1631,11 +1477,24 @@ void AliZDCv1::StepManager()
           AddHit(gAlice->CurrentTrack(), vol, hits);
         }
        } 
-       // (3) ZEM
-       if((vol[0]==3)) {
+       else if((vol[0]==3)) {  // (3) ZEM fibres
          if(ibe>fNbep) ibe=fNbep;
          out =  charge*charge*fTablep[ibeta][ialfa][ibe];
+        gMC->TrackPosition(s);
+         for(j=0; j<=2; 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];
+//      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]);
+//         printf("\n  xalic[0] = %f   xalic[1] = %f   xalic[2] = %f   z = %f  \n",
+//             xalic[0],xalic[1],xalic[2],z);
+        out = out*GuiEff;
         nphe = gRandom->Poisson(out);
+//         printf("    out*GuiEff = %f nphe = %d", out, nphe);
 //      printf("ZEM --- ibeta = %d, ialfa = %d, ibe = %d"
 //             "       -> out = %f, nphe = %d\n", ibeta, ialfa, ibe, out, nphe);
         hits[7] = 0;   
index 5293cf182aa1285ca40e74bfb17566b00da4d104..8b42df0ba64a9656cb6546881c0a65e296c1003b 100644 (file)
@@ -17,31 +17,17 @@ class AliZDCv1 : public AliZDC {
 public:
   AliZDCv1();
   AliZDCv1(const char *name, const char *title);
-  virtual      ~AliZDCv1() {}
+  virtual ~AliZDCv1() {}
   virtual void  CreateGeometry();
   virtual void  CreateBeamLine();
   virtual void  CreateZDC();
   virtual void  CreateMaterials();
-  Int_t         Digitize(Int_t Det, Int_t Quad, Int_t Light);
-  virtual void  SDigits2Digits();
-  virtual void  MakeBranch(Option_t* opt, const char *file=0);
-  virtual void  Hits2Digits(Int_t ntracks = 0);
   virtual Int_t IsVersion() const {return 1;}
   virtual void  DrawModule();
   virtual void  Init();
   virtual void  InitTables();
   virtual void  StepManager();
   
-//  // Digitization parameters setters and getters
-//  // PM gain
-//  void SetPMGain(Int_t Det, Int_t PMDet, Int_t PMGain)
-//       {fPMGain[Det][PMDet] = PMGain;}
-//  Float_t GetPMGain(Int_t Det, Int_t PMDet)
-//       {return fPMGain[Det][PMDet];}
-//  // Conversion factor from charge to ADC channels
-//  //   F = 1.6E-19 / Resolution [Coulomb/ch]
-//  void SetADCRes(Int_t ADCRes) {fADCRes =  ADCRes;}
-//  Float_t GetADCRes() {return fADCRes;}
  
 protected:
 
@@ -52,8 +38,8 @@ protected:
   Int_t   fMedSensZN;         // Sensitive medium for ZN
   Int_t   fMedSensZEM;        // Sensitive medium for EM ZDC
   Int_t   fMedSensGR;         // Other sensitive medium
-  Int_t   fMedSensPI;         // Beam pipe and magnet coils
-  Int_t   fMedSensTDI;        // TDI Cu shielding 
+//  Int_t   fMedSensPI;         // Beam pipe and magnet coils
+//  Int_t   fMedSensTDI;        // TDI Cu shielding 
   
   // Parameters for light tables
   Int_t   fNalfan;           // Number of Alfa (neutrons)
@@ -76,15 +62,12 @@ protected:
   // NB -> parameters used in CreateZDC() and in StepManager()
   // (other parameters are defined in CreateZDC())
   Float_t fPosZEM[3]; // Position of EM detector
+  Float_t fZEMLength; // ZEM length
   
   // Parameters for tracking studies
 //  Int_t fpLostIT, fpLostD1, fpLostTDI, fpDetected;
   
   
-//  // Parameters for conversion of light yield in ADC channels
-//  Float_t fPMGain[3][5];      // PM gain
-//  Float_t fADCRes;            // ADC conversion factor
-  
    ClassDef(AliZDCv1,1)  // Zero Degree Calorimeter version 1
 }; 
  
index 8bc9f2a0b18b90569254e0070e6dc526deb23c50..728de3c3a1134b7fd2e8f6589d6eb0b7c0214c29 100644 (file)
@@ -15,31 +15,17 @@ class AliZDCv2 : public AliZDC {
 public:
   AliZDCv2();
   AliZDCv2(const char *name, const char *title);
-  virtual      ~AliZDCv2() {}
+  virtual  ~AliZDCv2() {}
   virtual void  CreateGeometry();
   virtual void  CreateBeamLine();
   virtual void  CreateZDC();
   virtual void  CreateMaterials();
-  Int_t         Digitize(Int_t Det, Int_t Quad, Int_t Light);
-  virtual void  SDigits2Digits();
-  virtual void  Hits2Digits(Int_t ntracks = 0);
-  virtual void  MakeBranch(Option_t* opt, const char *file=0);
   virtual Int_t IsVersion() const {return 1;}
   virtual void  DrawModule();
   virtual void  Init();
   virtual void  InitTables();
   virtual void  StepManager();
   
-//  // Digitization parameters setters and getters
-//  // PM gain
-//  void SetPMGain(Int_t Det, Int_t PMDet, Int_t PMGain)
-//       {fPMGain[Det][PMDet] = PMGain;}
-//  Float_t GetPMGain(Int_t Det, Int_t PMDet)
-//       {return fPMGain[Det][PMDet];}
-//  // Conversion factor from charge to ADC channels
-//  //   F = 1.6E-19 / Resolution [Coulomb/ch]
-//  void SetADCRes(Int_t ADCRes) {fADCRes =  ADCRes;}
-//  Float_t GetADCRes() {return fADCRes;}
  
 protected:
 
@@ -74,15 +60,11 @@ protected:
   // NB -> parameters used in CreateZDC() and in StepManager()
   // (other parameters are defined in CreateZDC())
   Float_t fPosZEM[3]; // Position of EM detector
+  Float_t fZEMLength; // ZEM length
   
   // Parameters for tracking studies
 //  Int_t fpLostIT, fpLostD1, fpLostTDI, fpDetected;
   
-  
-//  // Parameters for conversion of light yield in ADC channels
-//  Float_t fPMGain[3][5];      // PM gain
-//  Float_t fADCRes;            // ADC conversion factor
-  
    ClassDef(AliZDCv2,1)  // Zero Degree Calorimeter version 1
 };