//
}
+
//______________________________________________________________________________
AliTOFClusterFinder::AliTOFClusterFinder(AliRunLoader* runLoader, AliTOFcalib *calib):
}
}
- // AliDebug(2,Form(" %2i %1i %2i %1i %2i ",dig[0],dig[1],dig[2],dig[3],dig[4]));
+ AliDebug(2,Form(" %2d %1d %2d %1d %2d ",dig[0],dig[1],dig[2],dig[3],dig[4]));
parTOF[0] = d->GetTdc(); //the TDC signal
parTOF[1] = d->GetToT(); //the ToT signal
GetClusterPars(dig, posClus,covClus);
AliTOFcluster *tofCluster = new AliTOFcluster(volIdClus,posClus[0],posClus[1],posClus[2],covClus[0],covClus[1],covClus[2],covClus[3],covClus[4],covClus[5],d->GetTracks(),dig,parTOF,status,ii);
InsertCluster(tofCluster);
-
+
}
- AliInfo(Form("Number of found clusters: %i for event: %i", fNumberOfTofClusters, iEvent));
+ AliInfo(Form("Number of found clusters: %d for event: %d", fNumberOfTofClusters, iEvent));
CalibrateRecPoint();
FillRecPoint();
}
}
- // AliDebug(2,Form(" %2i %1i %2i %1i %2i ",dig[0],dig[1],dig[2],dig[3],dig[4]));
+ // AliDebug(2,Form(" %2d %1d %2d %1d %2d ",dig[0],dig[1],dig[2],dig[3],dig[4]));
parTOF[0] = d->GetTdc(); //the TDC signal
parTOF[1] = d->GetToT(); //the ToT signal
Double_t covClus[6];
UShort_t volIdClus=GetClusterVolIndex(dig);
GetClusterPars(dig,posClus,covClus);
- AliTOFcluster *tofCluster = new AliTOFcluster(volIdClus,posClus[0],posClus[1],posClus[2],covClus[0],covClus[1],covClus[2],covClus[3],covClus[4],covClus[5],d->GetTracks(),dig,parTOF,status,ii);
+ AliTOFcluster *tofCluster = new AliTOFcluster(volIdClus,posClus[0],posClus[1],posClus[2],covClus[0],covClus[1],covClus[2],covClus[3],covClus[4],covClus[5],d->GetTracks(),dig,parTOF,status,ii);
InsertCluster(tofCluster);
}
- AliInfo(Form("Number of found clusters: %i", fNumberOfTofClusters));
+ AliInfo(Form("Number of found clusters: %d", fNumberOfTofClusters));
CalibrateRecPoint();
FillRecPoint();
Int_t inholes = 0;
- //const Int_t kDDL = fTOFGeometry->NDDL()*fTOFGeometry->NSectors();
const Int_t kDDL = AliDAQ::NumberOfDdls("TOF");
ResetRecpoint();
Int_t lab[3]={-1,-1,-1};
Bool_t status=kTRUE;
GetClusterPars(detectorIndex,posClus,covClus);
- AliTOFcluster *tofCluster = new AliTOFcluster(volIdClus,posClus[0],posClus[1],posClus[2],covClus[0],covClus[1],covClus[2],covClus[3],covClus[4],covClus[5],lab,detectorIndex,parTOF,status,-1);
+ AliTOFcluster *tofCluster = new AliTOFcluster(volIdClus,posClus[0],posClus[1],posClus[2],covClus[0],covClus[1],covClus[2],covClus[3],covClus[4],covClus[5],lab,detectorIndex,parTOF,status,-1);
InsertCluster(tofCluster);
if (fVerbose==2) {
if (fVerbose==2) ftxt.close();
- AliInfo(Form("Number of found clusters: %i", fNumberOfTofClusters));
+ AliInfo(Form("Number of found clusters: %d", fNumberOfTofClusters));
CalibrateRecPoint();
FillRecPoint();
Int_t inholes = 0;
- //const Int_t kDDL = fTOFGeometry->NDDL()*fTOFGeometry->NSectors();
const Int_t kDDL = AliDAQ::NumberOfDdls("TOF");
fRunLoader->GetEvent(iEvent);
- AliDebug(2,Form(" Event number %2i ", iEvent));
+ AliDebug(2,Form(" Event number %2d ", iEvent));
fTreeR = fTOFLoader->TreeR();
if (fVerbose==2) ftxt.close();
- AliInfo(Form("Number of found clusters: %i for event: %i", fNumberOfTofClusters, iEvent));
+ AliInfo(Form("Number of found clusters: %d for event: %d", fNumberOfTofClusters, iEvent));
CalibrateRecPoint();
FillRecPoint();
TStopwatch stopwatch;
stopwatch.Start();
- //const Int_t kDDL = fTOFGeometry->NDDL()*fTOFGeometry->NSectors();
const Int_t kDDL = AliTOFGeometry::NDDL()*AliTOFGeometry::NSectors();
fRunLoader->GetEvent(iEvent);
fRunLoader->GetEvent(iEvent);
- AliDebug(2,Form(" Event number %2i ", iEvent));
+ AliDebug(2,Form(" Event number %2d ", iEvent));
TClonesArray * clonesRawData;
fTOFLoader = fRunLoader->GetLoader("TOFLoader");
fTOFLoader->WriteDigits("OVERWRITE");
-
+
+ delete tofDigits;
+
AliDebug(1, Form("Execution time to read TOF raw data and to write TOF clusters : R:%.2fs C:%.2fs",
stopwatch.RealTime(),stopwatch.CpuTime()));
Int_t bufsize = 32000;
digitsTree->Branch("TOF", &tofDigits, bufsize);
- /// fRunLoader->GetEvent(iEvent);
-
- /// AliDebug(2,Form(" Event number %2i ", iEvent));
-
TClonesArray * clonesRawData;
Int_t dummy = -1;
digitsTree->Fill();
- AliDebug(1, Form("Got %d digits: ", tofDigits->GetEntries()));
+ delete tofDigits;
+
+ AliDebug(1, Form("Got %d digits: ", fDigits->GetEntries()));
AliDebug(1, Form("Execution time to read TOF raw data and fill TOF digit tree : R:%.2fs C:%.2fs",
stopwatch.RealTime(),stopwatch.CpuTime()));
GetClusterPars(detectorIndex,posClus,covClus);
new(lRecPoints[ii]) AliTOFcluster(volIdClus,posClus[0],posClus[1],posClus[2],covClus[0],covClus[1],covClus[2],covClus[3],covClus[4],covClus[5],trackLabels,detectorIndex, parTOF,status,digitIndex);
- AliDebug(2, Form(" %4i %4i %f %f %f %f %f %f %f %f %f %3i %3i %3i %2i %1i %2i %1i %2i %4i %3i %3i %4i %4i %1i %4i",
+ AliDebug(2, Form(" %4d %4d %f %f %f %f %f %f %f %f %f %3d %3d %3d %2d %1d %2d %1d %2d %4d %3d %3d %4d %4d %1d %4d",
ii, volIdClus, posClus[0], posClus[1], posClus[2],
fTofClusters[ii]->GetSigmaX2(),
fTofClusters[ii]->GetSigmaXY(),
//check the status, also unknown is fine!!!!!!!
- AliDebug(2, Form(" Status for channel %i = %i",index, (Int_t)status));
+ AliDebug(2, Form(" Status for channel %d = %d",index, (Int_t)status));
if((statusPulser & AliTOFChannelOnlineStatusArray::kTOFPulserBad)==(AliTOFChannelOnlineStatusArray::kTOFPulserBad)||(statusNoise & AliTOFChannelOnlineStatusArray::kTOFNoiseBad)==(AliTOFChannelOnlineStatusArray::kTOFNoiseBad)||(statusHW & AliTOFChannelOnlineStatusArray::kTOFHWBad)==(AliTOFChannelOnlineStatusArray::kTOFHWBad)){
- AliDebug(2, Form(" Bad Status for channel %i",index));
+ AliDebug(2, Form(" Bad Status for channel %d",index));
fTofClusters[ii]->SetStatus(kFALSE); //odd convention, to avoid conflict with calibration objects currently in the db (temporary solution).
}
else {
- AliDebug(2, Form(" Good Status for channel %i",index));
+ AliDebug(2, Form(" Good Status for channel %d",index));
}
// Get Rough channel online equalization
Double_t roughDelay=(Double_t)calDelay->GetDelay(index); // in ns
par[j]=(Double_t)calChannelOffline->GetSlewPar(j);
}
AliDebug(2,Form(" Calib Pars = %f, %f, %f, %f, %f, %f ",par[0],par[1],par[2],par[3],par[4],par[5]));
- AliDebug(2,Form(" The ToT and Time, uncorr (counts) = %i , %i", fTofClusters[ii]->GetToT(),fTofClusters[ii]->GetTDC()));
+ AliDebug(2,Form(" The ToT and Time, uncorr (counts) = %d , %d", fTofClusters[ii]->GetToT(),fTofClusters[ii]->GetTDC()));
tToT = (Double_t)(fTofClusters[ii]->GetToT()*AliTOFGeometry::ToTBinWidth());
tToT*=1.E-3; //ToT in ns
AliDebug(2,Form(" The ToT and Time, uncorr (ns)= %e, %e",fTofClusters[ii]->GetTDC()*AliTOFGeometry::TdcBinWidth()*1.E-3,tToT));
void AliTOFtrackerMI::MatchTracksMI(Bool_t mLastStep){
//Match ESD tracks to clusters in TOF
- const Float_t kTofOffset = 26; // time offset
- const Float_t kMinQuality = -6.; // minimal quality
- const Float_t kMaxQualityD = 1.; // max delta quality if cluster used
- const Float_t kForbiddenR = 0.1; // minimal PID according TPC
+ const Float_t kTofOffset = 26; // time offset
+ const Float_t kMinQuality = -6.; // minimal quality
+ const Float_t kMaxQualityD = 1.; // max delta quality if cluster used
+ const Float_t kForbiddenR = 0.1; // minimal PID according TPC
static const Double_t kMasses[]={
0.000511, 0.105658, 0.139570, 0.493677, 0.938272, 1.875613
Float_t scalefact=3.;
Double_t dphi=
scalefact*
- ((5*TMath::Sqrt(cov[0]) + 3.*fDy +10.*TMath::Abs(par[2]))/fR);
+ ((5*TMath::Sqrt(cov[0]) + 3.*fDy + 10.*TMath::Abs(par[2]))/fR);
Double_t dz=
scalefact*
- (5*TMath::Sqrt(cov[2]) + 3.*fDz +10.*TMath::Abs(par[3]));
+ (5*TMath::Sqrt(cov[2]) + 3.*fDz + 10.*TMath::Abs(par[3]));
Double_t phi=TMath::ATan2(par[0],x) + trackTOFin->GetAlpha();
if (phi<-TMath::Pi())phi+=2*TMath::Pi();
if (nfound>=1000) break;
index[nfound]=clind[icl];
AliTOFcluster *cluster = fClusters[clind[icl]];
- GetLinearDistances(trackTOFin,cluster, distances);
+ GetLinearDistances(trackTOFin, cluster, distances);
dist3D[nfound][0] = distances[4];
dist3D[nfound][1] = distances[1];
dist3D[nfound][2] = distances[2];
//
clusters[nfound] = cluster;
//
- //length and TOF updates
+ //length and TOF updates
trackTOFin->GetIntegratedTimes(times[nfound]);
length[nfound] = trackTOFin->GetIntegratedLength();
length[nfound]+=distances[4];