1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
19 //#include <TRandom.h>
21 //#include <Riostream.h>
23 #include "AliMUONSegmentationTrigger.h"
24 #include "AliMUONTriggerConstants.h"
25 #include "AliMUONConstants.h"
28 #include "AliMUONChamber.h"
31 ClassImp(AliMUONSegmentationTrigger)
33 //------------------------------------------------------------------
34 AliMUONSegmentationTrigger::AliMUONSegmentationTrigger()
35 : AliMUONSegmentationV0()
42 //------------------------------------------------------------------
43 void AliMUONSegmentationTrigger::Init(Int_t chamber)
45 // initialize Module geometry
46 AliMUON *pMUON = (AliMUON *) gAlice->GetModule("MUON");
47 AliMUONChamber* iChamber=&(pMUON->Chamber(chamber));
49 AliDebug(2,"Initialize Trigger Chamber Module Geometry");
51 Float_t zPos=iChamber->Z();
52 Float_t z1Pos=AliMUONConstants::DefaultChamberZ(10); //cp changed
59 Float_t z1PosPlus=z1Pos-dz/2.;
60 Float_t z1PosMinus=z1Pos+dz/2.;
62 Float_t z1pm=z1PosPlus/z1PosMinus;
63 Float_t z1mp=z1PosMinus/z1PosPlus;
65 AliDebug(2,Form("fZscale = %f ",fZscale));
66 // calculate yCmin and fYcmax
68 for (i=62; i>=0; i--) {
69 Int_t j=ModuleNumber(-AliMUONTriggerConstants::ModuleId(i)); // i == right, j == left
70 if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==5) { // start with middle chamber
71 if (AliMUONTriggerConstants::ModuleId(i)==51) { // special case (empty module)
72 fYcmin[i]=fYcmax[i]=fYcmin[j]=fYcmax[j]=0.;
74 y1Cmin[i]=y1Cmin[j]=-34;
75 y1Cmax[i]=y1Cmax[j]=34;
76 fYcmin[i]=fYcmin[j]=-34.;
78 fYcmax[i]=fYcmax[j]=34.;
80 } else if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==4) { // up
81 if (AliMUONTriggerConstants::ModuleId(i)!=41) {
82 y1Cmin[i]=y1Cmax[i+7]*z1pm;
83 y1Cmax[i]=y1Cmin[i]+68.;
85 fYcmax[i]=fYcmin[i]+68.;
87 y1Cmin[j]=y1Cmax[j+7]*z1mp;
88 y1Cmax[j]=y1Cmin[j]+68.;
90 fYcmax[j]=fYcmin[j]+68.;
92 y1Cmin[i]=y1Cmin[ModuleNumber(42)]+17;
93 y1Cmax[i]=y1Cmin[i]+51.;
95 fYcmax[i]=fYcmin[i]+51.;
97 y1Cmin[j]=y1Cmin[ModuleNumber(-42)]+17;
98 y1Cmax[j]=y1Cmin[j]+51.;
100 fYcmax[j]=fYcmin[j]+51.;
102 } else if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==3) {
103 y1Cmin[i]=y1Cmax[i+7]*z1mp;
104 y1Cmax[i]=y1Cmin[i]+68.;
106 fYcmax[i]=fYcmin[i]+68.;
108 y1Cmin[j]=y1Cmax[j+7]*z1pm;
109 y1Cmax[j]=y1Cmin[j]+68.;
111 fYcmax[j]=fYcmin[j]+68.;
112 } else if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==2) {
113 y1Cmin[i]=y1Cmax[i+7]*z1pm;
114 y1Cmax[i]=y1Cmin[i]+68.;
116 fYcmax[i]=fYcmin[i]+68.;
118 y1Cmin[j]=y1Cmax[j+7]*z1mp;
119 y1Cmax[j]=y1Cmin[j]+68.;
121 fYcmax[j]=fYcmin[j]+68.;
122 } else if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==1) {
123 y1Cmin[i]=y1Cmax[i+7]*z1mp;
124 y1Cmax[i]=y1Cmin[i]+68.;
126 fYcmax[i]=fYcmin[i]+68.;
128 y1Cmin[j]=y1Cmax[j+7]*z1pm;
129 y1Cmax[j]=y1Cmin[j]+68.;
131 fYcmax[j]=fYcmin[j]+68.;
135 for (i=0; i<63; i++) { // second loop (fill lower part)
136 Int_t j=ModuleNumber(-AliMUONTriggerConstants::ModuleId(i)); // i == right, j == left
137 if (TMath::Abs(Int_t(AliMUONTriggerConstants::ModuleId(i)/10))==6) {
138 fYcmin[i]=-fYcmax[i-14];
139 fYcmax[i]=-fYcmin[i-14];
140 fYcmin[j]=-fYcmax[j-14];
141 fYcmax[j]=-fYcmin[j-14];
142 } else if (TMath::Abs(Int_t(AliMUONTriggerConstants::ModuleId(i)/10))==7) {
143 fYcmin[i]=-fYcmax[i-28];
144 fYcmax[i]=-fYcmin[i-28];
145 fYcmin[j]=-fYcmax[j-28];
146 fYcmax[j]=-fYcmin[j-28];
147 } else if (TMath::Abs(Int_t(AliMUONTriggerConstants::ModuleId(i)/10))==8) {
148 fYcmin[i]=-fYcmax[i-42];
149 fYcmax[i]=-fYcmin[i-42];
150 fYcmin[j]=-fYcmax[j-42];
151 fYcmax[j]=-fYcmin[j-42];
152 } else if (TMath::Abs(Int_t(AliMUONTriggerConstants::ModuleId(i)/10))==9) {
153 fYcmin[i]=-fYcmax[i-56];
154 fYcmax[i]=-fYcmin[i-56];
155 fYcmin[j]=-fYcmax[j-56];
156 fYcmax[j]=-fYcmin[j-56];
163 // Set parent chamber number
164 fChamber=&(pMUON->Chamber(chamber));
167 //------------------------------------------------------------------
168 void AliMUONSegmentationTrigger::Draw(const char * /*opt*/) {}
169 //------------------------------------------------------------------
170 Int_t AliMUONSegmentationTrigger::ModuleNumber(Int_t imodule){
171 // returns module number (from 0 to 126) corresponding to module imodule
173 for (Int_t i=0; i<AliMUONTriggerConstants::Nmodule(); i++) {
174 if (AliMUONTriggerConstants::ModuleId(i)==imodule) {
181 //------------------------------------------------------------------
182 void AliMUONSegmentationTrigger::SetHit(Float_t xhit, Float_t yhit, Float_t /*zhit*/)
186 //------------------------------------------------------------------
187 Float_t AliMUONSegmentationTrigger::StripSizeX(Int_t imodule){
188 // Returns x-strip size for given module imodule
190 Int_t absimodule=TMath::Abs(imodule);
191 Int_t moduleNum=ModuleNumber(imodule);
192 if (absimodule==51) {
195 return TMath::Abs((fYcmax[moduleNum]-fYcmin[moduleNum])/
196 AliMUONTriggerConstants::NstripX(moduleNum));
200 //------------------------------------------------------------------
201 Float_t AliMUONSegmentationTrigger::StripSizeY(Int_t imodule, Int_t istrip){
202 // Returns y-strip size for given module imodule
204 Int_t absimodule=TMath::Abs(imodule);
205 Int_t moduleNum=ModuleNumber(imodule);
206 if (absimodule==51) {
208 } else if (TMath::Abs(imodule-10*Int_t(imodule/10.))==7) {
210 return AliMUONTriggerConstants::StripWidth(2);
213 return AliMUONTriggerConstants::StripWidth(1);
219 return TMath::Abs((AliMUONTriggerConstants::XcMax(moduleNum) -
220 AliMUONTriggerConstants::XcMin(moduleNum)) /
221 AliMUONTriggerConstants::NstripY(moduleNum));
225 //------------------------------------------------------------------
226 void AliMUONSegmentationTrigger::SetHit(Float_t xhit, Float_t yhit)
228 // Sets virtual hit position, needed for evaluating pad response
229 // outside the tracking program