2D Centrality files
[u/mrichter/AliRoot.git] / JETAN / centua1.C
1 void 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