//--------------------------------------------------------------------//
// //
-// AliTOFtrackerV1 Class //
+// AliTOFtrackerV1 Class //
// Task: Perform association of the ESD tracks to TOF Clusters //
// and Update ESD track with associated TOF Cluster parameters //
// //
#include <TFile.h>
#include <TH1F.h>
#include <TH2F.h>
+#include <TSeqCollection.h>
-#include "AliAlignObj.h"
#include "AliESDtrack.h"
#include "AliESDEvent.h"
#include "AliLog.h"
#include "AliTrackPointArray.h"
+#include "AliGeomManager.h"
+#include "AliCDBManager.h"
-#include "AliTOFcalib.h"
#include "AliTOFRecoParam.h"
+#include "AliTOFReconstructor.h"
#include "AliTOFcluster.h"
-#include "AliTOFGeometryV5.h"
+#include "AliTOFGeometry.h"
#include "AliTOFtrackerV1.h"
#include "AliTOFtrack.h"
+#include "AliTOFpidESD.h"
extern TROOT *gROOT;
//_____________________________________________________________________________
AliTOFtrackerV1::AliTOFtrackerV1():
fRecoParam(0x0),
- fGeom(0x0),
fPid(0x0),
fN(0),
fNseeds(0),
fnbadmatch(0),
fnunmatch(0),
fnmatch(0),
- fTracks(0x0),
- fSeeds(0x0),
+ fTracks(new TClonesArray("AliTOFtrack")),
+ fSeeds(new TClonesArray("AliESDtrack")),
fHDigClusMap(0x0),
fHDigNClus(0x0),
fHDigClusTime(0x0),
fHRecSigZVsPWin(0x0)
{
//AliTOFtrackerV1 main Ctor
-
- // Getting the geometry
- fGeom=new AliTOFGeometryV5();
- // Read the reconstruction parameters from the OCDB
- AliTOFcalib *calib = new AliTOFcalib(fGeom);
- fRecoParam = (AliTOFRecoParam*)calib->ReadRecParFromCDB("TOF/Calib",-1);
- if(fRecoParam->GetApplyPbPbCuts())fRecoParam=fRecoParam->GetPbPbparam();
- delete calib;
- Double_t parPID[2];
- parPID[0]=fRecoParam->GetTimeResolution();
- parPID[1]=fRecoParam->GetTimeNSigma();
- fPid=new AliTOFpidESD(parPID);
- InitCheckHists();
-
-}
-//_____________________________________________________________________________
-AliTOFtrackerV1::AliTOFtrackerV1(const AliTOFtrackerV1 &t):
- AliTracker(),
- fRecoParam(0x0),
- fGeom(0x0),
- fPid(0x0),
- fN(0),
- fNseeds(0),
- fNseedsTOF(0),
- fngoodmatch(0),
- fnbadmatch(0),
- fnunmatch(0),
- fnmatch(0),
- fTracks(0x0),
- fSeeds(0x0),
- fHDigClusMap(0x0),
- fHDigNClus(0x0),
- fHDigClusTime(0x0),
- fHDigClusToT(0x0),
- fHRecNClus(0x0),
- fHRecChi2(0x0),
- fHRecDistZ(0x0),
- fHRecSigYVsP(0x0),
- fHRecSigZVsP(0x0),
- fHRecSigYVsPWin(0x0),
- fHRecSigZVsPWin(0x0)
- {
- //AliTOFtrackerV1 copy Ctor
-
- fNseeds=t.fNseeds;
- fNseeds=t.fNseeds;
- fNseedsTOF=t.fNseedsTOF;
- fngoodmatch=t.fngoodmatch;
- fnbadmatch=t.fnbadmatch;
- fnunmatch=t.fnunmatch;
- fnmatch=t.fnmatch;
- fRecoParam=t.fRecoParam;
- fGeom=t.fGeom;
- fPid=t.fPid;
- fSeeds=t.fSeeds;
- fTracks=t.fTracks;
- fN=t.fN;
-}
-//_____________________________________________________________________________
-AliTOFtrackerV1& AliTOFtrackerV1::operator=(const AliTOFtrackerV1 &t)
-{
- //AliTOFtrackerV1 assignment operator
-
- this->fNseeds=t.fNseeds;
- this->fNseedsTOF=t.fNseedsTOF;
- this->fngoodmatch=t.fngoodmatch;
- this->fnbadmatch=t.fnbadmatch;
- this->fnunmatch=t.fnunmatch;
- this->fnmatch=t.fnmatch;
- this->fRecoParam = t.fRecoParam;
- this->fGeom = t.fGeom;
- this->fPid = t.fPid;
- this->fSeeds=t.fSeeds;
- this->fTracks=t.fTracks;
- this->fN=t.fN;
- return *this;
+ InitCheckHists();
}
//_____________________________________________________________________________
SaveCheckHists();
- delete fRecoParam;
- delete fGeom;
+ if(!(AliCDBManager::Instance()->GetCacheFlag())){
+ delete fRecoParam;
+ }
delete fPid;
delete fHDigClusMap;
delete fHDigNClus;
delete fHRecSigZVsP;
delete fHRecSigYVsPWin;
delete fHRecSigZVsPWin;
+ if (fTracks){
+ fTracks->Delete();
+ delete fTracks;
+ fTracks=0x0;
+ }
+ if (fSeeds){
+ fSeeds->Delete();
+ delete fSeeds;
+ fSeeds=0x0;
+ }
}
//_____________________________________________________________________________
Int_t AliTOFtrackerV1::PropagateBack(AliESDEvent* event) {
// Gets seeds from ESD event and Match with TOF Clusters
//
+ // initialize RecoParam for current event
+
+ AliInfo("Initializing params for TOF... ");
+
+ fRecoParam = AliTOFReconstructor::GetRecoParam(); // instantiate reco param from STEER...
+
+ if (fRecoParam == 0x0) {
+ AliFatal("No Reco Param found for TOF!!!");
+ }
+ //fRecoParam->Dump();
+ //if(fRecoParam->GetApplyPbPbCuts())fRecoParam=fRecoParam->GetPbPbparam();
+ //fRecoParam->PrintParameters();
+
+ Double_t parPID[2];
+ parPID[0]=fRecoParam->GetTimeResolution();
+ parPID[1]=fRecoParam->GetTimeNSigma();
+ fPid=new AliTOFpidESD(parPID);
//Initialise some counters
Int_t ntrk=event->GetNumberOfTracks();
fNseeds = ntrk;
- fSeeds= new TClonesArray("AliESDtrack",ntrk);
TClonesArray &aESDTrack = *fSeeds;
//Make TOF PID
fPid->MakePID(event,timeZero);
- if (fSeeds) {
- fSeeds->Delete();
- delete fSeeds;
- fSeeds = 0x0;
- }
- if (fTracks) {
- fTracks->Delete();
- delete fTracks;
- fTracks = 0x0;
- }
+ fSeeds->Clear();
+ fTracks->Clear();
return 0;
}
Int_t seedsTOF1=0;
Int_t seedsTOF2=0;
- fTracks= new TClonesArray("AliTOFtrack");
TClonesArray &aTOFTrack = *fTracks;
for (Int_t i=0; i<fNseeds; i++) {
AliESDtrack *t =(AliESDtrack*)fSeeds->UncheckedAt(i);
if ((t->GetStatus()&AliESDtrack::kTPCout)==0)continue;
- // TRD 'good' tracks, already propagated at 371 cm
+ // TRD 'good' tracks, already propagated at 372 cm
AliTOFtrack *track = new AliTOFtrack(*t); // New
Double_t x = track->GetX(); //New
if (((t->GetStatus()&AliESDtrack::kTRDout)!=0 ) &&
- ( x >= fGeom->RinTOF()) ){
+ ( x >= AliTOFGeometry::RinTOF()) ){
track->SetSeedIndex(i);
t->UpdateTrackParams(track,AliESDtrack::kTOFout);
new(aTOFTrack[fNseedsTOF]) AliTOFtrack(*track);
// Parameters regulating the reconstruction
- Float_t dY=fGeom->XPad();
- Float_t dZ=fGeom->ZPad();
+ Float_t dY=AliTOFGeometry::XPad();
+ Float_t dZ=AliTOFGeometry::ZPad();
- const Int_t ncmax = 100;
+ const Int_t kncmax = 100;
Float_t sensRadius = fRecoParam->GetSensRadius();
Float_t scaleFact = fRecoParam->GetWindowScaleFact();
Float_t dyMax=fRecoParam->GetWindowSizeMaxY();
//The matching loop
- AliTOFcalib *calib = new AliTOFcalib(fGeom);
-
for (Int_t iseed=0; iseed<fNseedsTOF; iseed++) {
AliTOFtrack *track =(AliTOFtrack*)fTracks->UncheckedAt(iseed);
Double_t cov[15]; trackTOFin->GetExternalCovariance(cov);
Double_t z = par[1];
- Double_t dz = scaleFact*3.*TMath::Sqrt(cov[2]+dZ*dZ/12);
+ Double_t dz = scaleFact*3.*TMath::Sqrt(cov[2]+dZ*dZ/12.);
Double_t dphi = scaleFact*3.*TMath::Sqrt(cov[0]+dY*dY/12.)/sensRadius;
Double_t phi=TMath::ATan2(par[0],x) + trackTOFin->GetAlpha();
// find the clusters inside the selected window
Int_t nc=0;
- AliTOFcluster *clusters[ncmax]; // pointers to the clusters in the window
- Int_t index[ncmax];//to keep track of the cluster index
+ AliTOFcluster *clusters[kncmax]; // pointers to the clusters in the window
+ Int_t index[kncmax];//to keep track of the cluster index
for (Int_t k=FindClusterIndex(z-dz); k<fN; k++) {
AliTOFcluster *c=fClusters[k];
- if(nc>ncmax)break;
+ // if(nc>kncmax)break; /* R+ fix (buffer overflow) */
+ if(nc>=kncmax)break; /* R+ fix (buffer overflow protection) */
if(c->GetZ() > z+dz) break;
if(c->IsUsed()) continue;
- if(!c->GetStatus()) continue; // skip bad channels as declared in OCDB
+ if(!c->GetStatus()) {
+ AliDebug(1,"Cluster in channel declared bad!");
+ continue; // skip bad channels as declared in OCDB
+ }
Float_t xyz[3]; c->GetGlobalXYZ(xyz);
Double_t clPhi=TMath::ATan2(xyz[1],xyz[0]);
Double_t dph=TMath::Abs(clPhi-phi);
nc++;
}
- //start propagation: go to the average TOF pad middle plane at ~378.5 cm
+ //start propagation: go to the average TOF pad middle plane at ~379.5 cm
Float_t xTOF = sensRadius;
- Double_t ymax = xTOF*TMath::Tan(0.5*fGeom->GetAlpha());
+ Double_t ymax = xTOF*TMath::Tan(0.5*AliTOFGeometry::GetAlpha());
Bool_t skip = kFALSE;
Double_t ysect = trackTOFin->GetYat(xTOF,skip);
if (skip) break;
if (ysect > ymax) {
- if (!trackTOFin->Rotate(fGeom->GetAlpha())) {
+ if (!trackTOFin->Rotate(AliTOFGeometry::GetAlpha())) {
break;
}
} else if (ysect <-ymax) {
- if (!trackTOFin->Rotate(fGeom->GetAlpha())) {
+ if (!trackTOFin->Rotate(-AliTOFGeometry::GetAlpha())) {
break;
}
}
AliTOFcluster *bestCluster=0;
Double_t bestChi2=maxChi2;
Int_t idclus=-1;
- for (Int_t i=0; i<nc; i++){
+ // for (Int_t i=0; i<nc; i++){ /* R+ fix (unsafe) */
+ for (Int_t i=0; i<nc && i<kncmax; i++){ /* R+ fix (buffer overflow protection) */
AliTOFcluster *c=clusters[i]; // one of the preselected clusters
Double_t chi2=trackTOFin->GetPredictedChi2((AliCluster3D*)c);
if (chi2 >= bestChi2) continue;
}
fnmatch++;
+
+ AliDebug(2, Form("%7i %7i %10i %10i %10i %10i %7i",
+ iseed,
+ fnmatch-1,
+ TMath::Abs(trackTOFin->GetLabel()),
+ bestCluster->GetLabel(0),
+ bestCluster->GetLabel(1),
+ bestCluster->GetLabel(2),
+ idclus)); // AdC
+
bestCluster->Use();
if (
(bestCluster->GetLabel(0)==TMath::Abs(trackTOFin->GetLabel()))
//Propagate the track to the best matched cluster
trackTOFin->PropagateTo(bestCluster);
+ // Fill the track residual histograms.
+ FillResiduals(trackTOFin,bestCluster,kFALSE);
//now take the local distance in Z from the pad center for time walk correction
- Float_t tiltangle = fGeom->GetAngles(bestCluster->GetDetInd(1),bestCluster->GetDetInd(2))*TMath::DegToRad();
+ Float_t tiltangle = AliTOFGeometry::GetAngles(bestCluster->GetDetInd(1),bestCluster->GetDetInd(2))*TMath::DegToRad();
Double_t dzTW=trackTOFin->GetZ()-bestCluster->GetZ(); // in cm
dzTW/=TMath::Cos(tiltangle);
delete trackTOFin;
// Store quantities to be used in the TOF Calibration
- Float_t tToT=fGeom->ToTBinWidth()*bestCluster->GetToT()*1E-3; // in ns
+ Float_t tToT=AliTOFGeometry::ToTBinWidth()*bestCluster->GetToT()*1E-3; // in ns
t->SetTOFsignalToT(tToT);
- Float_t rawTime=fGeom->TdcBinWidth()*bestCluster->GetTDCRAW()+32; // RAW time,in ps
+ Float_t rawTime=AliTOFGeometry::TdcBinWidth()*bestCluster->GetTDCRAW()+32; // RAW time,in ps
t->SetTOFsignalRaw(rawTime);
t->SetTOFsignalDz(dzTW);
AliDebug(2,Form(" Setting TOF raw time: %f z distance: %f time: %f = ",rawTime,dzTW));
ind[2]=bestCluster->GetDetInd(2);
ind[3]=bestCluster->GetDetInd(3);
ind[4]=bestCluster->GetDetInd(4);
- Int_t calindex = calib->GetIndex(ind);
+ Int_t calindex = AliTOFGeometry::GetIndex(ind);
t->SetTOFCalChannel(calindex);
// keep track of the track labels in the matched cluster
tlab[1]=bestCluster->GetLabel(1);
tlab[2]=bestCluster->GetLabel(2);
AliDebug(2,Form(" tdc time of the matched track %i = ",bestCluster->GetTDC()));
- Double_t tof=fGeom->TdcBinWidth()*bestCluster->GetTDC()+32; // in ps
+ Double_t tof=AliTOFGeometry::TdcBinWidth()*bestCluster->GetTDC()+32; // in ps
AliDebug(2,Form(" tof time of the matched track: %f = ",tof));
Double_t tofcorr=tof;
if(timeWalkCorr)tofcorr=CorrectTimeWalk(dzTW,tof);
// no longer there - all info is in the ESDs now
}
- delete calib;
}
//_________________________________________________________________________
//This function loads the TOF clusters
//--------------------------------------------------------------------
- Int_t npadX = fGeom->NpadX();
- Int_t npadZ = fGeom->NpadZ();
- Int_t nStripA = fGeom->NStripA();
- Int_t nStripB = fGeom->NStripB();
- Int_t nStripC = fGeom->NStripC();
+ Int_t npadX = AliTOFGeometry::NpadX();
+ Int_t npadZ = AliTOFGeometry::NpadZ();
+ Int_t nStripA = AliTOFGeometry::NStripA();
+ Int_t nStripB = AliTOFGeometry::NStripB();
+ Int_t nStripC = AliTOFGeometry::NStripC();
TBranch *branch=cTree->GetBranch("TOF");
if (!branch) {
return 1;
}
- TClonesArray dummy("AliTOFcluster",10000), *clusters=&dummy;
+ static TClonesArray dummy("AliTOFcluster",10000);
+ dummy.Clear();
+ TClonesArray *clusters=&dummy;
branch->SetAddress(&clusters);
cTree->GetEvent(0);
Int_t nc=clusters->GetEntriesFast();
fHDigNClus->Fill(nc);
+ AliInfo(Form("Number of clusters: %d",nc));
+
for (Int_t i=0; i<nc; i++) {
AliTOFcluster *c=(AliTOFcluster*)clusters->UncheckedAt(i);
fClusters[i]=new AliTOFcluster(*c); fN++;
Int_t ipadX = c->GetDetInd(4);
Int_t ipadZ = c->GetDetInd(3);
- Float_t time =(fGeom->TdcBinWidth()*c->GetTDC())*1E-3; // in ns
- Float_t tot = (fGeom->TdcBinWidth()*c->GetToT())*1E-3;//in ns
+ Float_t time =(AliTOFGeometry::TdcBinWidth()*c->GetTDC())*1E-3; // in ns
+ Float_t tot = (AliTOFGeometry::TdcBinWidth()*c->GetToT())*1E-3;//in ns
Int_t stripOffset = 0;
switch (iplate) {
Float_t phi=TMath::ATan2(xyz[1],xyz[0]);
Float_t phiangle = (Int_t(phi*TMath::RadToDeg()/20.)+0.5)*20.*TMath::DegToRad();
Float_t sinphi = TMath::Sin(phiangle), cosphi = TMath::Cos(phiangle);
- Float_t tiltangle = fGeom->GetAngles(cl->GetDetInd(1),cl->GetDetInd(2))*TMath::DegToRad();
+ Float_t tiltangle = AliTOFGeometry::GetAngles(cl->GetDetInd(1),cl->GetDetInd(2))*TMath::DegToRad();
Float_t sinth = TMath::Sin(tiltangle), costh = TMath::Cos(tiltangle);
- Float_t sigmay2 = fGeom->XPad()*fGeom->XPad()/12.;
- Float_t sigmaz2 = fGeom->ZPad()*fGeom->ZPad()/12.;
+ Float_t sigmay2 = AliTOFGeometry::XPad()*AliTOFGeometry::XPad()/12.;
+ Float_t sigmaz2 = AliTOFGeometry::ZPad()*AliTOFGeometry::ZPad()/12.;
Float_t cov[6];
cov[0] = sinphi*sinphi*sigmay2 + cosphi*cosphi*sinth*sinth*sigmaz2;
cov[1] = -sinphi*cosphi*sigmay2 + sinphi*cosphi*sinth*sinth*sigmaz2;
p.SetXYZ(xyz[0],xyz[1],xyz[2],cov);
// Detector numbering scheme
- Int_t nSector = fGeom->NSectors();
- Int_t nPlate = fGeom->NPlates();
- Int_t nStripA = fGeom->NStripA();
- Int_t nStripB = fGeom->NStripB();
- Int_t nStripC = fGeom->NStripC();
+ Int_t nSector = AliTOFGeometry::NSectors();
+ Int_t nPlate = AliTOFGeometry::NPlates();
+ Int_t nStripA = AliTOFGeometry::NStripA();
+ Int_t nStripB = AliTOFGeometry::NStripB();
+ Int_t nStripC = AliTOFGeometry::NStripC();
Int_t isector = cl->GetDetInd(0);
if (isector >= nSector)
//Init histos for Digits/Reco QA and Calibration
+ TDirectory *dir = gDirectory;
+ TFile *logFileTOF = 0;
+
+ TSeqCollection *list = gROOT->GetListOfFiles();
+ int n = list->GetEntries();
+ Bool_t isThere=kFALSE;
+ for(int i=0; i<n; i++) {
+ logFileTOF = (TFile*)list->At(i);
+ if (strstr(logFileTOF->GetName(), "TOFQA.root")){
+ isThere=kTRUE;
+ break;
+ }
+ }
+
+ if(!isThere)logFileTOF = new TFile( "TOFQA.root","RECREATE");
+ logFileTOF->cd();
+
//Digits "QA"
fHDigClusMap = new TH2F("TOFDig_ClusMap", "",182,0.5,182.5,864, 0.5,864.5);
fHDigNClus = new TH1F("TOFDig_NClus", "",200,0.5,200.5);
fHRecSigZVsP=new TH2F("TOFDig_SigZVsP", "",40,0.,4.,100, 0.,5.);
fHRecSigYVsPWin=new TH2F("TOFDig_SigYVsPWin", "",40,0.,4.,100, 0.,50.);
fHRecSigZVsPWin=new TH2F("TOFDig_SigZVsPWin", "",40,0.,4.,100, 0.,50.);
+
+ dir->cd();
+
}
//_________________________________________________________________________
}
}
+ if(!isThere) {
+ AliError(Form("File TOFQA.root not found!! not wring histograms...."));
+ return;
+ }
logFile->cd();
fHDigClusMap->Write(fHDigClusMap->GetName(), TObject::kOverwrite);
fHDigNClus->Write(fHDigNClus->GetName(), TObject::kOverwrite);
fHRecSigZVsPWin->Write(fHRecSigZVsPWin->GetName(), TObject::kOverwrite);
logFile->Flush();
- if(!isThere)logFileTOF = new TFile( "TOFQA.root","RECREATE");
- logFileTOF->cd();
- fHDigClusMap->Write(fHDigClusMap->GetName(), TObject::kOverwrite);
- fHDigNClus->Write(fHDigNClus->GetName(), TObject::kOverwrite);
- fHDigClusTime->Write(fHDigClusTime->GetName(), TObject::kOverwrite);
- fHDigClusToT->Write(fHDigClusToT->GetName(), TObject::kOverwrite);
- fHRecNClus->Write(fHRecNClus->GetName(), TObject::kOverwrite);
- fHRecChi2->Write(fHRecChi2->GetName(), TObject::kOverwrite);
- fHRecDistZ->Write(fHRecDistZ->GetName(), TObject::kOverwrite);
- fHRecSigYVsP->Write(fHRecSigYVsP->GetName(), TObject::kOverwrite);
- fHRecSigZVsP->Write(fHRecSigZVsP->GetName(), TObject::kOverwrite);
- fHRecSigYVsPWin->Write(fHRecSigYVsPWin->GetName(), TObject::kOverwrite);
- fHRecSigZVsPWin->Write(fHRecSigZVsPWin->GetName(), TObject::kOverwrite);
- logFileTOF->Flush();
-
dir->cd();
+
}
//_________________________________________________________________________
Float_t AliTOFtrackerV1::CorrectTimeWalk( Float_t dist, Float_t tof) {
//dummy, for the moment
Float_t tofcorr=0.;
- if(dist<fGeom->ZPad()*0.5){
+ if(dist<AliTOFGeometry::ZPad()*0.5){
tofcorr=tof;
//place here the actual correction
}else{
}
return 0.;
}
+//_________________________________________________________________________
+
+void AliTOFtrackerV1::FillClusterArray(TObjArray* arr) const
+{
+ //
+ // Returns the TOF cluster array
+ //
+
+ if (fN==0)
+ arr = 0x0;
+ else
+ for (Int_t i=0; i<fN; ++i) arr->Add(fClusters[i]);
+
+}