#include "AliITSRecPoint.h"
#include "AliAlignObj.h"
+using std::ios;
ClassImp(AliITSRecPoint)
//_____________________________________________________________
-AliITSRecPoint::AliITSRecPoint(): AliCluster(),
+AliITSRecPoint::AliITSRecPoint():
+AliCluster(),
fXloc(0),
fZloc(0),
fdEdX(0),
fChargeRatio(0),
fType(0),
fDeltaProb(0),
-fDriftTime(0.)
+fDriftTime(0.),
+fDriftSide(0)
{
// default constructor
}
fChargeRatio(0),
fType(0),
fDeltaProb(0),
-fDriftTime(0.)
+fDriftTime(0.),
+fDriftSide(0)
{
//standard constructor used in AliITSClusterFinderV2
case 4:
fdEdX=fQ*2.16;
SetSigmaYZ(hit[5]);
+ break;
case 5:
fdEdX=fQ*2.16;
hit[5]=-hit[5];
}
//_______________________________________________________________________
-AliITSRecPoint::AliITSRecPoint(const AliITSRecPoint& pt):AliCluster(pt),
-fXloc(pt.fXloc),
-fZloc(pt.fZloc),
-fdEdX(pt.fdEdX),
-fIndex(pt.fIndex),
-fQ(pt.fQ),
-fLayer(pt.fLayer),
-fNz(pt.fNz),
-fNy(pt.fNy),
-fChargeRatio(pt.fChargeRatio),
-fType(pt.fType),
-fDeltaProb(pt.fDeltaProb),
-fDriftTime(pt.fDriftTime)
+AliITSRecPoint::AliITSRecPoint(const AliITSRecPoint& pt):
+ AliCluster(pt),
+ fXloc(pt.fXloc),
+ fZloc(pt.fZloc),
+ fdEdX(pt.fdEdX),
+ fIndex(pt.fIndex),
+ fQ(pt.fQ),
+ fLayer(pt.fLayer),
+ fNz(pt.fNz),
+ fNy(pt.fNy),
+ fChargeRatio(pt.fChargeRatio),
+ fType(pt.fType),
+ fDeltaProb(pt.fDeltaProb),
+ fDriftTime(pt.fDriftTime),
+ fDriftSide(pt.fDriftSide)
{
//Copy constructor
*os << fXloc << " " << fZloc << " " << fdEdX << " ";
fmt = os->setf(ios::fixed); // every fixed
*os << fIndex <<" " << fQ << " "<<fLayer <<" "<<fNz<<" "<<fNy<<" ";
- *os << fChargeRatio<<" " << fType << " " << fDeltaProb << " " << fDriftTime;
+ *os << fChargeRatio<<" " << fType << " " << fDeltaProb << " " << fDriftTime<< " " << fDriftSide;
os->flags(fmt); // reset back to old formating.
return;
}
//----------------------------------------------------------------------
Int_t AliITSRecPoint::GetSPDclusterType() const {
//
-// returns the type of cluster according to conventional numbering
+// returns an Int_t with encoded information on cluster size
+// type <= 16: cluster type identifier according to conventional numbering
+// type > 16: Npixels+1000*Ny+1000000*Nz
//
Int_t type = -1;
else if(fNz == 2 && fNy == 2) type = 8;
else if(fNy == 2 && fNz == 3) type = 11;
else if(fNy == 3 && fNz == 2) type = 9;
- else if(fNz == 4) type = 15;
- else type = 17;
+ else type = 15;
break;
case 5 : if(fNy == 3 && fNz == 2) type = 10;
if(fNy == 2 && fNz == 3 ) type = 12;
if(fNy == 5) type = 16;
- else type = 17;
+ else type = fType+1000*fNy+1000000*fNz;
break;
case 6 : if(fNy ==3 && fNz == 2) type = 13;
if(fNy ==2 && fNz == 3) type = 14;
- else type = 17;
+ else type = fType+1000*fNy+1000000*fNz;
+ break;
+ default: type = fType+1000*fNy+1000000*fNz;
+ break;
}
return type;
}
+}
- }
+//----------------------------------------------------------------------
+Int_t AliITSRecPoint::GetSDDclusterType() const {
+// returns an Int_t with encoded information on cluster size
+// Byte1 = fNz Byte0=fNy, other two bytes empty for extra information
+// max. allowed cluster size = 255
+ Int_t typ=(fNz&0xFF)<<8;
+ typ+=fNy&0xFF;
+ if(fDriftSide==1) typ+=1<<16;
+ return typ;
+}
+//----------------------------------------------------------------------
+void AliITSRecPoint::DecodeSDDclusterType(Int_t cluType, Int_t &cluSizAn, Int_t& cluSizTb, Int_t &drSide){
+// Extract cluster sizes and drift side from cluster type
+ cluSizTb=cluType&0xFF;
+ cluSizAn=(cluType>>8)&0xFF;
+ drSide=(cluType>>16);
+ return;
+}
+//----------------------------------------------------------------------
+Int_t AliITSRecPoint::GetSSDclusterType() const {
+// returns an Int_t with encoded information on cluster size
+// Byte1 = fNz Byte0=fNy, other two bytes empty for extra information
+// max. allowed cluster size = 255
+ Int_t typ=(fNz&0xFF)<<8;
+ typ+=fNy&0xFF;
+ return typ;
+}
//----------------------------------------------------------------------
void AliITSRecPoint::Read(istream *is){
SetVolumeId(lab[0]);// fIsMisalinged = mis;
*is >> fXloc >> fZloc >> fdEdX;
*is >> fIndex >> fQ >> fLayer >> fNz >> fNy >> fChargeRatio >> fType;
- *is >> fDeltaProb >> fDriftTime;
+ *is >> fDeltaProb >> fDriftTime >> fDriftSide;
return;
}