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