|
| 1 | +// $Id: TpcResponseSimulator.y2009.C,v 1.13 2014/03/16 20:17:36 fisyak Exp $ |
| 2 | +// $Log: TpcResponseSimulator.y2009.C,v $ |
| 3 | +// Revision 1.13 2014/03/16 20:17:36 fisyak |
| 4 | +// New set of drift velocities for AuAu15 with fixed RHIC clock |
| 5 | +// |
| 6 | +// Revision 1.12 2013/02/06 14:34:17 fisyak |
| 7 | +// Set T0offset from Run 9 pp 200 GeV embedding |
| 8 | +// |
| 9 | +// Revision 1.11 2012/04/03 14:06:55 fisyak |
| 10 | +// Speed up using GetSaveL (__PAD_BLOCK__), sluggish shape histograms, Heed electron generation |
| 11 | +// |
| 12 | +// Revision 1.6 2012/01/18 13:57:57 fisyak |
| 13 | +// Adjust T0offset : Xianglei Zhu from Run 11 AuAu 27 & 19.6 GeV embedding |
| 14 | +// |
| 15 | +// Revision 1.5 2011/12/30 00:04:01 fisyak |
| 16 | +// Freeze parameters for y2011 base on TpcRS_2011_pp500LowLum_Q |
| 17 | +// |
| 18 | +// Revision 1.3 2011/10/11 19:09:23 fisyak |
| 19 | +// Add Yi Guo's tables for Run XI AuAu200 RFF dE/dx calibration |
| 20 | +// |
| 21 | +// Revision 1.2 2010/10/28 23:41:54 fisyak |
| 22 | +// extra t0 off set for Altro chip |
| 23 | +// |
| 24 | +// Revision 1.7 2010/06/14 23:36:08 fisyak |
| 25 | +// Freeze version V |
| 26 | +// |
| 27 | +// Revision 1.6 2010/05/24 21:39:53 fisyak |
| 28 | +// Fix bracket |
| 29 | +// |
| 30 | +// Revision 1.5 2010/05/24 16:07:20 fisyak |
| 31 | +// Add default dE/dx calibration tables, replace TpcAltroParameters and asic_thresholds_tpx by tpcAltroParams |
| 32 | +// |
| 33 | +// Revision 1.4 2010/04/19 15:05:58 fisyak |
| 34 | +// Final (2010_i) parameters for Run X |
| 35 | +// |
| 36 | +// Revision 1.3 2010/04/16 19:31:19 fisyak |
| 37 | +// Intermidiate version |
| 38 | +// |
| 39 | +// Revision 1.3 2010/04/04 23:14:33 fisyak |
| 40 | +// Add Row Correction |
| 41 | +// |
| 42 | +// Revision 1.2 2010/04/01 22:17:57 fisyak |
| 43 | +// Freeze version W |
| 44 | +// |
| 45 | +TDataSet *CreateTable() { |
| 46 | + // ----------------------------------------------------------------- |
| 47 | + // db/.const/StarDb/Calibrations/tpc/.TpcResponseSimulator/TpcResponseSimulator Allocated rows: 1 Used rows: 1 Row size: 124 bytes |
| 48 | + // Table: TpcResponseSimulator_st[0]--> TpcResponseSimulator_st[0] |
| 49 | + // ==================================================================== |
| 50 | + // ------ Test whether this table share library was loaded ------ |
| 51 | + if (!TClass::GetClass("St_TpcResponseSimulator")) return 0; |
| 52 | + TpcResponseSimulator_st row; |
| 53 | + St_TpcResponseSimulator *tableSet = new St_TpcResponseSimulator("TpcResponseSimulator",1); |
| 54 | + memset(&row,0,tableSet->GetRowSize()); |
| 55 | + row.I0 = 13.1;// eV, CH4 |
| 56 | + row.Cluster = 3.2; // average no. of electrons per primary |
| 57 | + row.W = 26.2;// eV |
| 58 | + row.OmegaTau = 3.02;// fit of data |
| 59 | + row.K3IP = 0.68;//(pads) for a/s = 2.5e-3 and h/s = 0.5 |
| 60 | + row.K3IR = 0.89;//(row) for a/s = 2.5e-3 and h/s = 0.5 |
| 61 | + row.K3OP = 0.55;//(pads) for a/s = 2.5e-3 and h/s = 1.0 |
| 62 | + row.K3OR = 0.61;//(row) for a/s = 2.5e-3 and h/s = 1.0 |
| 63 | + row.FanoFactor = 0.3; // |
| 64 | + row.AveragePedestal = 50.0;// |
| 65 | + row.AveragePedestalRMS = 1.4; // Old Tpc electronics |
| 66 | + row.AveragePedestalRMSX = 0.7; // New Tpx electronics |
| 67 | + row.tauIntegration = 2.5*74.6e-9;// secs |
| 68 | + row.tauF = 394.0e-9;// secs Tpc |
| 69 | + row.tauP = 775.0e-9;// secs Tpc |
| 70 | + row.tauXI = 60.0e-9;// secs Tpx Inner integration time |
| 71 | + row.tauXO = 74.6e-9;// secs Tpx Outer integration time |
| 72 | + row.tauCI = 0; |
| 73 | + row.tauCO = 0; |
| 74 | + row.SigmaJitterTI = 0;// i 0.25; // f 0.38;// 0.4317;// 0.25;//ad 0.0;// b for Tpx inner |
| 75 | + row.SigmaJitterTO = 0;// i 0.25; // f 0.38;// 0.4300;// E: 0.4801;//0.25;//ad 0.0;// b for Tpx outer |
| 76 | + row.SigmaJitterXI = 0.735*0.16;// 0.1027785; // P: 0.1353*1.05/1.10; //O: 0.1353*1.05;// N: 0.1353; // C:0.; |
| 77 | + row.SigmaJitterXO = 0.701*0.30;// 0.107525; // P: 0.1472*1.05/1.03; //O: 0.1472*1.05;// N: 0.1472; // C:0.; |
| 78 | + row.longitudinalDiffusion = 0.03624; // Magboltz // HD 0.03624*1.5; //HC 0.03624; // Magboltz |
| 79 | + // row.transverseDiffusion = 0.02218*TMath::Sqrt(1 + row.OmegaTau*row.OmegaTau) ; // Magboltz |
| 80 | + row.transverseDiffusion = 0.02218*TMath::Sqrt(1 + row.OmegaTau*row.OmegaTau)*0.62 ; // Magboltz / 1.62 |
| 81 | + row.NoElPerAdc = 335.; // No. of electrons per 1 ADC count |
| 82 | + row.OmegaTauScaleI = 1e3;// 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 |
| 83 | + row.OmegaTauScaleO = 1e3;// 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 |
| 84 | + // Inner_wire_to_plane_coupling ( 0.533 ) * Inner_wire_to_plane_couplingScale ( 0.843485 ) |
| 85 | + // Outer_wire_to_plane_coupling ( 0.512 ) * Outer_wire_to_plane_couplingScale ( 0.725267 ) |
| 86 | + row.SecRowCorIW[0] = row.SecRowCorIE[0] = 6.44083315017355607e-01 + 4.48625e-02; |
| 87 | + row.SecRowCorOW[0] = row.SecRowCorOE[0] = 1.12486177000493281e+00 - 1.45850e-01; |
| 88 | + /* SecRow3CGFTpcRS_2009_pp200_a - SecRow3CGFdaq_2009_pp200 */ |
| 89 | + /* SecRow3CGFTpcRS_2005_CuCu22_a : FitP->Draw("sigma:y>>s(45,0.5,45.5)","(i&&j)/dsigma**2","profg") |
| 90 | + SecRow3CGFdaq_2009_pp200 : Inner : 3.82589e-01, Outer: 2.97251e-01 |
| 91 | + SecRow3CGFTpcRS_2009_pp200_a: Inner : 3.40720e-01, Outer: 2.71217e-01 |
| 92 | + diff : 1.74024e-01, 1.21653e-01 |
| 93 | + SecRow3CGFTpcRS_2009_pp200_f: 3.66156e-01 2.85218e-01 |
| 94 | + SecRow3CGFTpcRS_2009_pp200_g: 3.45247e-01 2.26954e-01 |
| 95 | + */ |
| 96 | + const Double_t RowSigmaTrs[4] = {// |
| 97 | + 1.74024e-01, 0., |
| 98 | + 1.21653e-01, 0.}; |
| 99 | + Float_t *b = &row.SecRowSigIW[0]; |
| 100 | + for (Int_t i = 0; i < 8; i++) { |
| 101 | + b[i] = RowSigmaTrs[i%4]; |
| 102 | + } |
| 103 | + row.PolyaInner = 1.38; |
| 104 | + row.PolyaOuter = 1.38; |
| 105 | + // row.T0offset = 0.50; // From Lokesh Kumar for Run X |
| 106 | + // 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") |
| 107 | + // 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") |
| 108 | + // The corection has to be added M P |
| 109 | + //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 |
| 110 | + // row.T0offset = 0.25; // 02/06/13 Xianglei Zhu from Run 9 pp 200 GeV embedding |
| 111 | + row.T0offset = 0.25 - 0.56; // 03/08/14 Lukasz Fulek |
| 112 | + tableSet->AddAt(&row); |
| 113 | + // ----------------- end of code --------------- |
| 114 | + return (TDataSet *)tableSet; |
| 115 | +} |
0 commit comments