MUON/AliMUONchamber.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /*
  17 $Log$
  18 */
  19
  20 #include "AliMUONchamber.h"
  21 #include "TMath.h"
  22 ClassImp(AliMUONchamber)
  23     AliMUONchamber::AliMUONchamber()
  24 {
  25     fSegmentation = new TObjArray(2);
  26     fResponse=0;
  27     fnsec=1;
  28 }
  29
  30 void AliMUONchamber::Init()
  31 {
  32
  33     ((AliMUONsegmentation *) (*fSegmentation)[0])->Init(this);
  34     if (fnsec==2) {
  35         ((AliMUONsegmentation *) (*fSegmentation)[1])->Init(this);
  36     }
  37
  38 }
  39
  40 void AliMUONchamber::DisIntegration(Float_t eloss, Float_t xhit, Float_t yhit,
  41                                     Int_t& nnew,Float_t newclust[6][500])
  42 {
  43 //
  44 //  Generates pad hits (simulated cluster)
  45 //  using the segmentation and the response model
  46     Float_t dx, dy;
  47     //
  48     // Width of the integration area
  49     //
  50     dx=fResponse->SigmaIntegration()*fResponse->ChargeSpreadX();
  51     dy=fResponse->SigmaIntegration()*fResponse->ChargeSpreadY();
  52     //
  53     // Get pulse height from energy loss
  54     Float_t qtot = fResponse->IntPH(eloss);
  55     //
  56     // Loop Over Pads
  57
  58     Float_t qcheck=0, qp;
  59     nnew=0;
  60     for (Int_t i=1; i<=fnsec; i++) {
  61         qcheck=0;
  62         AliMUONsegmentation * segmentation=(AliMUONsegmentation *) (*fSegmentation)[i-1];
  63         for (segmentation->FirstPad(xhit, yhit, dx, dy);
  64              segmentation->MorePads();
  65              segmentation->NextPad())
  66         {
  67             qp=fResponse->IntXY(segmentation);
  68             qp=TMath::Abs(qp);
  69
  70 //
  71 //
  72             if (qp > 1.e-4) {
  73                 qcheck+=qp;
  74             //
  75             // --- store signal information
  76                 newclust[0][nnew]=qtot;
  77                 newclust[1][nnew]=segmentation->Ix();
  78                 newclust[2][nnew]=segmentation->Iy();
  79                 newclust[3][nnew]=qp * qtot;
  80                 newclust[4][nnew]=segmentation->ISector();
  81                 newclust[5][nnew]=(Float_t) i;
  82 //              printf("\n pad hit %d %d %f %f \n",nnew,i,newclust[1][nnew],newclust[2][nnew]);
  83                 nnew++;
  84
  85
  86             }
  87         } // Pad loop
  88 //      printf("\n check sum is %f %f %f %f %d \n",qcheck,qtot,xhit,yhit,fGid);
  89     } // Cathode plane loop
  90 }
  91
  92
  93
  94     void AliMUONchamber::InitGeo(Float_t)
  95 {
  96 //    sensitive gas gap
  97       fdGas= 0.5;
  98 //    3% radiation length of aluminum (X0=8.9 cm)
  99       fdAlu= 3.0/100*8.9;
 100 }