Int_t nsm=0;
Long64_t totdata[72]={0};
- for (Int_t i=0;i<nf; i++)
+ for (Int_t ii=0;ii<nf; ii++)
{
AliSurveyObj *so = new AliSurveyObj();
- const Char_t *nome=namefiles[i];
+ const Char_t *nome=namefiles[ii];
so->FillFromLocalFile(nome);
TObjArray *points = so->GetData();
Int_t nSurveyPoint=points->GetEntries();
- for(Int_t i=0;i<nSurveyPoint;i++){
- const char* pointName= ((AliSurveyPoint *) points->At(i))->GetPointName().Data();
+ for(Int_t jj=0;jj<nSurveyPoint;jj++){
+ const char* pointName= ((AliSurveyPoint *) points->At(jj))->GetPointName().Data();
nfm=atoi(&pointName[6]);
nsm=atoi(&pointName[2]);
- data[nsm*4+nfm][0][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(i))->GetX();
- data[nsm*4+nfm][2][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(i))->GetY();
- data[nsm*4+nfm][4][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(i))->GetZ();
- data[nsm*4+nfm][1][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(i))->GetPrecisionX();
- data[nsm*4+nfm][3][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(i))->GetPrecisionY();
- data[nsm*4+nfm][5][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(i))->GetPrecisionZ();
+ data[nsm*4+nfm][0][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(jj))->GetX();
+ data[nsm*4+nfm][2][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(jj))->GetY();
+ data[nsm*4+nfm][4][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(jj))->GetZ();
+ data[nsm*4+nfm][1][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(jj))->GetPrecisionX();
+ data[nsm*4+nfm][3][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(jj))->GetPrecisionY();
+ data[nsm*4+nfm][5][totdata[nsm*4+nfm]]=((AliSurveyPoint *) points->At(jj))->GetPrecisionZ();
totdata[nsm*4+nfm]=totdata[nsm*4+nfm]+1;
}
delete so;
}
-
for(Int_t i=0; i<72 ;i++){
Float_t numx=0, numy=0,numz=0, comodox=0, comodoy=0, comodoz=0,denx=0, deny=0, denz=0;
if(totdata[i]!=0){
if (fPackedAcquisition) {
- for (Int_t jj=0; jj<=localIndex; jj++) {
+ for (jj=0; jj<=localIndex; jj++) {
fIndex++;
buf[fIndex] = localBuffer[jj];
localBuffer[jj] = 0;
AliInfo(Form(" In the TOF DDL %2i, TRM %2i, TDC %2i, chain %1i, the maximum number of t.o.f. good measurements per channel is %2i",
nDDL, nTRM, iChain, nTDC, iCH, maxMeasuresPerChannelInTDC));
*/
- for (Int_t jj=0; jj<=localIndex; jj++) {
+ for (jj=0; jj<=localIndex; jj++) {
fIndex++;
buf[fIndex] = localBuffer[jj];
localBuffer[jj] = 0;
}
for (Int_t irun=0;irun<fNruns;irun++){
- Int_t i = ch[ich]+irun*fNChannels;
+ i = ch[ich]+irun*fNChannels;
//fTree->GetEntry(i);
fChain->GetEntry(i);
ntracksRun=nentries/3;
TProfile* hSlewingProf = new TProfile("hSlewingProf", "hSlewingProf",nusefulbins, binsProfile, "G"); // CHECK THE BUILD OPTION, PLEASE!!!!!!
for (Int_t irun=0;irun<fNruns;irun++){
- Int_t i = ch[ich]+irun*fNChannels;
+ i = ch[ich]+irun*fNChannels;
//fTree->GetEntry(i);
fChain->GetEntry(i);
ntracksRun=nentries/3;
Double_t yc=(c->GetPhi() - trackTOFin->GetAlpha())*c->GetR();
Double_t p[2]={yc, c->GetZ()};
- Double_t cov[3]= {dY*dY/12., 0., dZ*dZ/12.};
- if (trackTOFin->AliExternalTrackParam::GetPredictedChi2(p,cov) > maxChi2)continue;
+ Double_t cov2[3]= {dY*dY/12., 0., dZ*dZ/12.};
+ if (trackTOFin->AliExternalTrackParam::GetPredictedChi2(p,cov2) > maxChi2)continue;
clind[0][nc] = c->GetDetInd(0);
clind[1][nc] = c->GetDetInd(1);
// Some init
Int_t index[1000];
+ Float_t quality[1000];
Float_t dist3D[1000][6];
Double_t times[1000][6];
Float_t mintimedist[1000];
Float_t likelihood[1000];
Float_t length[1000];
AliTOFcluster *clusters[1000];
- Double_t tpcpid[5];
+ Double_t tpcpid[5];
dist3D[0][0]=1;
for (Int_t i=0; i<fNseedsTOF; i++) {
Double_t dph=TMath::Abs(c->GetPhi()-phi);
if (dph>TMath::Pi()) dph-=2.*TMath::Pi();
if (TMath::Abs(dph)>dphi) continue;
-
+
clind[0][nc] = c->GetDetInd(0);
clind[1][nc] = c->GetDetInd(1);
clind[2][nc] = c->GetDetInd(2);
clind[3][nc] = c->GetDetInd(3);
clind[4][nc] = c->GetDetInd(4);
- clind[5][nc] = k;
+ clind[5][nc] = k;
nc++;
}
//
//choose the best cluster
//
- Float_t quality[1000];
- Int_t index[1000];
+ //Float_t quality[1000];
+ //Int_t index[1000];
+ for (Int_t kk=0; kk<1000; kk++) quality[kk]=0;
//
AliTOFcluster * cgold=0;
Int_t igold =-1;