]>
Commit | Line | Data |
---|---|---|
5a15bc2d | 1 | void 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 |