---

readESD.cpp

Go to the documentation of this file.

#include <XANADOO/domain/interface/XANAEsdEvent.h>
#include <XANADOO/XANAClusters/interface/XANACluster.h>
#include <XANADOO/XANAClusters/interface/XANACaloRecHit.h>
#include <XANADOO/XANAClusters/interface/XANAEmCluster.h>
#include <XANADOO/XANAClusters/interface/XANASuperCluster.h>
#include <XANADOO/XANAClusters/interface/XANAHadCluster.h>
#include <XANADOO/XANATracks/interface/XANATrack.h>
#include <XANADOO/XANAMcInfo/interface/XANAGeneratorEvent.h>
#include <XANADOO/XANAMcInfo/interface/XANAGeantEvent.h>
#include <XANADOO/XANAMcInfo/interface/XANAGeneratorParticle.h>
#include <XANADOO/XANAElectronCandidate/interface/XANAElectronCandidate.h>
#include <XANADOO/XANAElectronCandidate/interface/XANAElectronTrack.h>
#include <XANADOO/XANAElectronCandidate/interface/XANAElectronSeed.h>
#include <XANADOO/XANAMuons/interface/XANAMuonCandidate.h>
#include <XANADOO/XANATriggerInfo/interface/XANATriggerInfo.h>
#include <XANADOO/XANAJets/interface/XANAJet.h>

#include "TTree.h"
#include "TBits.h"
#include "TObjectTable.h"
#include "TFile.h"
#include "TBranch.h"
#include "TClonesArray.h"
#include "TObjArray.h"
#include "TH1.h"
#include "TH2.h"
#include "TSystem.h"
#include "TChain.h"

#include <string>


using namespace std;

