MTCDetHit Class Reference

#include <MTCDetHit.h>

Inheritance diagram for MTCDetHit:

Collaboration diagram for MTCDetHit:

Public Member Functions
	MTCDetHit ()
	MTCDetHit (const MTCDetHit &hit)=default
MTCDetHit &	operator= (const MTCDetHit &hit)=default
	MTCDetHit (int detID, const std::vector< MTCDetPoint * > &, const std::vector< Float_t > &)
	~MTCDetHit () override=default
void	Print () const
Float_t	GetSignal () const
Float_t	GetTime () const
Int_t	GetChannelID () const
Int_t	GetLayer () const
Int_t	GetStationType () const
Int_t	GetSiPM () const
Int_t	GetSiPMChan () const
TVector3	GetChannelCoordinates () const
Float_t	GetEnergy () const
void	setInvalid ()
bool	isValid () const
Public Member Functions inherited from SHiP::DetectorHit
	DetectorHit ()
	DetectorHit (Int_t detID, Float_t digi)
	~DetectorHit () override=default
Double_t	GetDigi () const
Int_t	GetDetectorID () const
void	SetDigi (Float_t d)
void	SetDetectorID (Int_t detID)
void	Print (const Option_t *opt="") const override

Private Member Functions
	ClassDef (MTCDetHit, 6)

Private Attributes
Float_t	signals = 0
Float_t	time
Float_t	Xch = 0.0
Float_t	Ych = 0.0
Float_t	Zch = 0.0
Bool_t	flag {true}
	flag

Additional Inherited Members
Protected Member Functions inherited from SHiP::DetectorHit
	ClassDefOverride (SHiP::DetectorHit, 1)
Protected Attributes inherited from SHiP::DetectorHit
Float_t	fdigi
	digitized detector hit
Int_t	fDetectorID
	Detector unique identifier.

Detailed Description

Definition at line 14 of file MTCDetHit.h.

Constructor & Destructor Documentation

MTCDetHit() [1/3]

MTCDetHit::MTCDetHit

(

)

Default constructor

Definition at line 52 of file MTCDetHit.cxx.

: SHiP::DetectorHit() {}

MTCDetHit() [2/3]

MTCDetHit::MTCDetHit ( const MTCDetHit &  hit )

default

Copy constructor

MTCDetHit() [3/3]

MTCDetHit::MTCDetHit	(	int	detID,
		const std::vector< MTCDetPoint * > &	points,
		const std::vector< Float_t > &	weights
	)

Definition at line 56 of file MTCDetHit.cxx.

                                                        {
  auto* MTCDet =
      dynamic_cast<MTCDetector*>(gROOT->GetListOfGlobals()->FindObject("MTC"));
  if (!MTCDet) {
    LOG(fatal) << "MTCDetHit: MTCDetector not found";
    return;
  }
 
  fDetectorID = SiPMChan;  // Set the detector ID
  // Determine plane type once
  const int plane_type = GetStationType();
 
  Float_t total_light_yield = 0.0f;
  Float_t earliest_to_A = std::numeric_limits<Float_t>::max();
  Float_t earliest_to_B = std::numeric_limits<Float_t>::max();
  const size_t n = points.size();
  Float_t signal_sum = 0.0f;
  bool hit_flag = false;
 
  // Separate handling for scintillating mat (plane_type == 2)
  if (plane_type == 2) {
    Float_t x_temp = 0.0, y_temp = 0.0, z_temp = 0.0;
    for (auto* pt : points) {
      signal_sum += pt->GetEnergyLoss();
      // Track earliest arrival time
      Float_t arrival = pt->GetTime();
      earliest_to_B = std::min(earliest_to_B, arrival);
      x_temp += pt->GetX();
      y_temp += pt->GetY();
      z_temp += pt->GetZ();
    }
    time = gRandom->Gaus(earliest_to_B, time_res);
    // for scintillating tiles set simulated coordinates so far as the realistic
    // geometry is not yet done.
    Xch = x_temp / n;
    Ych = y_temp / n;
    Zch = z_temp / n;
    signals = signal_sum;
    return;
  }
 
  // Fiber hit processing
  total_light_yield = 0.0f;
 
  TVector3 sipmA, sipmB;
  MTCDet->GetSiPMPosition(SiPMChan, sipmA, sipmB);
  // Define only X and Z coordinates for Sci-Fi.
  Xch = sipmB.X();
  Zch = sipmB.Z();
 
  for (size_t i = 0; i < n; ++i) {
    auto* pt = points[i];
    Float_t energy = pt->GetEnergyLoss();
    Float_t weight = weights[i];
    Float_t signal = energy * weight;
 
    // Distance from deposit to SiPM
    TVector3 impact(pt->GetX(), pt->GetY(), pt->GetZ());
    const Float_t distance = (sipmB - impact).Mag();
 
    // Light yield before attenuation
    Float_t light_yield = signal * 1e6f * 0.16f;
    light_yield *= light_attenuation(distance);
    total_light_yield += light_yield;
 
    // Track earliest arrival time
    Float_t arrival = pt->GetTime() + distance * inv_signal_speed;
    earliest_to_A = std::min(earliest_to_A, arrival);
  }
 
  // Apply statistical smearing and saturation
  const Int_t smeared_light_yield = gRandom->Poisson(total_light_yield);
  const Float_t n_pixels = sipm_saturation(smeared_light_yield);
  signals = n_pixels_to_qdc(n_pixels);
 
  // Final hit decision and time
  hit_flag = (smeared_light_yield > n_photons_min);
  flag = hit_flag;
  time = gRandom->Gaus(earliest_to_A, time_res);
}

~MTCDetHit()

MTCDetHit::~MTCDetHit ( )

overridedefault

Destructor

Member Function Documentation

ClassDef()

MTCDetHit::ClassDef ( MTCDetHit  , 6    )

private

GetChannelCoordinates()

TVector3 MTCDetHit::GetChannelCoordinates ( ) const

inline

Definition at line 58 of file MTCDetHit.h.

{ return TVector3(Xch, Ych, Zch); }

GetChannelID()

Int_t MTCDetHit::GetChannelID ( ) const

inline

Definition at line 49 of file MTCDetHit.h.

{ return fDetectorID; }

GetEnergy()

Float_t MTCDetHit::GetEnergy

(

)

const

Definition at line 139 of file MTCDetHit.cxx.

                                   {
  // to be calculated from digis and calibration constants, missing!
  return signals;
}

GetLayer()

Int_t MTCDetHit::GetLayer ( ) const

inline

Definition at line 50 of file MTCDetHit.h.

                         {
    return static_cast<int>(fDetectorID / 1000000) % 100;
  }

GetSignal()

Float_t MTCDetHit::GetSignal ( ) const

inline

Definition at line 31 of file MTCDetHit.h.

{ return signals; };

GetSiPM()

Int_t MTCDetHit::GetSiPM ( ) const

inline

Definition at line 56 of file MTCDetHit.h.

{ return (static_cast<int>(fDetectorID / 1000) % 10); }

GetSiPMChan()

Int_t MTCDetHit::GetSiPMChan ( ) const

inline

Definition at line 57 of file MTCDetHit.h.

{ return (fDetectorID % 1000); }

GetStationType()

Int_t MTCDetHit::GetStationType ( ) const

inline

Definition at line 53 of file MTCDetHit.h.

                               {
    return static_cast<int>(fDetectorID / 100000) % 10;
  }

GetTime()

Float_t MTCDetHit::GetTime ( ) const

inline

Definition at line 32 of file MTCDetHit.h.

{ return time; };

isValid()

bool MTCDetHit::isValid ( ) const

inline

Definition at line 61 of file MTCDetHit.h.

{ return flag; }

operator=()

MTCDetHit & MTCDetHit::operator= ( const MTCDetHit &  hit )

default

Print()

void MTCDetHit::Print

(

)

const

Definition at line 145 of file MTCDetHit.cxx.

                            {
  std::cout << Form(
      "MTCDetHit: Detector ID %d, Layer %d, Station Type %d, SiPM %d, Channel "
      "%d, Signal %.2f, Time %.3f \n",
      fDetectorID, GetLayer(), GetStationType(), GetSiPM(), GetSiPMChan(),
      signals, time);
}

setInvalid()

void MTCDetHit::setInvalid ( )

inline

Definition at line 60 of file MTCDetHit.h.

{ flag = false; }

Member Data Documentation

flag

Bool_t MTCDetHit::flag {true}

private

flag

Definition at line 71 of file MTCDetHit.h.

signals

Float_t MTCDetHit::signals = 0

private

Definition at line 68 of file MTCDetHit.h.

time

Float_t MTCDetHit::time

private

Definition at line 69 of file MTCDetHit.h.

Xch

Float_t MTCDetHit::Xch = 0.0

private

Definition at line 70 of file MTCDetHit.h.

Ych

Float_t MTCDetHit::Ych = 0.0

private

Definition at line 70 of file MTCDetHit.h.

Zch

Float_t MTCDetHit::Zch = 0.0

private

Definition at line 70 of file MTCDetHit.h.

---

The documentation for this class was generated from the following files:

• SND/MTC/MTCDetHit.h
• SND/MTC/MTCDetHit.cxx