Skip to content

Commit

Permalink
Xianglei Zhu, 09/11/2021: FXT3p85_2018 and 09/14/2021: Isobar paramet…
Browse files Browse the repository at this point in the history 
…ers (#156)

* Xianglei Zhu, 09/11/2021: FXT3p85_2018 and 09/14/2021: Isobar parameters parameters

* Update StarDb/Calibrations/tpc/TpcResponseSimulator.20171220.000001.C

Co-authored-by: Dmitri Smirnov <dmixsmi@gmail.com>
  • Loading branch information
@fisyak @plexoos @jml985
2 people authored and jml985 committed Dec 7, 2021
1 parent 6faed27 commit a0a2bf9
Show file tree
Hide file tree
Showing 2 changed files with 91 additions and 6 deletions.
83 changes: 83 additions & 0 deletions StarDb/Calibrations/tpc/TpcResponseSimulator.20171220.000001.C
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,83 @@
TDataSet *CreateTable() {
// -----------------------------------------------------------------
// db/.const/StarDb/Calibrations/tpc/.TpcResponseSimulator/TpcResponseSimulator Allocated rows: 1 Used rows: 1 Row size: 124 bytes
// Table: TpcResponseSimulator_st[0]--> TpcResponseSimulator_st[0]
// ====================================================================
// ------ Test whether this table share library was loaded ------
if (!TClass::GetClass("St_TpcResponseSimulator")) return 0;
TpcResponseSimulator_st row;
St_TpcResponseSimulator *tableSet = new St_TpcResponseSimulator("TpcResponseSimulator",1);
memset(&row,0,tableSet->GetRowSize());
row.I0 = 13.1;// eV, CH4
row.Cluster = 3.2; // average no. of electrons per primary
row.W = 26.2;// eV
row.OmegaTau = 3.02;// fit of data
row.K3IP = 0.68;//(pads) for a/s = 2.5e-3 and h/s = 0.5
row.K3IR = 0.89;//(row) for a/s = 2.5e-3 and h/s = 0.5
row.K3OP = 0.55;//(pads) for a/s = 2.5e-3 and h/s = 1.0
row.K3OR = 0.61;//(row) for a/s = 2.5e-3 and h/s = 1.0
row.FanoFactor = 0.3; //
row.AveragePedestal = 50.0;//
row.AveragePedestalRMS = 1.4; // Old Tpc electronics or iTPC
row.AveragePedestalRMSX = 0.7; // New Tpx electronics
row.tauIntegration = 2.5*74.6e-9;// secs
row.tauF = 394.0e-9;// secs Tpc
row.tauP = 775.0e-9;// secs Tpc
row.tauXI = 60.0e-9;// secs Tpx Inner integration time
row.tauXO = 74.6e-9;// secs Tpx Outer integration time
row.tauCI = 0;
row.tauCO = 0;
row.SigmaJitterTI = 0.53;// 0.4317;// 0.25;//ad 0.0;// b for Tpx inner
row.SigmaJitterTO = 0.33;// 0.4300;// E: 0.4801;//0.25;//ad 0.0;// b for Tpx outer
row.SigmaJitterXI = 0;// 0.1027785; // P: 0.1353*1.05/1.10; //O: 0.1353*1.05;// N: 0.1353; // C:0.;
row.SigmaJitterXO = 0;// 0.107525; // P: 0.1472*1.05/1.03; //O: 0.1472*1.05;// N: 0.1472; // C:0.;
row.longitudinalDiffusion = 0.03624; // Magboltz // HD 0.03624*1.5; //HC 0.03624; // Magboltz
row.transverseDiffusion = 0.02218*TMath::Sqrt(1 + row.OmegaTau*row.OmegaTau) ; // Magboltz
row.NoElPerAdcI = 335./TMath::Exp(-0.81625); // No. of electrons per 1 ADC count for outer TPX
row.NoElPerAdcO = 335./TMath::Exp(-1.09801); // No. of electrons per 1 ADC count for outer TPX
row.NoElPerAdcX = 335./TMath::Exp(-0.49463);; // No. of electrons per 1 ADC count for iTPC
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
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
// Inner_wire_to_plane_coupling ( 0.533 ) * Inner_wire_to_plane_couplingScale ( 0.843485 )
// Outer_wire_to_plane_coupling ( 0.512 ) * Outer_wire_to_plane_couplingScale ( 0.725267 )
#if 1
row.SecRowCorIW[0] = row.SecRowCorIE[0] = 0.039-0.2+0.04;//- TMath::Log(0.533*0.843485) -5.84129e-01 + 4.52885e-01 + 3.09117e-02; // IRAKLI
row.SecRowCorOW[0] = row.SecRowCorOE[0] = 0.105-0.2+0.04;//- TMath::Log(0.512*0.725267) -5.47141e-01 + 5.23937e-01 + 1.19154e-02; // IRAKLI
#endif
// SecRow3CGFdaq_2011_pp500LowLum => Inner: 3.26428e-01 - -5.01720e-04*y; Outer: 2.68883e-01 + 1.23403e-04*y
// 3.22907e-01 2.72715e-01
// SecRow3CGFTpcRS_2011_pp500LowLum_f : 3.09711e-01 2.65342e-01
// diff : 9.13675e-02 6.29849e-02
// SecRow3CGFTpcRS_2011_pp500LowLum_g : 3.12857e-01 2.67379e-01
#if 0
const Double_t RowSigmaTrs[4] = {
9.13675e-02, 0, // Inner
6.29849e-02, 0}; // Outer
Float_t *b = &row.SecRowSigIW[0];
for (Int_t i = 0; i < 8; i++) {
b[i] = RowSigmaTrs[i%4];
}

/* Sigmas
Tpx inner = 0.395
Tpx outer = 0.314
iTpc = 0.322 Outer 0.316

row.PolyaInner = 1.38;
row.PolyaOuter = 1.38;
#endif
// row.T0offset = 0.50; // From Lokesh Kumar for Run X
// 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")
// 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")
// The corection has to be added M P
//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
row.T0offset = 0.50 -1.43663e-01;//g // 01/18/12 Xianglei Zhu from Run 11 AuAu 27 & 19.6 GeV embedding
// root.exe T0offset.C
// TI->FitSlicesY(); TI_1->Fit("pol2","er","",-100,100);
row.T0offsetI = 1.17437e-01 + 8.43584e-03;
// TO->FitSlicesY(); TO_1->Fit("pol2","er","",-100,100);
row.T0offsetO = -9.36725e-03 + 5.74947e-03;
tableSet->AddAt(&row);
// ----------------- end of code ---------------
return (TDataSet *)tableSet;
}
14 changes: 8 additions & 6 deletions StarDb/Calibrations/tpc/TpcResponseSimulator.20180531.071409.C
View file
Original file line number Diff line number Diff line change
Expand Up @@ -97,12 +97,14 @@ TDataSet *CreateTable() {
row.tauIntegration = 2.5*74.6e-9;// secs
row.tauF = 394.0e-9;// secs Tpc
row.tauP = 775.0e-9;// secs Tpc
//row.tauXI = 80.6e-9;// secs Tpx Inner integration time
//row.tauXO = 88.6e-9;// secs Tpx Outer integration time
row.tauXI = 60.0e-9;// secs Tpx Inner integration time
row.tauXO = 74.6e-9;// secs Tpx Outer integration time
row.tauCI = 0;
row.tauCO = 0;
row.SigmaJitterTI = 0;// 0.4317;// 0.25;//ad 0.0;// b for Tpx inner
row.SigmaJitterTO = 0;// 0.4300;// E: 0.4801;//0.25;//ad 0.0;// b for Tpx outer
row.SigmaJitterTI = 0.5300;// 0.4317;// 0.25;//ad 0.0;// b for Tpx inner
row.SigmaJitterTO = 0.3300;// 0.4300;// E: 0.4801;//0.25;//ad 0.0;// b for Tpx outer
row.SigmaJitterXI = 0;// 0.1027785; // P: 0.1353*1.05/1.10; //O: 0.1353*1.05;// N: 0.1353; // C:0.;
row.SigmaJitterXO = 0;// 0.107525; // P: 0.1472*1.05/1.03; //O: 0.1472*1.05;// N: 0.1472; // C:0.;
row.longitudinalDiffusion = 0.03624; // Magboltz // HD 0.03624*1.5; //HC 0.03624; // Magboltz
Expand All @@ -118,10 +120,10 @@ TDataSet *CreateTable() {
// row.SecRowCorOW[0] = row.SecRowCorOE[0] = 0.11; // IRAKLI : based on shift seen in the MuDst based simulation;
#if 1
const Double_t RowMuTrs[8] = {
4.45505e-02, 3.63901e-03, // Inner W
3.52547e-02, 3.85726e-03, // Outer W
5.28188e-02, 3.94541e-03, // Inner E
-9.11388e-03, 2.63198e-03}; // Outer E
4.45505e-02+2.e-01-9.531e-02-1.82321556793954592e-01+0.03243-0.02, 3.63901e-03-4.693e-03, // Inner W
3.52547e-02-9.1492e-02-1.222e-01+0.1858-0.02, 3.85726e-03-4.874e-03, // Outer W
5.28188e-02+2.e-01-9.531e-02-1.82321556793954592e-01-1.84571276527970035e-01+5.79418410380801580e-02, 3.94541e-03, // Inner E
-9.11388e-03-9.1492e-02-1.222e-01+9.01563105819495891e-02+0.09564, 2.63198e-03-1.972e-03}; // Outer E
Float_t *a = &row.SecRowCorIW[0];
for (Int_t i = 0; i < 8; i++) {
a[i] = RowMuTrs[i];
Expand Down

0 comments on commit a0a2bf9

Please sign in to comment.