Simulation/BOOST/EmcSim/include/EmcSim/BesEmcDigitizer.hh

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//      BOOST --- BESIII Object_Oreiented Simulation Tool                    //
//---------------------------------------------------------------------------//
// Descpirtion: EMC detector
// Author: Hemiao
// Created: Sep, 2004
// Comment:
//---------------------------------------------------------------------------//
// $Id:BesEmcDigitizer.hh
 
#ifndef BesEmcDigitizer_h
#define BesEmcDigitizer_h 1
 
#include "BesEmcDigi.hh"
#include "BesEmcHit.hh"
 
#include "G4VDigitizerModule.hh"
#include "GaudiKernel/NTuple.h"
#include "globals.hh"
#include <vector>
 
#include "EmcCalibConstSvc/IEmcCalibConstSvc.h"
#include "G4Svc/IG4Svc.h"
 
using namespace std;
 
// class BesEmcHit;
// class G4Svc;
// class IEmcCalibConstSvc;
 
class CrystalSingle {
  friend class BesEmcDigitizer;
 
public:
  CrystalSingle() { hitIndexes = new vector<G4int>; }
  CrystalSingle( const CrystalSingle& right ) {
    partId     = right.partId;
    nTheta     = right.nTheta;
    nPhi       = right.nPhi;
    edep       = right.edep;
    hitIndexes = right.hitIndexes;
  }
  const CrystalSingle& operator=( const CrystalSingle& right ) {
    partId     = right.partId;
    nTheta     = right.nTheta;
    nPhi       = right.nPhi;
    edep       = right.edep;
    hitIndexes = right.hitIndexes;
    return *this;
  }
  ~CrystalSingle() { delete hitIndexes; }
 
  G4int          GetPartId() { return partId; }
  G4int          GetNTheta() { return nTheta; }
  G4int          GetNPhi() { return nPhi; }
  G4double       GetEdep() { return edep; }
  vector<G4int>* GetHitIndexes() { return hitIndexes; }
 
  void SetPartId( G4int id ) { partId = id; }
  void SetNTheta( G4int theta ) { nTheta = theta; }
  void SetNPhi( G4int phi ) { nPhi = phi; }
  void SetEdep( G4double e ) { edep = e; }
  void AddEdep( G4double e ) { edep += e; }
 
private:
  G4int          partId;
  G4int          nTheta;
  G4int          nPhi;
  G4double       edep;
  vector<G4int>* hitIndexes;
};
 
class BesEmcDigitizer : public G4VDigitizerModule {
public:
  BesEmcDigitizer( G4String modName );
  ~BesEmcDigitizer();
 
public:
  virtual void Digitize();
  virtual void GroupHits( BesEmcHitsCollection* );
  // Add noise to 5x5 matrix around the seed
  virtual void AddNoise5x5( G4double coherentNoise );
  // Add noise to all crystals
  virtual void AddNoiseAll( G4double coherentNoise );
 
private: // memory size needed optimizing
  void Initialize();
  //  void Smear();
 
private:
  BesEmcDigitsCollection* m_besEmcDigitsCollection;
  vector<CrystalSingle*>* m_crystalGroup;
  G4double                m_energy;
 
private:
  // G4Svc* m_G4Svc; // mrli
  IG4Svc*              m_G4Svc;
  IEmcCalibConstSvc*   m_emcCalibConstSvc;
  NTuple::Tuple*       m_tupleEmc1; // digi level
  NTuple::Item<long>   m_partId;
  NTuple::Item<long>   m_nTheta;
  NTuple::Item<long>   m_nPhi;
  NTuple::Item<double> m_eDep;
  NTuple::Item<long>   m_nHits;
  NTuple::Item<double> m_adc;
  NTuple::Item<long>   m_tdc;
 
  NTuple::Tuple*       m_tupleEmc2; // event level
  NTuple::Item<double> m_eTot;
  NTuple::Item<long>   m_nDigi;
};
#endif