]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - EVE/alice-macros/histo2d.C
fixes for new CaloCells contents
[u/mrichter/AliRoot.git] / EVE / alice-macros / histo2d.C
index f030f2457b98e4ada82cbe6824b7a076783e8bd4..38fdec839be4d46832639b41ca06f225e90647f9 100644 (file)
@@ -1,15 +1,43 @@
+// Author: Stefano Carrazza 2010
+
 /**************************************************************************
  * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
  * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for          *
  * full copyright notice.                                                 *
  **************************************************************************/
 
+#if !defined(__CINT__) || defined(__MAKECINT__)
+#include <TGLViewer.h>
+#include <TGLWidget.h>
+#include <TH2.h>
+#include <TMath.h>
+#include <TTree.h>
+#include <TEveBrowser.h>
+#include <TEveCalo.h>
+#include <TEveCaloData.h>
+#include <TEveCaloLegoOverlay.h>
+#include <TEveLegoEventHandler.h>
+#include <TEveManager.h>
+#include <TEveScene.h>
+#include <TEveTrans.h>
+#include <TEveViewer.h>
+#include <TEveWindow.h>
+
+#include <AliESDEvent.h>
+#include <AliEveEventManager.h>
+#include <AliEveMultiView.h>
+#endif
+
 double pi = TMath::Pi();
 TEveViewer *g_histo2d_v   = 0;
 TEveScene  *g_histo2d_s   = 0;
 TEveScene  *g_histo2d_s2  = 0;
 TEveCaloLegoOverlay* g_histo2d_lego_overlay = 0;
 
+Double_t GetPhi(Double_t phi);
+TEveCaloLego* CreateHistoLego(TEveCaloData* data);
+TEveCalo3D* Create3DView(TEveCaloData* data);
+AliEveMultiView* CreateProjections(TEveCaloData* data, TEveCalo3D *calo3d);
 
 TEveCaloDataHist* histo2d()
 { 
@@ -21,19 +49,18 @@ TEveCaloDataHist* histo2d()
    TH2F *histopos = new TH2F("histopos","Histo 2d positive",100,-1.5,1.5,80,-pi,pi);
    TH2F *histoneg = new TH2F("histoneg","Histo 2d negative",100,-1.5,1.5,80,-pi,pi);
 
-   cout<<"Event: "<<AliEveEventManager::GetMaster()->GetEventId()
-       <<", Number of tracks: "<<esd->GetNumberOfTracks()<<endl;
+   Info("histo2d", "Event: %d, Number of tracks: %d\n", AliEveEventManager::GetMaster()->GetEventId(), esd->GetNumberOfTracks() );
 
    // Getting current tracks, filling histograms 
-   for (int n = 0; n < esd->GetNumberOfTracks(); ++n) {    
+   for ( int n = 0; n < esd->GetNumberOfTracks(); ++n ) {    
   
       if (esd->GetTrack(n)->GetSign() > 0) {
          histopos->Fill(esd->GetTrack(n)->Eta(),
-                       getphi(esd->GetTrack(n)->Phi()),
+                       GetPhi(esd->GetTrack(n)->Phi()),
                         fabs(esd->GetTrack(n)->Pt()));
       } else {
          histoneg->Fill(esd->GetTrack(n)->Eta(),
-                        getphi(esd->GetTrack(n)->Phi()),
+                        GetPhi(esd->GetTrack(n)->Phi()),
                         fabs(esd->GetTrack(n)->Pt()));
       }
    }
@@ -52,13 +79,13 @@ TEveCaloDataHist* histo2d()
    data->IncDenyDestroy();
 
    // Plotting the lego histogram in a new tab
-   Create_histo_lego(data);
+   CreateHistoLego(data);
    
    // Plotting the 3D histogram
-   TEveCalo3D *calo3d = Create_3D_view(data);
+   TEveCalo3D *calo3d = Create3DView(data);
 
    // Plotting projections RPhi and RhoZ
-   Create_projections(data, calo3d);
+   CreateProjections(data, calo3d);
    
    gEve->Redraw3D(kTRUE);
 
@@ -66,7 +93,7 @@ TEveCaloDataHist* histo2d()
 }
 
 //______________________________________________________________________________
-Double_t getphi(Double_t phi)
+Double_t GetPhi(Double_t phi)
 {
    if (phi > pi) {
       phi -= 2*pi;
@@ -75,12 +102,12 @@ Double_t getphi(Double_t phi)
 }
 
 //______________________________________________________________________________
-TEveCaloLego* Create_histo_lego(TEveCaloData* data){
+TEveCaloLego* CreateHistoLego(TEveCaloData* data){
 
    TGLViewer* glv;
 
    // Viewer initialization, tab creation
-   if (g_histo2d_v == 0) {
+   if ( g_histo2d_v == 0 ) {
       TEveWindowSlot *slot    = 0;
       TEveBrowser    *browser = gEve->GetBrowser();
 
@@ -119,10 +146,10 @@ TEveCaloLego* Create_histo_lego(TEveCaloData* data){
 }
 
 //______________________________________________________________________________
-TEveCalo3D* Create_3D_view(TEveCaloData* data){
+TEveCalo3D* Create3DView(TEveCaloData* data){
   
    //initialization
-   if (g_histo2d_s2 == 0) {
+   if ( g_histo2d_s2 == 0 ) {
       g_histo2d_s2 = gEve->SpawnNewScene("3D Histogram", "3D Histogram");
       gEve->GetDefaultViewer()->AddScene(g_histo2d_s2);
       g_histo2d_s2->SetElementName("3D Histogram Scene");
@@ -139,7 +166,7 @@ TEveCalo3D* Create_3D_view(TEveCaloData* data){
 }
 
 //______________________________________________________________________________
-AliEveMultiView* Create_projections(TEveCaloData* data, TEveCalo3D *calo3d){
+AliEveMultiView* CreateProjections(TEveCaloData* data, TEveCalo3D *calo3d){
 
    AliEveMultiView *al = AliEveMultiView::Instance();
    al->ImportEventRPhi(calo3d);