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 **************************************************************************/
18 #include "AliMUONSegmentationTrigger.h"
19 #include "AliMUONTriggerConstants.h"
25 #include "AliMUONChamber.h"
26 #include <Riostream.h>
28 ClassImp(AliMUONSegmentationTrigger)
30 //------------------------------------------------------------------
31 void AliMUONSegmentationTrigger::Init(Int_t chamber)
33 // initialize Module geometry
34 AliMUON *pMUON = (AliMUON *) gAlice->GetModule("MUON");
35 AliMUONChamber* iChamber=&(pMUON->Chamber(chamber));
37 if(pMUON->GetDebug()>1) printf("%s: Initialize Trigger Chamber Module Geometry\n",ClassName());
39 Float_t zPos=iChamber->Z();
47 Float_t z1PosPlus=z1Pos+dz/2.;
48 Float_t z1PosMinus=z1Pos-dz/2.;
50 Float_t z1pm=z1PosPlus/z1PosMinus;
51 Float_t z1mp=z1PosMinus/z1PosPlus;
53 if(pMUON->GetDebug()>1) printf("%s: fZscale = %f \n",ClassName(),fZscale);
55 // calculate yCmin and fYcmax
57 for (i=62; i>=0; i--) {
58 Int_t j=ModuleNumber(-AliMUONTriggerConstants::ModuleId(i)); // i == right, j == left
59 if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==5) { // start with middle chamber
60 if (AliMUONTriggerConstants::ModuleId(i)==51) { // special case (empty module)
61 fYcmin[i]=fYcmax[i]=fYcmin[j]=fYcmax[j]=0.;
63 y1Cmin[i]=y1Cmin[j]=-34;
64 y1Cmax[i]=y1Cmax[j]=34;
65 fYcmin[i]=fYcmin[j]=-34.;
67 fYcmax[i]=fYcmax[j]=34.;
69 } else if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==4) { // up
70 if (AliMUONTriggerConstants::ModuleId(i)!=41) {
71 y1Cmin[i]=y1Cmax[i+7]*z1pm;
72 y1Cmax[i]=y1Cmin[i]+68.;
74 fYcmax[i]=fYcmin[i]+68.;
76 y1Cmin[j]=y1Cmax[j+7]*z1mp;
77 y1Cmax[j]=y1Cmin[j]+68.;
79 fYcmax[j]=fYcmin[j]+68.;
81 y1Cmin[i]=y1Cmin[ModuleNumber(42)]+17;
82 y1Cmax[i]=y1Cmin[i]+51.;
84 fYcmax[i]=fYcmin[i]+51.;
86 y1Cmin[j]=y1Cmin[ModuleNumber(-42)]+17;
87 y1Cmax[j]=y1Cmin[j]+51.;
89 fYcmax[j]=fYcmin[j]+51.;
91 } else if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==3) {
92 y1Cmin[i]=y1Cmax[i+7]*z1mp;
93 y1Cmax[i]=y1Cmin[i]+68.;
95 fYcmax[i]=fYcmin[i]+68.;
97 y1Cmin[j]=y1Cmax[j+7]*z1pm;
98 y1Cmax[j]=y1Cmin[j]+68.;
100 fYcmax[j]=fYcmin[j]+68.;
101 } else if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==2) {
102 y1Cmin[i]=y1Cmax[i+7]*z1pm;
103 y1Cmax[i]=y1Cmin[i]+68.;
105 fYcmax[i]=fYcmin[i]+68.;
107 y1Cmin[j]=y1Cmax[j+7]*z1mp;
108 y1Cmax[j]=y1Cmin[j]+68.;
110 fYcmax[j]=fYcmin[j]+68.;
111 } else if (Int_t(AliMUONTriggerConstants::ModuleId(i)/10)==1) {
112 y1Cmin[i]=y1Cmax[i+7]*z1mp;
113 y1Cmax[i]=y1Cmin[i]+68.;
115 fYcmax[i]=fYcmin[i]+68.;
117 y1Cmin[j]=y1Cmax[j+7]*z1pm;
118 y1Cmax[j]=y1Cmin[j]+68.;
120 fYcmax[j]=fYcmin[j]+68.;
124 for (i=0; i<63; i++) { // second loop (fill lower part)
125 Int_t j=ModuleNumber(-AliMUONTriggerConstants::ModuleId(i)); // i == right, j == left
126 if (TMath::Abs(Int_t(AliMUONTriggerConstants::ModuleId(i)/10))==6) {
127 fYcmin[i]=-fYcmax[i-14];
128 fYcmax[i]=-fYcmin[i-14];
129 fYcmin[j]=-fYcmax[j-14];
130 fYcmax[j]=-fYcmin[j-14];
131 } else if (TMath::Abs(Int_t(AliMUONTriggerConstants::ModuleId(i)/10))==7) {
132 fYcmin[i]=-fYcmax[i-28];
133 fYcmax[i]=-fYcmin[i-28];
134 fYcmin[j]=-fYcmax[j-28];
135 fYcmax[j]=-fYcmin[j-28];
136 } else if (TMath::Abs(Int_t(AliMUONTriggerConstants::ModuleId(i)/10))==8) {
137 fYcmin[i]=-fYcmax[i-42];
138 fYcmax[i]=-fYcmin[i-42];
139 fYcmin[j]=-fYcmax[j-42];
140 fYcmax[j]=-fYcmin[j-42];
141 } else if (TMath::Abs(Int_t(AliMUONTriggerConstants::ModuleId(i)/10))==9) {
142 fYcmin[i]=-fYcmax[i-56];
143 fYcmax[i]=-fYcmin[i-56];
144 fYcmin[j]=-fYcmax[j-56];
145 fYcmax[j]=-fYcmin[j-56];
152 // Set parent chamber number
153 fChamber=&(pMUON->Chamber(chamber));
157 //------------------------------------------------------------------
158 void AliMUONSegmentationTrigger::Draw(const char * /*opt*/) const {}
159 //------------------------------------------------------------------
160 Int_t AliMUONSegmentationTrigger::ModuleNumber(Int_t imodule){
161 // returns module number (from 0 to 126) corresponding to module imodule
163 for (Int_t i=0; i<AliMUONTriggerConstants::Nmodule(); i++) {
164 if (AliMUONTriggerConstants::ModuleId(i)==imodule) {
171 //------------------------------------------------------------------
172 void AliMUONSegmentationTrigger::SetHit(Float_t xhit, Float_t yhit, Float_t /*zhit*/)
176 //------------------------------------------------------------------
177 Float_t AliMUONSegmentationTrigger::StripSizeX(Int_t imodule){
178 // Returns x-strip size for given module imodule
180 Int_t absimodule=TMath::Abs(imodule);
181 Int_t moduleNum=ModuleNumber(imodule);
182 if (absimodule==51) {
185 return TMath::Abs((fYcmax[moduleNum]-fYcmin[moduleNum])/
186 AliMUONTriggerConstants::NstripX(moduleNum));
190 //------------------------------------------------------------------
191 Float_t AliMUONSegmentationTrigger::StripSizeY(Int_t imodule){
192 // Returns y-strip size for given module imodule
194 Int_t absimodule=TMath::Abs(imodule);
195 Int_t moduleNum=ModuleNumber(imodule);
196 if (absimodule==51) {
199 return TMath::Abs((AliMUONTriggerConstants::XcMax(moduleNum) -
200 AliMUONTriggerConstants::XcMin(moduleNum)) /
201 AliMUONTriggerConstants::NstripY(moduleNum));
205 //------------------------------------------------------------------
206 void AliMUONSegmentationTrigger::SetHit(Float_t xhit, Float_t yhit)
208 // Sets virtual hit position, needed for evaluating pad response
209 // outside the tracking program