Classes
class	ScifiCluster
	Class of SciFi Tracker ######################################################. More...

Functions
	cm_to_channel (locpos, sipm_map=sipm_map, gaps_map=gaps_map, pitch=pitch, charr=charr, reverse=False, ch_max_num=ch_max_num, n_solid_sipms=n_solid_sipms)
	FUNCTIONS FOR SCIFICLASS ##################################################.

	channel_to_cm (channelpos, sipm_map=sipm_map, reverse=False, pitch=pitch)

	GetMatType (DetID)

	GetMatNum (DetID)

	GetMatLength ()

	GetMatQty ()

	GetSiPMs ()

	GetStationNum (DetID)

	global_to_local (DetID, globalpos)
	These functions assume that the detector is aligned with the beam ######.

	local_to_global (DetID, localpos)

	ly_loss_mean (distance, params=ly_loss_params)

	ly_attenuation (distance)

	cluster_width_mean (distance, params=cluster_width_mean_params)

	cluster_width_sigma (distance, params=cluster_width_sigma_params)

	cluster_width_random (distance, ly, percent=random_width_percent, cluster_width_min=cluster_width_min, cluster_width_max=cluster_width_max)

	approx_function (var, approx_data)

	edep_to_ly (energy, CDF_integral=CDF_integral, energy_range=energy_range, ly_linear_params=ly_linear_params, k_cdfs_corr=k_cdfs_corr, sigma_in_percent=sigma_in_percent, ly_CDF_params=ly_CDF_params, CDF_ly_params=CDF_ly_params, ly_CDF_landau_params=ly_CDF_landau_params)

	cluster_generator (amplitude, width, wmp, cluster_width_max=cluster_width_max, chpos_min=chpos_min, chpos_max=chpos_max)

	is_realistic (cluster, width)

	create_cluster (amplitude, width, wmp)

	fcheck_wall (cluster, wmp)

	weigthed_mean_pos (cluster)

Variables
int	sipm_hor_pos = +1
	GLOBAL VALUES AND PARAMETERS ############.

int	sipm_vert_pos = -1

bool	is_sipm_hor_reversed = True

bool	is_sipm_vert_reversed = True

float	hw = 19.528

float	pitch = 0.025

float	charr = 1.6

float	edge = 0.016

float	gap = 0.022

float	array = gap + 2.*charr

float	biggap = 0.038

int	ch_max_num = 128*12 - 1

int	n_solid_sipms = 12

tuple	sipm_map

dict	gaps_map

float	firstmatgap_startch = 521.5

float	firstmatgap_endch = 523.5

float	secmatgap_startch = 1046.5

float	secmatgap_endch = 1048.5

float	chpos_min = -0.5

float	chpos_max = 1535.5

int	cluster_width_min = 1

int	cluster_width_max = 10

float	ly_loss_params = 20.78, -0.26, 7.89, -3.06

float	ly_loss_sigma_params = 6.8482, -0.5757, 3.4458

float	cluster_width_mean_params = 0.6661, -0.006955, 2.163

float	cluster_width_sigma_params = 1.103, -0.0005751

float	random_width_percent = 0.01

float	ly_min = 4.5

int	ly_max = 104

float	energy_range = 0.18, 0.477

float	CDF_integral = 0.0185640424

float	ly_linear_params = 332.882, -58.7085

float	k_cdfs_corr = 0.993076766938

float	sigma_in_percent = 0.01

int	sigma_from_width = 1 / 4.

dict	ly_CDF_params

dict	ly_CDF_landau_params

dict	CDF_ly_params

	snd_geo = ConfigRegistry.loadpy("$SNDSW_ROOT/geometry/shipLHC_geom_config.py")
	LOAD THE SCIFI GEOMETRY ############.

	scifimat_width = snd_geo.Scifi.scifimat_width

	scifimat_length = snd_geo.Scifi.scifimat_length

	scifimat_gap = snd_geo.Scifi.scifimat_gap

	scifi_x = snd_geo.Scifi.xdim

	scifi_y = snd_geo.Scifi.ydim

	scifimat_z = snd_geo.Scifi.scifimat_z

	nmats = snd_geo.Scifi.nmats

	n_planes = snd_geo.Scifi.nscifi

	n_sipms = snd_geo.Scifi.nsipms

float	xshift = -27.8
	Need to automatize this part so it works for every geometry #########.

