const Float_t AliTOFGeometry::fgkToTBin = 48.8; // time-over-threshold bin width [ps]
const Float_t AliTOFGeometry::fgkBunchCrossingBin = fgkTdcBin * 1024; // bunch-crossing bin width [ps]
+const Float_t AliTOFGeometry::fgkSlewTOTMin = 10.; // min TOT for slewing correction [ns]
+const Float_t AliTOFGeometry::fgkSlewTOTMax = 16.; // max TOT for slewing correction [ns]
+
const Float_t AliTOFGeometry::fgkDeadTime = 25E+03; // Single channel dead time (ps)
const Float_t AliTOFGeometry::fgkMatchingWindow = fgkTdcBin*TMath::Power(2,13); // Matching window (ps)
}
//_____________________________________________________________________________
-Bool_t AliTOFGeometry::IsInsideThePad(TGeoHMatrix mat, const Float_t * const pos, Float_t *dist3d) const
+Bool_t AliTOFGeometry::IsInsideThePad(TGeoHMatrix *mat, const Float_t * const pos, Float_t *dist3d) const
{
//
// Returns true if space point with coor pos (x,y,z) [cm] falls inside
// from ALICE global reference system
// towards TOF pad reference system
Double_t posl[3] = {0., 0., 0.};
- mat.MasterToLocal(posg,posl);
+ mat->MasterToLocal(posg,posl);
Float_t xr = posl[0];
Float_t yr = posl[1];
Double_t padg[3] = {0., 0., 0.};
// from TOF pad local reference system
// towards ALICE global reference system
- TGeoHMatrix inverse = mat.Inverse();
+ TGeoHMatrix inverse = mat->Inverse();
inverse.MasterToLocal(padl,padg);
// returns the 3d distance
// This function returns the colume path of a given pad
//--------------------------------------------------------------------
Int_t sector = ind[0];
- Char_t string1[100];
- Char_t string2[100];
- Char_t string3[100];
+
+ const Int_t kSize = 100;
+
+ Char_t string1[kSize];
+ Char_t string2[kSize];
+ Char_t string3[kSize];
Int_t icopy=-1;
icopy=sector;
- sprintf(string1,"/ALIC_1/B077_1/BSEGMO%i_1/BTOF%i_1",icopy,icopy);
+ snprintf(string1,kSize,"/ALIC_1/B077_1/BSEGMO%i_1/BTOF%i_1",icopy,icopy);
Int_t iplate=ind[1];
Int_t istrip=ind[2];
if( iplate==3) icopy=istrip+NStripC()+NStripB()+NStripA();
if( iplate==4) icopy=istrip+NStripC()+2*NStripB()+NStripA();
icopy++;
- sprintf(string2,"FTOA_0/FLTA_0/FSTR_%i",icopy);
+ snprintf(string2,kSize,"FTOA_0/FLTA_0/FSTR_%i",icopy);
if(fHoles && (sector==13 || sector==14 || sector==15)){
- if(iplate<2) sprintf(string2,"FTOB_0/FLTB_0/FSTR_%i",icopy);
- if(iplate>2) sprintf(string2,"FTOC_0/FLTC_0/FSTR_%i",icopy);
+ if(iplate<2) snprintf(string2,kSize,"FTOB_0/FLTB_0/FSTR_%i",icopy);
+ if(iplate>2) snprintf(string2,kSize,"FTOC_0/FLTC_0/FSTR_%i",icopy);
}
Int_t padz = ind[3]+1;
Int_t padx = ind[4]+1;
- sprintf(string3,"FPCB_1/FSEN_1/FSEZ_%i/FPAD_%i",padz,padx);
- sprintf(path,"%s/%s/%s",string1,string2,string3);
+ snprintf(string3,kSize,"FPCB_1/FSEN_1/FSEZ_%i/FPAD_%i",padz,padx);
+ snprintf(path,2*kSize,"%s/%s/%s",string1,string2,string3);
}
//_____________________________________________________________________________
// This function returns the colume path of a given sector
//--------------------------------------------------------------------
- Char_t string[100];
+ const Int_t kSize = 100;
+
+ Char_t string[kSize];
Int_t icopy = sector;
- sprintf(string,"/ALIC_1/B077_1/BSEGMO%i_1/BTOF%i_1",icopy,icopy);
- sprintf(path,"%s",string);
+ snprintf(string,kSize,"/ALIC_1/B077_1/BSEGMO%i_1/BTOF%i_1",icopy,icopy);
+ snprintf(path,2*kSize,"%s",string);
}
//_____________________________________________________________________________
// This function returns the colume path of a given strip
//--------------------------------------------------------------------
- Char_t string1[100];
- Char_t string2[100];
- Char_t string3[100];
+ const Int_t kSize = 100;
+
+ Char_t string1[kSize];
+ Char_t string2[kSize];
+ Char_t string3[kSize];
Int_t icopy = sector;
- sprintf(string1,"/ALIC_1/B077_1/BSEGMO%i_1/BTOF%i_1",icopy,icopy);
+ snprintf(string1,kSize,"/ALIC_1/B077_1/BSEGMO%i_1/BTOF%i_1",icopy,icopy);
if(plate==0) icopy=strip;
if(plate==1) icopy=strip+NStripC();
if(plate==3) icopy=strip+NStripC()+NStripB()+NStripA();
if(plate==4) icopy=strip+NStripC()+2*NStripB()+NStripA();
icopy++;
- sprintf(string2,"FTOA_0/FLTA_0/FSTR_%i",icopy);
+ snprintf(string2,kSize,"FTOA_0/FLTA_0/FSTR_%i",icopy);
if(fHoles && (sector==13 || sector==14 || sector==15)){
- if(plate<2) sprintf(string2,"FTOB_0/FLTB_0/FSTR_%i",icopy);
- if(plate>2) sprintf(string2,"FTOC_0/FLTC_0/FSTR_%i",icopy);
+ if(plate<2) snprintf(string2,kSize,"FTOB_0/FLTB_0/FSTR_%i",icopy);
+ if(plate>2) snprintf(string2,kSize,"FTOC_0/FLTC_0/FSTR_%i",icopy);
}
- sprintf(string3,"FPCB_1/FSEN_1");
- sprintf(path,"%s/%s/%s",string1,string2,string3);
+ snprintf(string3,kSize,"FPCB_1/FSEN_1");
+ snprintf(path,2*kSize,"%s/%s/%s",string1,string2,string3);
}
//_____________________________________________________________________________
// Returns space point coor (x,y,z) (cm) for Detector
// Indices (iSect,iPlate,iStrip,iPadX,iPadZ)
//
- Char_t path[100];
- GetVolumePath(det,path );
+ Char_t path[200];
+ GetVolumePath(det,path);
if (!gGeoManager) {
printf("ERROR: no TGeo\n");
}
Translation(posLocal,step);
// FSTR reference frame -> FTOA/B/C = FLTA/B/C reference frame
- Double_t angles[6];
+ Double_t angles[6] = {0.,0.,0.,0.,0.,0.};
if (GetAngles(iplate,istrip) >0.) {
angles[0] = 90.;
angles[1] = 0.;
// FSTR reference frame -> FTOA/B/C = FLTA/B/C reference frame
- Double_t angles[6];
+ Double_t angles[6] = {0.,0.,0.,0.,0.,0.};
if (GetAngles(iplate,istrip) >0.) {
angles[0] = 90.;
angles[1] = 0.;
Translation(posLocal,step);
// FSTR reference frame -> FTOA/B/C = FLTA/B/C reference frame
- Double_t angles[6];
+ Double_t angles[6] = {0.,0.,0.,0.,0.,0.};
if (GetAngles(iplate,istrip) >0.) {
angles[0] = 90.;
angles[1] = 0.;
}
//_____________________________________________________________________________
-void AliTOFGeometry::DetToSectorRF(Int_t vol[5], Double_t **coord)
+void AliTOFGeometry::DetToSectorRF(Int_t vol[5], Double_t coord[4][3])
{
//
// Returns the local coordinates (x, y, z) in sector reference frame
if (!gGeoManager) printf("ERROR: no TGeo\n");
// ALICE -> TOF Sector
- Char_t path1[100]="";
+ Char_t path1[200];
GetVolumePath(vol[0],path1);
gGeoManager->cd(path1);
TGeoHMatrix aliceToSector;
//TGeoHMatrix sectorToALICE = aliceToSector.Inverse();
// ALICE -> TOF Pad
- Char_t path2[100]="";
+ Char_t path2[200];
GetVolumePath(vol,path2);
gGeoManager->cd(path2);
TGeoHMatrix aliceToPad;
for (Int_t jj=0; jj<4; jj++) padToSector.MasterToLocal(&cornerPad[jj][0], &coord[jj][0]);
- delete cornerPad;
+ delete [] cornerPad;
//sectorToPad.LocalToMaster(cornerPad, coord);
}
//-------------------------------------------------------------------------
+void AliTOFGeometry::GetStripAndModule(Int_t iStripPerSM, Int_t &iplate, Int_t &istrip)
+{
+ //
+ // Convert the serial number of the TOF strip number iStripPerSM [0,90]
+ // in module number iplate [0,4] and strip number istrip [0,14/18].
+ //
+
+ if (iStripPerSM<0 || iStripPerSM>=kNStripC+kNStripB+kNStripA+kNStripB+kNStripC) {
+ iplate = -1;
+ istrip = -1;
+ }
+ else if (iStripPerSM<kNStripC) {
+ iplate = 0;
+ istrip = iStripPerSM;
+ }
+ else if (iStripPerSM>=kNStripC && iStripPerSM<kNStripC+kNStripB) {
+ iplate = 1;
+ istrip = iStripPerSM-kNStripC;
+ }
+ else if (iStripPerSM>=kNStripC+kNStripB && iStripPerSM<kNStripC+kNStripB+kNStripA) {
+ iplate = 2;
+ istrip = iStripPerSM-kNStripC-kNStripB;
+ }
+ else if (iStripPerSM>=kNStripC+kNStripB+kNStripA && iStripPerSM<kNStripC+kNStripB+kNStripA+kNStripB) {
+ iplate = 3;
+ istrip = iStripPerSM-kNStripC-kNStripB-kNStripA;
+ }
+ else if (iStripPerSM>=kNStripC+kNStripB+kNStripA+kNStripB && iStripPerSM<kNStripC+kNStripB+kNStripA+kNStripB+kNStripC) {
+ iplate = 4;
+ istrip = iStripPerSM-kNStripC-kNStripB-kNStripA-kNStripB;
+ }
+
+
+}
+//-------------------------------------------------------------------------
+
Int_t AliTOFGeometry::GetStripNumberPerSM(Int_t iplate, Int_t istrip)
{
//
l2t.MasterToLocal(padCentreL,padCentreT);
- Char_t path[100];
+ Char_t path[200];
// pad centre coordinates in its ref. frame
Double_t padCentreL2[3] = {0., 0., 0.};
// pad centre coordinates in the ALICE global ref. frame
l2g.MasterToLocal(padCentreL2,padCentreG);
- Char_t path2[100];
+ Char_t path2[200];
// strip centre coordinates in its ref. frame
Double_t stripCentreL[3] = {0., 0., 0.};
// strip centre coordinates in the ALICE global ref. frame
for (Int_t ii=0; ii<3; ii++) differenceT[ii] = (Float_t)(posT[ii]-padCentreT[ii]);
}
+//-------------------------------------------------------------------------
+
+Int_t AliTOFGeometry::GetTOFsupermodule(const Int_t index)
+{
+ // Return the TOF supermodule where TOF channel index is located
+
+ if (index<0 || index>=NPadXSector()*NSectors()) return -1;
+ else return index/NpadXStrip()/NStripXSector();
+
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff