ROOT logo
void centua1(char* bin)
{
  gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libJETAN");
  
  char dir[100];
  //sprintf(dir, "/data1/jgcn/partonic_events/unquenched/%s", bin);
  sprintf(dir, "/home/guest/alice/data/cent1_nq/%s", bin);

  AliJetESDReaderHeader *krh = new AliJetESDReaderHeader(); // hijing
  //AliJetKineReaderHeader *krh = new AliJetKineReaderHeader(); // pythia
  krh->SetComment("");
  krh->SetDirectory(dir);
  krh->SetPattern("miniesd"); // hijing
  krh->SetFirstEvent(0);
  krh->SetLastEvent(100000);
  krh->SetPtCut(1.0);
  //krh->SetFastSimTPC(kFALSE); // pythia
  //krh->SetFastSimEMCAL(kFALSE); // pythia
 
  // define reader and set its header
  AliJetESDReader *kr = new AliJetESDReader(); // hijing
  //AliJetKineReader *kr = new AliJetKineReader(); // pythia
  kr->SetReaderHeader(krh);
 
  // define jet header
  AliUA1JetHeader *jh=new AliUA1JetHeader();
  jh->SetComment("UA1 jet code with radius 1");
  jh->SetMode(1);
  jh->SetRadius(0.4);
  jh->SetMinCellEt(0.);
  jh->SetEtSeed(4.);
  jh->SetLegoNbinPhi(420.);
  jh->SetLegoNbinEta(120.);
  jh->SetLegoEtaMin(-0.9);
  jh->SetLegoEtaMax(+0.9);
    
  // define jet finder. Set its header and reader
  AliUA1JetFinder *jf = new AliUA1JetFinder();
  jf->SetJetHeader(jh);
  jf->SetJetReader(kr);
  jf->SetPlotMode(kTRUE);
  jf->SetOutputFile("jets.root");
  // do the job
  jf->Run();
}
 
 centua1.C:1
 centua1.C:2
 centua1.C:3
 centua1.C:4
 centua1.C:5
 centua1.C:6
 centua1.C:7
 centua1.C:8
 centua1.C:9
 centua1.C:10
 centua1.C:11
 centua1.C:12
 centua1.C:13
 centua1.C:14
 centua1.C:15
 centua1.C:16
 centua1.C:17
 centua1.C:18
 centua1.C:19
 centua1.C:20
 centua1.C:21
 centua1.C:22
 centua1.C:23
 centua1.C:24
 centua1.C:25
 centua1.C:26
 centua1.C:27
 centua1.C:28
 centua1.C:29
 centua1.C:30
 centua1.C:31
 centua1.C:32
 centua1.C:33
 centua1.C:34
 centua1.C:35
 centua1.C:36
 centua1.C:37
 centua1.C:38
 centua1.C:39
 centua1.C:40
 centua1.C:41
 centua1.C:42
 centua1.C:43
 centua1.C:44
 centua1.C:45
 centua1.C:46
 centua1.C:47