float	yshift = +35.3

Function Documentation

◆ approx_function()

ScifiCluster.approx_function	(	var,
		approx_data
	)

    This universal function substitutes the parameters to the function. 
    The parameters and the function are in the dictionary

Definition at line 466 of file ScifiCluster.py.

def approx_function(var, approx_data):
 
    """
    This universal function substitutes the parameters to the function. 
    The parameters and the function are in the dictionary
    """
 
    for (left, right), (params, func) in approx_data.items():
        if left <= var <= right:
            return func(var, *params)
    return False
 

◆ channel_to_cm()

ScifiCluster.channel_to_cm	(	channelpos,
		sipm_map = `sipm_map`,
		reverse = `False`,
		pitch = `pitch`
	)

    It converts a particle position measured channels to a position measured 
    in cm. The SiPM map is used. The position is in the scifi module frame.

Definition at line 304 of file ScifiCluster.py.

def channel_to_cm(channelpos, sipm_map=sipm_map, reverse=False, pitch=pitch):
 
    """
    It converts a particle position measured channels to a position measured 
    in cm. The SiPM map is used. The position is in the scifi module frame.
    """
 
    if reverse is True:
        for left, _, ch_range in sipm_map:
            if not ch_range is False:
                ch_start, ch_end = ch_range
                if ch_start <= channelpos <= ch_end:
                    return -(left + (channelpos - ch_start) * pitch)
    if reverse is False:
        for left, _, ch_range in sipm_map:
            if not ch_range is False:
                ch_start, ch_end = ch_range
                if ch_start <= channelpos <= ch_end:
                    return left + (channelpos - ch_start) * pitch
 

◆ cluster_generator()

ScifiCluster.cluster_generator	(	amplitude,
		width,
		wmp,
		cluster_width_max = `cluster_width_max`,
		chpos_min = `chpos_min`,
		chpos_max = `chpos_max`
	)

    It generates an event cluster with given weighted mean position in channels, width and 
    amplitude. 

    If right side of the cluster can be out of range, the maximum of the right side will be 
    right channel.

    At first an array [0, 0, 0, ... ] is generated which corresponds to the channels. 
    Next the cluster generated in the array. 
    Final array will be like [0, 0, ..., 1, 2, 5, 1, 0, ...], 
    [0, 17, 0, ...] or etc.

Definition at line 513 of file ScifiCluster.py.

    chpos_min=chpos_min, chpos_max=chpos_max):
 
    """
    It generates an event cluster with given weighted mean position in channels, width and 
    amplitude. 
 
    If right side of the cluster can be out of range, the maximum of the right side will be 
    right channel.
 
    At first an array [0, 0, 0, ... ] is generated which corresponds to the channels. 
    Next the cluster generated in the array. 
    Final array will be like [0, 0, ..., 1, 2, 5, 1, 0, ...], 
    [0, 17, 0, ...] or etc.
    """
 
    if int(ceil(wmp + 0.5 + cluster_width_max)) < chpos_max:
        max_fired_ch = int(ceil(wmp + 0.5 + cluster_width_max))
    else:
        max_fired_ch = int(ceil(chpos_max))
    cluster = [0 for _ in range(max_fired_ch)]     
    mean = wmp
    if width == 1:
        if wmp == chpos_min:    
            fired_channel = 0
        elif wmp == chpos_max: 
            fired_channel = int(chpos_max)
        else: 
            fired_channel = int(wmp + 0.5)
        cluster[fired_channel] += amplitude 
    else:
        sigma = width * sigma_from_width
        for _ in range(amplitude):
            fired_channel = int(round(random.gauss(mean, sigma)))
            while not 0 <= fired_channel < len(cluster):
                fired_channel = int(round(random.gauss(mean, sigma)))
            cluster[fired_channel] += 1
   
    return cluster
 

◆ cluster_width_mean()

ScifiCluster.cluster_width_mean	(	distance,
		params = `cluster_width_mean_params`
	)

    It return a mean cluster width depending on the distance to SiPMs

Definition at line 428 of file ScifiCluster.py.

def cluster_width_mean(distance, params=cluster_width_mean_params):
 
    """
    It return a mean cluster width depending on the distance to SiPMs
    """
 
    A, B, C = params
    return exp(A + B*distance) + C
 

◆ cluster_width_random()

