HLT/JET/macros/tasks/ConfigJetAnalysisHLTMC.C

   1 /** ---------------------------------------------------------------------
   2  *  @file   ConfigJetAnalysisHLTMC.C
   3  *  @author Jochen Thaeder <thaeder@kip.uni-heidelberg.de>
   4  *  @brief  Run HLT cone finder in analysis framework, on Kinematics
   5  *
   6  *  ---------------------------------------------------------------------
   7  */
   8
   9 AliJetFinder*  ConfigJetAnalysis() {
  10
  11   printf("ConfigJetAnalysis() -- HLT \n");
  12
  13   // ---------------------------------------------------------------------
  14   // -- Defaults
  15   // ---------------------------------------------------------------------
  16
  17   TString comment       = "HLT Fast Fixed Seeded Cone finder on MC";
  18   AliHLTJETBase::JetAlgorithmType_t algorithm = AliHLTJETBase::kFFSCSquareCell;
  19
  20   Bool_t  leading       = kFALSE;
  21   Float_t coneRadius    =  0.4;
  22   Float_t trackCutMinPt =  1.0;
  23   Float_t seedCutMinPt  =  5.0;
  24   Float_t jetCutMinEt   = 15.0;
  25   Bool_t  useMC         = kTRUE;
  26
  27   // -- Jet Track Cuts
  28   // ---------------------------------------------------------------------
  29   AliHLTJETTrackCuts *trackCuts = new AliHLTJETTrackCuts();
  30   trackCuts->SetChargedOnly( kTRUE );
  31   trackCuts->SetMinPt( trackCutMinPt );
  32
  33   // -- Jet Seed Cuts
  34   // ---------------------------------------------------------------------
  35   AliHLTJETConeSeedCuts *seedCuts = new AliHLTJETConeSeedCuts();
  36   seedCuts->SetMinPt( seedCutMinPt );
  37
  38   // -- Jet Jet Cuts
  39   // ---------------------------------------------------------------------
  40   AliHLTJETJetCuts *jetCuts = new AliHLTJETJetCuts();
  41   jetCuts->SetMinEt( jetCutMinEt );
  42
  43   // -- Jet Reader Header
  44   // ---------------------------------------------------------------------
  45   AliHLTJETReaderHeader *jetReaderHeader = new AliHLTJETReaderHeader();
  46
  47   // Set Algorithm
  48   jetReaderHeader->SetJetAlgorithm(algorithm);
  49
  50   // Set prt to track cuts
  51   jetReaderHeader->SetTrackCuts( trackCuts );
  52   jetReaderHeader->SetSeedCuts( seedCuts );
  53
  54   // Set Eta min/max and Phi min/max
  55   jetReaderHeader->SetFiducialEta( -0.9, 0.9) ;
  56   jetReaderHeader->SetFiducialPhi(  0.0, TMath::TwoPi() ) ;
  57
  58   // Set grid binning
  59   jetReaderHeader->SetGridEtaBinning( 0.05 );
  60   jetReaderHeader->SetGridPhiBinning( 0.05 );
  61
  62   // Set cone radius
  63   jetReaderHeader->SetConeRadius(coneRadius);
  64
  65   // Use Kinematics
  66   jetReaderHeader->SetUseMC(useMC);
  67
  68   // -- Jet Reader
  69   // ---------------------------------------------------------------------
  70   AliHLTJETReader *jetReader = new AliHLTJETReader();
  71   jetReader->SetReaderHeader(jetReaderHeader);
  72
  73   // ---------------------------------------------------------------------
  74   // -- Jet Container
  75   // ---------------------------------------------------------------------
  76   AliHLTJets *jets = new AliHLTJets();
  77   jets->SetComment(comment);
  78
  79   // ---------------------------------------------------------------------
  80   // -- Jet Header
  81   // ---------------------------------------------------------------------
  82   AliHLTJETConeHeader *jetHeader = new AliHLTJETConeHeader();
  83   jetHeader->SetJetCuts(jetCuts);
  84   jetHeader->SetUseLeading(leading);
  85
  86   // ---------------------------------------------------------------------
  87   // -- Jet Finder
  88   // ---------------------------------------------------------------------
  89   AliHLTJETConeFinder *jetFinder = new AliHLTJETConeFinder();
  90   jetFinder->SetJetHeader(jetHeader);
  91   jetFinder->SetJetReader(jetReader);
  92   jetFinder->SetOutputJets(jets);
  93
  94   return jetFinder;
  95 }