Changing once more (hopefully we get it correct this time...) the logic to trig the...
[u/mrichter/AliRoot.git] / ITS / AliITSsimulationFastPoints.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 /* $Id$ */
17 //////////////////////////////////////////////////////////
18 // implements fast simulation                           //
19 //                                                      //
20 //                                                      //
21 //////////////////////////////////////////////////////////
22
23
24 #include <TRandom.h>
25
26 #include "AliITS.h"
27 #include "AliITShit.h"
28 #include "AliITSRecPoint.h"
29 #include "AliITSmodule.h"
30 #include "AliITSgeom.h"
31 #include "AliRun.h"
32 #include "AliITSsimulationFastPoints.h"
33
34
35 ClassImp(AliITSsimulationFastPoints)
36
37 AliITSsimulationFastPoints::AliITSsimulationFastPoints()
38 {
39   //constructor
40   fSigmaRPhi[0] = fSigmaRPhi[1] = 12e-4;
41   fSigmaRPhi[2] = fSigmaRPhi[3] = 38e-4;
42   fSigmaRPhi[4] = fSigmaRPhi[5] = 20e-4;
43   fSigmaZ[0] = fSigmaZ[1] = 120e-4;        // resolution for 425 micron pixels
44   fSigmaZ[2] = fSigmaZ[3] = 28e-4;
45   fSigmaZ[4] = fSigmaZ[5] = 830e-4;
46   fSigmaDe[0] = fSigmaDe[1] = 0.72e-6;
47   fSigmaDe[2] = fSigmaDe[3] = 0.90e-6;
48   fSigmaDe[4] = fSigmaDe[5] =  5e-6;
49   fThrDe[0] = fThrDe[1] = 7.2e-6;
50   fThrDe[2] = fThrDe[3] = 2.70e-6;
51   fThrDe[4] = fThrDe[5] = 10e-6;
52 }
53
54 //-------------------------------------------------------------
55 void AliITSsimulationFastPoints::CreateFastRecPoints(Int_t module, TClonesArray* recp){
56     // Fast points simulator
57     AliITS *aliITS  = (AliITS*)gAlice->GetModule("ITS");
58
59     CreateFastRecPoints((AliITSmodule *)(aliITS->GetModule(module)),
60                         module,gRandom,recp);
61 }
62 //-------------------------------------------------------------
63 void AliITSsimulationFastPoints::CreateFastRecPoints(AliITSmodule *mod,
64                                                      Int_t module,
65                                                      TRandom *random,
66                                                      TClonesArray* recp) {
67   // Fast points simulator 
68
69   TClonesArray &pt=*recp;
70   AliITS *aliITS  = (AliITS*)gAlice->GetModule("ITS");
71   AliITSgeom *gm = aliITS->GetITSgeom();
72   const Float_t kdEdXtoQ = 1.0e+6;  // GeV->KeV
73
74   Int_t lay,lad,det;
75   gm->GetModuleId(module,lay,lad,det);
76   Int_t ind=(lad-1)*gm->GetNdetectors(lay)+(det-1);
77   Int_t lyr=(lay-1);
78
79
80   Int_t ihit,flag,numofhits;
81   Float_t locals[3];
82   Float_t globals[3];
83   Double_t sigmarphi=0., sigmaz=0., sigmade=0., thrde=0.;
84   Float_t deltaXl,deltaZl,deltaDe;
85
86   Int_t hitlay, hitlad, hitdet, hitstatus;
87   Float_t hitpx, hitpy, hitpz, hitdestep;
88
89   Int_t   hitstatus1, hittrack1;
90   Float_t hitx1, hity1, hitz1;
91   Float_t hitdestep1;
92   Float_t xMg,yMg,zMg;
93   Int_t irecp=0;
94   numofhits = mod->GetNhits();
95   //printf("numofhits %d \n",numofhits);
96   for(ihit=0;ihit<numofhits;ihit++){
97     AliITShit *hit=mod->GetHit(ihit);
98     hit->GetPositionG(hitx1,hity1,hitz1);
99     hitstatus1 = hit->GetTrackStatus();
100     hitdestep1 = hit->GetIonization();
101     hittrack1 = hit->GetTrack();
102     
103     mod->MedianHit(module,hitx1,hity1,hitz1,hitstatus1,xMg,yMg,zMg,flag);
104     if (flag!=1) {
105       hitdestep = hit->GetIonization();
106       
107       if (hitdestep > 0) {
108         hit->GetDetectorID(hitlay,hitlad,hitdet);
109         hit->GetMomentumG(hitpx,hitpy,hitpz);            
110         hitstatus = hitstatus1;
111                 // Transform to the module local frame
112         globals[0] = xMg; 
113         globals[1] = yMg;
114         globals[2] = zMg;
115         gm->GtoL(hitlay,hitlad,hitdet,globals,locals);
116         // Retrieve sigma values for position and energy, and energy
117         // threshold
118         sigmarphi = SigmaRPhi(hitlay);
119         sigmaz = SigmaZ(hitlay);
120         sigmade = SigmaDe(hitlay);
121         thrde = ThrDe(hitlay);
122         deltaXl = random->Gaus(0,sigmarphi);
123         deltaZl = random->Gaus(0,sigmaz);
124         deltaDe = random->Gaus(0,sigmade);
125         
126         // Apply energy threshold and trasform back to global reference
127         // system
128         
129         if ( (hitdestep+deltaDe) > thrde ){
130           locals[0] += deltaXl;
131           locals[2] += deltaZl;
132           Int_t lab[4] = {hit->GetTrack(),-3,-3,ind};
133           Float_t q=kdEdXtoQ*(hitdestep+deltaDe);
134           if(hitlay<3) q=1.; // SPD binary readout
135           Float_t hitv[5] = {locals[0],locals[2],sigmarphi*sigmarphi,sigmaz*sigmaz,q};
136           Int_t info[3] = {0,0,lyr};
137           AliITSRecPoint rp(lab,hitv,info,kTRUE);
138
139           new (pt[irecp]) AliITSRecPoint(rp);
140           irecp++;
141         } // end if ( (hitdestep+deltaDe)
142       } // end if (hitdestep > 0)
143     } // end if (flag!=1)
144   } // end for ihit
145 }
146 //_______________________________________________________________________
147 void AliITSsimulationFastPoints::SetSigmaRPhi(Double_t  srphi[6])
148 {
149   // set sigmas in rphi
150
151     Int_t i;
152     for (i=0; i<6; i++) {
153         fSigmaRPhi[i]=srphi[i];
154     }
155 }
156 //_______________________________________________________________________
157 void AliITSsimulationFastPoints::SetSigmaZ(Double_t  sz[6])
158 {
159   // set sigmas in z
160
161     Int_t i;
162     for (i=0; i<6; i++) {
163         fSigmaZ[i]=sz[i];
164     }
165 }
166 //_______________________________________________________________________
167 void AliITSsimulationFastPoints::SetSigmaDe(Double_t  sde[6])
168 {
169   // set sigmas in energy
170
171     Int_t i;
172     for (i=0; i<6; i++) {
173         fSigmaDe[i]=sde[i];
174     }
175 }
176 //_______________________________________________________________________
177 void AliITSsimulationFastPoints::SetThrDe(Double_t  thrde[6])
178 {
179   // set energy thersholds
180
181     Int_t i;
182     for (i=0; i<6; i++) {
183         fThrDe[i]=thrde[i];
184     }
185 }
186