ScifiCluster.cluster_width_random	(	distance,
		ly,
		percent = `random_width_percent`,
		cluster_width_min = `cluster_width_min`,
		cluster_width_max = `cluster_width_max`
	)

    It generates a cluster. The cluster have 'ly' photoelectrons. 
    The cluster width depends on the distance to SiPM

Definition at line 446 of file ScifiCluster.py.

                         cluster_width_min=cluster_width_min, cluster_width_max=cluster_width_max):
    
    """
    It generates a cluster. The cluster have 'ly' photoelectrons. 
    The cluster width depends on the distance to SiPM
    """
 
    mean = cluster_width_mean(distance)
    sigma = cluster_width_sigma(distance)
    if random.random() <= percent:
        return random.randint(cluster_width_min, cluster_width_max)
    random_width = int(round(random.gauss(mean - 1, sigma))) + 1
 
    # Generate again if the width < minimal width and the light yield < the width
    while random_width < 1 or ly < random_width:
        random_width = int(round(random.gauss(mean - 1, sigma))) + 1
 
    return random_width
 

◆ cluster_width_sigma()

ScifiCluster.cluster_width_sigma	(	distance,
		params = `cluster_width_sigma_params`
	)

    It return a standard deviation of the mean cluster width depending on the distance to SiPMs

Definition at line 437 of file ScifiCluster.py.

def cluster_width_sigma(distance, params=cluster_width_sigma_params):
 
    """
    It return a standard deviation of the mean cluster width depending on the distance to SiPMs
    """
 
    A, B = params
    return A + B*distance
 

◆ cm_to_channel()

ScifiCluster.cm_to_channel	(	locpos,
		sipm_map = `sipm_map`,
		gaps_map = `gaps_map`,
		pitch = `pitch`,
		charr = `charr`,
		reverse = `False`,
		ch_max_num = `ch_max_num`,
		n_solid_sipms = `n_solid_sipms`
	)

FUNCTIONS FOR SCIFICLASS ##################################################.

    It converts a particle position (an event) measured in cm to a position measured 
    in channels. The SiPM map is used. The position is in the scifi module frame.

Definition at line 276 of file ScifiCluster.py.

                  reverse=False, ch_max_num=ch_max_num, n_solid_sipms=n_solid_sipms):
    
    """
    It converts a particle position (an event) measured in cm to a position measured 
    in channels. The SiPM map is used. The position is in the scifi module frame.
    """
    
    if reverse is True:
        for left, right, ch_range in sipm_map:
            if left <= locpos <= right:
                if not ch_range is False:
                    ch_start = ch_range[0]
                    return int(round(ch_max_num - ((locpos-left) / pitch + ch_start)))
                else:
                    mapindex = int(round(n_solid_sipms - (locpos + hw) / charr))
                    return gaps_map.get(mapindex, False)
 
    elif reverse is False:
        for left, right, ch_range in sipm_map:
            if left <= locpos <= right:
                if not ch_range is False:
                    ch_start = ch_range[0]
                    return int(round(locpos-left) / pitch + ch_start)
                else:
                    mapindex2 = int(round(locpos + hw) / charr)
                    return gaps_map.get( mapindex2, False)
 

◆ create_cluster()

ScifiCluster.create_cluster	(	amplitude,
		width,
		wmp
	)

    The final function for creating a signal cluster

Definition at line 568 of file ScifiCluster.py.

def create_cluster(amplitude, width, wmp):
 
    """
    The final function for creating a signal cluster
 
    """
    if wmp is False: 
        return False
    if wmp is None: 
        return False
 
    if not ((chpos_min < wmp) and (wmp< chpos_max)):
        return False
 
    if not ly_min <= amplitude >= width: return False
    shifted_wmp = wmp + 0.5 # For right counting
    cluster = cluster_generator(amplitude, width, shifted_wmp)
 
    # Generate again if it doesn't look like real cluster
    while is_realistic(cluster, width) is False:
        cluster = cluster_generator(amplitude, width, shifted_wmp)
 
    return cluster
 

◆ edep_to_ly()

