]> git.uio.no Git - u/mrichter/AliRoot.git/blame - JETAN/centua1.C
Dixed a crash when no proper files/data to display (Mikolaj)
[u/mrichter/AliRoot.git] / JETAN / centua1.C
CommitLineData
8d3b4eb1 1void centua1(char* bin)
2{
3 gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libJETAN");
4
5 char dir[100];
6 //sprintf(dir, "/data1/jgcn/partonic_events/unquenched/%s", bin);
7 sprintf(dir, "/home/guest/alice/data/cent1_nq/%s", bin);
8
9 AliJetESDReaderHeader *krh = new AliJetESDReaderHeader(); // hijing
10 //AliJetKineReaderHeader *krh = new AliJetKineReaderHeader(); // pythia
11 krh->SetComment("");
12 krh->SetDirectory(dir);
13 krh->SetPattern("miniesd"); // hijing
14 krh->SetFirstEvent(0);
15 krh->SetLastEvent(100000);
16 krh->SetPtCut(1.0);
17 //krh->SetFastSimTPC(kFALSE); // pythia
18 //krh->SetFastSimEMCAL(kFALSE); // pythia
19
20 // define reader and set its header
21 AliJetESDReader *kr = new AliJetESDReader(); // hijing
22 //AliJetKineReader *kr = new AliJetKineReader(); // pythia
23 kr->SetReaderHeader(krh);
24
25 // define jet header
26 AliUA1JetHeader *jh=new AliUA1JetHeader();
27 jh->SetComment("UA1 jet code with radius 1");
28 jh->SetMode(1);
29 jh->SetRadius(0.4);
30 jh->SetMinCellEt(0.);
31 jh->SetEtSeed(4.);
32 jh->SetLegoNbinPhi(420.);
33 jh->SetLegoNbinEta(120.);
34 jh->SetLegoEtaMin(-0.9);
35 jh->SetLegoEtaMax(+0.9);
36
37 // define jet finder. Set its header and reader
38 AliUA1JetFinder *jf = new AliUA1JetFinder();
39 jf->SetJetHeader(jh);
40 jf->SetJetReader(kr);
41 jf->SetPlotMode(kTRUE);
42 jf->SetOutputFile("jets.root");
43 // do the job
44 jf->Run();
45}
46