---

StandardHistogramMaker.cxx

Go to the documentation of this file.

// ================================================
// StandardHistogramMaker class Implementation
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
// Namespace ATLF::
//


#include "AtlfastAlgs/StandardHistogramMaker.h"
#include "AtlfastEvent/ReconstructedParticle.h"
#include "AtlfastEvent/CollectionDefs.h"

// standard library includes
#include <math.h> 

// Gaudi includes
#include "GaudiKernel/DataSvc.h"
#include "AIDA/IHistogram1D.h"
#include "GaudiKernel/IHistogramSvc.h"
//htp killed for Gaudi v4
//include "GaudiKernel/H1D.h"

namespace Atlfast {

//--------------------------------
// Constructors and destructors
//-------------------------------

StandardHistogramMaker::StandardHistogramMaker 
  ( const std::string& name, ISvcLocator* pSvcLocator ) 
    : Algorithm( name, pSvcLocator )
{




  // Default paths for entities in the TES
  m_electronLocation         = "/Event/Electrons";
  m_photonLocation           = "/Event/Photons";
  m_isolatedElectronLocation = "/Event/IsolatedElectrons";
  m_isolatedPhotonLocation   = "/Event/IsolatedPhotons";
  m_cellLocation             = "/Event/AtlfastCells";
  m_clusterLocation          = "/Event/Clusters";
  m_mcTruthLocation          = "/Event/McEventCollection";


  // Defaults for which sets of histograms to do
  m_doElectronHistograms         = true;
  m_doPhotonHistograms           = true;
  m_doIsolatedElectronHistograms = true;
  m_doIsolatedPhotonHistograms   = true;
  m_doCellHistograms             = true;
  m_doClusterHistograms          = true;


  // Default number at which to start the histograms
  m_histStart                    = 0;

  declareProperty( "ElectronLocation", m_electronLocation ) ;
  declareProperty( "PhotonLocation",   m_photonLocation );
  declareProperty( "IsolatedElectronLocation", m_isolatedElectronLocation );
  declareProperty( "IsolatedPhotonLocation", m_isolatedPhotonLocation );
  declareProperty( "CellLocation", m_cellLocation );
  declareProperty( "ClusterLocation", m_clusterLocation );
  declareProperty( "McTruthLocation", m_mcTruthLocation );

  declareProperty( "DoElectronHistograms",m_doElectronHistograms );
  declareProperty( "DoPhotonHistograms", m_doPhotonHistograms );
  declareProperty( "DoIsolatedElectronHistograms",m_doIsolatedElectronHistograms );
  declareProperty( "DoIsolatedPhotonHistograms",m_doIsolatedPhotonHistograms );
  declareProperty( "DoCellHistograms",m_doCellHistograms );
  declareProperty( "DoClusterHistograms",m_doClusterHistograms );

  declareProperty( "HistogramStart", m_histStart );


}
StandardHistogramMaker::~StandardHistogramMaker(){}

//---------------------------------
// initialise() 
//---------------------------------

StatusCode StandardHistogramMaker::initialize()
{
  // .. put any initialisation in here...
  MsgStream log(messageService(), name());
  log << MSG::DEBUG << "in initialize()" << endreq;
  
  m_nHist = m_histStart;
  if ( m_doElectronHistograms ) bookElectronHistograms();
  if ( m_doPhotonHistograms ) bookPhotonHistograms();
  if ( m_doIsolatedElectronHistograms ) bookIsolatedElectronHistograms();
  if ( m_doIsolatedPhotonHistograms ) bookIsolatedPhotonHistograms();
  if ( m_doCellHistograms ) bookCellHistograms();
  if ( m_doClusterHistograms ) bookClusterHistograms();

  log << MSG::DEBUG << "Booked histograms OK" << endreq;
  //  m_tesIO = new TesIO(eventDataService());
  m_tesIO = new TesIO();
  return StatusCode::SUCCESS ;
}



//---------------------------------
// finalise() 
//---------------------------------

StatusCode StandardHistogramMaker::finalize()
{
  // .. put any finalisation in here...
  return StatusCode::SUCCESS ;
}




//----------------------------------------------
// execute() method called once per event
//----------------------------------------------

StatusCode StandardHistogramMaker::execute( )
{

  //................................
  // make a message logging stream

  MsgStream log( messageService(), name() ) ;
  log << MSG::DEBUG << "Executing " << endreq;

  if ( m_doElectronHistograms ) fillElectronHistograms();
  if ( m_doPhotonHistograms ) fillPhotonHistograms();
  if ( m_doIsolatedElectronHistograms ) fillIsolatedElectronHistograms();
  if ( m_doIsolatedPhotonHistograms ) fillIsolatedPhotonHistograms();
  if ( m_doCellHistograms ) fillCellHistograms();
  if ( m_doClusterHistograms ) fillClusterHistograms();
    
  return StatusCode::SUCCESS ;
}




void StandardHistogramMaker::bookElectronHistograms () {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "Booking electron histograms" << endreq;

  // first increment the histogram count
  m_nHist++;
  // then book the histogram
  m_h_electronMultiplicity = 
    histoSvc()->book("/stat/simple/",m_nHist,"Electron Multiplicity",20,0.0,20.0);
  // or be a C++ guru and do both at once
  m_h_electronPt = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Electron Pt",100,0.0,1000.0);

  //MC truth electrons
  

  return;
}



void StandardHistogramMaker::bookPhotonHistograms () {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "Booking photon histograms" << endreq;

  m_h_photonMultiplicity = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Photon Multiplicity",100,0.0,100.0);
  // or be a C++ guru and do both at once
  m_h_photonPt = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Photon Pt",100,0.0,1000.0);
  return;
}



void StandardHistogramMaker::bookIsolatedElectronHistograms () {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "Booking isolated electron histograms" << endreq;

  m_h_isolatedElectronMultiplicity = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Isolated Electron Multiplicity",10,0.0,10.0);
  // or be a C++ guru and do both at once
  m_h_isolatedElectronPt = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Isolated Electron Pt",100,0.0,1000.0);
  return;
}



void StandardHistogramMaker::bookIsolatedPhotonHistograms () {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "Booking isolated photon histograms" << endreq;

  m_h_isolatedPhotonMultiplicity = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Isolated Photon Multiplicity",20,0.0,20.0);
  // or be a C++ guru and do both at once
  m_h_isolatedPhotonPt = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Isolated Photon Pt",100,0.0,1000.0);
  return;
}



void StandardHistogramMaker::bookCellHistograms () {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "Booking Cell histograms" << endreq;

  m_h_cellMultiplicity = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Cell Multiplicity",200,0.0,4000.0);
  // or be a C++ guru and do both at once
  m_h_cellPt = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Cell Pt",100,0.0,1000.0);
  return;
}



void StandardHistogramMaker::bookClusterHistograms () {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "Booking Cluster histograms" << endreq;

  m_h_clusterMultiplicity = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Cluster Multiplicity",10,0.0,10.0);
  // or be a C++ guru and do both at once
  m_h_clusterPt = 
    histoSvc()->book("/stat/simple/",++m_nHist,"Cluster Pt",100,0.0,1000.0);
  return;
}





void StandardHistogramMaker::fillElectronHistograms() {

  MsgStream log( messageService(), name() ) ;
  log << MSG::DEBUG << "Fillin Electron Histograms " << endreq;
  std::vector<ReconstructedParticle*> amyElectrons;
  if(!m_tesIO->
     copy<ReconstructedParticleCollection >(amyElectrons, m_electronLocation)){
    log << MSG::DEBUG << "no electrons found in the TES. This is probably normal" << endreq; 
    // fill any histograms that have meaning in the absence of electrons
    m_h_electronMultiplicity->fill(0.0,1.0);
  } else {
    // fill all electron histograms
    m_h_electronMultiplicity->fill(float(amyElectrons.size()),1.0);
    std::vector<ReconstructedParticle*>::const_iterator iter  = amyElectrons.begin();
    for (; iter != amyElectrons.end(); ++iter) {
      m_h_electronPt->fill(float( (*iter)->pT()), 1.0);
    }
  }
  return;
}




void StandardHistogramMaker::fillPhotonHistograms() {

  MsgStream log( messageService(), name() ) ;
  log << MSG::DEBUG << "Fillins Photon Histograms " << endreq;
  std::vector<ReconstructedParticle*> amyPhotons;
  if(!m_tesIO->copy<ReconstructedParticleCollection >
     (amyPhotons, m_photonLocation)){
    log << MSG::DEBUG << "no photons found in the TES. This is probably normal" << endreq; 
    // fill any histograms that have meaning in the absence of photons
    m_h_photonMultiplicity->fill(0.0,1.0);
  } else {
    // fill all photon histograms
    m_h_photonMultiplicity->fill(float(amyPhotons.size()),1.0);
    std::vector<ReconstructedParticle*>::const_iterator iter  = amyPhotons.begin();
    for (; iter != amyPhotons.end(); ++iter) {
      m_h_photonPt->fill(float( (*iter)->pT()), 1.0);
    }
  }
  return;
}

void StandardHistogramMaker::fillIsolatedElectronHistograms() {

  MsgStream log( messageService(), name() ) ;
  log << MSG::DEBUG << "Filling Isolated Electron Histograms " << endreq;
  std::vector<ReconstructedParticle*> amyElectrons;
  if(!m_tesIO->copy<ReconstructedParticleCollection >
     (amyElectrons, m_isolatedElectronLocation)){
    log << MSG::DEBUG << "no electrons found in the TES. This is probably normal" << endreq; 
    // fill any histograms that have meaning in the absence of electrons
    m_h_isolatedElectronMultiplicity->fill(0.0,1.0);
  } else {
    // fill all electron histograms
    m_h_isolatedElectronMultiplicity->fill(float(amyElectrons.size()),1.0);
    std::vector<ReconstructedParticle*>::const_iterator iter  = amyElectrons.begin();
    for (; iter != amyElectrons.end(); ++iter) {
      m_h_isolatedElectronPt->fill(float( (*iter)->pT()), 1.0);
    }
  }
  return;
}



void StandardHistogramMaker::fillIsolatedPhotonHistograms() {

  MsgStream log( messageService(), name() ) ;
  log << MSG::DEBUG << "Filling isolated Photon Histograms " << endreq;
  std::vector<ReconstructedParticle*> amyPhotons;
  if(!m_tesIO->copy<ReconstructedParticleCollection >
     (amyPhotons, m_isolatedPhotonLocation)){
    log << MSG::DEBUG << "no photons found in the TES. This is probably normal" << endreq; 
    // fill any histograms that have meaning in the absence of photons
    m_h_isolatedPhotonMultiplicity->fill(0.0,1.0);
  } else {
    // fill all photon histograms
    m_h_isolatedPhotonMultiplicity->fill(float(amyPhotons.size()),1.0);
    std::vector<ReconstructedParticle*>::const_iterator iter  = amyPhotons.begin();
    for (; iter != amyPhotons.end(); ++iter) {
      m_h_isolatedPhotonPt->fill(float( (*iter)->pT()), 1.0);
    }
  }
  return;
}



void StandardHistogramMaker::fillCellHistograms() {

  MsgStream log( messageService(), name() ) ;
  log << MSG::DEBUG << "Filling Cell Histograms, currently empty " << endreq;

  return;
}
 

void StandardHistogramMaker::fillClusterHistograms() {

  MsgStream log( messageService(), name() ) ;
  log << MSG::DEBUG << "Filling Cluster Histograms, currently empty " << endreq;

  return;
}




} // end of namespace bracket

---