#include "AliTPCtracker.h"
#endif
-struct GoodTrack {
+struct GoodTrackTPC {
Int_t fEventN; //event number
Int_t lab;
Int_t code;
Float_t px,py,pz;
Float_t x,y,z;
};
-Int_t good_tracks(GoodTrack *gt, Int_t max, Int_t eventn=1);
Int_t AliTPCComparison(Int_t eventn=1) {
+ Int_t good_tracks_tpc(GoodTrackTPC *gt, Int_t max, Int_t eventn=1);
+
cerr<<"Doing comparison...\n";
Int_t i;
gBenchmark->Start("AliTPCComparison");
// Load clusters
eventptr[0]=0;
eventptr[1]=0;
+ char tname[100]; //Y.B.
for (Int_t event=0;event<eventn; event++){
cf->cd();
tracker = new AliTPCtracker(digp,event);
tracker->LoadInnerSectors();
tracker->LoadOuterSectors();
- char tname[100];
+ //Y.B. char tname[100];
if (eventn==-1) {
sprintf(tname,"TreeT_TPC");
}
}
eventptr[event+1] = nentr; //store end of the event
delete tracker;
+ delete tracktree; //Thanks to Mariana Bondila
}
tf->Close();
+ delete digp; //Thanks to Mariana Bondila
cf->Close();
- //MI change for a moment
-
/////////////////////////////////////////////////////////////////////////
- // GoodTrack gt[15000];
- GoodTrack gt[45000];
+ const Int_t MAX=45000;
+ GoodTrackTPC gt[MAX];
Int_t ngood=0;
ifstream in("good_tracks_tpc");
gt[ngood].x >>gt[ngood].y >>gt[ngood].z) {
ngood++;
cerr<<ngood<<'\r';
- if (ngood==45000) {
+ if (ngood==MAX) {
cerr<<"Too many good tracks !\n";
break;
}
if (!in.eof()) cerr<<"Read error (good_tracks_tpc) !\n";
} else {
cerr<<"Marking good tracks (this will take a while)...\n";
- ngood=good_tracks(gt,45000,eventn); //mi change
+ ngood=good_tracks_tpc(gt,MAX,eventn); //mi change
ofstream out("good_tracks_tpc");
if (out) {
for (Int_t ngd=0; ngd<ngood; ngd++)
} else cerr<<"Can not open file (good_tracks_tpc) !\n";
out.close();
- /*
cerr<<"Preparing tracks for matching with the ITS...\n";
tarray.Sort();
tf=TFile::Open("AliTPCtracks.root","recreate");
- tracktree=new TTree("TPCf","Tree with TPC tracks");
+ tracktree=new TTree(tname,"Tree with \"forward\" TPC tracks");
tracktree->Branch("tracks","AliTPCtrack",&iotrack,32000,0);
for (i=0; i<nentr; i++) {
iotrack=(AliTPCtrack*)tarray.UncheckedAt(i);
}
tracktree->Write();
tf->Close();
- */
}
cerr<<"Number of good tracks : "<<ngood<<endl;
}
-Int_t good_tracks(GoodTrack *gt, Int_t max, Int_t eventn) {
+Int_t good_tracks_tpc(GoodTrackTPC *gt, Int_t max, Int_t eventn) {
//eventn - number of events in file
- TFile *file=TFile::Open("rfio:galice.root");
+ TFile *file=TFile::Open("galice.root");
if (!file->IsOpen()) {cerr<<"Can't open galice.root !\n"; exit(4);}
if (!(gAlice=(AliRun*)file->Get("gAlice"))) {
AliTPCClustersArray *ca=new AliTPCClustersArray;
ca->Setup(digp);
ca->SetClusterType("AliTPCcluster");
- ca->ConnectTree("Segment Tree");
+ //ca->ConnectTree("Segment Tree");
+ ca->ConnectTree("TreeC_TPC_0");
Int_t nrows=Int_t(ca->GetTree()->GetEntries());
for (Int_t n=0; n<nrows; n++) {
AliSegmentID *s=ca->LoadEntry(n);
}
}
delete[] count;
+ delete TD; //Thanks to Mariana Bondila
break;
default:
cerr<<"Invalid TPC version !\n";
delete[] good;
}
delete gAlice; gAlice=0;
-
file->Close();
gBenchmark->Stop("AliTPCComparison");
gBenchmark->Show("AliTPCComparison");
/*
$Log$
+Revision 1.13 2001/07/23 12:01:30 hristov
+Initialisation added
+
Revision 1.12 2001/07/20 14:32:44 kowal2
Processing of many events possible now
#include "AliTPCParam.h"
#include "AliTPCClustersRow.h"
-
-AliTPCtracker::AliTPCtracker(const AliTPCParam *par):
-fkNIS(par->GetNInnerSector()/2),
-fkNOS(par->GetNOuterSector()/2)
-{;
-//MI change only provisore - need change in the ITS code which depend on it
-}
-
-
//_____________________________________________________________________________
AliTPCtracker::AliTPCtracker(const AliTPCParam *par, Int_t eventn):
fkNIS(par->GetNInnerSector()/2),
}
delete seed;
+
+ delete seedTree; //Thanks to Mariana Bondila
+
savedir->cd();
return 0;
}
in->cd();
- TTree *bckTree=(TTree*)in->Get("ITSb");
+ TTree *bckTree=(TTree*)in->Get("TreeT_ITSb_0");
if (!bckTree) {
cerr<<"AliTPCtracker::PropagateBack() ";
cerr<<"can't get a tree with back propagated ITS tracks !\n";
AliTPCtrack *bckTrack=new AliTPCtrack;
bckTree->SetBranchAddress("tracks",&bckTrack);
- TTree *tpcTree=(TTree*)in->Get("TPCf");
+ TTree *tpcTree=(TTree*)in->Get("TreeT_TPC_0");
if (!tpcTree) {
cerr<<"AliTPCtracker::PropagateBack() ";
cerr<<"can't get a tree with TPC tracks !\n";
}
out->cd();
- TTree backTree("TPCb","Tree with back propagated TPC tracks");
+ TTree backTree("TreeT_TPCb_0","Tree with back propagated TPC tracks");
AliTPCtrack *otrack=0;
backTree.Branch("tracks","AliTPCtrack",&otrack,32000,0);
delete bckTrack;
delete tpcTrack;
+ delete bckTree; //Thanks to Mariana Bondila
+ delete tpcTree; //Thanks to Mariana Bondila
+
return 0;
}