* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.3 2000/07/03 11:54:57 morsch
-AliMUONSegmentation and AliMUONHitMap have been replaced by AliSegmentation and AliHitMap in STEER
-The methods GetPadIxy and GetPadXxy of AliMUONSegmentation have changed name to GetPadI and GetPadC.
+/* $Id$ */
-Revision 1.2 2000/06/15 07:58:48 morsch
-Code from MUON-dev joined
-
-Revision 1.1.2.2 2000/06/12 07:58:06 morsch
-include TMath.h
-
-Revision 1.1.2.1 2000/06/09 22:01:09 morsch
-Code from AliMUONHitMap.h
-Most coding rule violations corrected.
-
-*/
+// -----------------------
+// Class AliMUONHitMapA1
+// -----------------------
+// Implements AliHitMap as a 1-dim array
+// ...
#include "AliMUONHitMapA1.h"
-#include "AliSegmentation.h"
-#include "AliMUONResponse.h"
#include "AliMUONDigit.h"
+#include "AliLog.h"
+
#include <TObjArray.h>
#include <TMath.h>
ClassImp(AliMUONHitMapA1)
-
-AliMUONHitMapA1::AliMUONHitMapA1(AliSegmentation *seg, TObjArray *dig)
+AliMUONHitMapA1::AliMUONHitMapA1()
+ : AliHitMap(),
+ fNpx(0),
+ fNpy(0),
+ fDigits(0),
+ fMaxIndex(0),
+ fHitMap(0)
{
-// Constructor
- fSegmentation = seg;
- fNpx = fSegmentation->Npx();
- fNpy = fSegmentation->Npy();
- fMaxIndex=2*(fNpx+1)*2*(fNpy+1)+2*fNpy;
-
- fHitMap = new Int_t[fMaxIndex];
- fDigits = dig;
- fNdigits = fDigits->GetEntriesFast();
- Clear();
+/// Default constructor
}
-AliMUONHitMapA1::AliMUONHitMapA1(const AliMUONHitMapA1 & hitMap)
+//________________________________________________________________________________
+AliMUONHitMapA1::AliMUONHitMapA1(Int_t npx, Int_t npy, TObjArray* dig)
+ : AliHitMap(),
+ fNpx(npx),
+ fNpy(npy),
+ fDigits(dig),
+ fMaxIndex(2*(fNpx+1)*2*(fNpy+1)+2*fNpy),
+ fHitMap(new Int_t[fMaxIndex])
{
-// Dummy copy constructor
- ;
+/// Standard constructor
+ Clear();
}
-
-
+//_________________________________
AliMUONHitMapA1::~AliMUONHitMapA1()
{
-// Destructor
-// if (fDigits) delete fDigits;
+/// Destructor
if (fHitMap) delete[] fHitMap;
}
-
+//______________________________________
void AliMUONHitMapA1::Clear(const char *)
{
-// Clear hitmap
+/// Clear hitmap
memset(fHitMap,0,sizeof(int)*fMaxIndex);
}
+//___________________________________________________
+Bool_t AliMUONHitMapA1::ValidateHit(Int_t ix, Int_t iy)
+{
+/// Check if pad coordinates are within boundaries
-Int_t AliMUONHitMapA1::CheckedIndex(Int_t ix, Int_t iy)
+// printf("\n Validate %d %d %d %d", ix, iy, fNpx, fNpy);
+
+ return (TMath::Abs(ix) <= fNpx && TMath::Abs(iy) <= fNpy);
+}
+//_________________________________________________________
+Int_t AliMUONHitMapA1::CheckedIndex(Int_t ix, Int_t iy) const
{
-// Return checked indices ix, iy
+/// Return checked indices ix, iy
Int_t index=2*fNpy*(ix+fNpx)+(iy+fNpy);
- if (index > fMaxIndex) {
- printf("\n \n \n Try to read/write outside array !!!! \n \n %d %d %d %d %d %d",ix,iy, fMaxIndex, index, fNpx, fNpy);
+ if ( index < 0 || index >= fMaxIndex ) {
+ AliWarning(Form("index outside array ix %d iy %d MaxIndex %d index %d Npx %d Npy %d",
+ ix, iy, fMaxIndex, index, fNpx, fNpy));
return fMaxIndex-1;
} else {
return index;
}
}
-
-
+//_____________________________
void AliMUONHitMapA1::FillHits()
{
-// Fill hits from digits list
+/// Fill hits from digits list
Int_t ndigits = fDigits->GetEntriesFast();
//printf("\n Filling hits into HitMap\n");
//printf("FindRawClusters -- ndigits %d \n",ndigits);
if (!ndigits) return;
AliMUONDigit *dig;
- for (Int_t ndig=0; ndig<fNdigits; ndig++) {
+ for (Int_t ndig=0; ndig<ndigits; ndig++) {
dig = (AliMUONDigit*)fDigits->UncheckedAt(ndig);
- SetHit(dig->fPadX,dig->fPadY,ndig);
+ SetHit(dig->PadX(),dig->PadY(),ndig);
}
}
-
-
+//___________________________________________________________
void AliMUONHitMapA1::SetHit(Int_t ix, Int_t iy, Int_t idigit)
{
-// Assign digit to hit cell ix,iy
+/// Assign digit to hit cell ix,iy
+
// fHitMap[kMaxNpady*(ix+fNpx)+(iy+fNpy)]=idigit+1;
fHitMap[CheckedIndex(ix, iy)]=idigit+1;
}
-
+//_______________________________________________
void AliMUONHitMapA1::DeleteHit(Int_t ix, Int_t iy)
{
-// Delete hit at cell ix,iy
+/// Delete hit at cell ix,iy
+
// fHitMap[kMaxNpady*(ix+fNpx)+(iy+fNpy)]=0;
fHitMap[CheckedIndex(ix, iy)]=0;
}
-
+//_____________________________________________
void AliMUONHitMapA1::FlagHit(Int_t ix, Int_t iy)
{
-// Flag hit as used
+/// Flag hit as used
fHitMap[CheckedIndex(ix, iy)]=
-TMath::Abs(fHitMap[CheckedIndex(ix, iy)]);
}
-
-Int_t AliMUONHitMapA1::GetHitIndex(Int_t ix, Int_t iy)
+//________________________________________________________
+Int_t AliMUONHitMapA1::GetHitIndex(Int_t ix, Int_t iy) const
{
-// Get absolute value of contents of hit cell ix,iy
+/// Get absolute value of contents of hit cell ix,iy
return TMath::Abs(fHitMap[CheckedIndex(ix, iy)])-1;
}
-
-TObject* AliMUONHitMapA1::GetHit(Int_t ix, Int_t iy)
+//_______________________________________________________
+TObject* AliMUONHitMapA1::GetHit(Int_t ix, Int_t iy) const
{
- // Get pointer to object at hit cell ix, iy
- // Force crash if index does not exist ! (Manu)
+/// Get pointer to object at hit cell ix, iy
+/// Force crash if index does not exist ! (Manu)
Int_t index=GetHitIndex(ix,iy);
return (index <0) ? 0 : fDigits->UncheckedAt(GetHitIndex(ix,iy));
}
-
+//_________________________________________________
FlagType AliMUONHitMapA1::TestHit(Int_t ix, Int_t iy)
{
-// Check if hit cell is empty, used or unused
-//
- Int_t inf=fHitMap[CheckedIndex(ix, iy)];
+/// Check if hit cell is empty, used or unused
+
+ Int_t index = CheckedIndex(ix, iy);
+ if (index<0 || index >= fMaxIndex) return kEmpty;
+
+ Int_t inf=fHitMap[index];
if (inf < 0) {
return kUsed;
} else if (inf == 0) {
return kUnused;
}
}
-
-AliMUONHitMapA1 & AliMUONHitMapA1::operator = (const AliMUONHitMapA1 & rhs)
-{
-// Dummy assignment operator
- return *this;
-}
-
-