sndScifiHit Class Reference

#include <sndScifiHit.h>

Inheritance diagram for sndScifiHit:

Collaboration diagram for sndScifiHit:

Public Member Functions
	sndScifiHit ()
	sndScifiHit (const sndScifiHit &hit)=default
sndScifiHit &	operator= (const sndScifiHit &hit)=default
	sndScifiHit (Int_t detID)
	sndScifiHit (int detID, std::vector< ScifiPoint * >, std::vector< Float_t >)
virtual	~sndScifiHit ()
void	Print ()
Float_t	GetEnergy ()
void	setInvalid ()
bool	isValid () const
Int_t	GetStation ()
bool	isVertical ()
Int_t	GetMat ()
Int_t	GetSiPM ()
Int_t	GetSiPMChan ()
Int_t	GetChannelID ()
Public Member Functions inherited from SndlhcHit
	SndlhcHit ()
	SndlhcHit (Int_t detID, Int_t nSiPMs=1, Int_t nSides=0)
virtual	~SndlhcHit ()
Int_t	GetDetectorID () const
Float_t	GetSignal (Int_t nChannel=0)
Float_t	GetTime (Int_t nChannel=0)
Int_t	GetnSiPMs () const
Int_t	GetnSides () const
void	SetDigi (Float_t s, Float_t t, Int_t i=0)
void	SetDetectorID (Int_t detID)
void	SetDaqID (Int_t i, Int_t k, Int_t board_id, Int_t tofpet_id, Int_t tofpet_channel)
Int_t	GetBoardID (Int_t i)
Int_t	GetTofpetID (Int_t i)
Int_t	Getchannel (Int_t i)
Int_t	GetRawHitIndex (Int_t i=0)
virtual void	Print (const Option_t *opt="") const
virtual void	GetPosition (TVector3 L, TVector3 R) const
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Private Member Functions
Float_t	ly_loss (Float_t distance)
Float_t	sipm_saturation (Float_t ly, Float_t nphe_max)
Float_t	npix_to_qdc (Float_t npix)
	ClassDef (sndScifiHit, 4)

Private Attributes
Float_t	flag
	flag

Additional Inherited Members
Protected Member Functions inherited from SndlhcHit
	ClassDef (SndlhcHit, 2)
	encodes rawhitindex*100000+(board_id * 1000) + (tofpet_id * 100) + tofpet_channel
Protected Attributes inherited from SndlhcHit
Int_t	fDetectorID
	Detector unique identifier.
Int_t	nSiPMs
Int_t	nSides
	number of SiPMs per side
Float_t	signals [16]
	number of sides
Float_t	times [16]
	SiPM signal.
Int_t	fDaqID [16]
	SiPM time.

Detailed Description

Definition at line 9 of file sndScifiHit.h.

Constructor & Destructor Documentation

sndScifiHit() [1/4]

sndScifiHit::sndScifiHit

(

)

Default constructor

Definition at line 18 of file sndScifiHit.cxx.

  : SndlhcHit()
{
 flag = true;
}

sndScifiHit() [2/4]

sndScifiHit::sndScifiHit ( const sndScifiHit &  hit )

default

Copy constructor

sndScifiHit() [3/4]

sndScifiHit::sndScifiHit ( Int_t  detID )

explicit

Definition at line 24 of file sndScifiHit.cxx.

  : SndlhcHit(detID)
{
 flag = true;
}

sndScifiHit() [4/4]

sndScifiHit::sndScifiHit	(	int	detID,
		std::vector< ScifiPoint * >	V,
		std::vector< Float_t >	W
	)

Definition at line 30 of file sndScifiHit.cxx.

{
     Scifi* ScifiDet = dynamic_cast<Scifi*> (gROOT->GetListOfGlobals()->FindObject("Scifi") );
     Float_t nphe_min = ScifiDet->GetConfParF("Scifi/nphe_min");
     Float_t nphe_max = ScifiDet->GetConfParF("Scifi/nphe_max");
     Float_t timeResol = ScifiDet->GetConfParF("Scifi/timeResol");
     Float_t signalSpeed = ScifiDet->GetConfParF("Scifi/signalSpeed");
     
     nSides   = 1;
     nSiPMs   = 1;
     for (unsigned int j=0; j<16; ++j){
        signals[j] = -1;
        times[j]    =-1;
     }
 
     fDetectorID  = SiPMChan;
     Float_t ly_total = 0;
     Float_t earliestToA   = 1E20;
     for( int i = 0; i <V.size();i++) {
        
        Double_t signal = V[i]->GetEnergyLoss()*W[i];
 
    // Find distances from MCPoint centre to ends of fibre
        TVector3 a, b;
        TVector3 impact(V[i]->GetX(),V[i]->GetY() ,V[i]->GetZ() );
        ScifiDet->GetSiPMPosition(V[i]->GetDetectorID(),a, b);
         
    Double_t distance;
 
    bool verticalHit = int(fDetectorID/100000)%10 == 1;
 
    // Calculate distance from energy deposit to SiPM.
    // Vertical
    if (verticalHit) {
      distance = (b - impact).Mag();
      // Horizontal
    } else {
      distance = (impact - a).Mag();
    }
 
    // convert energy deposit to light yield (here Np.e. == avg. N fired pixels)
        Float_t ly = signal*1E+6*0.16; //0.16 p.e per 1 keV
        // account for the light attenuation in the fibers
        ly*= ly_loss(distance);
    ly_total+= ly;
    
    // for the timing, find earliest light to arrive at SiPM and smear with time resolution
    Float_t arrival_time = V[i]->GetTime() + distance/signalSpeed;
    if (arrival_time < earliestToA){earliestToA = arrival_time;}
    
     }
     // smear the total light yield using Poisson distribution
     ly_total= gRandom->Poisson(ly_total);
 
     // account for limited SiPM dyn. range
     Float_t Npix = sipm_saturation(ly_total,nphe_max);
     // convert Npix to QDC
     signals[0] = npix_to_qdc(Npix);
     
     if (ly_total > nphe_min){   // nominal threshold at 3.5 p.e.
            flag=true;
      }else{
            flag=false;
      }
      times[0] = gRandom->Gaus(earliestToA, timeResol);
 
     LOG(DEBUG) << "signal created";
}

~sndScifiHit()

sndScifiHit::~sndScifiHit ( )

virtual

Destructor

Definition at line 100 of file sndScifiHit.cxx.

{ }

Member Function Documentation

ClassDef()

sndScifiHit::ClassDef ( sndScifiHit  , 4    )

private

GetChannelID()

Int_t sndScifiHit::GetChannelID ( )

inline

Definition at line 39 of file sndScifiHit.h.

{return fDetectorID;}

GetEnergy()

Float_t sndScifiHit::GetEnergy

(

)

Definition at line 104 of file sndScifiHit.cxx.

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

GetMat()

Int_t sndScifiHit::GetMat ( )

inline

Definition at line 36 of file sndScifiHit.h.

{return ( int(fDetectorID/10000)%10);}

GetSiPM()

Int_t sndScifiHit::GetSiPM ( )

inline

Definition at line 37 of file sndScifiHit.h.

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

GetSiPMChan()

Int_t sndScifiHit::GetSiPMChan ( )

inline

Definition at line 38 of file sndScifiHit.h.

{return ( fDetectorID%1000);}

GetStation()

Int_t sndScifiHit::GetStation ( )

inline

Definition at line 31 of file sndScifiHit.h.

{return floor(fDetectorID/1000000);}

isValid()

bool sndScifiHit::isValid ( ) const

inline

Definition at line 30 of file sndScifiHit.h.

{return flag;}

isVertical()

bool sndScifiHit::isVertical ( )

inline

Definition at line 32 of file sndScifiHit.h.

                     {
        if ( int(fDetectorID/100000)%10 == 1){return true;}
        else{return (false);}
    }

ly_loss()

Float_t sndScifiHit::ly_loss ( Float_t  distance )

private

Definition at line 110 of file sndScifiHit.cxx.

                                            {
//  It returns the light yield attenuation depending on the distance to SiPM
    return TMath::Exp(-(distance-ly_loss_params[0])/ly_loss_params[1]);
}

npix_to_qdc()

Float_t sndScifiHit::npix_to_qdc ( Float_t  npix )

private

Definition at line 121 of file sndScifiHit.cxx.

                                            {
//  It returns QDC per channel after Gaussian smearing of the parameters
        Float_t A = gRandom->Gaus(npix_to_qdc_params[0], npix_to_qdc_params[2]);
        Float_t B = gRandom->Gaus(npix_to_qdc_params[1], npix_to_qdc_params[3]);
        return A*npix + B;
}

operator=()

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

default

Print()

void sndScifiHit::Print

(

)

Output to screen

Definition at line 129 of file sndScifiHit.cxx.

{
  std::cout << "-I- sndScifiHit: Scifi hit " << " in station " << GetStation();
  if ( isVertical()) {     std::cout << " vertical plane ";}
  else {     std::cout << " horizontal plane ";}
  std::cout << "SiPM nr "<<GetSiPM()<< " SiPM channel "<<GetSiPMChan()<<std::endl;
}

setInvalid()

void sndScifiHit::setInvalid ( )

inline

Definition at line 29 of file sndScifiHit.h.

{flag = false;}

sipm_saturation()

Float_t sndScifiHit::sipm_saturation ( Float_t  ly, Float_t  nphe_max  )

private

Definition at line 115 of file sndScifiHit.cxx.

                                                                {
//  It returns the number of fired pixels per channel
        Float_t factor = 1 - TMath::Exp(-ly/nphe_max);
    return nphe_max*factor;
}

Member Data Documentation

flag

Float_t sndScifiHit::flag

private

flag

Definition at line 48 of file sndScifiHit.h.

---

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

• shipLHC/sndScifiHit.h
• shipLHC/sndScifiHit.cxx