Made more robust
[u/mrichter/AliRoot.git] / EMCAL / testJetFinderBg.C
CommitLineData
5a15bc2d 1void testJetFinderBg(Int_t evNumber1=0, Int_t evNumber2=0)
2{
3//*-- Author: Andreas Morsch (CERN)
4// Dynamically link some shared libs
5 if (gClassTable->GetID("AliRun") < 0) {
6 gROOT->LoadMacro("../macros/loadlibs.C");
7 loadlibs();
8 }
9//
10// Create and configure JetFinder
11//
12 AliEMCALJetFinder* jetFinder =
13 new AliEMCALJetFinder("UA1 Jet Finder", "Test");
14//
15// Debug
16 jetFinder->SetDebug(1);
17//
18// Input and fast simulation
19 jetFinder->SetPtCut(2.);
20 jetFinder->SetIncludeK0andN(0);
21// jetFinder->SetMomentumSmearing(1);
22// jetFinder->SetEfficiencySim(1);
23 jetFinder->SetHadronCorrection(0);
24// jetFinder->SetHadronCorrector(AliEMCALHadronCorrectionv0::Instance());
25 jetFinder->SetSamplingFraction(12.9);
26//
27// Parameters for jet finding
28 jetFinder->SetConeRadius(0.5);
29 jetFinder->SetEtSeed(4.);
30 jetFinder->SetMinJetEt(10.);
31 jetFinder->SetMinCellEt(1.);
32
33//.............
34// |
35// |
36// V
37// This part will go into class
38//
39// Open background file
40//
41//
42 printf("\n Opening Background File !!\n");
43
44 TFile* fileB = new TFile("bg.root");
45 gAlice = (AliRun*)(fileB->Get("gAlice"));
46 Int_t nparticles = gAlice->GetEvent(0);
47 // Read and save background event
48 jetFinder->FillFromHits();
49 jetFinder->FillFromTracks();
50 jetFinder->SaveBackgroundEvent();
51 delete gAlice;
52
53
54//
55// Open signal file
56//
57 printf("\n Opening Signal file !!\n");
58 TFile* fileS = new TFile("galice.root", "update");
59 gAlice = (AliRun*)(fileS->Get("gAlice"));
60
61//
62// Loop over signal events
63//
64
65 Int_t nhit=0;
66 for (Int_t nev = evNumber1; nev<= evNumber2; nev++) {
67 Int_t nparticles = gAlice->GetEvent(nev);
68 if (nev < evNumber1) continue;
69 if (nparticles <= 0) return;
70// Load background
71 jetFinder->InitFromBackground();
72// ECAL information
73 jetFinder->FillFromHits(1);
74// ^ preserves info from background
75// jetFinder->FillFromDigits();
76// TPC information
77 jetFinder->FillFromTracks(1, 0);
78// ^ preserves info from hit
79
80// TPC information from Hits associated to the EMCAL
81// jetFinder->FillFromHitFlaggedTracks(1);
82//
83 jetFinder->Find();
84
85//
86// Look at results
87 printf("\n Test Jets: %d\n", jetFinder->Njets());
88 Int_t njet = jetFinder->Njets();
89 for (Int_t nj=0; nj<njet; nj++)
90 {
91 printf("\n Jet Energy %5d %8.2f \n",
92 nj, jetFinder->JetEnergy(nj));
93 printf("\n Jet Phi %5d %8.2f %8.2f \n",
94 nj, jetFinder->JetPhiL(nj), jetFinder->JetPhiW(nj));
95 printf("\n Jet Eta %5d %8.2f %8.2f \n",
96 nj, jetFinder->JetEtaL(nj), jetFinder->JetEtaW(nj));
97 }
98// TCanvas *c1 = new TCanvas("c1","Canvas 1",400,10,600,700);
99// (jetFinder->GetLego())->Draw();
100 } // event loop
101
102// ^
103// |
104// |
105// |
106//.............
107 fileB->Close();
108 fileS->Close();
109}
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125