1 void testJetFinderBg(Int_t evNumber1=0, Int_t evNumber2=0)
3 //*-- Author: Andreas Morsch (CERN)
4 // Dynamically link some shared libs
5 if (gClassTable->GetID("AliRun") < 0) {
6 gROOT->LoadMacro("../macros/loadlibs.C");
10 // Create and configure JetFinder
12 AliEMCALJetFinder* jetFinder =
13 new AliEMCALJetFinder("UA1 Jet Finder", "Test");
16 jetFinder->SetDebug(1);
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);
27 // Parameters for jet finding
28 jetFinder->SetConeRadius(0.5);
29 jetFinder->SetEtSeed(4.);
30 jetFinder->SetMinJetEt(10.);
31 jetFinder->SetMinCellEt(1.);
37 // This part will go into class
39 // Open background file
42 printf("\n Opening Background File !!\n");
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();
57 printf("\n Opening Signal file !!\n");
58 TFile* fileS = new TFile("galice.root", "update");
59 gAlice = (AliRun*)(fileS->Get("gAlice"));
62 // Loop over signal events
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;
71 jetFinder->InitFromBackground();
73 jetFinder->FillFromHits(1);
74 // ^ preserves info from background
75 // jetFinder->FillFromDigits();
77 jetFinder->FillFromTracks(1, 0);
78 // ^ preserves info from hit
80 // TPC information from Hits associated to the EMCAL
81 // jetFinder->FillFromHitFlaggedTracks(1);
87 printf("\n Test Jets: %d\n", jetFinder->Njets());
88 Int_t njet = jetFinder->Njets();
89 for (Int_t nj=0; nj<njet; nj++)
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));
98 // TCanvas *c1 = new TCanvas("c1","Canvas 1",400,10,600,700);
99 // (jetFinder->GetLego())->Draw();