First commit.

[u/mrichter/AliRoot.git] / EMCAL / AnaJets.C

diff --git a/EMCAL/AnaJets.C b/EMCAL/AnaJets.C
new file mode 100644 (file)
index 0000000..becb68e
--- /dev/null
+++ b/EMCAL/AnaJets.C
@@ -0,0 +1,64 @@
+void AnaJets(Int_t evNumber1=0, Int_t evNumber2=0) 
+{
+//*-- Author: Andreas Morsch (CERN)
+
+    if (gClassTable->GetID("AliRun") < 0) {
+       gROOT->LoadMacro("../macros/loadlibs.C");
+       loadlibs();
+    }
+// Connect the Root Galice file containing Geometry, Kine and Hits
+
+    TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
+    
+    if (!file) {
+       printf("\n Creating galice.root \n");
+       file = new TFile("galice.root");
+    } else {
+       printf("\n galice.root found in file list");
+    }
+// Get AliRun object from file or create it if not on file
+    if (!gAlice) {
+       gAlice = (AliRun*)(file->Get("gAlice"));
+       if (gAlice) printf("AliRun object found on file\n");
+       if (!gAlice) {
+           printf("\n create new gAlice object");
+           gAlice = new AliRun("gAlice","Alice test program");
+       }
+    }
+// Book histos    
+    TH1F *eH   = new TH1F("eH","Energy", 150,  0.0, 150.);
+    TH1F *etaH = new TH1F("eEta","Eta",  180, -0.9, 0.9);
+    TH1F *phiH = new TH1F("ePhi","Phi",   62, -3.1, 3.1);
+ 
+
+    TClonesArray* jets = new TClonesArray("AliEMCALJet",10000);
+
+    for (int nev=0; nev<= evNumber2; nev++) {
+       printf("\n Event .............%d", nev);
+       Int_t nparticles = gAlice->GetEvent(nev);
+       Int_t nbytes     = 0;
+       AliEMCAL *pEMCAL  = (AliEMCAL*) gAlice->GetModule("EMCAL");
+       if (pEMCAL) {
+           TTree *TR = gAlice->TreeR();
+           Int_t nent=TR->GetEntries();
+           TR->SetBranchAddress("Jets", &jets);
+           nbytes += TR->GetEntry(0);
+           Int_t nJet = jets->GetEntries();
+           printf("\n Number of Jets %d", nJet);
+           AliEMCALJet  *mJet;
+           for (Int_t ij=0; ij < nJet; ij++) {
+               mJet = (AliEMCALJet*)jets->UncheckedAt(ij);
+               printf("\n Jet:%d E %f phi %f eta %f\n", ij, 
+                      mJet->Energy(), mJet->Phi(), mJet->Eta());
+               etaH->Fill(mJet->Eta());
+               phiH->Fill(mJet->Phi());
+               eH  ->Fill(mJet->Energy());             
+          } // jet
+       } // ?EMCAL
+   } // event
+    TCanvas *c1 = new TCanvas("c1","Canvas 1",400,10,600,700);
+    c1->Divide(2,2);
+    c1->cd(1); eH->Draw();
+    c1->cd(2); etaH->Draw();
+    c1->cd(3); phiH->Draw();
+}