TString inPidTPCFile  = "pid_tpccombi.root";

  TFile *inFile = TFile::Open(inPidTPCFile,"READ");
  TTree *lhe=(TTree *) inFile->Get("cbmsim") ;

...

  TClonesArray* cluster_array=new TClonesArray("PndEmcCluster");
  lhe->SetBranchAddress("EmcCluster",&cluster_array);

...

// loop over  events
    NTevents=lhe->GetEntriesFast();
    for (Int_t j=0; j< NTevents ; j++){
       lhe->GetEntry(j);   // kinematics

// loop over clusters
       Int_t nclusters= cluster_array->GetEntriesFast();
       for (Int_t nc = 0; nc < nclusters; nc++) {
           PndEmcCluster *cl = (PndEmcCluster *)cluster_array->At(nc);
// access to energy and moments of cluster
           Double_t Energy = cl->energy(); 
           Double_t Z20 = cl->Z20(); 
           Double_t Z53 = cl->Z53(); 
           Double_t LatMom = cl->LatMom(); 

// (hypothetical) loop over crystals in this cluster
           Int_t ncrystals= crystal_array->GetEntriesFast();
           for (Int_t nx = 0; nx < ncrystals; nx++) {

              ...etc...
           }
        }
     }