dc/ddd/Pixel3DDigitizerAlgorithm_8h_source.html

#ifndef _SimTracker_SiPhase2Digitizer_Pixel3DDigitizerAlgorithm_h

#define _SimTracker_SiPhase2Digitizer_Pixel3DDigitizerAlgorithm_h


//-------------------------------------------------------------

// class Pixel3DDigitizerAlgorithm

//

// Specialization of the tracker digitizer for the 3D pixel

// sensors placed at PXB-Layer1 (and possibly Layer2), and

// at PXF-Disk1 and Disk2

//

// Authors: Jordi Duarte-Campderros (CERN/IFCA)

//          Clara Lasaosa Garcia (IFCA)

//--------------------------------------------------------------


#include "SimTracker/SiPhase2Digitizer/plugins/Phase2TrackerDigitizerAlgorithm.h"


// Data formats

#include "DataFormats/GeometryVector/interface/LocalPoint.h"


// system

#include <functional>


class Pixel3DDigitizerAlgorithm : public Phase2TrackerDigitizerAlgorithm {

public:

  Pixel3DDigitizerAlgorithm(const edm::ParameterSet& conf);

  ~Pixel3DDigitizerAlgorithm() override;


  // initialization that cannot be done in the constructor

  void init(const edm::EventSetup& es) override;


  // run the algorithm to digitize a single det

  void accumulateSimHits(const std::vector<PSimHit>::const_iterator inputBegin,

                         const std::vector<PSimHit>::const_iterator inputEnd,

                         const size_t inputBeginGlobalIndex,

                         const unsigned int tofBin,

                         const Phase2TrackerGeomDetUnit* pixdet,

                         const GlobalVector& bfield) override;


  // overload the mother (drift3d maybe)

  std::vector<DigitizerUtility::SignalPoint> drift(

      const PSimHit& hit,

      const Phase2TrackerGeomDetUnit* pixdet,

      const GlobalVector& bfield,

      const std::vector<DigitizerUtility::EnergyDepositUnit>& ionization_points,

      bool diffusion_activated);

  // New diffusion function: check implementation

  std::vector<DigitizerUtility::EnergyDepositUnit> diffusion(const LocalPoint& pos,

                                                             const float& ncarriers,

                                                             const std::function<LocalVector(float, float)>& u_drift,

                                                             const std::pair<float, float> pitches,

                                                             const float& thickness);

  // Specific for 3D-pixel

  void induce_signal(const PSimHit& hit,

                     const size_t hitIndex,

                     const unsigned int tofBin,

                     const Phase2TrackerGeomDetUnit* pixdet,

                     const std::vector<DigitizerUtility::SignalPoint>& collection_points);


private:

  // Raidus of Column np and ohmic

  float _np_column_radius;

  float _ohm_column_radius;


  // Check if a carrier is inside the column: The point should

  // be described in the pixel cell frame

  const bool _is_inside_n_column(const LocalPoint& p) const;

  const bool _is_inside_ohmic_column(const LocalPoint& p, const std::pair<float, float>& pitch) const;

};

#endif