1 void centua1(char* bin)
3 gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libJETAN");
6 //sprintf(dir, "/data1/jgcn/partonic_events/unquenched/%s", bin);
7 sprintf(dir, "/home/guest/alice/data/cent1_nq/%s", bin);
9 AliJetESDReaderHeader *krh = new AliJetESDReaderHeader(); // hijing
10 //AliJetKineReaderHeader *krh = new AliJetKineReaderHeader(); // pythia
12 krh->SetDirectory(dir);
13 krh->SetPattern("miniesd"); // hijing
14 krh->SetFirstEvent(0);
15 krh->SetLastEvent(100000);
17 //krh->SetFastSimTPC(kFALSE); // pythia
18 //krh->SetFastSimEMCAL(kFALSE); // pythia
20 // define reader and set its header
21 AliJetESDReader *kr = new AliJetESDReader(); // hijing
22 //AliJetKineReader *kr = new AliJetKineReader(); // pythia
23 kr->SetReaderHeader(krh);
26 AliUA1JetHeader *jh=new AliUA1JetHeader();
27 jh->SetComment("UA1 jet code with radius 1");
32 jh->SetLegoNbinPhi(420.);
33 jh->SetLegoNbinEta(120.);
34 jh->SetLegoEtaMin(-0.9);
35 jh->SetLegoEtaMax(+0.9);
37 // define jet finder. Set its header and reader
38 AliUA1JetFinder *jf = new AliUA1JetFinder();
41 jf->SetPlotMode(kTRUE);
42 jf->SetOutputFile("jets.root");