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 "AliMUONChamberTrigger.h"
19 #include "AliMUONResponseTrigger.h"
20 #include "AliMUONHit.h"
22 #include "AliMUONSegmentation.h"
23 #include "AliMUONGeometrySegmentation.h"
24 #include "AliMUONGeometryTransformer.h"
28 /// \class AliMUONChamberTrigger
30 /// Implementation of AliMUONChamber for the trigger
32 /// This class is to be deprecated.
36 ClassImp(AliMUONChamberTrigger)
39 //-------------------------------------------
41 AliMUONChamberTrigger::AliMUONChamberTrigger()
45 // Default constructor
48 //-------------------------------------------
50 AliMUONChamberTrigger::AliMUONChamberTrigger(const AliMUONChamberTrigger& right)
51 : AliMUONChamber(right)
53 /// Protected copy constructor (not implemented)
55 AliFatal("Copy constructor not provided.");
58 //-------------------------------------------
60 AliMUONChamberTrigger::AliMUONChamberTrigger(Int_t id,
61 const AliMUONGeometryTransformer* kGeometryTransformer)
63 fkGeomTransformer(kGeometryTransformer)
65 // Constructor using chamber id
68 //-------------------------------------------
69 AliMUONChamberTrigger&
70 AliMUONChamberTrigger::operator=(const AliMUONChamberTrigger& right)
72 /// Protected assignement operator (not implemented)
74 // check assignement to self
75 if (this == &right) return *this;
77 AliFatal("Assignement operator not provided.");
82 //-------------------------------------------
83 void AliMUONChamberTrigger::DisIntegration(AliMUONHit* hit,
85 Float_t newclust[6][500])
88 // Generates pad hits (simulated cluster)
89 // using the segmentation and the response model
92 Float_t tof = hit->Age();
93 Float_t xhit = hit->X();
94 Float_t yhit = hit->Y();
95 Float_t zhit = hit->Z();
96 Int_t id = hit->DetElemId();
99 if (tof<75*TMath::Power(10,-9)) {
107 for (Int_t i = 1; i <= 2; i++) {
109 AliMUONGeometrySegmentation* segmentation=
110 fMUON->GetSegmentation()->GetModuleSegmentation(fId, i-1);
113 // Find the module & strip Id. which has fired
116 segmentation->GetPadI(id,xhit,yhit,0,ix,iy);
117 // treatment of GEANT hits w/o corresponding strip (due to the fact that
118 // geometry & segmentation are computed in a very slightly different way)
122 fkGeomTransformer->Global2Local(id,xhit,yhit,0,lx,ly,lz);
123 AliWarning(Form("AliMUONChamberTrigger hit w/o strip %i-%d %e %e "
124 "local %e %e %e ix,iy=%d,%d\n",id,i-1,xhit,yhit,lx,ly,lz,ix,iy));
127 segmentation->SetPad(id,ix,iy);
128 if (xhit<0) ix = -ix;
129 // printf(" fId id fnsec xhit yhit zhit ix iy %i %i %i %f %f %f %i %i \n",fId,i,id,xhit,yhit,zhit,ix,iy);
130 // if (ix < 0 || ix > 10000) return;
131 // if (iy < 0 || iy > 10000) return;
133 // --- store signal information for this strip
134 newclust[0][nnew]=1.; // total charge
135 newclust[1][nnew]=ix; // ix-position of pad
136 newclust[2][nnew]=iy; // iy-position of pad
137 newclust[3][nnew]=twentyNano; // time of flight
138 newclust[4][nnew]=segmentation->ISector(); // sector id
139 newclust[5][nnew]=(Float_t) i; // counter
142 // cluster-size if AliMUONResponseTriggerV1, nothing if AliMUONResponseTrigger
143 if (((AliMUONResponseTrigger*) fResponse)->SetGenerCluster()) {
146 segmentation->SetHit(id,xhit,yhit,zhit);
147 // get the list of nearest neighbours
148 Int_t nList, xList[10], yList[10];
149 segmentation->Neighbours(id,ix,iy,&nList,xList,yList);
152 for (Int_t j=0; j<nList; j++){ // loop over neighbours
153 if (xList[j]!=0) { // existing neighbour
154 if (j==0||j==5||qp!=0) { // built-up cluster-size
156 // neighbour real coordinates (just for checks here)
158 segmentation->GetPadC(id,xList[j],yList[j],x,y,z);
159 // set pad (fx fy & fix fiy are the current pad coord. & Id.)
160 segmentation->SetPad(id,xList[j],yList[j]);
161 // get the chamber (i.e. current strip) response
162 qp=fResponse->IntXY(id,segmentation);
165 // --- store signal information for neighbours
166 newclust[0][nnew]=qp; // total charge
167 newclust[1][nnew]=segmentation->Ix(); // ix-pos. of pad
168 newclust[2][nnew]=segmentation->Iy(); // iy-pos. of pad
169 newclust[3][nnew]=twentyNano; // time of flight
170 newclust[4][nnew]=segmentation->ISector(); // sector id
171 newclust[5][nnew]=(Float_t) i; // counter
174 } // built-up cluster-size
175 } // existing neighbour
176 } // loop over neighbours
177 } // endif hit w/o strip
178 } // loop over planes
179 } // if AliMUONResponseTriggerV1