Moving to the new VMC naming convention
[u/mrichter/AliRoot.git] / EVGEN / AliGenDoubleScan.cxx
CommitLineData
df725edf 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
88cb7938 16/* $Id$ */
df725edf 17
675e9664 18// As AliGenScan, generation of particles on a 3-dim grid
19// but here double hits with a predefined distance are generated.
20// The second particle is generated at a constant distance but with random phi.
21// Generator can be used to evaluate double hit resolutions.
22// Author: andreas.morsch@cern.ch
23
df725edf 24#include "AliGenDoubleScan.h"
df725edf 25#include "AliRun.h"
f87cfe57 26
df725edf 27 ClassImp(AliGenDoubleScan)
28
29 AliGenDoubleScan::AliGenDoubleScan()
30 :AliGenScan(-1)
31{
32}
33
34AliGenDoubleScan::AliGenDoubleScan(Int_t npart)
35 :AliGenScan(npart)
36{
f87cfe57 37// Constructor
8b31bfac 38 fName = "Double Scan";
39 fTitle= "Particle Generator for two correlated particles on a grid";
df725edf 40}
41
42//____________________________________________________________
43AliGenDoubleScan::~AliGenDoubleScan()
f87cfe57 44{
45// Destructor
46}
df725edf 47
48//____________________________________________________________
49void AliGenDoubleScan::Generate()
50{
f87cfe57 51 //
52 // Generate one trigger
53 //
df725edf 54
f87cfe57 55 Float_t polar[3]= {0,0,0};
56 //
57 Float_t origin[3];
58 Float_t p[3];
59 Int_t nt;
60 Float_t pmom, theta, phi;
61 //
62 Float_t random[6];
63 Float_t dx,dy,dz;
64
65 //
833149de 66 if (fNx > 0) {
38f1bd58 67 dx=(fXCmax-fXCmin)/fNx;
f87cfe57 68 } else {
69 dx=1e10;
70 }
df725edf 71
f87cfe57 72 if (fNy > 0) {
38f1bd58 73 dy=(fYCmax-fYCmin)/fNy;
f87cfe57 74 } else {
75 dy=1e10;
76 }
77
78 if (fNz > 0) {
df725edf 79 dz=(fZmax-fZmin)/fNz;
f87cfe57 80 } else {
81 dz=1e10;
82 }
83 for (Int_t ix=0; ix<fNx; ix++) {
df725edf 84 for (Int_t iy=0; iy<fNy; iy++) {
85 for (Int_t iz=0; iz<fNz; iz++){
65fb704d 86 Rndm(random,6);
38f1bd58 87 origin[0]=fXCmin+ix*dx+2*(random[0]-0.5)*fOsigma[0];
88 origin[1]=fYCmin+iy*dy+2*(random[1]-0.5)*fOsigma[1];
df725edf 89 origin[2]=fZmin+iz*dz+2*(random[2]-0.5)*fOsigma[2];
90 pmom=fPMin+random[3]*(fPMax-fPMin);
91 theta=fThetaMin+random[4]*(fThetaMax-fThetaMin);
92 phi=fPhiMin+random[5]*(fPhiMax-fPhiMin);
93 p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
94 p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
95 p[2] = pmom*TMath::Cos(theta);
642f15cf 96 PushTrack(fTrackIt,-1,fIpart,p,origin,polar,0,kPPrimary,nt);
df725edf 97//
98// Generate 2nd particle at distance fDistance from the first
99//
65fb704d 100 Rndm(random,6);
df725edf 101 Float_t phi2=2.*TMath::Pi()*random[0];
102 Float_t dx =fDistance*TMath::Sin(phi2);
103 Float_t dy =fDistance*TMath::Cos(phi2);
104 origin[0]=origin[0]+dx;
105 origin[1]=origin[1]+dy;
106 pmom=fPMin+random[1]*(fPMax-fPMin);
107 theta=fThetaMin+random[2]*(fThetaMax-fThetaMin);
108 phi=fPhiMin+random[3]*(fPhiMax-fPhiMin);
109 p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
110 p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
111 p[2] = pmom*TMath::Cos(theta);
642f15cf 112 PushTrack(fTrackIt,-1,fIpart,p,origin,polar,0,kPPrimary,nt);
df725edf 113 }
114 }
115 }
116}
117
118
119
120
121
122
123
124
125
126
127
128
129