Adding AliTPCTracklet to the repository (M.Mager)
[u/mrichter/AliRoot.git] / TPC / AliTPCPropagateBack.C
CommitLineData
c2b0eaa4 1/****************************************************************************
2 * Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch *
3 ****************************************************************************/
4
b9de75e1 5#ifndef __CINT__
ddae8318 6 #include <Riostream.h>
b9de75e1 7 #include "AliTPCtracker.h"
8
9 #include "TFile.h"
10 #include "TStopwatch.h"
11#endif
12
13Int_t AliTPCPropagateBack() {
14 cerr<<"Propagating tracks back through the TPC...\n";
15
16 TFile *in=TFile::Open("AliTPCtracks.root");
17 if (!in->IsOpen()) {
18 cerr<<"Can't open AliTPCtracks.root !\n";
19 return 1;
20 }
21 TFile *out=TFile::Open("AliTPCBackTracks.root","new");
22 if (!out->IsOpen()) {
23 cerr<<"Delete old AliTPCBackTracks.root !\n"; return 2;
24 }
25 TFile *file=TFile::Open("AliTPCclusters.root");
26 if (!file->IsOpen()) {
27 cerr<<"Can't open AliTPCclusters.root !\n";return 3;
28 }
024a7fe9 29 AliTPCParam *param=(AliTPCParam*)file->Get("75x40_100x60_150x60");
30 if (!param) {cerr<<"Can't get TPC parameters !\n"; return 4;}
b9de75e1 31
32 TStopwatch timer;
33 AliTPCtracker *tracker=new AliTPCtracker(param);
34 Int_t rc=tracker->PropagateBack(in,out);
35 delete tracker;
36 timer.Stop(); timer.Print();
37
38 file->Close();
39 in->Close();
40 out->Close();
41
42 return rc;
43}