]>
Commit | Line | Data |
---|---|---|
5b312f44 | 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->SetPtCut(.0); | |
23 | jrh->SetReadSignalOnly(kFALSE); | |
24 | jrh->SetReadBkgdOnly(kFALSE); | |
25 | // Detector options: 0 = Charged particles only (MomentumArray) | |
26 | // 1 = Charged particles only (UnitArray) | |
27 | // 2 = Neutral cells only (UnitArray) | |
28 | // 3 = Charged particles + neutral cells (UnitArray) | |
29 | jrh->SetDetector(0); | |
30 | jrh->SetDebug(0); | |
31 | jrh->SetFiducialEta(-0.9,0.9); | |
32 | jrh->SetFiducialPhi(0,2*TMath::Pi()); | |
33 | ||
34 | // Define reader and set its header | |
35 | AliJetESDReader *er = new AliJetESDReader(); | |
36 | er->SetReaderHeader(jrh); | |
37 | er->SetTPCGrid(grid); | |
38 | er->SetEMCalGrid(grid2); | |
39 | ||
40 | // Define jet header | |
41 | AliUA1JetHeaderV1 *jh=new AliUA1JetHeaderV1(); | |
42 | jh->SetComment("UA1 jet code with default parameters"); | |
43 | jh->BackgMode(0); | |
44 | jh->SetRadius(1.0); | |
45 | jh->SetEtSeed(2.); | |
46 | jh->SetLegoNbinPhi(420.); | |
47 | jh->SetLegoNbinEta(120.); | |
48 | jh->SetLegoEtaMin(-0.9); | |
49 | jh->SetLegoEtaMax(+0.9); | |
50 | jh->SetMinJetEt(10.); | |
51 | ||
52 | // Define jet finder. Set its header and reader | |
53 | jetFinder = new AliUA1JetFinderV2(); | |
54 | jetFinder->SetJetHeader(jh); | |
55 | jetFinder->SetJetReader(er); | |
56 | jetFinder->SetPlotMode(kTRUE); | |
57 | ||
58 | return jetFinder; | |
59 | } |