---

CellMaker.h

Go to the documentation of this file.

// Cell Maker.h
//
// HTP tidied up and put into namespace

#ifndef _Atlfast_CellMaker_H
#define _Atlfast_CellMaker_H

// Normal STL and physical vectors

#include "CLHEP/Vector/LorentzVector.h"
#include "CLHEP/Vector/ThreeVector.h"
#include <vector>
#include <string>

// Gaudi
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/DataObject.h"

#include "HepMC/GenEvent.h"
#include "HepMC/GenParticle.h"

//Atlfast
#include "AtlfastCode/Cell.h"
#include "AtlfastCode/Calorimeter.h"
#include "AtlfastCode/MagField.h"
#include "AtlfastCode/HepMC_helper.h"
#include "AtlfastCode/TesIO.h"
#include "AtlfastCode/CommonData.h"

//Transient event store paths

#define   DEFAULT_inputLocation      "/Event/McEventCollection"
#define   DEFAULT_commonDataLocation "/Event/AtlfastCommonData"
#define   DEFAULT_outputLocation     "/Event/AtlfastCells" 

//Steering parameters needed to set up a calorimeter object
// it is clearly very ugly to have to read all of this in 
// for CellMaker - we should instead hand CellMaker the right


// sort of Calorimeter object. But right now we do it this way
// because we don't yet have Gaudi tools.
//
#define   DEFAULT_etaCoverage 5.0
#define   DEFAULT_minETCell 0.0
// Starting granularities. Will be adjusted to give integral numbers of bins
#define   DEFAULT_granBarrelEta  0.1
#define   DEFAULT_granBarrelPhi  0.1
#define   DEFAULT_granForwardEta 0.2
#define   DEFAULT_granForwardPhi 0.2


namespace Atlfast {

  class CellMaker : public Algorithm {

  public:
    //Gaudi style constructor and execution methods
    CellMaker(const std::string& name, ISvcLocator* pSvcLocator);
    ~CellMaker();
    
    StatusCode          initialize();
    StatusCode          execute();
    //    StatusCode          finalize() {return StatusCode::SUCCESS;}
    StatusCode          finalize();

  private:
    Calorimeter*                         m_calorimeter;  
    MagField*                            m_magField;     
    ICellSelector*                       m_cellSelector;
    HepMC_helper::IMCselector*          m_mcSelector;

    // parameters needed to construct a calorimeter object...
    double    m_etaCoverage;      
 // entire eta coverage
    double    m_minETCell;         // minimum Et of cell to write out
    double    m_barrelForwardEta;  // boundary between forward and barrel
    bool      m_fieldOn;           // field bending calculted if true


    //...............Cell Granularities for...............................
    double    m_granBarrelEta;     // barrel eta
    double    m_granBarrelPhi;     // barrel phi
    double    m_granForwardEta;    // forward eta
    double    m_granForwardPhi;    // forward phi



    //Paths in the transient event store to get/put things
    //Currently set in member variables to be overwritten
    //by job options service
    std::string m_inputLocation;
    std::string m_commonDataLocation;
    std::string m_outputLocation;

    TesIO* m_tesIO;
    bool m_notInit;
  };

} // end of namespace bracket
#endif

---