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 //-----------------------------------------------------------------------------
19 /// \class AliMUONVDigit
21 /// This is the base class of a MUON digit that most client code should deal with.
22 /// There should be no reason to have to use a concrete class in most cases.
24 /// All digits have basic features, like :
26 /// - a way to identify it : detection element, electronics card and
27 /// channel, cathode. Note that some static methods exists to compact
28 /// those 4 informations into a single 4 bytes integer (stored in the
29 /// fUniqueID data member present in all TObjects).
33 /// - a set of boolean methods to indicate whether the digit has been calibrated, etc...
35 /// In addition, if HasMCInformation is true, the digit store also the list
36 /// of MC tracks that contributed to its charge
38 /// Also, if HasGeometryInformation is true, the digit knows the position and
39 /// the (half) dimensions (in cm) of the pad it corresponds to.
41 /// \author Laurent Aphecetche, Subatech
42 //-----------------------------------------------------------------------------
44 #include "AliMUONVDigit.h"
46 #include <Riostream.h>
50 ClassImp(AliMUONVDigit)
53 //_____________________________________________________________________________
54 AliMUONVDigit::AliMUONVDigit(Int_t detElemId, Int_t eCardId,
55 Int_t eCardChannel, Int_t cathode)
58 /// Normal constructor for trigger digits
59 SetUniqueID(BuildUniqueID(detElemId,eCardId,eCardChannel,cathode));
62 //_____________________________________________________________________________
63 AliMUONVDigit::AliMUONVDigit()
69 //_____________________________________________________________________________
70 AliMUONVDigit::~AliMUONVDigit()
75 //_____________________________________________________________________________
77 AliMUONVDigit::IsEqual(const TObject* object) const
79 /// Whether we're equal to object.
80 /// WARNING : only based on our identifiers (de,manu,channel,cathode), not our
81 /// content (i.e. charge, status...)
83 const AliMUONVDigit* d = static_cast<const AliMUONVDigit*>(object);
85 return ( DetElemId() == d->DetElemId() &&
86 Cathode() == d->Cathode() &&
87 ManuId() == d->ManuId() &&
88 ManuChannel() == d->ManuChannel() );
91 //_____________________________________________________________________________
93 AliMUONVDigit::Compare(const TObject* object) const
95 /// Compare two digits, trying to get as complete an order as possible.
96 /// We sort by DE, then by charge, then by manu, etc...
98 const AliMUONVDigit* d = static_cast<const AliMUONVDigit*>(object);
100 if ( DetElemId() > d->DetElemId() )
104 else if ( DetElemId() < d->DetElemId() )
110 if ( Charge() > d->Charge() )
114 else if ( Charge() < d->Charge() )
120 if ( ManuId() > d->ManuId() )
124 else if ( ManuId() < d->ManuId() )
130 if ( ManuChannel() > d->ManuChannel() )
134 else if ( ManuChannel() < d->ManuChannel() )
144 //_____________________________________________________________________________
146 AliMUONVDigit::BuildUniqueID(Int_t detElemId, Int_t manuId,
147 Int_t manuChannel, Int_t cathode)
149 /// Build a single integer with id information
150 return ( ( detElemId ) | ( manuId << 12 ) | ( manuChannel << 24 )
151 | ( cathode << 30 ) );
154 //_____________________________________________________________________________
156 AliMUONVDigit::DetElemId(UInt_t uniqueID)
158 /// Return detection element id part of the uniqueID
159 return uniqueID & 0xFFF;
162 //_____________________________________________________________________________
164 AliMUONVDigit::ManuChannel(UInt_t uniqueID)
166 /// Return manuChannel part of the uniqueID
167 return ( uniqueID & 0x3F000000 ) >> 24;
170 //_____________________________________________________________________________
172 AliMUONVDigit::ManuId(UInt_t uniqueID)
174 /// Return manuId part of the uniqueID
175 return ( uniqueID & 0xFFF000 ) >> 12;
178 //_____________________________________________________________________________
180 AliMUONVDigit::Cathode(UInt_t uniqueID)
182 /// Return the cathode part of the uniqueID
183 return ( uniqueID & 0x40000000 ) >> 30;
186 //_____________________________________________________________________________
188 AliMUONVDigit::DecodeUniqueID(UInt_t uniqueID,
189 Int_t& detElemId, Int_t& manuId,
190 Int_t& manuChannel, Int_t& cathode)
192 /// Unpack uniqueID into 4 elements
193 detElemId = DetElemId(uniqueID);
194 manuId = ManuId(uniqueID);
195 manuChannel = ManuChannel(uniqueID);
196 cathode = Cathode(uniqueID);
199 //_____________________________________________________________________________
201 AliMUONVDigit::GetName() const
203 /// Return the name of this digit, composed of its id parts.
204 return Form("DE%04d-%04d-%02d-%d",
205 DetElemId(),ManuId(),ManuChannel(),Cathode());
208 //_____________________________________________________________________________
210 AliMUONVDigit::Print(Option_t* opt) const
213 /// If opt=="tracks", info on tracks are printed too.
215 cout << Form("<%s>: ID %12u DE %4d Cath %d (Ix,Iy)=(%3d,%3d) (Manu,Channel)=(%4d,%2d)"
217 ClassName(),GetUniqueID(),
218 DetElemId(),Cathode(),PadX(),PadY(),ManuId(),ManuChannel(),Charge());
228 if ( IsCalibrated() )
246 cout << Form(" ADC=%4d StatusMap=%04x",ADC(),StatusMap());
248 TString options(opt);
250 if ( options.Contains("tracks") && HasMCInformation() )
252 cout << " Hit " << setw(3) << Hit();
253 Int_t ntracks = Ntracks();
256 cout << " Tracks : " << setw(2) << ntracks;
257 for ( Int_t i = 0; i < ntracks; ++i )
259 cout << " Track(" << i << ")=" << setw(3) << Track(i)
260 << " Charge(" << i << ")=" << setw(5) << TrackCharge(i);
265 cout << " no track info.";