Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

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

Merged
merged 2 commits into from
Nov 2, 2021
Merged

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

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
7807df9
Xianglei Zhu, 09/11/2021: FXT3p85_2018 and 09/14/2021: Isobar paramet…
fisyak Sep 27, 2021
89aee7c
Update StarDb/Calibrations/tpc/TpcResponseSimulator.20171220.000001.C
plexoos Oct 9, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
83 changes: 83 additions & 0 deletions StarDb/Calibrations/tpc/TpcResponseSimulator.20171220.000001.C
View file
Open in desktop
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
Open in desktop
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