Adding extra check for GPU_FORCE_64BIT_PTR env var
[u/mrichter/AliRoot.git] / JETAN / ConfigJetAnalysisUA1.C
CommitLineData
5b312f44 1AliJetFinder* 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}