TestEMCALHit.C

Go to the documentation of this file.

#if !defined(__CINT__) || defined(__MAKECINT__)

//Root include files 
#include <Riostream.h>
#include <TClonesArray.h>
#include <TGeoManager.h>

//AliRoot include files 
#include "AliRun.h"
#include "AliRunLoader.h"
#include "AliEMCALLoader.h"
#include "AliEMCAL.h"
#include "AliEMCALHit.h"
#include "AliEMCALGeometry.h"

#endif

void TestEMCALHit()
{  
  // Getting EMCAL Detector and Geometry.
  
  AliRunLoader *rl = AliRunLoader::Open("galice.root",AliConfig::GetDefaultEventFolderName(),"read");
  
  if (rl == 0x0)
    cout<<"Can not instatiate the Run Loader"<<endl;
  
  rl->LoadgAlice();//Needed to get geometry
  
  AliEMCALLoader *emcalLoader = dynamic_cast<AliEMCALLoader*>
    (rl->GetDetectorLoader("EMCAL"));
  
  TGeoManager::Import("geometry.root");
  
  AliRun * alirun   = rl->GetAliRun(); // Needed to get Geometry
  
  AliEMCALGeometry * geom ;
  if(alirun)
  {
    AliEMCAL * emcal  = (AliEMCAL*)alirun->GetDetector("EMCAL");
    geom = emcal->GetGeometry();
  }
  
  if (geom == 0) cout<<"Did not get geometry from EMCALLoader"<<endl;
  //else   geom->PrintGeometry();
  
  //Load Hits
  rl->LoadHits("EMCAL");
  
  //Get maximum number of events
  Int_t maxevent =  rl->GetNumberOfEvents();
  cout<<"Number of events "<<maxevent<<endl;
  //maxevent = 8000 ;
  
  AliEMCALHit * hit;
  TClonesArray *hits = 0;
  
  for (Int_t iEvent=0; iEvent<maxevent; iEvent++)
  {
    //cout <<  " ======> Event " << iEvent <<endl ;  
    // Load Event
    rl->GetEvent(iEvent);
    Float_t elos=-1;
    Float_t time  = -1 ;
    Int_t id      = -1 ;
    Int_t iSupMod =  0 ;
    Int_t iTower  =  0 ;
    Int_t iIphi   =  0 ;
    Int_t iIeta   =  0 ;
    Int_t iphi    =  0 ;
    Int_t ieta    =  0 ;
    
    cout <<  " ======> Event " << iEvent << endl; 
    
    // Get hits from the list      
    // Hits are stored in different branches in the hits Tree, 
    // first get the branch and then access the hits in the branch
    TTree *treeH = emcalLoader->TreeH();
    if (!treeH)
    {
      printf("TreeH not available!\n");
      continue;
    }
    
    // TreeH exists, get the branch
    Int_t nTrack = treeH->GetEntries();  // TreeH has array of hits for every primary
    TBranch * branchH = treeH->GetBranch("EMCAL");
    branchH->SetAddress(&hits);
    
    for (Int_t iTrack = 0; iTrack < nTrack; iTrack++) 
    {
      branchH->GetEntry(iTrack);
      
      // Now get the hits in this branch
      Int_t nHit = hits->GetEntriesFast();
      for(Int_t ihit = 0; ihit < nHit;ihit++)
      {
        hit = static_cast<AliEMCALHit *>(hits->At(ihit));//hits->At(ihit) ;
        
        if(!hit)
        {
          cout<<"Hit pointer 0x0"<<endl;
          continue;
        }
        
        id   = hit->GetId() ;    // cell (hit) label
        elos = hit->GetEnergy(); // amplitude in cell (hit)
        time = hit->GetTime();   // time of creation of hit after collision
        
        cout<<"Hit ID "<<id<<" ELoss "<<elos;
        
        // Geometry methods  
        if(geom)
        {
          geom->GetCellIndex(id,iSupMod,iTower,iIphi,iIeta); 
          
          // Gives SuperModule and Tower numbers
          geom->GetCellPhiEtaIndexInSModule(iSupMod,iTower,
                                            iIphi, iIeta,iphi,ieta);
          //Gives label of cell in eta-phi position per each supermodule
          // cout<< "SModule "<<iSupMod<<"; Tower "<<iTower <<"; Eta "<<iIeta
          //<<"; Phi "<<iIphi<<"; Cell Eta "<<ieta<<"; Cell Phi "<<iphi<<endl;
          cout<< ";  SModule "<<iSupMod<<"; Cell Eta "<<ieta<<"; Cell Phi "<<iphi<<endl;
        }//geom?
      }//hit loop
    }// track loop
  }//event loop
}