EVGEN/AliGenDoubleScan.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 Revision 1.6  2000/12/21 16:24:06  morsch
  19 Coding convention clean-up
  20
  21 Revision 1.5  2000/12/06 15:11:38  morsch
  22 Correct double declared data members.
  23
  24 Revision 1.4  2000/11/30 07:12:50  alibrary
  25 Introducing new Rndm and QA classes
  26
  27 Revision 1.3  2000/10/02 21:28:06  fca
  28 Removal of useless dependecies via forward declarations
  29
  30 Revision 1.2  2000/06/09 20:37:51  morsch
  31 All coding rule violations except RS3 corrected
  32
  33 Revision 1.1  2000/02/23 16:25:14  morsch
  34 First commit of this file
  35
  36 */
  37
  38 // As AliGenScan,  generation of particles on a 3-dim grid
  39 // but here double hits with a predefined distance are generated.
  40 // The second particle is generated at a constant distance but with random phi.
  41 // Generator can be used to evaluate double hit resolutions.
  42 // Author: andreas.morsch@cern.ch
  43
  44 #include "AliGenDoubleScan.h"
  45 #include "AliRun.h"
  46
  47  ClassImp(AliGenDoubleScan)
  48
  49  AliGenDoubleScan::AliGenDoubleScan()
  50          :AliGenScan(-1)
  51 {
  52 }
  53
  54 AliGenDoubleScan::AliGenDoubleScan(Int_t npart)
  55     :AliGenScan(npart)
  56 {
  57 // Constructor
  58 }
  59
  60 //____________________________________________________________
  61 AliGenDoubleScan::~AliGenDoubleScan()
  62 {
  63 // Destructor
  64 }
  65
  66 //____________________________________________________________
  67 void AliGenDoubleScan::Generate()
  68 {
  69     //
  70     // Generate one trigger
  71     //
  72
  73     Float_t polar[3]= {0,0,0};
  74     //
  75     Float_t origin[3];
  76     Float_t p[3];
  77     Int_t nt;
  78     Float_t pmom, theta, phi;
  79     //
  80     Float_t random[6];
  81     Float_t dx,dy,dz;
  82
  83     //
  84     if (fNy > 0) {
  85         dx=(fXCmax-fXCmin)/fNx;
  86     } else {
  87         dx=1e10;
  88     }
  89
  90     if (fNy > 0) {
  91         dy=(fYCmax-fYCmin)/fNy;
  92     } else {
  93         dy=1e10;
  94     }
  95
  96     if (fNz > 0) {
  97       dz=(fZmax-fZmin)/fNz;
  98     } else {
  99         dz=1e10;
 100     }
 101     for (Int_t ix=0; ix<fNx; ix++) {
 102       for (Int_t iy=0; iy<fNy; iy++) {
 103           for (Int_t iz=0; iz<fNz; iz++){
 104               Rndm(random,6);
 105               origin[0]=fXCmin+ix*dx+2*(random[0]-0.5)*fOsigma[0];
 106               origin[1]=fYCmin+iy*dy+2*(random[1]-0.5)*fOsigma[1];
 107               origin[2]=fZmin+iz*dz+2*(random[2]-0.5)*fOsigma[2];
 108               pmom=fPMin+random[3]*(fPMax-fPMin);
 109               theta=fThetaMin+random[4]*(fThetaMax-fThetaMin);
 110               phi=fPhiMin+random[5]*(fPhiMax-fPhiMin);
 111               p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
 112               p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
 113               p[2] = pmom*TMath::Cos(theta);
 114               SetTrack(fTrackIt,-1,fIpart,p,origin,polar,0,kPPrimary,nt);
 115 //
 116 // Generate 2nd particle at distance fDistance from  the first
 117 //
 118               Rndm(random,6);
 119               Float_t phi2=2.*TMath::Pi()*random[0];
 120               Float_t dx  =fDistance*TMath::Sin(phi2);
 121               Float_t dy  =fDistance*TMath::Cos(phi2);
 122               origin[0]=origin[0]+dx;
 123               origin[1]=origin[1]+dy;
 124               pmom=fPMin+random[1]*(fPMax-fPMin);
 125               theta=fThetaMin+random[2]*(fThetaMax-fThetaMin);
 126               phi=fPhiMin+random[3]*(fPhiMax-fPhiMin);
 127               p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
 128               p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
 129               p[2] = pmom*TMath::Cos(theta);
 130               SetTrack(fTrackIt,-1,fIpart,p,origin,polar,0,kPPrimary,nt);
 131           }
 132       }
 133   }
 134 }
 135
 136
 137
 138
 139
 140
 141
 142
 143
 144
 145
 146
 147