JETAN/ConfigJetAnalysisSISCone.C

   1 AliJetFinder*  ConfigJetAnalysis()
   2 {
   3   //
   4   // Configuration goes here
   5   //
   6   printf("ConfigJetAnalysis() \n");
   7   // Define the grids
   8   AliJetGrid *grid = new AliJetGrid(419,119,0.,2*TMath::Pi(),-0.9,0.9);
   9   grid->SetGridType(1);
  10   grid->InitParams(80.*TMath::Pi()/180,190.*TMath::Pi()/180,-0.7,0.7);
  11   grid->SetMatrixIndexes();
  12   grid->SetIndexIJ();
  13   AliJetGrid *grid2 = new AliJetGrid(131,95,80.*TMath::Pi()/180.,190.*TMath::Pi()/180.,-0.7,0.7);
  14   grid2->SetGridType(0);
  15   grid2->SetMatrixIndexes();
  16   grid2->SetIndexIJ();
  17
  18   // Define ESD reader header
  19   AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader();
  20   jrh->SetComment("Testing");
  21   jrh->SetReadSignalOnly(kFALSE);
  22   // Detector options: 0 = Charged particles only (MomentumArray)
  23   //                   1 = Charged particles only (UnitArray)
  24   //                   2 = Neutral cells only (UnitArray)
  25   //                   3 = Charged particles + neutral cells (UnitArray)
  26   jrh->SetDetector(0);
  27   jrh->SetDebug(0);
  28   jrh->SetFiducialEta(-0.9,0.9);
  29   jrh->SetFiducialPhi(0,2*TMath::Pi());
  30
  31   // Define reader and set its header
  32   AliJetESDReader *er = new AliJetESDReader();
  33   er->SetReaderHeader(jrh);
  34   er->SetTPCGrid(grid);
  35   er->SetEMCalGrid(grid2);
  36
  37   // Define jet finder header
  38   AliSISConeJetHeader * jh = new AliSISConeJetHeader();
  39   //siscone parameters
  40   jh->SetConeRadius(1.0);                   // cone radius
  41   jh->SetOverlapThreshold(0.75);            // overlap parameter, between 0 and 1 excluded!! 0.75 value is advised
  42   jh->SetMinJetPt(20);                      // Ptmin of jets
  43   //do you want to subtract BG (0 = no, 1 = yes)
  44   jh->SetBGMode(0);
  45   //to determine jets area (for BG subtraction)
  46   jh->SetAreaTypeNumber(4);                 // from 1 to 4 : 1 = active_area, 2 = active_area_explicit_ghosts, 3 = one_ghost_passive_area, 4 = passive_area
  47
  48   // Define jet finder
  49   AliSISConeJetFinder *jetFinder = new AliSISConeJetFinder();
  50   jetFinder->SetJetHeader(jh);
  51   jetFinder->SetJetReader(er);
  52
  53   return jetFinder;
  54 }
  55
  56
  57
  58
  59