#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 "AliTOFGeometry.h"
#include "AliTOFtrackerV1.h"
#include "AliTOFtrack.h"
+#include "AliTOFpidESD.h"
extern TROOT *gROOT;
fnbadmatch(0),
fnunmatch(0),
fnmatch(0),
- fTracks(0x0),
- fSeeds(0x0),
+ fTracks(new TClonesArray("AliTOFtrack")),
+ fSeeds(new TClonesArray("AliESDtrack")),
fHDigClusMap(0x0),
fHDigNClus(0x0),
fHDigClusTime(0x0),
fnbadmatch(0),
fnunmatch(0),
fnmatch(0),
- fTracks(0x0),
- fSeeds(0x0),
+ fTracks(new TClonesArray("AliTOFtrack")),
+ fSeeds(new TClonesArray("AliESDtrack")),
fHDigClusMap(0x0),
fHDigNClus(0x0),
fHDigClusTime(0x0),
SaveCheckHists();
- delete fRecoParam;
+ 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) {
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
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();
// 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;
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*AliTOFGeometry::GetAlpha());
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++;
}
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]);
+
+}