Class moved to STEERBase.
[u/mrichter/AliRoot.git] / JETAN / testua1.C
1 void testua1()
2 {
3   gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libJETAN");
4
5   // define reader header
6   AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader();
7   jrh->SetComment("testing");
8   //jrh->SetDirectory("rfio:///castor/cern.ch/user/p/phristov/cent1_nq/104-125GeV");
9   jrh->SetDirectory("/home/morsch/jets_104_125");
10   jrh->SetPattern("00165");
11   jrh->SetFirstEvent(0);
12   jrh->SetLastEvent(1000);
13   jrh->SetPtCut(2.);
14   jrh->SetReadSignalOnly(kFALSE);
15
16   AliJetKineReaderHeader *krh = new AliJetKineReaderHeader();
17   krh->SetComment("testing");
18   krh->SetDirectory("/home/morsch/AliRoot/newio/FASTSIM/ctest/100/uq");
19   krh->SetFirstEvent(0);
20   krh->SetLastEvent(1000);
21   krh->SetPtCut(2.);
22   krh->SetFastSimTPC(kTRUE);
23   
24
25   // define reader and set its header
26   AliJetESDReader *jr = new AliJetESDReader();
27   jr->SetReaderHeader(jrh);
28
29   // define reader and set its header
30   AliJetKineReader *kr = new AliJetKineReader();
31   kr->SetReaderHeader(krh);
32
33   // define jet header
34   AliUA1JetHeader *jh=new AliUA1JetHeader();
35   jh->SetComment("UA1 jet code with default parameters");
36   jh->SetMode(0);
37   jh->SetRadius(0.4);
38   jh->SetMinCellEt(0.);
39   jh->SetEtSeed(4.);
40   jh->SetNbinPhi(420.);
41   jh->SetNbinEta(120.);
42   jh->SetEtaMin(-0.9);
43   jh->SetEtaMax(+0.9);  
44   
45
46   // define jet finder. Set its header and reader
47   AliUA1JetFinder *jf = new AliUA1JetFinder();
48   jf->SetJetHeader(jh);
49   jf->SetJetReader(jr);
50   jf->SetPlotMode(kTRUE);
51   jf->SetOutputFile("jets.root");
52   // do the job
53   jf->Run();
54 }