Class moved to STEERBase.
[u/mrichter/AliRoot.git] / JETAN / testua1.C
CommitLineData
99e5fe42 1void 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}