]>
Commit | Line | Data |
---|---|---|
99e5fe42 | 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 | } |