tof_v20f_mcbm.digibdf.par 20.8 KB
Newer Older
Norbert Herrmann's avatar
Norbert Herrmann committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
####################################################################################################
# Class:   CbmTofDigiBdfPar                                                                                                         
# Context: TestDefaultContext                                                                                                       
####################################################################################################
[CbmTofDigiBdfPar]                                                                                                                  
// ----------------------------------------------------------------------------------------//
// --------------------------- General switches -------------------------------------------//                                     
// Digi type to use:                                                                                                                
// 1 = expanded digis (faster),                                                                                                     
// 0 = reDigitized digis (probably closer to what could be transported in CBMnet                                                    
UseExpDigi: Int_t           1
// Switch to use only primary tracks
// 0 = all tracks, 1 = only primary
UseOnlyPrim: Int_t          0
// Switch method to try avoiding counting all gaps/TofPoint from same track
// with different time/charge
UseOneGapPerTrk: Int_t      1                                      
// ----------------------------------------------------------------------------------------//
// ------------------------------ Cluster Model -------------------------------------------//  
//******************************//                                                                                                  
//   Cluster Model selector     //                                                                                                  
//******************************//                                                                                                  
// 0 -> Directly Cluster Size                                                                                                       
// ------- Convert TofPoints in input TClonesArray to Tof Digis using directly the                                                  
// ------- cluster size distribution from beamtime. The charge is encoded as ToT in ns and                                          
// ------- shared among strips using a gauss distribution centered on the TofPoints.                                                
// ------- Works only for strip RPC !!!!!                                          
// ------- Quite slow due to the integration over a TF1 for charge.                                                                                           
// 1 -> A simple "Flat disc" cluster approximation                                                                                  
// ------- Convert TofPoints in input TClonesArray to Tof Digis using an approximation of the                                       
// ------- cluster charge distribution by a "Flat disc". The charge is encoded as ToT in ns.                                        
// 2 -> A 2D Gauss cluster approximation                                                                                            
// ------- Convert TofPoints in input TClonesArray to Tof Digis using an approximation of the                                       
// ------- cluster charge distribution by a 2D symetric Gauss function, with the approx. that                                       
// ------- sigmaX = sigmaY = cluster size/2.                  
// ------- Really slow due to the integration over a TF2.                                                                                      
ClusterModel: Int_t         1                                                                                                       
// Width of the electronic channels gain distribution                                                                               
// A value of 0 or smaller set all channels to gain 1                                                                               
FeeGainSigma: Double_t      0.01                                                                                                    
// "Charge" Threshold in ToT[ps], cut channels with too small values                                                                
//FeeTotThr: Double_t         0.01                                                             
FeeTotThr: Double_t         0.1                                                                                                         
// Time resolution of the digitizer channels in [ps]                                                                                
TimeResElec: Double_t       0.04                                                                                                       
// Time resolution of the digitizer channels in [ps]                                                                                
TimeResStart: Double_t      0.04                                                                                                   
// Propagation speed of the signal in the readout electrodes, in [cm/ns]]                                                           
SignalPropSpeed: Double_t  18.0                                                                                                     
// ----------------------------------------------------------------------------------------//                                       
// ------- Geometry variables, text to be generated in the CreateGeometry macros ----------//                                       
// Number of Super Module types (1)                                                                                                 
// !!!! this parameter HAS to match the value used in the geometry !!!!                                                             
NbSmTypes: Int_t            10                                                                                                       
// Number of SMs for each SM type                                                                                                   
// The number of entries should match (1), the number of types                                                                      
// !!!! this parameter HAS to xmatch the value used in the geometry !!!!                                                             
//  0    1    2    3    4    5    6    7    8    9                                                                                  
NbSm: Int_t \                                                                                                                       
Norbert Herrmann's avatar
Norbert Herrmann committed
60
    5    0    0    0    0    1    1    0    0    1                                                
Norbert Herrmann's avatar
Norbert Herrmann committed
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
// Number of RPCs for each SM type (2)                                                                                              
// The number of entries should match (1), the number of types                                                                      
// !!!! this parameter HAS to match the value used in the geometry !!!!                                                             
//  0    1    2    3    4    5    6    7    8    9                                                                                  
NbRpc: Int_t \                                                                                                                      
    5    3    1    1    1    1    2    4    6    2                                                
// Number of gaps per RPC for each SM type                                                                                      
// 1 line per SM type, the number of entries for each line should match                                                             
// the number of RPCs defined in (2)                                                                                                
// !!!! this parameter HAS to match the value used in the geometry !!!!                                                             
//  0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15                                                    
NbGaps000: Int_t \                                                                                                                    
    8    8    8    8    8    
NbGaps001: Int_t \                                                                                                                   
    8    8    8    8    8
NbGaps002: Int_t \                                                                                                                   
    8    8    8    8    8						
NbGaps003: Int_t \                                                                                                                    
    8    8    8    8    8    
