]>
Commit | Line | Data |
---|---|---|
6aebab11 | 1 | AliJetFinder* ConfigJetAnalysis() |
2 | { | |
3 | // | |
4 | // Configuration goes here | |
5 | // | |
6 | printf("ConfigJetAnalysis() \n"); | |
7 | ||
8 | // Define the grids | |
9 | AliJetGrid *grid = new AliJetGrid(419,119,0.,2*TMath::Pi(),-0.9,0.9); | |
10 | grid->SetGridType(1); | |
11 | grid->InitParams(80.*TMath::Pi()/180,190.*TMath::Pi()/180,-0.7,0.7); | |
12 | grid->SetMatrixIndexes(); | |
13 | grid->SetIndexIJ(); | |
14 | AliJetGrid *grid2 = new AliJetGrid(131,95,80.*TMath::Pi()/180.,190.*TMath::Pi()/180.,-0.7,0.7); | |
15 | grid2->SetGridType(0); | |
16 | grid2->SetMatrixIndexes(); | |
17 | grid2->SetIndexIJ(); | |
18 | ||
19 | // Define ESD reader header | |
20 | AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader(); | |
21 | jrh->SetComment("Testing"); | |
22 | jrh->SetReadSignalOnly(kFALSE); | |
23 | // Detector options: 0 = Charged particles only (MomentumArray) | |
24 | // 1 = Charged particles only (UnitArray) | |
25 | // 2 = Neutral cells only (UnitArray) | |
26 | // 3 = Charged particles + neutral cells (UnitArray) | |
27 | jrh->SetDetector(0); | |
28 | jrh->SetDebug(0); | |
29 | jrh->SetFiducialEta(-0.9,0.9); | |
30 | jrh->SetFiducialPhi(0,2*TMath::Pi()); | |
31 | ||
32 | // Define reader and set its header | |
33 | AliJetESDReader *er = new AliJetESDReader(); | |
34 | er->SetReaderHeader(jrh); | |
35 | er->SetTPCGrid(grid); | |
36 | er->SetEMCalGrid(grid2); | |
37 | ||
38 | // Define jet header | |
39 | AliFastJetHeaderV1 *jh=new AliFastJetHeaderV1(); | |
40 | jh->SetComment("Fast jet code with default parameters"); | |
41 | jh->SetDebug(0); | |
42 | jh->SetBGMode(0); | |
43 | jh->SetRparam(1.0); // setup parameters | |
44 | jh->SetPtMin(20.); | |
45 | ||
46 | // Define jet finder. Set its header and reader | |
47 | jetFinder = new AliFastJetFinder(); | |
48 | jetFinder->SetJetHeader(jh); | |
49 | jetFinder->SetJetReader(er); | |
50 | jetFinder->SetPlotMode(kTRUE); | |
51 | ||
52 | return jetFinder; | |
53 | } |