+++ /dev/null
-#ifndef ALIHLTMUONBASICTYPES_H
-#define ALIHLTMUONBASICTYPES_H
-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-////////////////////////////////////////////////////////////////////////////////
-
-
-#ifndef NULL
-# define NULL 0x0
-#endif
-
-
-typedef char Char;
-typedef short Short;
-typedef int Int;
-typedef long Long;
-typedef unsigned char UChar;
-typedef unsigned short UShort;
-typedef unsigned int UInt;
-typedef unsigned long ULong;
-typedef float Float;
-typedef double Double;
-
-typedef char* NullString;
-
-
-#endif // ALIHLTMUONBASICTYPES_H
+++ /dev/null
-#ifndef ALIHLTMUONCOREADCSTREAM_H
-#define ALIHLTMUONCOREADCSTREAM_H
-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-////////////////////////////////////////////////////////////////////////////////
-
-#include "AliHLTMUONBasicTypes.h"
-
-
-class AliHLTMUONCoreADCStream
-{
-};
-
-
-#endif // ALIHLTMUONCOREADCSTREAM_H
+++ /dev/null
-#ifndef ALIHLTMUONCORECLUSTER_H
-#define ALIHLTMUONCORECLUSTER_H
-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-////////////////////////////////////////////////////////////////////////////////
-
-#include "AliHLTMUONCorePoint.h"
-
-
-typedef AliHLTMUONCorePoint AliHLTMUONCoreClusterPoint;
-
-
-#endif // ALIHLTMUONCORECLUSTER_H
+++ /dev/null
-#ifndef ALIHLTMUONCOREEVENTID_H
-#define ALIHLTMUONCOREEVENTID_H
-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-////////////////////////////////////////////////////////////////////////////////
-
-#include "AliHLTMUONBasicTypes.h"
-#include "AliHLTMUONUtils.h"
-
-
-// Must correct this definition!!!!
-struct AliHLTMUONCoreEventID
-{
- UInt fBunch;
- UInt fTimeStamp;
-
-
- AliHLTMUONCoreEventID(UInt bunch = 0, UInt timestamp = 0)
- : fBunch(bunch), fTimeStamp(timestamp)
- {};
-
-
- bool operator == (const AliHLTMUONCoreEventID& rhs)
- {
- return fBunch == rhs.fBunch and fTimeStamp == rhs.fTimeStamp;
- };
-
- bool operator != (const AliHLTMUONCoreEventID& rhs)
- {
- return not (*this == rhs);
- };
-};
-
-
-#endif // ALIHLTMUONCOREEVENTID_H
+++ /dev/null
-#ifndef ALIHLTMUONCOREPOINT_H
-#define ALIHLTMUONCOREPOINT_H
-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-// A 2D point structure using floats.
-// These are used to store impact points on the trigger chambers and
-// cluster centroids.
-//
-////////////////////////////////////////////////////////////////////////////////
-
-#include "AliHLTMUONBasicTypes.h"
-
-class AliHLTMUONCorePoint
-{
-public:
-
- AliHLTMUONCorePoint() : fX(0.0), fY(0.0) {}
-
- AliHLTMUONCorePoint(Float x, Float y) : fX(x), fY(y) {}
-
- Float X() const { return fX; }
- Float& X() { return fX; }
- void X(Float x) { fX = x; }
- Float Y() const { return fY; }
- Float& Y() { return fY; }
- void Y(Float y) { fY = y; }
-
-private:
-
- Float fX, fY; // X and Y coordinate.
-};
-
-#endif // ALIHLTMUONCOREPOINT_H
+++ /dev/null
-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-// The region of interest object is used to encode/decode and work with boundary
-// box type regions of interest. The 32 bit ROI codes are used to communicate
-// regions of interest between different parts of the dHLT system. This is more
-// efficient than sending 20 byte long region of interest objects.
-//
-////////////////////////////////////////////////////////////////////////////////
-
-#include "AliHLTMUONCoreRegionOfInterest.h"
-#include <math.h>
-
-
-Float AliHLTMUONCoreRegionOfInterest::fgPlaneScale[AliHLTMUONCoreRegionOfInterest::kNumberOfTrackingChambers]
- = { 102.0f,
- 104.0f,
- 130.0f,
- 132.0f,
- 184.0f,
- 188.0f,
- 238.0f,
- 244.0f,
- 270.0f,
- 275.0f
- };
-
-
-void AliHLTMUONCoreRegionOfInterest::CreateToContain(
- const AliHLTMUONCoreClusterPoint& point, AliHLTMUONCoreChamberID chamber
- )
-{
-// Creates a region of interest around the given point for the
-// specified chamber.
-
- Assert( 0 <= chamber && chamber < (AliHLTMUONCoreChamberID)kNumberOfTrackingChambers );
- fChamber = chamber;
- fLeft = fRight = point.X();
- fBottom = fTop = point.Y();
-}
-
-
-void AliHLTMUONCoreRegionOfInterest::ExpandToContain(const AliHLTMUONCoreClusterPoint& point)
-{
-// Extends the region of interest to contain the specified point.
-
- if (point.X() < fLeft)
- fLeft = point.X();
- else
- if (point.X() > fRight) fRight = point.X();
- if (point.Y() < fBottom)
- fBottom = point.Y();
- else
- if (point.Y() > fTop) fTop = point.Y();
-}
-
-
-void AliHLTMUONCoreRegionOfInterest::CreateToContain(
- const AliHLTMUONCoreClusterPoint* points, UInt count,
- AliHLTMUONCoreChamberID chamber
- )
-{
-// Creates a region of interest around all the given points and for the
-// specified chamber.
-
- Assert( 0 <= chamber && chamber < (AliHLTMUONCoreChamberID)kNumberOfTrackingChambers );
- Assert( count > 0 );
-
- CreateToContain(points[0], chamber);
- for (UInt i = 1; i < count; i++)
- ExpandToContain(points[i]);
-}
-
-
-bool AliHLTMUONCoreRegionOfInterest::InBounds()
-{
-// Checks if the region of interest is within the boundaries imposed on
-// the specific chamber plane. This boundary is aproximately the square
-// box around the chamber's detection region.
-
- Assert( 0 <= fChamber && fChamber < (AliHLTMUONCoreChamberID)kNumberOfTrackingChambers );
- register Float bound = fgPlaneScale[fChamber];
- return -bound <= fLeft
- && fRight <= bound
- && -bound <= fBottom
- && fTop <= bound;
-}
-
-
-#define MAX_LEVELS 13
-#define GRID_SIZE (2 << MAX_LEVELS)
-
-// The following numbers specifiy the total indices for lower levels.
-// Computed with: x_level = sum from 0 to level ( ( 2^level - 1 )^2 )
-#define MAX_INDICES 357848422
-#define INDICES_TO_LEVEL_13 89445733
-#define INDICES_TO_LEVEL_12 22353252
-#define INDICES_TO_LEVEL_11 5584227
-#define INDICES_TO_LEVEL_10 1394018
-#define INDICES_TO_LEVEL_9 347489
-#define INDICES_TO_LEVEL_8 86368
-#define INDICES_TO_LEVEL_7 21343
-#define INDICES_TO_LEVEL_6 5214
-#define INDICES_TO_LEVEL_5 1245
-#define INDICES_TO_LEVEL_4 284
-#define INDICES_TO_LEVEL_3 59
-#define INDICES_TO_LEVEL_2 10
-#define INDICES_TO_LEVEL_1 1
-
-
-UInt AliHLTMUONCoreRegionOfInterest::fgIndexOffsets[13]
- = { INDICES_TO_LEVEL_1,
- INDICES_TO_LEVEL_2,
- INDICES_TO_LEVEL_3,
- INDICES_TO_LEVEL_4,
- INDICES_TO_LEVEL_5,
- INDICES_TO_LEVEL_6,
- INDICES_TO_LEVEL_7,
- INDICES_TO_LEVEL_8,
- INDICES_TO_LEVEL_9,
- INDICES_TO_LEVEL_10,
- INDICES_TO_LEVEL_11,
- INDICES_TO_LEVEL_12,
- INDICES_TO_LEVEL_13
- };
-
-
-inline void AliHLTMUONCoreRegionOfInterest::ConvertToGrid(
- register UInt& l, register UInt& r, register UInt& b, register UInt& t
- ) const
-{
-// Converts the internal region of interest boundary box, which is
-// specified in as floats, into a regular integer grid.
-// l = left boundary, r = right boundary, b = bottom boundary,
-// t = top boundary.
-
- Assert( 0 <= fChamber && fChamber < (AliHLTMUONCoreChamberID)kNumberOfTrackingChambers );
-
- register Float scale = fgPlaneScale[fChamber];
- l = (UInt) floor( (fLeft / scale + 1.0f) * 0.5f * GRID_SIZE );
- r = (UInt) ceil( (fRight / scale + 1.0f) * 0.5f * GRID_SIZE );
- b = (UInt) floor( (fBottom / scale + 1.0f) * 0.5f * GRID_SIZE );
- t = (UInt) ceil( (fTop / scale + 1.0f) * 0.5f * GRID_SIZE );
-
- /*
- // Not required with proper coding. This case is handled properly by
- // encoding as a global boundary box.
- Assert( l <= GRID_SIZE );
- Assert( r <= GRID_SIZE );
- Assert( b <= GRID_SIZE );
- Assert( t <= GRID_SIZE );
- */
-}
-
-
-inline void AliHLTMUONCoreRegionOfInterest::ConvertBackFromGrid(
- register UInt l, register UInt r, register UInt b, register UInt t
- )
-{
-// Performs the inverse conversion of the method ConvertToGrid.
-// That is converts from a regular integer grid back to the internal
-// floating point boundary box specification.
-
- Assert( l <= GRID_SIZE );
- Assert( r <= GRID_SIZE );
- Assert( b <= GRID_SIZE );
- Assert( t <= GRID_SIZE );
-
- Assert( 0 <= fChamber && fChamber < (AliHLTMUONCoreChamberID)kNumberOfTrackingChambers );
- register Float scale = fgPlaneScale[fChamber];
- fLeft = ((Float)l / (Float)GRID_SIZE - 0.5f) * 2.0f * scale;
- fRight = ((Float)r / (Float)GRID_SIZE - 0.5f) * 2.0f * scale;
- fBottom = ((Float)b / (Float)GRID_SIZE - 0.5f) * 2.0f * scale;
- fTop = ((Float)t / (Float)GRID_SIZE - 0.5f) * 2.0f * scale;
-}
-
-
-AliHLTMUONCoreROI AliHLTMUONCoreRegionOfInterest::Encode() const
-{
-// Encodes the region of interest into a 32 bit code.
-
- UInt l, r, b, t, maxwidth, index;
-
- ConvertToGrid(l, r, b, t);
-
- // Work out which resolution level the location of the ROI
- // boundary box needs to be coded at. The higher the level
- // the smaller the box size: i.e. box_size = scale / 2^(level)
- // We use a binary search type method for finding the level.
- // More specificaly we search for which n,
- // (l * 2^(-n) + 2) * 2^n >= r and (b * 2^(-n) + 2) * 2^n >= t
-
- if ( ((l >> 6) + 2) << 6 >= r && ((b >> 6) + 2) << 6 >= t )
- {
- if ( ((l >> 2) + 2) << 2 >= r && ((b >> 2) + 2) << 2 >= t )
- {
- if ( l + 2 >= r && b + 2 >= t )
- {
- index = INDICES_TO_LEVEL_13;
- maxwidth = (2 << 13) - 1;
- }
- else
- {
- if ( ((l >> 1) + 2) << 1 >= r && ((b >> 1) + 2) << 1 >= t )
- {
- index = INDICES_TO_LEVEL_12;
- maxwidth = (2 << 12) - 1;
- l >>= 1; b >>= 1;
- }
- else
- {
- index = INDICES_TO_LEVEL_11;
- maxwidth = (2 << 11) - 1;
- l >>= 2; b >>= 2;
- }
- }
- }
- else
- {
- if ( ((l >> 4) + 2) << 4 >= r && ((b >> 4) + 2) << 4 >= t )
- {
- if ( ((l >> 3) + 2) << 3 >= r && ((b >> 3) + 2) << 3 >= t )
- {
- index = INDICES_TO_LEVEL_10;
- maxwidth = (2 << 10) - 1;
- l >>= 3; b >>= 3;
- }
- else
- {
- index = INDICES_TO_LEVEL_9;
- maxwidth = (2 << 9) - 1;
- l >>= 4; b >>= 4;
- }
- }
- else
- {
- if ( ((l >> 5) + 2) << 5 >= r && ((b >> 5) + 2) << 5 >= t )
- {
- index = INDICES_TO_LEVEL_8;
- maxwidth = (2 << 8) - 1;
- l >>= 5; b >>= 5;
- }
- else
- {
- index = INDICES_TO_LEVEL_7;
- maxwidth = (2 << 7) - 1;
- l >>= 6; b >>= 6;
- }
- }
- }
- }
- else
- {
- if ( ((l >> 10) + 2) << 10 >= r && ((b >> 10) + 2) << 10 >= t )
- {
- if ( ((l >> 8) + 2) << 8 >= r && ((b >> 8) + 2) << 8 >= t )
- {
- if ( ((l >> 7) + 2) << 7 >= r && ((b >> 7) + 2) << 7 >= t )
- {
- index = INDICES_TO_LEVEL_6;
- maxwidth = (2 << 6) - 1;
- l >>= 7; b >>= 7;
- }
- else
- {
- index = INDICES_TO_LEVEL_5;
- maxwidth = (2 << 5) - 1;
- l >>= 8; b >>= 8;
- }
- }
- else
- {
- if ( ((l >> 9) + 2) << 9 >= r && ((b >> 9) + 2) << 9 >= t )
- {
- index = INDICES_TO_LEVEL_4;
- maxwidth = (2 << 4) - 1;
- l >>= 9; b >>= 9;
- }
- else
- {
- index = INDICES_TO_LEVEL_3;
- maxwidth = (2 << 3) - 1;
- l >>= 10; b >>= 10;
- }
- }
- }
- else
- {
- if ( ((l >> 12) + 2) << 12 >= r && ((b >> 12) + 2) << 12 >= t )
- {
- if ( ((l >> 11) + 2) << 11 >= r && ((b >> 11) + 2) << 11 >= t )
- {
- index = INDICES_TO_LEVEL_2;
- maxwidth = (2 << 2) - 1;
- l >>= 11; b >>= 11;
- }
- else
- {
- index = INDICES_TO_LEVEL_1;
- maxwidth = (2 << 1) - 1;
- l >>= 12; b >>= 12;
- }
- }
- else
- {
- index = 0;
- maxwidth = 0;
- l = 0; b = 0;
- }
- }
- }
-
- // Make sure the ROI boundary box does not go outside the
- // global region of interest.
- if ( l > maxwidth - 1 ) l = maxwidth - 1;
- if ( b > maxwidth - 1 ) b = maxwidth - 1;
-
- return MAX_INDICES * fChamber + b * maxwidth + l + index;
-}
-
-
-AliHLTMUONCoreROI AliHLTMUONCoreRegionOfInterest::Encode(UChar& level, UInt& l, UInt& b) const
-{
-// Encodes the region of interest into a 32 bit code, and returns the
-// hierarchal level the region was encoded at and the left and right
-// grid coordinate of the bottom left corner of the region boundary box.
-
- UInt r, t, maxwidth, index;
-
- ConvertToGrid(l, r, b, t);
-
- // Work out which resolution level the location of the ROI
- // boundary box needs to be coded at. The higher the level
- // the smaller the box size: i.e. box_size = scale / 2^(level)
- // We use a binary search type method for finding the level.
- // More specificaly we search for which n,
- // (l * 2^(-n) + 2) * 2^n >= r and (b * 2^(-n) + 2) * 2^n >= t
-
- if ( ((l >> 6) + 2) << 6 >= r && ((b >> 6) + 2) << 6 >= t )
- {
- if ( ((l >> 2) + 2) << 2 >= r && ((b >> 2) + 2) << 2 >= t )
- {
- if ( l + 2 >= r && b + 2 >= t )
- {
- level = 13;
- index = INDICES_TO_LEVEL_13;
- maxwidth = (2 << 13) - 1;
- }
- else
- {
- if ( ((l >> 1) + 2) << 1 >= r && ((b >> 1) + 2) << 1 >= t )
- {
- level = 12;
- index = INDICES_TO_LEVEL_12;
- maxwidth = (2 << 12) - 1;
- l >>= 1; b >>= 1;
- }
- else
- {
- level = 11;
- index = INDICES_TO_LEVEL_11;
- maxwidth = (2 << 11) - 1;
- l >>= 2; b >>= 2;
- }
- }
- }
- else
- {
- if ( ((l >> 4) + 2) << 4 >= r && ((b >> 4) + 2) << 4 >= t )
- {
- if ( ((l >> 3) + 2) << 3 >= r && ((b >> 3) + 2) << 3 >= t )
- {
- level = 10;
- index = INDICES_TO_LEVEL_10;
- maxwidth = (2 << 10) - 1;
- l >>= 3; b >>= 3;
- }
- else
- {
- level = 9;
- index = INDICES_TO_LEVEL_9;
- maxwidth = (2 << 9) - 1;
- l >>= 4; b >>= 4;
- }
- }
- else
- {
- if ( ((l >> 5) + 2) << 5 >= r && ((b >> 5) + 2) << 5 >= t )
- {
- level = 8;
- index = INDICES_TO_LEVEL_8;
- maxwidth = (2 << 8) - 1;
- l >>= 5; b >>= 5;
- }
- else
- {
- level = 7;
- index = INDICES_TO_LEVEL_7;
- maxwidth = (2 << 7) - 1;
- l >>= 6; b >>= 6;
- }
- }
- }
- }
- else
- {
- if ( ((l >> 10) + 2) << 10 >= r && ((b >> 10) + 2) << 10 >= t )
- {
- if ( ((l >> 8) + 2) << 8 >= r && ((b >> 8) + 2) << 8 >= t )
- {
- if ( ((l >> 7) + 2) << 7 >= r && ((b >> 7) + 2) << 7 >= t )
- {
- level = 6;
- index = INDICES_TO_LEVEL_6;
- maxwidth = (2 << 6) - 1;
- l >>= 7; b >>= 7;
- }
- else
- {
- level = 5;
- index = INDICES_TO_LEVEL_5;
- maxwidth = (2 << 5) - 1;
- l >>= 8; b >>= 8;
- }
- }
- else
- {
- if ( ((l >> 9) + 2) << 9 >= r && ((b >> 9) + 2) << 9 >= t )
- {
- level = 4;
- index = INDICES_TO_LEVEL_4;
- maxwidth = (2 << 4) - 1;
- l >>= 9; b >>= 9;
- }
- else
- {
- level = 3;
- index = INDICES_TO_LEVEL_3;
- maxwidth = (2 << 3) - 1;
- l >>= 10; b >>= 10;
- }
- }
- }
- else
- {
- if ( ((l >> 12) + 2) << 12 >= r && ((b >> 12) + 2) << 12 >= t )
- {
- if ( ((l >> 11) + 2) << 11 >= r && ((b >> 11) + 2) << 11 >= t )
- {
- level = 2;
- index = INDICES_TO_LEVEL_2;
- maxwidth = (2 << 2) - 1;
- l >>= 11; b >>= 11;
- }
- else
- {
- level = 1;
- index = INDICES_TO_LEVEL_1;
- maxwidth = (2 << 1) - 1;
- l >>= 12; b >>= 12;
- }
- }
- else
- {
- level = 0;
- index = 0;
- maxwidth = 0;
- l = 0; b = 0;
- }
- }
- }
-
- // Make sure the ROI boundary box does not go outside the
- // global region of interest.
- if ( l > maxwidth - 1 ) l = maxwidth - 1;
- if ( b > maxwidth - 1 ) b = maxwidth - 1;
-
- return MAX_INDICES * fChamber + b * maxwidth + l + index;
-}
-
-
-void AliHLTMUONCoreRegionOfInterest::Decode(AliHLTMUONCoreROI code)
-{
-// Decodes a 32 bit region of interest code into this region of interest object.
-
- UInt l, r, b, t;
- UChar colevel;
- DecodeBits(code, fChamber, colevel, l, b);
-
- // Complete decoding of bottom left corner.
- b = b << colevel;
- l = l << colevel;
-
- // Decode top right corner of boundry box.
- r = l + (2 << colevel);
- t = b + (2 << colevel);
-
- ConvertBackFromGrid(l, r, b, t);
-}
-
-
-void AliHLTMUONCoreRegionOfInterest::Decode(
- AliHLTMUONCoreROI code, AliHLTMUONCoreChamberID& chamber,
- UChar& level, UInt& l, UInt& b
- )
-{
-// Decodes the 32 bit region of interest code into the chamber number,
-// hierarchal level, left and right grid coordinates of the region
-// boundary box.
-
- UChar colevel;
- DecodeBits(code, chamber, colevel, l, b);
- level = MAX_LEVELS - colevel;
-}
-
-
-void AliHLTMUONCoreRegionOfInterest::DecodeBits(
- AliHLTMUONCoreROI code, AliHLTMUONCoreChamberID& chamber,
- UChar& colevel, UInt& l, UInt& b
- )
-{
- // First decode the chamber number and the remainder
- // contains the location index.
- chamber = (AliHLTMUONCoreChamberID)(code / MAX_INDICES);
- UInt index = code % MAX_INDICES;
-
- UInt width;
-
- // Binary search which colevel the index is for.
- if (index >= INDICES_TO_LEVEL_7)
- {
- if (index >= INDICES_TO_LEVEL_11)
- {
- if (index >= INDICES_TO_LEVEL_13)
- {
- index -= INDICES_TO_LEVEL_13;
- width = (2 << 13) - 1;
- colevel = 0;
- }
- else
- {
- if (index >= INDICES_TO_LEVEL_12)
- {
- index -= INDICES_TO_LEVEL_12;
- width = (2 << 12) - 1;
- colevel = 1;
- }
- else
- {
- index -= INDICES_TO_LEVEL_11;
- width = (2 << 11) - 1;
- colevel = 2;
- }
- }
- }
- else
- {
- if (index >= INDICES_TO_LEVEL_9)
- {
- if (index >= INDICES_TO_LEVEL_10)
- {
- index -= INDICES_TO_LEVEL_10;
- width = (2 << 10) - 1;
- colevel = 3;
- }
- else
- {
- index -= INDICES_TO_LEVEL_9;
- width = (2 << 9) - 1;
- colevel = 4;
- }
- }
- else
- {
- if (index >= INDICES_TO_LEVEL_8)
- {
- index -= INDICES_TO_LEVEL_8;
- width = (2 << 8) - 1;
- colevel = 5;
- }
- else
- {
- index -= INDICES_TO_LEVEL_7;
- width = (2 << 7) - 1;
- colevel = 6;
- }
- }
- }
- }
- else
- {
- if (index >= INDICES_TO_LEVEL_3)
- {
- if (index >= INDICES_TO_LEVEL_5)
- {
- if (index >= INDICES_TO_LEVEL_6)
- {
- index -= INDICES_TO_LEVEL_6;
- width = (2 << 6) - 1;
- colevel = 7;
- }
- else
- {
- index -= INDICES_TO_LEVEL_5;
- width = (2 << 5) - 1;
- colevel = 8;
- }
- }
- else
- {
- if (index >= INDICES_TO_LEVEL_4)
- {
- index -= INDICES_TO_LEVEL_4;
- width = (2 << 4) - 1;
- colevel = 9;
- }
- else
- {
- index -= INDICES_TO_LEVEL_3;
- width = (2 << 3) - 1;
- colevel = 10;
- }
- }
- }
- else
- {
- if (index >= INDICES_TO_LEVEL_1)
- {
- if (index >= INDICES_TO_LEVEL_2)
- {
- index -= INDICES_TO_LEVEL_2;
- width = (2 << 2) - 1;
- colevel = 11;
- }
- else
- {
- index -= INDICES_TO_LEVEL_1;
- width = (2 << 1) - 1;
- colevel = 12;
- }
- }
- else
- {
- width = 1;
- colevel = 13;
- }
- }
- }
-
- // Can now decode the x, y position of the bottom left corner
- // of the ROI boundary box.
- b = (index / width);
- l = (index % width);
-}
-
-
-AliHLTMUONCoreChamberID AliHLTMUONCoreRegionOfInterest::DecodeChamber(AliHLTMUONCoreROI code)
-{
-// Decodes the chamber number of the region of interest 32 bit code.
-
- return (AliHLTMUONCoreChamberID)(code / MAX_INDICES);
-}
-
+++ /dev/null
-#ifndef ALIHLTMUONCOREREGIONOFINTEREST_H
-#define ALIHLTMUONCOREREGIONOFINTEREST_H
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * See cxx source for full Copyright notice */
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-// The region of interest object is used to encode/decode and work with boundary
-// box type regions of interest. The 32 bit ROI codes are used to communicate
-// regions of interest between different parts of the dHLT system. This is more
-// efficient than sending 20 byte long region of interest objects.
-//
-////////////////////////////////////////////////////////////////////////////////
-
-#include "AliHLTMUONBasicTypes.h"
-#include "AliHLTMUONUtils.h"
-#include "AliHLTMUONCoreCluster.h"
-
-
-typedef UInt AliHLTMUONCoreROI;
-
-enum
-{
- kInvalidROI = 0xFFFFFFFF
-};
-
-
-/* Identification numbers specifying the tracking chambers of the dimuon
- spectrometer.
- */
-enum AliHLTMUONCoreChamberID
-{
- kUnknownChamber = -1,
- kChamber1 = 0,
- kChamber2 = 1,
- kChamber3 = 2,
- kChamber4 = 3,
- kChamber5 = 4,
- kChamber6 = 5,
- kChamber7 = 6,
- kChamber8 = 7,
- kChamber9 = 8,
- kChamber10 = 9
-};
-
-
-class AliHLTMUONCoreRegionOfInterest
-{
-public:
-
- AliHLTMUONCoreRegionOfInterest()
- : fChamber(kChamber1), fLeft(0.0), fRight(0.0), fBottom(0.0), fTop(0.0)
- {
- fChamber = kChamber1;
- fLeft = fRight = fBottom = fTop = 0.0;
- }
-
- /* This constructor decodes the ROI bit pattern into a region of
- interest object.
- */
- AliHLTMUONCoreRegionOfInterest(const AliHLTMUONCoreROI& code)
- : fChamber(kChamber1), fLeft(0.0), fRight(0.0), fBottom(0.0), fTop(0.0)
- {
- Decode(code);
- }
-
- /* Creates a region of interest around the given point for the
- specified chamber.
- */
- AliHLTMUONCoreRegionOfInterest(
- const AliHLTMUONCoreClusterPoint& point0,
- AliHLTMUONCoreChamberID chamber0
- )
- : fChamber(kChamber1), fLeft(0.0), fRight(0.0), fBottom(0.0), fTop(0.0)
- {
- CreateToContain(point0, chamber0);
- }
-
- /* Creates a region of interest around all the given points and for the
- specified chamber.
- */
- AliHLTMUONCoreRegionOfInterest(
- const AliHLTMUONCoreClusterPoint* points0, UInt count0,
- AliHLTMUONCoreChamberID chamber0
- )
- : fChamber(kChamber1), fLeft(0.0), fRight(0.0), fBottom(0.0), fTop(0.0)
- {
- CreateToContain(points0, count0, chamber0);
- }
-
- /* Creates a region of interest with the specified boundaries and for
- the specified chamber.
- */
- AliHLTMUONCoreRegionOfInterest(
- Float left, Float right, Float bottom, Float top,
- AliHLTMUONCoreChamberID chamber
- );
-
- /* Checks if the point is contained in this region of interest.
- */
- bool Contains(const AliHLTMUONCoreClusterPoint& point) const;
-
- /* Checks if the point is contained in this region of interest and the
- chamber number corresponds to this region object.
- */
- bool Contains(
- const AliHLTMUONCoreClusterPoint& point,
- AliHLTMUONCoreChamberID chamber
- ) const;
-
- /* Checks if the specified region of interest is contained in this
- region of interest object.
- */
- bool Contains(const AliHLTMUONCoreRegionOfInterest& roi) const;
-
- /* Creates a region of interest around the given point for the
- specified chamber.
- */
- void CreateToContain(
- const AliHLTMUONCoreClusterPoint& point,
- AliHLTMUONCoreChamberID chamber
- );
-
- /* Extends the region of interest to contain the specified point.
- */
- void ExpandToContain(const AliHLTMUONCoreClusterPoint& point);
-
- /* Creates a region of interest around all the given points and for the
- specified chamber.
- */
- void CreateToContain(
- const AliHLTMUONCoreClusterPoint* points, UInt count,
- AliHLTMUONCoreChamberID chamber
- );
-
- /* Checks if the region of interest is within the boundaries imposed on
- the specific chamber plane. This boundary is aproximately the square
- box around the chamber's detection region.
- */
- bool InBounds();
-
- /* Encodes the region of interest into a 32 bit code.
- */
- AliHLTMUONCoreROI Encode() const;
-
- /* Encodes the region of interest into a 32 bit code, and returns the
- hierarchal level the region was encoded at and the left and right
- grid coordinate of the bottom left corner of the region boundary box.
- */
- AliHLTMUONCoreROI Encode(UChar& level, UInt& l, UInt& b) const;
-
- /* Decodes a 32 bit region of interest code into this region of interest
- object.
- */
- void Decode(AliHLTMUONCoreROI code);
-
- /* Decodes the chamber number of the region of interest 32 bit code.
- */
- static AliHLTMUONCoreChamberID DecodeChamber(AliHLTMUONCoreROI code);
-
- /* Decodes the 32 bit region of interest code into the chamber number,
- hierarchal level, left and right grid coordinates of the region
- boundary box.
- */
- static void Decode(
- AliHLTMUONCoreROI code, AliHLTMUONCoreChamberID& chamber,
- UChar& level, UInt& l, UInt& b
- );
-
- /* Returns the chamber number of the region of interest.
- */
- AliHLTMUONCoreChamberID Chamber() const { return fChamber; };
-
- /* Returns the left hand boundary of the region of interest.
- */
- Float Left() const { return fLeft; };
-
- /* Returns the right hand boundary of the region of interest.
- */
- Float Right() const { return fRight; };
-
- /* Returns the bottom boundary of the region of interest.
- */
- Float Bottom() const { return fBottom; };
-
- /* Returns the top boundary of the region of interest.
- */
- Float Top() const { return fTop; };
-
-
- /* Typecast operator for implicit typecasing to 32 bit ROI codes.
- */
- operator AliHLTMUONCoreROI () const { return Encode(); };
-
-
-private:
-
- /* Converts the internal region of interest boundary box, which is
- specified in as floats, into a regular integer grid.
- l = left boundary, r = right boundary, b = bottom boundary,
- t = top boundary.
- */
- inline void ConvertToGrid(UInt& l, UInt& r, UInt& b, UInt& t) const;
-
- /* Performs the inverse conversion of the method ConvertToGrid.
- That is converts from a regular integer grid back to the internal
- floating point boundary box specification.
- */
- inline void ConvertBackFromGrid(register UInt l, register UInt r, register UInt b, register UInt t);
-
- /* Internal method for decoding 32 bit region codes. This method is
- called by the Decode methods.
- */
- static void DecodeBits(
- AliHLTMUONCoreROI code, AliHLTMUONCoreChamberID& chamber,
- UChar& colevel, UInt& l, UInt& b
- );
-
- enum {kNumberOfTrackingChambers = 10}; // Number of tracking chambers.
-
- // Boundary box scale numbers for each chamber. These are the boundary
- // boxes around the chambers detection surface.
- static Float fgPlaneScale[kNumberOfTrackingChambers]; // scale numbers.
-
- static UInt fgIndexOffsets[13]; // Offset numbers used in the encoding and decoding process.
-
-
- AliHLTMUONCoreChamberID fChamber; // Specifies the chamber the region of interest is on.
- Float fLeft; // Left boundary of boundary box.
- Float fRight; // Right boundary of boundary box.
- Float fBottom; // Bottom boundary of boundary box.
- Float fTop; // Top boundary of boundary box.
-};
-
-
-//-----------------------------------------------------------------------------
-// Inline methods:
-
-
-inline AliHLTMUONCoreRegionOfInterest::AliHLTMUONCoreRegionOfInterest(
- Float left, Float right, Float bottom, Float top,
- AliHLTMUONCoreChamberID chamber
- )
- : fChamber(kChamber1), fLeft(0.0), fRight(0.0), fBottom(0.0), fTop(0.0)
-{
-// Creates a region of interest with the specified boundaries and for
-// the specified chamber.
-
- Assert( 0 <= chamber && chamber < (AliHLTMUONCoreChamberID)kNumberOfTrackingChambers );
- fChamber = chamber;
- fLeft = left;
- fRight = right;
- fBottom = bottom;
- fTop = top;
-}
-
-
-inline bool AliHLTMUONCoreRegionOfInterest::Contains(const AliHLTMUONCoreClusterPoint& point) const
-{
-// Checks if the point is contained in this region of interest.
-
- return fLeft <= point.X()
- && point.X() <= fRight
- && fBottom <= point.Y()
- && point.Y() <= fTop;
-}
-
-
-inline bool AliHLTMUONCoreRegionOfInterest::Contains(
- const AliHLTMUONCoreClusterPoint& point,
- AliHLTMUONCoreChamberID chamber
- ) const
-{
-// Checks if the point is contained in this region of interest and the
-// chamber number corresponds to this region object.
-
- return fLeft <= point.X()
- && point.X() <= fRight
- && fBottom <= point.Y()
- && point.Y() <= fTop
- && fChamber == chamber;
-}
-
-
-inline bool AliHLTMUONCoreRegionOfInterest::Contains(
- const AliHLTMUONCoreRegionOfInterest& roi
- ) const
-{
-// Checks if the specified region of interest is contained in this
-// region of interest object.
-
- return fChamber == roi.fChamber
- && fLeft <= roi.fLeft
- && fRight >= roi.fRight
- && fBottom <= roi.fBottom
- && fTop >= roi.fTop;
-}
-
-
-#endif // ALIHLTMUONCOREREGIONOFINTEREST_H
-
+++ /dev/null
-#ifndef ALIHLTMUONCORETRACK_H
-#define ALIHLTMUONCORETRACK_H
-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-////////////////////////////////////////////////////////////////////////////////
-
-#include "AliHLTMUONCorePoint.h"
-#include "AliHLTMUONCoreRegionOfInterest.h"
-#include "AliHLTMUONCoreTriggerRecord.h"
-
-
-typedef UInt AliHLTMUONCoreTrackID;
-
-
-struct AliHLTMUONCoreTrack
-{
-
- AliHLTMUONCoreTriggerRecordID fTriggerid;
- AliHLTMUONCoreParticleSign fSign;
- Float fP; // momentum.
- Float fPt; // transverse momentum.
- AliHLTMUONCorePoint fPoint[10]; // Computed track coordinates on the 10 tracking chambers.
- AliHLTMUONCoreROI fRegion[10]; // Regions of interest from which clusters were used to compute this track.
-
-
- AliHLTMUONCoreTrack() : fTriggerid(0), fSign(kUnknownSign), fP(0.0), fPt(0.0) {}
-
-};
-
-
-#endif // ALIHLTMUONCORETRACK_H
+++ /dev/null
-#ifndef ALIHLTMUONCORETRIGGERRECORD_H
-#define ALIHLTMUONCORETRIGGERRECORD_H
-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
-
-/* $Id$ */
-
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-////////////////////////////////////////////////////////////////////////////////
-
-
-#include "AliHLTMUONBasicTypes.h"
-#include "AliHLTMUONCorePoint.h"
-
-
-/* The sign of the particle as given by L0.
- */
-enum AliHLTMUONCoreParticleSign
-{
- kSignMinus = -1,
- kUnknownSign = 0,
- kSignPlus = 1
-};
-
-
-typedef UInt AliHLTMUONCoreTriggerRecordID;
-
-
-/* Data structure containing information about L0 validated trigger hit.
- */
-struct AliHLTMUONCoreTriggerRecord
-{
- AliHLTMUONCoreParticleSign fSign; // The sign of the particle.
- Float fPt; // Transverse momentum of the particle.
- AliHLTMUONCorePoint fStation1impact; // Impact point of particle on trigger station 1.
- AliHLTMUONCorePoint fStation2impact; // Impact point of particle on trigger station 2.
-
-
- /* Default constructor.
- Sets the fSign to UnknownSign, fPt to -1 and the impact points are set to zero.
- */
- AliHLTMUONCoreTriggerRecord()
- : fSign(kUnknownSign), fPt(-1.0), fStation1impact(), fStation2impact() {}
-
- /* Creates a trigger record with the specifed particle sign, pt and impact points.
- The impactpoint1 corresponds to trigger station 1 and simmilarly impactpoint2
- corresponds to station 2.
- */
- AliHLTMUONCoreTriggerRecord(
- const AliHLTMUONCoreParticleSign sign, const Float pt,
- const AliHLTMUONCorePoint impactpoint1, const AliHLTMUONCorePoint impactpoint2
- )
- : fSign(sign), fPt(pt), fStation1impact(impactpoint1), fStation2impact(impactpoint2)
- {}
-};
-
-
-#endif // ALIHLTMUONCORETRIGGERRECORD_H
+++ /dev/null
-////////////////////////////////////////////////////////////////////////////////
-//
-// Author: Artur Szostak
-// Email: artur@alice.phy.uct.ac.za | artursz@iafrica.com
-//
-////////////////////////////////////////////////////////////////////////////////
-
-#ifndef dHLT_DECISION_RECORD_HPP
-#define dHLT_DECISION_RECORD_HPP
-
-#include "AliHLTMUONBasicTypes.h"
-#include "AliHLTMUONCoreTrack.h"
-
-namespace dHLT
-{
-
-
-enum DecisionType
-{
- UnknownDecision = 0,
- EventAccepted,
- EventRejected
-};
-
-
-struct DecisionRecord
-{
- DecisionType type;
-};
-
-
-} // dHLT
-
-#endif // dHLT_DECISION_RECORD_HPP
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff