Skip to content

Commit cd30cca

Browse files
plexoosfisyak
authored and
genevb
committed
Add local DB parameters for TpcRS 26p5GeV_fixedTarget_2018 (star-bnl#107, star-bnl#88)
New T0 parameters are tuned by Yue-Hang Leung Co-authored-by: Yuri Fisyak <fisyak@bnl.gov>
1 parent 6e991fa commit cd30cca

1 file changed

Lines changed: 86 additions & 0 deletions

File tree

‎StarDb/Calibrations/tpc/TpcResponseSimulator.20180605.071409.C‎

Lines changed: 86 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,86 @@
1+
TDataSet *CreateTable() {
2+
// -----------------------------------------------------------------
3+
// db/.const/StarDb/Calibrations/tpc/.TpcResponseSimulator/TpcResponseSimulator Allocated rows: 1 Used rows: 1 Row size: 124 bytes
4+
// Table: TpcResponseSimulator_st[0]--> TpcResponseSimulator_st[0]
5+
// ====================================================================
6+
// ------ Test whether this table share library was loaded ------
7+
if (!TClass::GetClass("St_TpcResponseSimulator")) return 0;
8+
TpcResponseSimulator_st row;
9+
St_TpcResponseSimulator *tableSet = new St_TpcResponseSimulator("TpcResponseSimulator",1);
10+
memset(&row,0,tableSet->GetRowSize());
11+
row.I0 = 13.1;// eV, CH4
12+
row.Cluster = 3.2; // average no. of electrons per primary
13+
row.W = 26.2;// eV
14+
row.OmegaTau = 3.02;// fit of data
15+
row.K3IP = 0.68;//(pads) for a/s = 2.5e-3 and h/s = 0.5
16+
row.K3IR = 0.89;//(row) for a/s = 2.5e-3 and h/s = 0.5
17+
row.K3OP = 0.55;//(pads) for a/s = 2.5e-3 and h/s = 1.0
18+
row.K3OR = 0.61;//(row) for a/s = 2.5e-3 and h/s = 1.0
19+
row.FanoFactor = 0.3; //
20+
row.AveragePedestal = 50.0;//
21+
row.AveragePedestalRMS = 1.4; // Old Tpc electronics or iTPC
22+
row.AveragePedestalRMSX = 0.7; // New Tpx electronics
23+
row.tauIntegration = 2.5*74.6e-9;// secs
24+
row.tauF = 394.0e-9;// secs Tpc
25+
row.tauP = 775.0e-9;// secs Tpc
26+
row.tauXI = 60.0e-9;// secs Tpx Inner integration time
27+
row.tauXO = 74.6e-9;// secs Tpx Outer integration time
28+
row.tauCI = 0;
29+
row.tauCO = 0;
30+
row.SigmaJitterTI = 0;// 0.4317;// 0.25;//ad 0.0;// b for Tpx inner
31+
row.SigmaJitterTO = 0;// 0.4300;// E: 0.4801;//0.25;//ad 0.0;// b for Tpx outer
32+
row.SigmaJitterXI = 0;// 0.1027785; // P: 0.1353*1.05/1.10; //O: 0.1353*1.05;// N: 0.1353; // C:0.;
33+
row.SigmaJitterXO = 0;// 0.107525; // P: 0.1472*1.05/1.03; //O: 0.1472*1.05;// N: 0.1472; // C:0.;
34+
row.longitudinalDiffusion = 0.03624; // Magboltz // HD 0.03624*1.5; //HC 0.03624; // Magboltz
35+
row.transverseDiffusion = 0.02218*TMath::Sqrt(1 + row.OmegaTau*row.OmegaTau) ; // Magboltz
36+
row.NoElPerAdcI = 335./TMath::Exp(-0.81625); // No. of electrons per 1 ADC count for outer TPX
37+
row.NoElPerAdcO = 335./TMath::Exp(-1.09801); // No. of electrons per 1 ADC count for outer TPX
38+
row.NoElPerAdcX = 335./TMath::Exp(-0.49463);; // No. of electrons per 1 ADC count for iTPC
39+
row.OmegaTauScaleI = 2.145*1.515;// HC 1.;// 2.145*1.515; //i; 2.145*1.4; //h 2.145; //ad 2.145*1.25; //b effective reduction of OmegaTau near Inner sector anode wire
40+
row.OmegaTauScaleO = 1.8 *1.201; //HC 1.;// 1.8 *1.201; //i 1.8 *1.1; //h 1.8; //ad 1.8 *1.25; //b effective reduction of OmegaTau near Outer sector anode wire
41+
// Inner_wire_to_plane_coupling ( 0.533 ) * Inner_wire_to_plane_couplingScale ( 0.843485 )
42+
// Outer_wire_to_plane_coupling ( 0.512 ) * Outer_wire_to_plane_couplingScale ( 0.725267 )
43+
44+
// row.SecRowCorOW[0] = row.SecRowCorOE[0] = 0.11; // IRAKLI : based on shift seen in the MuDst based simulation;
45+
row.SecRowCorIW[0] = 0.013;
46+
row.SecRowCorOW[0] = -0.013;
47+
row.SecRowCorIE[0] = -0.035;
48+
row.SecRowCorOE[0] = -0.013;
49+
// row.SecRowCorIW[0] = row.SecRowCorIE[0] = 6.99114715017355337e-01;//- TMath::Log(0.533*0.843485) -5.84129e-01 + 4.52885e-01 + 3.09117e-02;
50+
// row.SecRowCorOW[0] = row.SecRowCorOE[0] = 9.79357470004933006e-01;//- TMath::Log(0.512*0.725267) -5.47141e-01 + 5.23937e-01 + 1.19154e-02;
51+
// SecRow3CGFdaq_2011_pp500LowLum => Inner: 3.26428e-01 - -5.01720e-04*y; Outer: 2.68883e-01 + 1.23403e-04*y
52+
// 3.22907e-01 2.72715e-01
53+
// SecRow3CGFTpcRS_2011_pp500LowLum_f : 3.09711e-01 2.65342e-01
54+
// diff : 9.13675e-02 6.29849e-02
55+
// SecRow3CGFTpcRS_2011_pp500LowLum_g : 3.12857e-01 2.67379e-01
56+
const Double_t RowSigmaTrs[8] = {
57+
-5.00e-03, 0, // Inner W
58+
-5.00e-03, 0, // Outer W
59+
-1.76e-02, 0, // Inner E
60+
-1.76e-02, 0}; // Outer E
61+
Float_t *b = &row.SecRowSigIW[0];
62+
for (Int_t i = 0; i < 8; i++) {
63+
b[i] = RowSigmaTrs[i%8];
64+
}
65+
/* Sigmas
66+
Tpx inner = 0.395
67+
Tpx outer = 0.314
68+
iTpc = 0.322 Outer 0.316 */
69+
row.PolyaInner = 1.38;
70+
row.PolyaOuter = 1.38;
71+
// row.T0offset = 0.50; // From Lokesh Kumar for Run X
72+
// TpcT->Draw("fMcHit.mMcl_t+0.165*Frequency-fRcHit.mMcl_t/64:fMcHit.mPosition.mX3>>T(210,-210,210,100,-2,3)","fNoMcHit==1&&fNoRcHit==1&&fRcHit.mQuality>90","colz")
73+
// TpcT->Draw("fMcHit.mPosition.mX3-fRcHit.mPosition.mX3:fMcHit.mPosition.mX3>>Z(210,-210,210,100,-2,3)","fNoMcHit==1&&fNoRcHit==1&&fRcHit.mQuality>90","colz")
74+
// The corection has to be added M P
75+
//row.T0offset = 0.50 + 1.65431e-01 - 3.45247e-01 -1.54583e+00 -2.90686e-03+ 1.54353e+00 + 0.0191135 -1.20938e-03 ; //E
76+
// row.T0offset = 0.50 -1.43663e-01;//g // 01/18/12 Xianglei Zhu from Run 11 AuAu 27 & 19.6 GeV embedding
77+
row.T0offset = 0.50 -1.43663e-01 -0.0839244;
78+
// root.exe T0offset.C
79+
// TI->FitSlicesY(); TI_1->Fit("pol2","er","",-100,100);
80+
row.T0offsetI = 1.17437e-01 + 8.43584e-03 + 7.537344e-02 - 0.73383303057411764 + 0.60752;
81+
// TO->FitSlicesY(); TO_1->Fit("pol2","er","",-100,100);
82+
row.T0offsetO = -9.36725e-03 + 5.74947e-03 + 1.140144e-01 - 0.73383303057411764 + 0.60752;
83+
tableSet->AddAt(&row);
84+
// ----------------- end of code ---------------
85+
return (TDataSet *)tableSet;
86+
}

0 commit comments

Comments
 (0)
⚡