ScifiCluster.edep_to_ly	(	energy,
		CDF_integral = `CDF_integral`,
		energy_range = `energy_range`,
		ly_linear_params = `ly_linear_params`,
		k_cdfs_corr = `k_cdfs_corr`,
		sigma_in_percent = `sigma_in_percent`,
		ly_CDF_params = `ly_CDF_params`,
		CDF_ly_params = `CDF_ly_params`,
		ly_CDF_landau_params = `ly_CDF_landau_params`
	)

    It returns the light yield calculated from the energy deposit. The calculations are based 
    on the equality of the cumulative distribution functions (CDF) :
    energy => CDF(energy) => CDF(light yield) => ly

    The linear converting range 0.18 MeV < dE < 0.477 MeV corresponds 4.5 < LY < 104 ph.e.

    If energy more then 0.477 MeV the light yield calculated randomly (uniformly in the range) 
    according to the distribution

    Also a little randomness is added to the CDF value with a normal distribution and 
    standard deviation with 'sigma_in_percent' (in percent of the whole range 0 - max CDF)

Definition at line 478 of file ScifiCluster.py.

    ly_CDF_params=ly_CDF_params, CDF_ly_params=CDF_ly_params, ly_CDF_landau_params=ly_CDF_landau_params):
    
    """
    It returns the light yield calculated from the energy deposit. The calculations are based 
    on the equality of the cumulative distribution functions (CDF) :
    energy => CDF(energy) => CDF(light yield) => ly
 
    The linear converting range 0.18 MeV < dE < 0.477 MeV corresponds 4.5 < LY < 104 ph.e.
 
    If energy more then 0.477 MeV the light yield calculated randomly (uniformly in the range) 
    according to the distribution
 
    Also a little randomness is added to the CDF value with a normal distribution and 
    standard deviation with 'sigma_in_percent' (in percent of the whole range 0 - max CDF)
    """
 
    e_min, e_max = energy_range
    A, C = ly_linear_params
    if e_min < energy < e_max:
        ly_lin = A * energy + C
        cdf_raw = approx_function(ly_lin, ly_CDF_landau_params) * k_cdfs_corr
    elif e_max <= energy:
        cdf_raw = CDF_integral * np.random.uniform(0., 1.0) 
    else:
        return 0.    
    cdf_random = random.gauss(cdf_raw, sigma_in_percent * CDF_integral)
 
    # Generate again while the CDF value is not in the range
    while cdf_random < 0 or cdf_random > CDF_integral:
        cdf_random = random.gauss(cdf_raw, sigma_in_percent * CDF_integral)
 
    return approx_function(cdf_random, CDF_ly_params)
 

◆ fcheck_wall()

ScifiCluster.fcheck_wall	(	cluster,
		wmp
	)

    A function that prevents clusters from propagating to the gaps between the mats

Definition at line 592 of file ScifiCluster.py.

def fcheck_wall(cluster, wmp):
 
    """
    A function that prevents clusters from propagating to the gaps between the mats
    """
 
    if (wmp + 0.5)<= firstmatgap_startch:                                               #1st mat
        for chann in cluster:
            if cluster[(chann + 0.5) > firstmatgap_startch] != 0.: cluster[chann] = 0. 
            return cluster
    elif firstmatgap_endch <= (wmp + 0.5)<= secmatgap_startch:                          #2nd mat
        for chann in cluster:
            if (cluster[(chann + 0.5) < firstmatgap_endch] != 0.) or (cluster[(chann + 0.5) > secmatgap_startch] != 0.): cluster[chann] = 0. 
            return cluster
    elif (wmp + 0.5) >= secmatgap_endch:                                                #3rd mat
        for chann in cluster:
            if cluster[(chann + 0.5)< secmatgap_endch] != 0.: cluster[chann] = 0.
            return cluster
    else:
        return False                                                                    #"False" is returned if wmp lies in the gaps
 

◆ GetMatLength()

ScifiCluster.GetMatLength ( )

    It returns a length of a scifi mat. Identical for vertical and horizontal.

Definition at line 343 of file ScifiCluster.py.

def GetMatLength():
 
    """
    It returns a length of a scifi mat. Identical for vertical and horizontal.
    """
 
    return scifimat_length
 

◆ GetMatNum()

ScifiCluster.GetMatNum ( DetID )

    It returns an id (number) of a scifi module. In current version one plane has 3 vertical
    and 3 horizontal scifi assemblies.

Definition at line 334 of file ScifiCluster.py.

def GetMatNum(DetID):
 
    """
    It returns an id (number) of a scifi module. In current version one plane has 3 vertical
    and 3 horizontal scifi assemblies. 
    """
 
    return int(DetID % 100 % 10)
 

