JETAN/testua1.C

   1 void testua1()
   2 {
   3   gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libJETAN");
   4
   5   // define reader header
   6   AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader();
   7   jrh->SetComment("testing");
   8   //jrh->SetDirectory("rfio:///castor/cern.ch/user/p/phristov/cent1_nq/104-125GeV");
   9   jrh->SetDirectory("/home/morsch/jets_104_125");
  10   jrh->SetPattern("00165");
  11   jrh->SetFirstEvent(0);
  12   jrh->SetLastEvent(1000);
  13   jrh->SetPtCut(2.);
  14   jrh->SetReadSignalOnly(kFALSE);
  15
  16   AliJetKineReaderHeader *krh = new AliJetKineReaderHeader();
  17   krh->SetComment("testing");
  18   krh->SetDirectory("/home/morsch/AliRoot/newio/FASTSIM/ctest/100/uq");
  19   krh->SetFirstEvent(0);
  20   krh->SetLastEvent(1000);
  21   krh->SetPtCut(2.);
  22   krh->SetFastSimTPC(kTRUE);
  23
  24
  25   // define reader and set its header
  26   AliJetESDReader *jr = new AliJetESDReader();
  27   jr->SetReaderHeader(jrh);
  28
  29   // define reader and set its header
  30   AliJetKineReader *kr = new AliJetKineReader();
  31   kr->SetReaderHeader(krh);
  32
  33   // define jet header
  34   AliUA1JetHeader *jh=new AliUA1JetHeader();
  35   jh->SetComment("UA1 jet code with default parameters");
  36   jh->SetMode(0);
  37   jh->SetRadius(0.4);
  38   jh->SetMinCellEt(0.);
  39   jh->SetEtSeed(4.);
  40   jh->SetNbinPhi(420.);
  41   jh->SetNbinEta(120.);
  42   jh->SetEtaMin(-0.9);
  43   jh->SetEtaMax(+0.9);
  44
  45
  46   // define jet finder. Set its header and reader
  47   AliUA1JetFinder *jf = new AliUA1JetFinder();
  48   jf->SetJetHeader(jh);
  49   jf->SetJetReader(jr);
  50   jf->SetPlotMode(kTRUE);
  51   jf->SetOutputFile("jets.root");
  52   // do the job
  53   jf->Run();
  54 }