int main(int argc, char** argv)
{
  //   Connect to file generated
  TFile *f;
    //  if (argc==1) f = new TFile("SE04_E25_Pall_PT30.root");
  if (argc==1) f = new TFile("h300eemm.root");
  else f = new TFile(*++argv);
 
    //   Read Tree named "T" in memory. Tree pointer is assigned the same name
  TTree *T = (TTree*)f->Get("ESD");
   
// version TChain
//   TChain *T = new TChain("ESD");
//   T->Add("file1.root");
//   T->Add("file2.root");
 
// example of how to suppress reading of branches
// has to be put before setting the branch addresses!
// T->SetBranchStatus("GeantEvent*",0);

// Reconstructed event
  XANAEsdEvent *event = new XANAEsdEvent();  
  T->SetBranchAddress("Event.",&event);
  // MC information (generator and GEANT event)
  XANAGeneratorEvent *generatorEvent = new XANAGeneratorEvent();
  T->SetBranchAddress("GeneratorEvent.",&generatorEvent);
    
  XANAGeantEvent *geantEvent = new XANAGeantEvent();
  T->SetBranchAddress("GeantEvent.",&geantEvent);
  
  Short_t nc = 0, nh = 0, nhc = 0, ns = 0, npc = 0,  nt = 0, nv = 0, nvt = 0;  
  Short_t ne = 0, nm=0, nelt = 0;
  
  TObjArray *clusters = 0;
  TObjArray *hadClusters = 0;
  TObjArray *superClusters = 0;
  TObjArray *preshClusters = 0;
  TRefArray *hits = 0;
  TObjArray *tracks = 0;
  TObjArray *vertices = 0;
  TObjArray *electrons = 0;
  TObjArray *electronTracks = 0;
  TObjArray *muons = 0;

  Short_t mc_np = 0, mc_nsp = 0, mc_nssp = 0, mc_npp = 0, mc_nspp = 0, mc_ngt = 0, mc_ngv = 0;  
  
  //   Start main loop on all events
  Int_t nb = 0; 
  Int_t nevent = (Int_t)T->GetEntries();
  cout << "================== Total number of events: " << nevent << " ================" << endl;    
 
  for (Int_t i=0; i<nevent; i++) {
    nb += T->GetEntry(i); 

     cout << "================ Event: " << i << " ================" << endl;    
    
    // Reconstructed event
    
    cout << "Reconstructed:" ;
    
    nc  = event->getNumberOfEmClusters();      clusters = event->getEmClusters(); 
    nh  = event->getNumberOfEmClusterHits();          hits = event->getEmClusterHits();   
    nhc = event->getNumberOfHadClusters();   hadClusters = event->getHadClusters();   
    ns  = event->getNumberOfSuperClusters(); superClusters = event->getSuperClusters(); 
    npc = event->getNumberOfPreshClusters(); preshClusters = event->getPreshClusters(); 
    nt  = event->getNumberOfTracks();        tracks = event->getTracks();
    nv  = event->getNumberOfVertices();      vertices = event->getVertices();   
    nvt = event->getNumberOfVertexTracks();  
    ne  = event->getNumberOfElectronCandidates(); electrons = event->getElectronCandidates();
    nm  = event->getNumberOfMuonCandidates(); muons = event->getMuonCandidates();
    nelt    = event->getNumberOfElectronTracks();
    electronTracks = event->getElectronTracks();
 
    XANATriggerInfo trigger(event->getTriggerInfo());
    
    cout << "  clust, hits, hadClu, supers, presho, tracks, vertex, vTrack, nElec, nMuons" << endl;
    printf("  \t \t %d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\t%d\n", nc, nh, nhc, ns, npc, nt, nv, nvt, ne, nm);
    cout << " Number of electron tracks: "<<nelt <<endl;
    
    // electrons
    for (Int_t ie=0; ie<ne; ie++) {
      XANAElectronCandidate *elec = (XANAElectronCandidate*)((*electrons)[ie]);
      cout << "new electron candidate : " << elec << endl;
      cout << "  4-momentum at vertex = " << elec->getMomentumAtVertex() << endl;
      cout << "  Pt at vertex = " << elec->getPt() << endl;
      cout << "  eta at vertex = " << elec->getEta() << endl;
      cout << "  phi at vertex = " << elec->getPhi() << endl;
      cout << "  track momentum at calo = " << elec->getTrackMomentumAtCalo() << endl;
      cout << "  super cluster position = " << elec->getSuperClusterPosition() << endl;
      cout << "  super cluster (eta,phi) = " << "("<<elec->getSuperClusterPosition().pseudoRapidity()<<
                                          ","<<elec->getSuperClusterPosition().phi()<<")"<< endl;
      cout << "  track position at calo = " << elec->getTrackPositionAtCalo() << endl;
      cout << "  Esuper/p = " << elec->getESuperClusterOverP() << endl;
      cout << "  Eseed/p = " << elec->getESeedClusterOverP() << endl;
      cout << "  deta (seed-track) at calo = " << elec->getDeltaEtaSeedClusterAtCalo() << endl;
      cout << "  dphi (seed-track) at calo = " << elec->getDeltaPhiSeedClusterAtCalo() << endl;
      cout << "  deta (super-track) at calo = " << elec->getDeltaEtaSuperClusterAtCalo() << endl;
      cout << "  dphi (super-track) at calo = " << elec->getDeltaPhiSuperClusterAtCalo() << endl;
      cout << "  deta (super-track) at vtx = " << elec->getDeltaEtaSuperClusterAtVtx() << endl;
      cout << "  dphi (super-track) at vtx = " << elec->getDeltaPhiSuperClusterAtVtx() << endl;
      cout << "  track iso = " << elec->getTrackIsolation() << endl;
      cout << "  calo iso = " << elec->getCaloIsolation() << endl;
      cout << "  had/em = " << elec->getHadronicOverEm() << endl;
      cout << "  algo = " << elec->getAlgoName() << endl;
      cout << " constituant seed cluster : " << endl;
      cout << "  seed " << elec->getSeedCluster() << endl;
      cout << " constituant super cluster : " << endl;
      cout << "  super " << elec->getSuperCluster() << endl;
      if (elec->getSuperCluster()) {
      cout << "  E = " << elec->getSuperCluster()->getEnergy() << endl ;
      cout << "  (x,y,z) = " << elec->getSuperCluster()->getPosition() << endl ;
      }
      cout << " constituant brem clusters : " << endl;
      cout << "  brems " << elec->getBremClusters() << endl;
      if (elec->getBremClusters()) {
//      for (Int_t ibr=0; ibr<elec->getBremClusters()->GetEntries(); ibr++)
//        cout << "new brem cluster with energy " << ((XANACluster *)(elec->getBremClusters())[i])->getEnergy() << endl;
      }
      cout << " constituant electron track : " << endl;
      cout << "  eletk " << elec->getElectronTrack() << endl;
      if (elec->getElectronTrack()) {
      cout << "  q = " << elec->getElectronTrack()->getCharge() << endl ;
      cout << "  chi2 = " << elec->getElectronTrack()->getChi2OverDof() << endl ;
      cout << "  n hits = " << elec->getElectronTrack()->getNumberOfHits() << endl ;
      cout << "  lost hits = " << elec->getElectronTrack()->getNumberOfLostHits() << endl ;
      cout << "  long impact param = " << elec->getElectronTrack()->getLongImpactParameter() << endl ;
      cout << "  trans impact param = " << elec->getElectronTrack()->getTransImpactParameter() << endl ;
      cout << "  3D impact param = " << elec->getElectronTrack()->getImpactParameter() << endl ;
      cout << "  momentum at vertex = " << elec->getElectronTrack()->getMomentumAtVertex() << endl ;
      cout << "  momentum at first hit = " << elec->getElectronTrack()->getMomentumAtFirstPoint() << endl ;
      cout << "  momentum at last hit = " << elec->getElectronTrack()->getMomentumAtLastPoint() << endl ;
      cout << "  position at first hit = " << elec->getElectronTrack()->getPositionAtFirstPoint() << endl ;
      cout << "  position at last hit = " << elec->getElectronTrack()->getPositionAtLastPoint() << endl ;
      cout << "  track algo = " << elec->getElectronTrack()->getAlgoName() << endl;
      cout << "  vertex = " << elec->getElectronTrack()->getVertex() << endl;
      cout << "  super = " << elec->getElectronTrack()->getSuperCluster() << endl;
//      cout << "  super cluster E = " << elec->getElectronTrack()->getSuperCluster()->getEnergy() << endl;
      cout << "  elec = " << elec->getElectronTrack()->getElectronCandidate() << endl;
      cout << "  seed = " << elec->getElectronTrack()->getElectronSeed() << endl;
      if (elec->getElectronTrack()->getElectronSeed()) {
      cout << "  seed n hits = " << elec->getElectronTrack()->getElectronSeed()->getNumberOfHits() << endl;
      cout << "  seed algo = " << elec->getElectronTrack()->getElectronSeed()->getAlgoName() << endl;
      cout << "  seed super cluster = " << elec->getElectronTrack()->getElectronSeed()->getSuperCluster() << endl;
      cout << "  seed electron track = " << elec->getElectronTrack()->getElectronSeed()->getElectronTrack() << endl;
      if (event->getEleTrackHits()) cout << " Nr of entries for electron track hits = " << event->getEleTrackHits()->GetEntries() << endl;
      cout << " Number of electron track hits = "<<event->getNumberOfEleTrackHits() <<  endl;
      }
      }
    }
    cout << endl;
 

    // hits
    Short_t nchits=event->getNumberOfCaloRecHits();
    Short_t ntrhits=event->getNumberOfAllTrackHits();
    Short_t nelerechits=0;
    if (event->getEleTrackHits()) nelerechits=event->getEleTrackHits()->GetEntries();
    cout << nchits<< " CaloRecHits, " << ntrhits << " TrackHits, "<<nelerechits << " ElectronTrackHits. "<< endl;

     // jets 
    cout << " Jets: " << event->getNumberOfJets() <<endl;
    TObjArray *a =event->getJets();
    for (int i=0;i<event->getNumberOfJets();i++) {
      cout <<" Jet "<< i <<" has "<<((XANAJet* )a->At(i))->getNumberOfClusters()<<" HadClusters "<<endl;
    }

   // Monte Carlo information
   
    // generator event
    mc_np   = generatorEvent->getNumberOfParticles();
    mc_nsp  = generatorEvent->getNumberOfSignalParticles();
    mc_nssp = generatorEvent->getNumberOfStableSignalParticles();
    mc_npp  = generatorEvent->getNumberOfPileupParticles();
    mc_nspp = generatorEvent->getNumberOfStablePileupParticles();
    
    // GEANT event
    mc_ngt = geantEvent->getNumberOfGeantTracks();
    mc_ngv = geantEvent->getNumberOfGeantVertices();
    
    cout <<endl;
    cout  << "MC Info: ";
    cout << "parti, signl, stSig, pileup, stPil, geantTr, geantVe" << endl; 
    printf("\t   %d\t %d\t%d\t%d\t%d\t%d\t%d\n", mc_np, mc_nsp, mc_nssp, mc_npp, mc_nspp, mc_ngt, mc_ngv);  
       
    cout << endl;
    cout << "Trigger info:  ";
    cout << "level1 decision " << trigger.getGlobalL1Decision() << " hlt decision : " << trigger.getHltDecision() << endl; 
    cout << "Level 1 detailed response " << endl;
    trigger.getGlobalL1Response().Print();
    cout << "HLT detailed response " << endl;
    trigger.getHltResponse().Print();
   
    //    gObjectTable->Print();

    event->clear(); 
    generatorEvent->clear();
    geantEvent->clear();
  }
  
  cout << endl;
  cout << nevent << " events and " << nb << " bytes read " << endl;  
  cout << endl;
  
  //  f->Close();  
 
}

---