Split: HLT/JET -> HLTANALYSIS/JET
[u/mrichter/AliRoot.git] / HLTANALYSIS / JET / macros / tasks / ConfigJetAnalysisHLT.C
CommitLineData
1f9fec4a 1/** ---------------------------------------------------------------------
2 * @file ConfigJetAnalysisHLT.C
79fe198d 3 * @author Jochen Thaeder <jochen@thaeder.de>
1f9fec4a 4 * @brief Run HLT cone finder in analysis framework, on ESDs
5 *
6 * ---------------------------------------------------------------------
7 */
8
9AliJetFinder* 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
79fe198d 51 jetReaderHeader->SetTrackCuts(trackCuts);
52 jetReaderHeader->SetSeedCuts(seedCuts);
1f9fec4a 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}