Correction.

[u/mrichter/AliRoot.git] / JETAN / testua1.C

void testua1()
{
  gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libJETAN");

  // define reader header
  AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader();
  jrh->SetComment("testing");
  //jrh->SetDirectory("rfio:///castor/cern.ch/user/p/phristov/cent1_nq/104-125GeV");
  jrh->SetDirectory("/home/morsch/jets_104_125");
  jrh->SetPattern("00165");
  jrh->SetFirstEvent(0);
  jrh->SetLastEvent(1000);
  jrh->SetPtCut(2.);
  jrh->SetReadSignalOnly(kFALSE);

  AliJetKineReaderHeader *krh = new AliJetKineReaderHeader();
  krh->SetComment("testing");
  krh->SetDirectory("/home/morsch/AliRoot/newio/FASTSIM/ctest/100/uq");
  krh->SetFirstEvent(0);
  krh->SetLastEvent(1000);
  krh->SetPtCut(2.);
  krh->SetFastSimTPC(kTRUE);
  

  // define reader and set its header
  AliJetESDReader *jr = new AliJetESDReader();
  jr->SetReaderHeader(jrh);

  // define reader and set its header
  AliJetKineReader *kr = new AliJetKineReader();
  kr->SetReaderHeader(krh);

  // define jet header
  AliUA1JetHeader *jh=new AliUA1JetHeader();
  jh->SetComment("UA1 jet code with default parameters");
  jh->SetMode(0);
  jh->SetRadius(0.4);
  jh->SetMinCellEt(0.);
  jh->SetEtSeed(4.);
  jh->SetNbinPhi(420.);
  jh->SetNbinEta(120.);
  jh->SetEtaMin(-0.9);
  jh->SetEtaMax(+0.9);  
  

  // define jet finder. Set its header and reader
  AliUA1JetFinder *jf = new AliUA1JetFinder();
  jf->SetJetHeader(jh);
  jf->SetJetReader(jr);
  jf->SetPlotMode(kTRUE);
  jf->SetOutputFile("jets.root");
  // do the job
  jf->Run();
}