NbGaps004: Int_t \                                                                                                                    
    8    8    8    8    8    
NbGaps005: Int_t \                                                                                                                    
    8    8    8    8    8   
NbGaps006: Int_t \                                                                                                                    
    8    8    8    8    8   
NbGaps007: Int_t \                                                                                                                    
    8    8    8    8    8    8
NbGaps008: Int_t \                                                                                                                    
    8    8    8    8    8    8  
NbGaps009: Int_t \                                                                                                                    
    8    8    8    8    8 
// Gap Size in mm per RPC for each SM type, including the glass type
// 1 line per SM type, the number of entries for each line should match                                                             
// the number of RPCs defined in (2)                                                                                                
// !!!! this parameter HAS to match the value used in the geometry !!!!                                                             
//  0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15                                                    
GapSize000: Double_t \                                                                                                                    
 1.25 1.25 1.25 1.25 1.25
GapSize001: Double_t \                                                                                                                
 1.25 1.25 1.25 1.25 1.25
GapSize002: Double_t \                                                                                                            
 1.25 1.25 1.25 1.25 1.25
GapSize003: Double_t \                                                                                                                
 1.25 1.25 1.25 1.25 1.25 
GapSize004: Double_t \                                                                                                            
 1.25 1.25 1.25 1.25 1.25 
GapSize005: Double_t \                                                                                                              
 1.25 1.25 1.25 1.25 1.25 
GapSize006: Double_t \                                                                                                              
 1.25 1.25 1.25 1.25 1.25 
GapSize007: Double_t \                                                                                                              
 1.25 1.25 1.25 1.25 1.25 1.25
GapSize008: Double_t \                                                                                                              
 1.25 1.25 1.25 1.25 1.25 1.25 
GapSize009: Double_t \                                                                                                              
 1.25 1.25 1.25 1.25 1.25
// Signal velocity in cm/ns per RPC for each SM type, including the glass type
// 1 line per SM type, the number of entries for each line should match                                                             
// the number of RPCs defined in (2)                                                                                                
// !!!! this parameter HAS to match the value used in the geometry !!!!                                                             
//  0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15                                                    
SigVel000: Double_t \
16.30 16.43 16.38 16.50 16.58 \
16.16 16.21 16.22 16.20 16.20 \
16.15 16.20 16.30 16.20 16.13 \
Norbert Herrmann's avatar
Norbert Herrmann committed
126
16.27 16.50 16.40 16.49 16.60 \
Norbert Herrmann's avatar
Norbert Herrmann committed
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
16.25 16.25 16.19 16.35 16.27 
SigVel001: Double_t \
 17.75 18.18 18.07 16.2 16.2 15.9 16.2 16.2 16.2 
SigVel002: Double_t \
 18.0 18.0 18.0 
SigVel003: Double_t \
 18.26 
SigVel004: Double_t \
 13.49 
SigVel005: Double_t \
 18. 18. 18. 18. 18.
SigVel006: Double_t \
 17.55 14.6
SigVel007: Double_t \
 15. 15. 15. 15. 
SigVel008: Double_t \
 18. 18. 18. 18. 18. 18.
SigVel009: Double_t \
 17.77 17.75 16. 16. 16.
TrkStation000: Int_t \
    1    0    1    0    1 \
    3    2    3    2    3 \
    5    4    5    4    5 \
    1    0    1    0    1 \
    5    4    5    4    5 
TrkStation001: Int_t \
    0    0    0    0    0 
TrkStation002: Int_t \
    0    0    0    0    0  
TrkStation003: Int_t \
    0    0    0    0    0 
TrkStation004: Int_t \
    0    0    0    0    0 
TrkStation005: Int_t \
    0    0    0    0    0 
TrkStation006: Int_t \
    0    0    0    0    0 
TrkStation007: Int_t \
    0    0    0    0    0 
TrkStation008: Int_t \
    0    0    0    0    0     
TrkStation009: Int_t \
    0    0    0    0    0         
// Number of channels per RPC for each SM type                                                                                      
// 1 line per SM type, the number of entries for each line should match                                                             
// the number of RPCs defined in (2)                                                                                                
// !!!! this parameter HAS to match the value used in the geometry !!!!                                                             
//  0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15                                                    
NbCh000: Int_t \                                                                                                                    
   32  32  32  32  32                                                                                                            
NbCh001: Int_t \                                                                                                                    
   32  32  32                                                                                                             
NbCh002: Int_t \                                                                                                                    
    1    1                                                                                            
NbCh003: Int_t \                                                                                                                    
   32    
NbCh004: Int_t \                                                                                                                    
Norbert Herrmann's avatar
Norbert Herrmann committed
184
   32  32  32  32       
Norbert Herrmann's avatar
Norbert Herrmann committed
185
186
187
NbCh005: Int_t \                                                                                                                    
   8   16   16  
