calorimeter_step_hit_processor.h

Go to the documentation of this file.

/* Author(s) :   Francois Mauger <mauger@lpccaen.in2p3.fr>
 * Creation date: 2010-05-26
 * Last modified: 2014-04-23
 *
 * License:
 *
 * Description:
 *
 *  A basic processor of simulated step hits in some plastic scintillator
 *  volume.
 *
 *  This algorithm clusterizes step calorimeter hits (from a
 *  simulation program) and build new hits that cumulate the
 *  physics of original MC step hits (particle type, geometry ID,
 *  energy deposit, timing...)
 *
 *  CAUTION: this processor can manage only one geometry mapping category
 *  so care should be taken to attach only one geometry model (with the proper
 *  mapping category) to this processor.
 *
 *  QUESTION:
 *  Should delta rays energy deposit along the track of an alpha particle
 *  be clusterized within the parent alpha clusterized hit (quenching effect) ?
 *
 */

#ifndef MCTOOLS_CALORIMETER_STEP_HIT_PROCESSOR_H
#define MCTOOLS_CALORIMETER_STEP_HIT_PROCESSOR_H 1

// Standard library:
#include <string>

// Third party:
// - Boost:
#include <boost/cstdint.hpp>
// - Bayeux/datatools:
#include <datatools/ioutils.h>
// - Bayeux/geomtools:
#include <geomtools/i_shape_3d.h>
#include <geomtools/smart_id_locator.h>
#include <geomtools/box.h>
#include <geomtools/sphere.h>
#include <geomtools/placement.h>

// This project:
#include <mctools/base_step_hit_processor.h>

namespace mctools {

  class calorimeter_step_hit_processor
    : public base_step_hit_processor
  {
  public:

    virtual bool locate_calorimeter_block(const geomtools::vector_3d & position_,
                                          geomtools::geom_id & gid_) const;

    const std::string & get_mapping_category() const;

    void set_mapping_category(const std::string & sc_);

    const std::vector<int> & get_mapping_category_any_addresses() const;

    calorimeter_step_hit_processor();

    virtual ~calorimeter_step_hit_processor();

    virtual void initialize(const ::datatools::properties & config_,
                            ::datatools::service_manager & service_mgr_);

    virtual void process(const ::mctools::base_step_hit_processor::step_hit_ptr_collection_type & base_step_hits_,
                         ::mctools::simulated_data::hit_handle_collection_type & handle_hits_);

    virtual void process(const ::mctools::base_step_hit_processor::step_hit_ptr_collection_type & base_step_hits_,
                         ::mctools::simulated_data::hit_collection_type & plain_hits_);

    bool match_scintillation_hit(const base_step_hit & scintillation_hit_,
                                 const base_step_hit & step_hit_) const;

    void merge_scintillation_hits(simulated_data::hit_handle_collection_type & scintillation_hits_);

    void tree_dump(std::ostream & out_ = std::clog,
                   const std::string & title_ = "" ,
                   const std::string & indent_ = "",
                   bool inherit_ = false) const;

    static void init_ocd(datatools::object_configuration_description &);

  protected:

    void _init(const ::datatools::properties & config_,
               ::datatools::service_manager & service_mgr_);

    void _process(const base_step_hit_processor::step_hit_ptr_collection_type & base_step_hits_,
                  simulated_data::hit_handle_collection_type * calo_hits_,
                  simulated_data::hit_collection_type * plain_calo_hits_);

  private:

    double _scintillation_cluster_time_range_;

    double _scintillation_cluster_space_range_;

    std::string _mapping_category_;

    std::vector<int> _mapping_category_any_addresses_;

    // internals:
    const geomtools::mapping * _mapping_; 

    const geomtools::id_mgr::categories_by_name_col_type * _categories_; 

    uint32_t _calo_block_type_; 

    geomtools::smart_id_locator _calo_block_locator_;

    bool _alpha_quenching_; 

    // Registration macro :
    MCTOOLS_STEP_HIT_PROCESSOR_REGISTRATION_INTERFACE(calorimeter_step_hit_processor)

  };

} // end of namespace mctools

#include <datatools/ocd_macros.h>
DOCD_CLASS_DECLARATION(mctools::calorimeter_step_hit_processor)

#endif // MCTOOLS_CALORIMETER_STEP_HIT_PROCESSOR_H

// Local Variables: --
// mode: c++ --
// c-file-style: "gnu" --
// tab-width: 2 --
// End: --