]> git.uio.no Git - u/mrichter/AliRoot.git/blame - JETAN/ConfigJetAnalysisSISCone.C
Adding new macro, loadFromOCDB.C, for loading magnetic field, mapping and
[u/mrichter/AliRoot.git] / JETAN / ConfigJetAnalysisSISCone.C
CommitLineData
5b312f44 1AliJetFinder* ConfigJetAnalysis()
2{
3 //
4 // Configuration goes here
5 //
6 printf("ConfigJetAnalysis() \n");
7 // Define the grids
8 AliJetGrid *grid = new AliJetGrid(419,119,0.,2*TMath::Pi(),-0.9,0.9);
9 grid->SetGridType(1);
10 grid->InitParams(80.*TMath::Pi()/180,190.*TMath::Pi()/180,-0.7,0.7);
11 grid->SetMatrixIndexes();
12 grid->SetIndexIJ();
13 AliJetGrid *grid2 = new AliJetGrid(131,95,80.*TMath::Pi()/180.,190.*TMath::Pi()/180.,-0.7,0.7);
14 grid2->SetGridType(0);
15 grid2->SetMatrixIndexes();
16 grid2->SetIndexIJ();
17
18 // Define ESD reader header
19 AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader();
20 jrh->SetComment("Testing");
21 jrh->SetReadSignalOnly(kFALSE);
22 // Detector options: 0 = Charged particles only (MomentumArray)
23 // 1 = Charged particles only (UnitArray)
24 // 2 = Neutral cells only (UnitArray)
25 // 3 = Charged particles + neutral cells (UnitArray)
26 jrh->SetDetector(0);
27 jrh->SetDebug(0);
28 jrh->SetFiducialEta(-0.9,0.9);
29 jrh->SetFiducialPhi(0,2*TMath::Pi());
30
31 // Define reader and set its header
32 AliJetESDReader *er = new AliJetESDReader();
33 er->SetReaderHeader(jrh);
34 er->SetTPCGrid(grid);
35 er->SetEMCalGrid(grid2);
36
37 // Define jet finder header
38 AliSISConeJetHeader * jh = new AliSISConeJetHeader();
39 //siscone parameters
40 jh->SetConeRadius(1.0); // cone radius
41 jh->SetOverlapThreshold(0.75); // overlap parameter, between 0 and 1 excluded!! 0.75 value is advised
42 jh->SetMinJetPt(20); // Ptmin of jets
43 //do you want to subtract BG (0 = no, 1 = yes)
44 jh->SetBGMode(0);
45 //to determine jets area (for BG subtraction)
46 jh->SetAreaTypeNumber(4); // from 1 to 4 : 1 = active_area, 2 = active_area_explicit_ghosts, 3 = one_ghost_passive_area, 4 = passive_area
47
48 // Define jet finder
49 AliSISConeJetFinder *jetFinder = new AliSISConeJetFinder();
50 jetFinder->SetJetHeader(jh);
51 jetFinder->SetJetReader(er);
52
53 return jetFinder;
54}
55
56
57
58
59