fNy(0),
fChargeRatio(0),
fType(0),
-fDeltaProb(0)
+fDeltaProb(0),
+fDriftTime(0.)
{
// default constructor
}
fNy(info[0]),
fChargeRatio(0),
fType(0),
-fDeltaProb(0)
+fDeltaProb(0),
+fDriftTime(0.)
{
//standard constructor used in AliITSClusterFinderV2
Double_t lxyz[3] = {0, 0, 0};
GetTracking2LocalMatrix()->LocalToMaster(txyz,lxyz);
fXloc = lxyz[0]; fZloc = lxyz[2];
+ if(fLayer==4) hit[5]=-hit[5];
+ if( (fLayer==4) || (fLayer==5) ) SetSigmaYZ(hit[5]);
}
else {
switch (fLayer) {
fdEdX=fQ*1e-6;
break;
case 4:
+ fdEdX=fQ*2.16;
+ SetSigmaYZ(hit[5]);
case 5:
fdEdX=fQ*2.16;
+ hit[5]=-hit[5];
+ SetSigmaYZ(hit[5]);
break;
default:
AliError(Form("Wrong ITS layer %d (0 -> 5)",fLayer));
fNy(pt.fNy),
fChargeRatio(pt.fChargeRatio),
fType(pt.fType),
-fDeltaProb(pt.fDeltaProb)
+fDeltaProb(pt.fDeltaProb),
+fDriftTime(pt.fDriftTime)
{
//Copy constructor
*os << fXloc << " " << fZloc << " " << fdEdX << " ";
fmt = os->setf(ios::fixed); // every fixed
*os << fIndex <<" " << fQ << " "<<fLayer <<" "<<fNz<<" "<<fNy<<" ";
- *os << fChargeRatio<<" " << fType << " " << fDeltaProb;
+ *os << fChargeRatio<<" " << fType << " " << fDeltaProb << " " << fDriftTime;
os->flags(fmt); // reset back to old formating.
return;
}
+
+//----------------------------------------------------------------------
+Int_t AliITSRecPoint::GetNpixels() const {
+//
+// returns the number of pixels used for the SPD clusters
+//
+
+ if(fLayer > 1) return -1;
+ else return fType;
+
+}
+
+//----------------------------------------------------------------------
+Int_t AliITSRecPoint::GetSPDclusterType() const {
+//
+// 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;
+ if(fLayer > 1) return type;
+ else {
+
+ switch (fType) {
+ case 1 : type = 1 ;break;
+ case 2 : if(fNy == 2) type = 2;
+ else type = 3;
+ break;
+ case 3 : if(fNy == 3) type = 4;
+ else if(fNz == 3) type = 6;
+ else type = 5;
+ break;
+ case 4 : if(fNz == 1) type = 7;
+ 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 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 = fType+1000*fNy+1000000*fNz;
+ break;
+ case 6 : if(fNy ==3 && fNz == 2) type = 13;
+ if(fNy ==2 && fNz == 3) type = 14;
+ 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;
+ return typ;
+}
+
+//----------------------------------------------------------------------
+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;
+ *is >> fDeltaProb >> fDriftTime;
return;
}