* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.2 2000/06/15 07:58:48 morsch
-Code from MUON-dev joined
+/* $Id$ */
-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.
-
-*/
+#include <TObjArray.h>
+#include <TMath.h>
#include "AliMUONHitMapA1.h"
-#include "AliSegmentation.h"
-#include "AliMUONResponse.h"
+#include "AliMUONGeometrySegmentation.h"
#include "AliMUONDigit.h"
-
-#include <TObjArray.h>
-#include <TMath.h>
+#include "AliLog.h"
ClassImp(AliMUONHitMapA1)
+AliMUONHitMapA1::AliMUONHitMapA1()
+ : AliHitMap()
+{
+ // Default constructor
+ fNpx = 0;
+ fNpy = 0;
+ fMaxIndex = 0;
+
+ fHitMap = 0;
+ fDigits = 0;
+}
-AliMUONHitMapA1::AliMUONHitMapA1(AliSegmentation *seg, TObjArray *dig)
+//________________________________________________________________________________
+AliMUONHitMapA1::AliMUONHitMapA1(Int_t idDE, AliMUONGeometrySegmentation* seg, TObjArray* dig)
+ : AliHitMap()
{
-// Constructor
- fSegmentation = seg;
- fNpx = fSegmentation->Npx();
- fNpy = fSegmentation->Npy();
+// Constructor with new segmentation
+ fNpx = seg->Npx(idDE)+1;
+ fNpy = seg->Npy(idDE)+1;
fMaxIndex=2*(fNpx+1)*2*(fNpy+1)+2*fNpy;
fHitMap = new Int_t[fMaxIndex];
fDigits = dig;
- fNdigits = fDigits->GetEntriesFast();
Clear();
}
-
+//______________________________________________________________
AliMUONHitMapA1::AliMUONHitMapA1(const AliMUONHitMapA1 & hitMap)
+ : AliHitMap(hitMap)
{
-// Dummy copy constructor
- ;
-}
+// Protected copy constructor
-
+ AliFatal("Not implemented.");
+}
+//_________________________________
AliMUONHitMapA1::~AliMUONHitMapA1()
{
// Destructor
-// if (fDigits) delete fDigits;
if (fHitMap) delete[] fHitMap;
}
-
-void AliMUONHitMapA1::Clear()
+//______________________________________
+void AliMUONHitMapA1::Clear(const char *)
{
// Clear hitmap
memset(fHitMap,0,sizeof(int)*fMaxIndex);
}
-
-Int_t AliMUONHitMapA1::CheckedIndex(Int_t ix, Int_t iy)
+//___________________________________________________
+Bool_t AliMUONHitMapA1::ValidateHit(Int_t ix, Int_t iy)
+{
+ //
+ // Check if pad coordinates are within boundaries
+ //
+// 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
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 >= 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);
return fMaxIndex-1;
} else {
return index;
}
}
-
-
+//_____________________________
void AliMUONHitMapA1::FillHits()
{
// Fill hits from digits list
//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
// 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
// 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
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
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)
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)];
+ 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;
+// Protected assignement operator
+
+ if (this == &rhs) return *this;
+
+ AliFatal( "Not implemented.");
+
+ return *this;
}
+
+
+