update test macros
[u/mrichter/AliRoot.git] / HLT / JET / macros / tasks / ConfigJetAnalysisHLT.C
1 /** ---------------------------------------------------------------------
2  *  @file   ConfigJetAnalysisHLT.C
3  *  @author Jochen Thaeder <jochen@thaeder.de>
4  *  @brief  Run HLT cone finder in analysis framework, on ESDs
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 ESD";
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         = kFALSE;
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 }