Protections for coverity: DIVIDE_BY_ZERO
[u/mrichter/AliRoot.git] / JETAN / ConfigJetAnalysisUA1.C
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 }