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 **************************************************************************/
16 ////////////////////////////////////
18 /// \class AliMUONTrackHitPattern
20 /// This class propagates tracks to trigger chambers
21 /// searching for fired strips.
23 /// To each track, a hit pattern for trigger chambers is set.
24 /// The hit pattern is a UShort_t with 8 bits used:
28 /// ----------- ------------
29 /// chamber: 11 12 13 14 | 11 12 13 14
30 /// cathode: bending | non-bending
32 /// The main method is:
35 /// \author Diego Stocco
37 ////////////////////////////////////
40 #include "AliMUONTrackHitPattern.h"
41 #include "AliMUONRecData.h"
42 #include "AliMUONTrack.h"
43 #include "AliMUONTrackParam.h"
44 #include "AliMUONTrackExtrap.h"
45 #include "AliMUONConstants.h"
46 #include "AliMUONGeometryTransformer.h"
47 #include "AliMUONDigit.h"
48 #include "AliMUONLocalTrigger.h"
49 #include "AliMUONTriggerCrateStore.h"
50 #include "AliMUONLocalTriggerBoard.h"
51 #include "AliMUONTriggerCircuit.h"
52 #include "AliMUONDigitMaker.h"
55 #include "AliMpVSegmentation.h"
56 #include "AliMpSegmentation.h"
59 #include "AliTracker.h"
62 #include <Riostream.h>
63 #include <TClonesArray.h>
69 ClassImp(AliMUONTrackHitPattern) // Class implementation in ROOT context
73 //______________________________________________________________________________
74 AliMUONTrackHitPattern::AliMUONTrackHitPattern(AliMUONRecData *data)
77 fTransformer(new AliMUONGeometryTransformer(kTRUE)),
78 fCrateManager(new AliMUONTriggerCrateStore()),
79 fDigitMaker(new AliMUONDigitMaker())
81 /// Default constructor
84 const AliMagF* kField = AliTracker::GetFieldMap();
85 if (!kField) AliFatal("No field available");
86 AliMUONTrackExtrap::SetField(kField);
88 // Geometry transformer
89 fTransformer->ReadGeometryData("volpath.dat", "geometry.root");
91 // Crate manager to retrieve local boards
92 fCrateManager->ReadFromFile();
95 fDigitMaker->SetCrateManager(fCrateManager);
97 for(Int_t ch=0; ch<4; ch++){
98 fTriggerDigitsList[ch].Clear();
103 //______________________________________________________________________________
104 AliMUONTrackHitPattern::~AliMUONTrackHitPattern(void)
107 for(Int_t ch=0; ch<4; ch++){
108 fTriggerDigitsList[ch].Delete();
110 delete fCrateManager;
114 //______________________________________________________________________________
115 void AliMUONTrackHitPattern::GetHitPattern(TClonesArray *recTracksPtr)
119 /// Loops on reco tracks, extrapolates them to trigger chambers
120 /// and searches for matching digits
123 const Int_t kMask[2][4]={{0x80, 0x40, 0x20, 0x10},
124 {0x08, 0x04, 0x02, 0x01}};
128 Int_t nRecTracks = (Int_t)recTracksPtr->GetEntriesFast();
129 for(Int_t iTrack=0; iTrack<nRecTracks; iTrack++){
131 AliMUONTrack *muonTrack = (AliMUONTrack*) recTracksPtr->At(iTrack);
132 AliMUONTrackParam *trackParam = (AliMUONTrackParam*) ((muonTrack->GetTrackParamAtHit())->Last());
133 for(Int_t ch=0; ch<4; ch++){
134 AliMUONTrackExtrap::ExtrapToZCov(trackParam, AliMUONConstants::DefaultChamberZ(10+ch));
135 FindPadMatchingTrack(trackParam, isMatch, ch);
136 for(Int_t cath=0; cath<2; cath++){
137 if(isMatch[cath]) pattern |= kMask[cath][ch];
140 muonTrack->SetHitsPatternInTrigCh(pattern);
146 //______________________________________________________________________________
147 void AliMUONTrackHitPattern::FindPadMatchingTrack(AliMUONTrackParam *trackParam,
148 Bool_t isMatch[2], Int_t iChamber)
151 /// Given track position, searches for matching digits.
154 Float_t minMatchDist[2];
156 for(Int_t cath=0; cath<2; cath++){
157 isMatch[cath]=kFALSE;
158 minMatchDist[cath]=9999.;
161 Int_t ndigits = (Int_t)fTriggerDigitsList[iChamber].GetEntries();
162 AliMUONDigit * mDigit = 0x0;
163 for(Int_t idigit=0; idigit<ndigits; idigit++) { // digit loop
164 mDigit = (AliMUONDigit*)fTriggerDigitsList[iChamber].At(idigit);
165 Int_t currDetElemId = mDigit->DetElemId();
167 Int_t cathode = mDigit->Cathode();
168 Int_t ix = mDigit->PadX();
169 Int_t iy = mDigit->PadY();
170 Float_t xpad, ypad, zpad;
171 const AliMpVSegmentation* seg = AliMpSegmentation::Instance()
172 ->GetMpSegmentation(currDetElemId,AliMp::GetCathodType(cathode));
174 AliMpPad pad = seg->PadByIndices(AliMpIntPair(ix,iy),kTRUE);
175 Float_t xlocal1 = pad.Position().X();
176 Float_t ylocal1 = pad.Position().Y();
177 Float_t dpx = pad.Dimensions().X();
178 Float_t dpy = pad.Dimensions().Y();
179 fTransformer->Local2Global(currDetElemId, xlocal1, ylocal1, 0, xpad, ypad, zpad);
180 Float_t matchDist = MinDistanceFromPad(xpad, ypad, zpad, dpx, dpy, trackParam);
181 if(matchDist>minMatchDist[cathode])continue;
182 isMatch[cathode] = kTRUE;
183 minMatchDist[cathode] = matchDist;
190 //______________________________________________________________________________
191 Float_t AliMUONTrackHitPattern::MinDistanceFromPad(Float_t xPad, Float_t yPad, Float_t zPad,
192 Float_t dpx, Float_t dpy, AliMUONTrackParam *trackParam)
195 /// Decides if the digit belongs to the track.
197 Float_t xTrackAtPad = trackParam->GetNonBendingCoor();
198 Float_t yTrackAtPad = trackParam->GetBendingCoor();
200 Float_t sigmaX, sigmaY, sigmaMS;
201 GetPosUncertainty(trackParam, zPad, sigmaX, sigmaY, sigmaMS);
203 Float_t maxDistX = 3.*(sigmaX + sigmaMS); // in cm
204 Float_t maxDistY = 3.*(sigmaY + sigmaMS); // in cm
206 Float_t deltaX = TMath::Abs(xPad-xTrackAtPad)-dpx;
207 Float_t deltaY = TMath::Abs(yPad-yTrackAtPad)-dpy;
209 Float_t matchDist = 99999.;
210 if(deltaX<=maxDistX && deltaY<=maxDistY) matchDist = TMath::Max(deltaX, deltaY);
215 //______________________________________________________________________________
216 void AliMUONTrackHitPattern::GetPosUncertainty(AliMUONTrackParam *trackParam, Float_t zChamber,
217 Float_t &sigmaX, Float_t &sigmaY, Float_t &sigmaMS)
220 /// Returns uncertainties on extrapolated position.
221 /// Takes into account Branson plane corrections in the iron wall.
224 const Float_t kAlpha = 0.1123; // GeV/c
226 // Find a better way to get such parameters ???
227 const Float_t kZFilterIn = 1471.; // From STRUCT/SHILConst2.h
228 const Float_t kZFilterOut = kZFilterIn + 120.; // From STRUCT/SHILConst2.h
230 const Float_t kZBranson = - (kZFilterIn + (kZFilterOut - kZFilterIn)*2./3. ); // - sign because distance are positive
231 Float_t zDistFromWall = TMath::Abs(zChamber - kZBranson);
232 Float_t zDistFromLastTrackCh = TMath::Abs(zChamber - AliMUONConstants::DefaultChamberZ(9));
234 TMatrixD *covParam = trackParam->GetCovariances();
236 sigmaX = (*covParam)(0,0);
237 sigmaY = (*covParam)(2,2);
239 // If covariance matrix is not extrapolated, use "reasonable" errors
240 // (To be removed as soon as covariance matrix is correctly propagated).
241 if (sigmaX==0.)sigmaX = 0.003 * zDistFromLastTrackCh;
242 if (sigmaY==0.)sigmaY = 0.004 * zDistFromLastTrackCh;
244 Float_t p = trackParam->P();
245 Float_t thetaMS = kAlpha/p;
246 sigmaMS = zDistFromWall * TMath::Tan(thetaMS);
252 //____________________________________________________________________
253 Bool_t AliMUONTrackHitPattern::TriggerDigits()
256 /// make (S)Digit for trigger
259 for(Int_t ch=0; ch<4; ++ch) fTriggerDigitsList[ch].Clear();
267 AliMUONLocalTrigger *locTrg = 0x0;
269 fMUONData->SetTreeAddress("RC,TC");
270 fMUONData->GetTrigger();
272 TClonesArray *localTrigger = fMUONData->LocalTrigger();
273 Int_t nLocTrig = (Int_t) localTrigger->GetEntriesFast();
275 for(Int_t iLoc=0; iLoc<nLocTrig; iLoc++) {
276 locTrg = (AliMUONLocalTrigger*)localTrigger->UncheckedAt(iLoc);
278 TArrayS xyPattern[2];
282 xyPattern[0].AddAt(locTrg->GetX1Pattern(),0);
283 xyPattern[0].AddAt(locTrg->GetX2Pattern(),1);
284 xyPattern[0].AddAt(locTrg->GetX3Pattern(),2);
285 xyPattern[0].AddAt(locTrg->GetX4Pattern(),3);
287 xyPattern[1].AddAt(locTrg->GetY1Pattern(),0);
288 xyPattern[1].AddAt(locTrg->GetY2Pattern(),1);
289 xyPattern[1].AddAt(locTrg->GetY3Pattern(),2);
290 xyPattern[1].AddAt(locTrg->GetY4Pattern(),3);
292 for(Int_t cath=0; cath<2; cath++){
293 for(Int_t ch=0; ch<4; ch++){
294 if(xyPattern[cath][ch]==0) continue;
298 nBoard = locTrg->LoCircuit();
299 fDigitMaker->TriggerDigits(nBoard, xyPattern, digitList);
302 } // loop on localTriggers
304 for (Int_t iEntry = 0; iEntry < digitList.GetEntries(); ++iEntry) {
305 AliMUONDigit* digit = (AliMUONDigit*)digitList.At(iEntry);
306 Int_t detElemId = digit->DetElemId();
307 Int_t iChamber = detElemId/100 - 11; //FIXEME should be given by mapping
308 fTriggerDigitsList[iChamber].Add(digit);