JETAN/ConfigJetAnalysisUA1.C

   1 AliJetFinder*  ConfigJetAnalysis()
   2 {
   3     //
   4     // Configuration goes here
   5     //
   6     printf("ConfigJetAnalysis() \n");
   7
   8     // Define the grids
   9     AliJetGrid *grid = new AliJetGrid(419,119,0.,2*TMath::Pi(),-0.9,0.9);
  10     grid->SetGridType(1);
  11     grid->InitParams(80.*TMath::Pi()/180,190.*TMath::Pi()/180,-0.7,0.7);
  12     grid->SetMatrixIndexes();
  13     grid->SetIndexIJ();
  14     AliJetGrid *grid2 = new AliJetGrid(131,95,80.*TMath::Pi()/180.,190.*TMath::Pi()/180.,-0.7,0.7);
  15     grid2->SetGridType(0);
  16     grid2->SetMatrixIndexes();
  17     grid2->SetIndexIJ();
  18
  19     // Define ESD reader header
  20     AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader();
  21     jrh->SetComment("Testing");
  22     jrh->SetPtCut(.0);
  23     jrh->SetReadSignalOnly(kFALSE);
  24     jrh->SetReadBkgdOnly(kFALSE);
  25     // Detector options: 0 = Charged particles only (MomentumArray)
  26     //                   1 = Charged particles only (UnitArray)
  27     //                   2 = Neutral cells only (UnitArray)
  28     //                   3 = Charged particles + neutral cells (UnitArray)
  29     jrh->SetDetector(0);
  30     jrh->SetDebug(0);
  31     jrh->SetFiducialEta(-0.9,0.9);
  32     jrh->SetFiducialPhi(0,2*TMath::Pi());
  33
  34     // Define reader and set its header
  35     AliJetESDReader *er = new AliJetESDReader();
  36     er->SetReaderHeader(jrh);
  37     er->SetTPCGrid(grid);
  38     er->SetEMCalGrid(grid2);
  39
  40     // Define jet header
  41     AliUA1JetHeaderV1 *jh=new AliUA1JetHeaderV1();
  42     jh->SetComment("UA1 jet code with default parameters");
  43     jh->BackgMode(0);
  44     jh->SetRadius(1.0);
  45     jh->SetEtSeed(2.);
  46     jh->SetLegoNbinPhi(420.);
  47     jh->SetLegoNbinEta(120.);
  48     jh->SetLegoEtaMin(-0.9);
  49     jh->SetLegoEtaMax(+0.9);
  50     jh->SetMinJetEt(10.);
  51
  52     // Define jet finder. Set its header and reader
  53     jetFinder = new AliUA1JetFinderV2();
  54     jetFinder->SetJetHeader(jh);
  55     jetFinder->SetJetReader(er);
  56     jetFinder->SetPlotMode(kTRUE);
  57
  58     return jetFinder;
  59 }