JETAN/ConfigJetAnalysisMC.C

   1 AliJetFinder*  ConfigJetAnalysis()
   2 {
   3     //
   4     // Configuration goes here
   5     //
   6     printf("ConfigJetAnalysis() \n");
   7     Printf("For MC events");
   8     AliJetKineReaderHeader *jrh = new AliJetKineReaderHeader();
   9     jrh->SetComment("MC full Kinematics");
  10     jrh->SetFastSimTPC(kFALSE);
  11     jrh->SetFastSimEMCAL(kFALSE);
  12     jrh->SetPtCut(0.);
  13     jrh->SetFiducialEta(-2.1,2.1); // to take all MC particles default is 0.9
  14
  15     // Define reader and set its header
  16     AliJetKineReader *er = new AliJetKineReader();
  17     er->SetReaderHeader(jrh);
  18
  19
  20     // Define jet header
  21     AliUA1JetHeaderV1 *jh=new AliUA1JetHeaderV1();
  22     jh->SetComment("UA1 jet code with  parameters similar to PYCELL during sim.");
  23     jh->BackgMode(0);
  24     jh->SetRadius(1.0);
  25     jh->SetEtSeed(4.);
  26     jh->SetLegoNbinPhi(432);
  27     jh->SetLegoNbinEta(274);
  28     jh->SetLegoEtaMin(-2);
  29     jh->SetLegoEtaMax(+2);
  30     jh->SetJetEtaMax(0.5);
  31     jh->SetJetEtaMin(-0.5);
  32     jh->SetMinJetEt(10.);
  33
  34
  35     // Define jet finder. Set its header and reader
  36     jetFinder = new AliUA1JetFinderV1();
  37     jetFinder->SetJetHeader(jh);
  38     jetFinder->SetJetReader(er);
  39     return jetFinder;
  40 }