prod/LHC09a1/ConfigJetAnalysis.C

   1 AliJetFinder*  ConfigJetAnalysis()
   2 {
   3     //
   4     // Configuration goes here
   5     //
   6     printf("ConfigJetAnalysis() \n");
   7     printf("Reading AODs");
   8     AliJetAODReaderHeader *jrh = new AliJetAODReaderHeader();
   9     jrh->SetComment("AOD Reader");
  10     jrh->SetPtCut(0.);
  11     jrh->SetTestFilterMask(1<<1); // hard cuts...
  12     // Define reader and set its header
  13     AliJetAODReader *er = new AliJetAODReader();
  14     er->SetReaderHeader(jrh);
  15
  16     // Define jet header
  17     AliUA1JetHeaderV1 *jh=new AliUA1JetHeaderV1();
  18     jh->SetComment("UA1 jet code with default parameters");
  19     jh->BackgMode(0);
  20     jh->SetRadius(0.4);
  21     jh->SetEtSeed(2.);
  22     jh->SetLegoNbinPhi(420.);
  23     jh->SetLegoNbinEta(120.);
  24     jh->SetLegoEtaMin(-1.9);
  25     jh->SetLegoEtaMax(+1.9);
  26     jh->SetMinJetEt(5.);
  27     jh->SetJetEtaMin(-1.5);
  28     jh->SetJetEtaMax(1.5);
  29
  30
  31
  32     // Define jet finder. Set its header and reader
  33     jetFinder = new AliUA1JetFinderV1();
  34     jetFinder->SetJetHeader(jh);
  35     jetFinder->SetJetReader(er);
  36     /*
  37     jetFinder->SetPlotMode(kTRUE);
  38     jetFinder->SetOutputFile("jets.root");
  39     */
  40     //
  41     return jetFinder;
  42 }