◆ GetMatQty()

ScifiCluster.GetMatQty ( )

    It returns a number of scifi mats in a plane. Identical for vertical and horizontal fiber planes.

Definition at line 351 of file ScifiCluster.py.

def GetMatQty():
 
    """
    It returns a number of scifi mats in a plane. Identical for vertical and horizontal fiber planes.
    """
    
    return n_planes
 

◆ GetMatType()

ScifiCluster.GetMatType ( DetID )

    It returns a type of a scifi module.
    1 - vertical scifi assembly
    0 - horizontal scifi assembly

Definition at line 324 of file ScifiCluster.py.

def GetMatType(DetID):
    
    """
    It returns a type of a scifi module.
    1 - vertical scifi assembly
    0 - horizontal scifi assembly
    """
    return int((DetID % 100) // 10)
 
 

◆ GetSiPMs()

ScifiCluster.GetSiPMs ( )

    It returns a number of SiPMs per plane. Identical for vertical and horizontal fiber planes.

Definition at line 359 of file ScifiCluster.py.

def GetSiPMs():
    
    """
    It returns a number of SiPMs per plane. Identical for vertical and horizontal fiber planes.
    """
    
    return n_sipms
 

◆ GetStationNum()

ScifiCluster.GetStationNum ( DetID )

    It returns an id of a plane. In current version the detector has 5 TT stations (numbering starts
    at 1.

Definition at line 367 of file ScifiCluster.py.

def GetStationNum(DetID):
 
    """
    It returns an id of a plane. In current version the detector has 5 TT stations (numbering starts
    at 1.
    """
    return int(DetID // 100)
 

◆ global_to_local()

ScifiCluster.global_to_local	(	DetID,
		globalpos
	)

These functions assume that the detector is aligned with the beam ######.

    It returns the local coordinates in one scifi assembly frame from global coordinates.

Definition at line 379 of file ScifiCluster.py.

def global_to_local(DetID, globalpos):
 
    """
    It returns the local coordinates in one scifi assembly frame from global coordinates.
    """
    
    global xshift, yshift
 
    if GetMatType(DetID) == 0:
        return globalpos - yshift
    
    if GetMatType(DetID) == 1:
        return globalpos - xshift
 
 

◆ is_realistic()

ScifiCluster.is_realistic	(	cluster,
		width
	)

    It returns TRUE if cluster is realistic: it doesn't have a gap between numders, like
    [..., 0, 1, 2, 0, 0, 5, 6, ...], and it doens't have the light yield less then width.

Definition at line 553 of file ScifiCluster.py.

def is_realistic(cluster, width):
 
    """
    It returns TRUE if cluster is realistic: it doesn't have a gap between numders, like
    [..., 0, 1, 2, 0, 0, 5, 6, ...], and it doens't have the light yield less then width.
    """
 
    cluster_only_values = [(channel, value) for channel, value in enumerate(cluster) if value > 0]
    first_channel, _ = cluster_only_values[0]
    last_channel, _ = cluster_only_values[-1]
    if len(cluster_only_values) != width or (last_channel - first_channel + 1) != width:
        return False
    else:
        return True
    

◆ local_to_global()

ScifiCluster.local_to_global	(	DetID,
		localpos
	)

    It returns the global coordinates from the local coordinates in one scifi assembly frame.

Definition at line 394 of file ScifiCluster.py.

def local_to_global(DetID, localpos):
 
    """
    It returns the global coordinates from the local coordinates in one scifi assembly frame.
    """
 
    global xshift, yshift
 
    if GetMatType(DetID) == 0:
        return localpos + yshift
 
    if GetMatType(DetID) == 1:
        return localpos + xshift
 

◆ ly_attenuation()

ScifiCluster.ly_attenuation ( distance )

    It return the light yield losses in percent depending on the distance to SiPMs

Definition at line 419 of file ScifiCluster.py.

def ly_attenuation(distance):
 
    """
    It return the light yield losses in percent depending on the distance to SiPMs
    """
 
    res_ly_loss = ly_loss_mean(distance) / ly_loss_mean(0.) 
    return res_ly_loss
 

◆ ly_loss_mean()

ScifiCluster.ly_loss_mean	(	distance,
		params = `ly_loss_params`
	)

    It returns the light yield depending on the distance to SiPMs

Definition at line 410 of file ScifiCluster.py.

def ly_loss_mean(distance, params=ly_loss_params):
 
    """
    It returns the light yield depending on the distance to SiPMs
    """
 
    A1, k1, A2, k2 = params
    return A1 * exp(k1 * distance / 100.) + A2 * exp(k2 * distance / 100.)
 

◆ weigthed_mean_pos()

ScifiCluster.weigthed_mean_pos ( cluster )

    Calculate the weighted mean position of the cluster

Definition at line 613 of file ScifiCluster.py.

def weigthed_mean_pos(cluster):
 
    """
    Calculate the weighted mean position of the cluster
    """
 
    if cluster is False:
        return False
    sumup = 0.
    sumdown = 0.
    wmp = 0.
    for channel, value in enumerate(cluster):
        if value > 0:
            sumup += value * channel
            sumdown += value
    if sumdown != 0:
        wmp = sumup / sumdown - 0.5
    return wmp
 

Variable Documentation

◆ array

float ScifiCluster.array = gap + 2.*charr

Definition at line 43 of file ScifiCluster.py.

◆ biggap

float ScifiCluster.biggap = 0.038

Definition at line 44 of file ScifiCluster.py.

◆ CDF_integral

float ScifiCluster.CDF_integral = 0.0185640424

Definition at line 171 of file ScifiCluster.py.

◆ CDF_ly_params

dict ScifiCluster.CDF_ly_params

Initial value:

=  {
    (0., 0.0006): (
        (89., 4.152, 0.0001574, -1.326e+04, 4.3), 
        lambda cdf, *p: p[0] * sqrt(p[1]*(cdf+p[2])) * (1-exp(p[3]*cdf)) + p[4]
    ),
    (0.0006, 0.012): (
        (158, 1.035, 0.24, 217),
        lambda cdf, *p: p[0] * log(sqrt(cdf)+p[1]) + p[2]*exp(p[3]*cdf)
    ),
    (0.012, 0.018561405): (
        (9.36, 335.984, -18100, -400, 15),
        lambda cdf, *p: p[0] * log((p[1]-p[2]*cdf)/(p[1]+p[2]*cdf)) + p[3]*cdf + p[4]
    ),
    (0.018561405, 0.0185640424): (
        (9.36, 335.984, -18100, -400, 15),
        lambda cdf, *p: (p[0] * log((p[1]-p[2]*0.018561405)/(p[1]+p[2]*0.018561405)) 
            + p[3]*0.018561405 + p[4])
    )
}

Definition at line 222 of file ScifiCluster.py.

◆ ch_max_num

int ScifiCluster.ch_max_num = 128*12 - 1

Definition at line 45 of file ScifiCluster.py.

◆ charr

float ScifiCluster.charr = 1.6

Definition at line 40 of file ScifiCluster.py.

◆ chpos_max

float ScifiCluster.chpos_max = 1535.5

Definition at line 145 of file ScifiCluster.py.

◆ chpos_min

float ScifiCluster.chpos_min = -0.5

Definition at line 144 of file ScifiCluster.py.

◆ cluster_width_max

int ScifiCluster.cluster_width_max = 10

Definition at line 147 of file ScifiCluster.py.

◆ cluster_width_mean_params

float ScifiCluster.cluster_width_mean_params = 0.6661, -0.006955, 2.163

Definition at line 155 of file ScifiCluster.py.

◆ cluster_width_min

int ScifiCluster.cluster_width_min = 1

Definition at line 146 of file ScifiCluster.py.

◆ cluster_width_sigma_params

float ScifiCluster.cluster_width_sigma_params = 1.103, -0.0005751

Definition at line 156 of file ScifiCluster.py.

◆ edge

float ScifiCluster.edge = 0.016

Definition at line 41 of file ScifiCluster.py.

◆ energy_range

float ScifiCluster.energy_range = 0.18, 0.477

Definition at line 168 of file ScifiCluster.py.

◆ firstmatgap_endch

float ScifiCluster.firstmatgap_endch = 523.5

Definition at line 139 of file ScifiCluster.py.

◆ firstmatgap_startch

float ScifiCluster.firstmatgap_startch = 521.5

Definition at line 138 of file ScifiCluster.py.

◆ gap

float ScifiCluster.gap = 0.022

Definition at line 42 of file ScifiCluster.py.

◆ gaps_map

dict ScifiCluster.gaps_map

Initial value:

=  {
: -0.5, 
: 63.5,
: 127.5, 
: 191.5,
: 255.5, 
: 319.5,
: 383.5,
: 447.5,
 
: 511.5,
: 575.5, 
: 639.5,
: 703.5, 
: 767.5,
: 831.5, 
: 895.5,
: 959.5,
 
: 1023.5, 
: 1087.5,
: 1151.5, 
: 1215.5,
: 1279.5, 
: 1343.5,
: 1407.5, 
: 1471.5,
: 1535.5
}

Definition at line 105 of file ScifiCluster.py.

◆ hw

float ScifiCluster.hw = 19.528

Definition at line 38 of file ScifiCluster.py.

◆ is_sipm_hor_reversed

bool ScifiCluster.is_sipm_hor_reversed = True

Definition at line 23 of file ScifiCluster.py.

◆ is_sipm_vert_reversed

bool ScifiCluster.is_sipm_vert_reversed = True

Definition at line 24 of file ScifiCluster.py.

◆ k_cdfs_corr

float ScifiCluster.k_cdfs_corr = 0.993076766938

Definition at line 178 of file ScifiCluster.py.

◆ ly_CDF_landau_params

dict ScifiCluster.ly_CDF_landau_params

Initial value:

=  {
    (4.5, 15): (
        (0.001038, -0.000378, 3.53e-05), 
        lambda ly, *p: p[0] + p[1]*ly + p[2]*ly**2
    ),
    (15, 40): (
        (-0.001986, -0.0003014, 7.031e-05, -2.067e-06, 1.892e-08),
        lambda ly, *p: p[0] + p[1]*ly + p[2]*ly**2 + p[3]*ly**3 + p[4]*ly**4
    ),
    (40, 104): (
        (-0.007149, 0.001056, -1.779e-05, 1.41e-07, -4.29e-10),
        lambda ly, *p: p[0] + p[1]*ly + p[2]*ly**2 + p[3]*ly**3 + p[4]*ly**4
    )
}

Definition at line 206 of file ScifiCluster.py.

◆ ly_CDF_params

dict ScifiCluster.ly_CDF_params

Initial value:

=  {
    (4.5, 13): (
        (-13.2, 1.976), 
        lambda ly, *p: exp(p[0] + p[1]*sqrt(ly))
    ),
    (13, 104): (
        (0.0108183, -0.179752, -19.2, 0.00772965),
        lambda ly, *p: p[0]*(1 - exp(p[1]*(ly+p[2]))) / (1 + exp(p[1]*(ly+p[2]))) + p[3]
    )
}

Definition at line 192 of file ScifiCluster.py.

◆ ly_linear_params

float ScifiCluster.ly_linear_params = 332.882, -58.7085

Definition at line 175 of file ScifiCluster.py.

◆ ly_loss_params

float ScifiCluster.ly_loss_params = 20.78, -0.26, 7.89, -3.06

Definition at line 153 of file ScifiCluster.py.

◆ ly_loss_sigma_params

float ScifiCluster.ly_loss_sigma_params = 6.8482, -0.5757, 3.4458

Definition at line 154 of file ScifiCluster.py.

◆ ly_max

int ScifiCluster.ly_max = 104

Definition at line 164 of file ScifiCluster.py.

◆ ly_min

float ScifiCluster.ly_min = 4.5

Definition at line 163 of file ScifiCluster.py.

◆ n_planes

ScifiCluster.n_planes = snd_geo.Scifi.nscifi

Definition at line 259 of file ScifiCluster.py.

◆ n_sipms

ScifiCluster.n_sipms = snd_geo.Scifi.nsipms

Definition at line 260 of file ScifiCluster.py.

◆ n_solid_sipms

int ScifiCluster.n_solid_sipms = 12

Definition at line 46 of file ScifiCluster.py.

◆ nmats

ScifiCluster.nmats = snd_geo.Scifi.nmats

Definition at line 258 of file ScifiCluster.py.

◆ pitch

float ScifiCluster.pitch = 0.025

Definition at line 39 of file ScifiCluster.py.

◆ random_width_percent

float ScifiCluster.random_width_percent = 0.01

Definition at line 160 of file ScifiCluster.py.

◆ scifi_x

ScifiCluster.scifi_x = snd_geo.Scifi.xdim

Definition at line 254 of file ScifiCluster.py.

◆ scifi_y

ScifiCluster.scifi_y = snd_geo.Scifi.ydim

Definition at line 255 of file ScifiCluster.py.

◆ scifimat_gap

ScifiCluster.scifimat_gap = snd_geo.Scifi.scifimat_gap

Definition at line 252 of file ScifiCluster.py.

◆ scifimat_length

ScifiCluster.scifimat_length = snd_geo.Scifi.scifimat_length

Definition at line 251 of file ScifiCluster.py.

◆ scifimat_width

ScifiCluster.scifimat_width = snd_geo.Scifi.scifimat_width

Definition at line 250 of file ScifiCluster.py.

◆ scifimat_z

ScifiCluster.scifimat_z = snd_geo.Scifi.scifimat_z

Definition at line 256 of file ScifiCluster.py.

◆ secmatgap_endch

float ScifiCluster.secmatgap_endch = 1048.5

Definition at line 141 of file ScifiCluster.py.

◆ secmatgap_startch

float ScifiCluster.secmatgap_startch = 1046.5

Definition at line 140 of file ScifiCluster.py.

◆ sigma_from_width

int ScifiCluster.sigma_from_width = 1 / 4.

Definition at line 185 of file ScifiCluster.py.

◆ sigma_in_percent

float ScifiCluster.sigma_in_percent = 0.01

Definition at line 182 of file ScifiCluster.py.

◆ sipm_hor_pos

int ScifiCluster.sipm_hor_pos = +1

GLOBAL VALUES AND PARAMETERS ############.

Definition at line 18 of file ScifiCluster.py.

◆ sipm_map

tuple ScifiCluster.sipm_map

Definition at line 50 of file ScifiCluster.py.

◆ sipm_vert_pos

int ScifiCluster.sipm_vert_pos = -1

Definition at line 19 of file ScifiCluster.py.

◆ snd_geo

ScifiCluster.snd_geo = ConfigRegistry.loadpy("$SNDSW_ROOT/geometry/shipLHC_geom_config.py")

LOAD THE SCIFI GEOMETRY ############.

Definition at line 248 of file ScifiCluster.py.

◆ xshift

float ScifiCluster.xshift = -27.8

Need to automatize this part so it works for every geometry #########.

Definition at line 269 of file ScifiCluster.py.

◆ yshift

float ScifiCluster.yshift = +35.3

Definition at line 270 of file ScifiCluster.py.

Classes

Functions

Variables

Function Documentation

◆ approx_function()

◆ channel_to_cm()

◆ cluster_generator()

◆ cluster_width_mean()

◆ cluster_width_random()

◆ cluster_width_sigma()

◆ cm_to_channel()

◆ create_cluster()

◆ edep_to_ly()

◆ fcheck_wall()

◆ GetMatLength()

◆ GetMatNum()

◆ GetMatQty()

◆ GetMatType()

◆ GetSiPMs()

◆ GetStationNum()

◆ global_to_local()

◆ is_realistic()

◆ local_to_global()

◆ ly_attenuation()

◆ ly_loss_mean()

◆ weigthed_mean_pos()

Variable Documentation

◆ array

◆ biggap

◆ CDF_integral

◆ CDF_ly_params

◆ ch_max_num

◆ charr

◆ chpos_max

◆ chpos_min

◆ cluster_width_max

◆ cluster_width_mean_params

◆ cluster_width_min

◆ cluster_width_sigma_params

◆ edge

◆ energy_range

◆ firstmatgap_endch

◆ firstmatgap_startch

◆ gap

◆ gaps_map

◆ hw

◆ is_sipm_hor_reversed

◆ is_sipm_vert_reversed

◆ k_cdfs_corr

◆ ly_CDF_landau_params

◆ ly_CDF_params

◆ ly_linear_params

◆ ly_loss_params

◆ ly_loss_sigma_params

◆ ly_max

◆ ly_min

◆ n_planes

◆ n_sipms

◆ n_solid_sipms

◆ nmats

◆ pitch

◆ random_width_percent

◆ scifi_x

◆ scifi_y

◆ scifimat_gap

◆ scifimat_length

◆ scifimat_width

◆ scifimat_z

◆ secmatgap_endch

◆ secmatgap_startch

◆ sigma_from_width

◆ sigma_in_percent

◆ sipm_hor_pos

◆ sipm_map

◆ sipm_vert_pos

◆ snd_geo

◆ xshift

◆ yshift