NbCh006: Int_t \                                                                                                                    
Norbert Herrmann's avatar
Norbert Herrmann committed
188
   32   32   
Norbert Herrmann's avatar
Norbert Herrmann committed
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
NbCh007: Int_t \                                                                                                                    
   40    40   40   40    
NbCh008: Int_t \                                                                                                                    
    8    8    8    8    8    8    
NbCh009: Int_t \                                                                                                                    
   32   32   32   32   32  
// Type of readout channel for each rpc for each SM type                                                                            
// 0 = strip (readout on both end), 1 = pad (readout on 1 side)
// 1 line per SM type, the number of entries for each line should match
// the number of RPCs defined in (2)
// !!!! this parameter HAS to match the value used in the geometry !!!!
//  0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15
ChType000: Int_t \
    0    0    0    0    0
ChType001: Int_t \
    0    0    0    0    0
ChType002: Int_t \
    0    0    0    0    0 
ChType003: Int_t \
    0    0    0    0    0    
ChType004: Int_t \
    0    0    0    0    0    
ChType005: Int_t \
    0    0    0    0    0    
ChType006: Int_t \
    0    0    0    0    0    
ChType007: Int_t \
    0    0    0    0    0    0    
ChType008: Int_t \
    0    0    0    0    0    0
ChType009: Int_t \
    0    0    0    0    0  
// Orientation of the channels/rpc
// reused to store tracking station number 
// 1 line per SM type, the number of entries for each line should match
// the number of RPCs defined in (2)
// !!!! this parameter HAS to match the value used in the geometry !!!!
//  0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15
ChOrientation000: Int_t \
    0    0    0    0    0
ChOrientation001: Int_t \
    0    0    0    0    0
ChOrientation002: Int_t \
    0    0    0    0    0
ChOrientation003: Int_t \
    0    0    0    0    0    
ChOrientation004: Int_t \
    0    0    0    0    0    
ChOrientation005: Int_t \
    0    0    0    0    0  
ChOrientation006: Int_t \
    0    0    0    0    0  
ChOrientation007: Int_t \
    0    0    0    0    0    0
ChOrientation008: Int_t \
    0    0    0    0    0    0
ChOrientation009: Int_t \
    0    0    0    0    0    
// ----------------------------------------------------------------------------------------//
// ------------------------- Beamtime variables -------------------------------------------//
// Name of the root file holding the values/histograms extracted from beamtime data
BeamInputFile: Text_t ./test_bdf_input.root
BeamCalibFile: Text_t ./CALIB
// Model to use for the cluster radius distribution
// (not used for Cluster Model 0)
// -1 = Fixed value at 0.0002 to get a cluster size as close to 1 as possible,
//      Combined with a MaxTimeDistClust at clustering of 0.0001ns, should emulate hit producer
//  0 = fixed value from beam data mean, 
//  1 = Landau distrib. with fixed param, 
//  2 = Landau distrib. with params extracted from fit on beam data
ClusterRadiusModel: Int_t    1
// Mapping between geometry SM type and SM types in the beamtime file
// The number of entries should match (1), the number of types in geometry
//  0    1    2    3    4    5    6    7    8    9
SmTypeInpMapp: Int_t \
    0    0    0    0    0    0    0    0    0    0
// "Most Probable Value" setting of the Landau distribution used to modelize the cluster radius
// Used if ClusterRadiusModel = 1 and as default in case of fit failure if ClusterRadiusModel = 2
// The number of entries should match (1), the number of types in geometry
RadiusLandauMpv: Double_t \
 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06 0.06
// "Sigma (scale)" setting of the Landau distribution used to modelize the cluster radius
// Used if ClusterRadiusModel = 1 and as default in case of fit failure if ClusterRadiusModel = 2
// Empirical best: 1/3 of MPV
// The number of entries should match (1), the number of types in geometry
RadLandauSigma: Double_t \
 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02
// ----------------------------------------------------------------------------------------//
// ------------------------- Simple clusterizer -------------------------------------------//
// Switch what to use for the multiplicity flag in TofHit
// 0 = number of different TofPoints generating the cluster 
// 1 = number of different TrackId generating the cluster
// Both should be equivalent when using 
MulUseTrkId: Int_t          1
// Maximal Time distance between current Digi pair (L+R) and the last Digi pair (L+R) in cluster
// in [ns]
// If -1 a values of 5*0.08 = 0.4 will be used instead (5*Nominal System Resolution)
MaxTimeDistClust: Double_t  -1.0
//MaxTimeDistClust: Double_t  0.0001
// Maximal distance along the channel direction (Y for vertical strips) between current 
// Digi pair (L+R) and the last Digi pair (L+R) in cluster
//  in [cm along channel] 
// "Nominal" Res. (80ps) * Prop. speed ~ 1.44
// sqrt(2)*Fee Res. (40ps) * Prop. speed ~ 1
MaxSpaceDistClust: Double_t 2.5
####################################################################################################