Skip to content

Guannan's tune up for 7p7GeV_2021, https://www.star.bnl.gov/~xgn1992/… #505

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

Merged
starsdong merged 2 commits into from
Mar 22, 2023
Merged

Guannan's tune up for 7p7GeV_2021, https://www.star.bnl.gov/~xgn1992/… #505

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
85c44e2
Guannan's tune up for 7p7GeV_2021, https://www.star.bnl.gov/~xgn1992/…
fisyak Mar 19, 2023
122024d
Clean up tpcRDOMasks, move deactivate for them in MySQL from 2023-03-…
fisyak Mar 21, 2023
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
227 changes: 153 additions & 74 deletions StRoot/StTpcRSMaker/StTpcRSMaker.cxx
View file
Open in desktop

Large diffs are not rendered by default.

10 changes: 6 additions & 4 deletions StRoot/StTpcRSMaker/StTpcRSMaker.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -55,8 +55,8 @@ class StTpcRSMaker : public StMaker {
TF1F *GetChargeFraction(Int_t io = 0, Int_t sector = 20) {return (TF1F *) mChargeFraction[io][sector-1];}
TF1F *GetPadResponseFunction(Int_t io = 0, Int_t sector = 20){return (TF1F *) mPadResponseFunction[io][sector-1];}
TF1F *GetPolya(Int_t io = 0) {return (TF1F *) mPolya[io];}
TF1F *GetTimeShape0(Int_t io = 0) {return fgTimeShape0[io];}
TF1F *GetTimeShape3(Int_t io = 0) {return fgTimeShape3[io];}
TF1 *GetTimeShape0(Int_t io = 0) {return fgTimeShape0[io];}
TF1 *GetTimeShape3(Int_t io = 0) {return fgTimeShape3[io];}
TF1 *GetHeed() {return mHeed;}
Double_t GetNoPrimaryClusters(Double_t betaGamma, Int_t charge);
virtual void Print(Option_t *option="") const;
Expand All @@ -70,6 +70,8 @@ class StTpcRSMaker : public StMaker {
static Double_t shapeEI3(Double_t *x, Double_t *par=0);
static Double_t shapeEI3_I(Double_t *x, Double_t *par=0);
static Double_t fei(Double_t t, Double_t t0, Double_t T);
static Double_t feiFunc(Double_t *x, Double_t *par) {return fei(x[0],par[0],par[1]);}
TF1 *Fei();
static Double_t polya(Double_t *x, Double_t *par);
SignalSum_t *GetSignalSum(Int_t sector);
SignalSum_t *ResetSignalSum(Int_t sector);
Expand All @@ -93,8 +95,8 @@ class StTpcRSMaker : public StMaker {
void GenerateSignal(HitPoint_t &TrackSegmentHits, Int_t sector, Int_t rowMin, Int_t rowMax, Double_t sigmaJitterT, Double_t sigmaJitterX);
Double_t dEdxCorrection(HitPoint_t &TrackSegmentHits);
#endif
static TF1F *fgTimeShape3[2]; //!
static TF1F *fgTimeShape0[2]; //!
static TF1 *fgTimeShape3[2]; //!
static TF1 *fgTimeShape0[2]; //!
Char_t beg[1]; //!
TTree *fTree; //!
SignalSum_t *m_SignalSum; //!
Expand Down
28 changes: 28 additions & 0 deletions StRoot/StTpcRSMaker/TpcRS_tuneup.doc
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,28 @@
y2009 TpcRS
TpcRS_pp200Lo trans diffusion from y2008
TpcRS_pp200Lo_GT adjust gain, OmegaTauScaleX = 1e3, 0.62*trasf.diffusion,
TpcRS_pp200Lo_GnoT -"- , trans diff. back to Mogboltz
TpcRS_pp200Lo_GT2 tras diff * 0.92
TpcRS_pp200Lo_GT3 + adjust Jitter
GT4N0 + adjust Jitter
GT5N0 + adjust Jitter
GT6N0 + adjust Jitter
GT7N0 + adjust Jitter
TpcRS_pp200Lo_GT7N0x2 + tauXx/2
TpcRS_pp200Lo_GT7N0X2 + tauXx*2

01/25/2023 ================================================================================
static const Double_t xmaxD = 100;// XC11 , XC15
// static const Double_t xmaxD = 20;// XC12
// static const Double_t xmaxD = 50;// XC13
// static const Double_t xmaxD = 35;// XC14
TpcRS_pp200Lo_Y: static Double_t tmax = 6e-6; // Y
TpcRS_pp200Lo_Y10.5 : row.longitudinalDiffusion = 0.03624*0.5;
TpcRS_pp200Lo_Z: static Double_t tmax = 2e-6; // Z
TpcRS_pp200Lo_Zx10: row.tauX? = 10e-9
TpcRS_pp200Lo_Zx10B: T0offset and SigmaJitterTX : 0.25 => 0.50
TpcRS_pp200Lo_Zx200: row.tauXI = 200e-9, restore T0offset

The best description is TpcRS_pp200Lo_Z
TpcRS_pp200Lo_ZZ: restore row.tauX?, static Double_t tmax = 1e-6; // ZZ
TpcRS_pp200Lo_ZZ1.5: static Double_t tmax = 1.5e-6; // ZZ1.5
0 StRoot/StTpcRSMaker/TF1F.cxx → StRoot/StarRoot/TF1F.cxx
View file
Open in desktop
File renamed without changes.
0 StRoot/StTpcRSMaker/TF1F.h → StRoot/StarRoot/TF1F.h
View file
Open in desktop
File renamed without changes.
115 changes: 115 additions & 0 deletions StarDb/Calibrations/tpc/TpcResponseSimulator.y2008.C
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,115 @@
// $Id: TpcResponseSimulator.y2009.C,v 1.13 2014/03/16 20:17:36 fisyak Exp $
// $Log: TpcResponseSimulator.y2009.C,v $
// Revision 1.13 2014/03/16 20:17:36 fisyak
// New set of drift velocities for AuAu15 with fixed RHIC clock
//
// Revision 1.12 2013/02/06 14:34:17 fisyak
// Set T0offset from Run 9 pp 200 GeV embedding
//
// Revision 1.11 2012/04/03 14:06:55 fisyak
// Speed up using GetSaveL (__PAD_BLOCK__), sluggish shape histograms, Heed electron generation
//
// Revision 1.6 2012/01/18 13:57:57 fisyak
// Adjust T0offset : Xianglei Zhu from Run 11 AuAu 27 & 19.6 GeV embedding
//
// Revision 1.5 2011/12/30 00:04:01 fisyak
// Freeze parameters for y2011 base on TpcRS_2011_pp500LowLum_Q
//
// Revision 1.3 2011/10/11 19:09:23 fisyak
// Add Yi Guo's tables for Run XI AuAu200 RFF dE/dx calibration
//
// Revision 1.2 2010/10/28 23:41:54 fisyak
// extra t0 off set for Altro chip
//
// Revision 1.7 2010/06/14 23:36:08 fisyak
// Freeze version V
//
// Revision 1.6 2010/05/24 21:39:53 fisyak
// Fix bracket
//
// Revision 1.5 2010/05/24 16:07:20 fisyak
// Add default dE/dx calibration tables, replace TpcAltroParameters and asic_thresholds_tpx by tpcAltroParams
//
// Revision 1.4 2010/04/19 15:05:58 fisyak
// Final (2010_i) parameters for Run X
//
// Revision 1.3 2010/04/16 19:31:19 fisyak
// Intermidiate version
//
// Revision 1.3 2010/04/04 23:14:33 fisyak
// Add Row Correction
//
// Revision 1.2 2010/04/01 22:17:57 fisyak
// Freeze version W
//
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
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;// i 0.25; // f 0.38;// 0.4317;// 0.25;//ad 0.0;// b for Tpx inner
row.SigmaJitterTO = 0;// i 0.25; // f 0.38;// 0.4300;// E: 0.4801;//0.25;//ad 0.0;// b for Tpx outer
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.;
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.;
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.transverseDiffusion = 0.02218*TMath::Sqrt(1 + row.OmegaTau*row.OmegaTau)*0.62 ; // Magboltz / 1.62
row.NoElPerAdc = 335.; // No. of electrons per 1 ADC count
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
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
// 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 )
row.SecRowCorIW[0] = row.SecRowCorIE[0] = 6.44083315017355607e-01 + 4.48625e-02;
row.SecRowCorOW[0] = row.SecRowCorOE[0] = 1.12486177000493281e+00 - 1.45850e-01;
/* SecRow3CGFTpcRS_2009_pp200_a - SecRow3CGFdaq_2009_pp200 */
/* SecRow3CGFTpcRS_2005_CuCu22_a : FitP->Draw("sigma:y>>s(45,0.5,45.5)","(i&&j)/dsigma**2","profg")
SecRow3CGFdaq_2009_pp200 : Inner : 3.82589e-01, Outer: 2.97251e-01
SecRow3CGFTpcRS_2009_pp200_a: Inner : 3.40720e-01, Outer: 2.71217e-01
diff : 1.74024e-01, 1.21653e-01
SecRow3CGFTpcRS_2009_pp200_f: 3.66156e-01 2.85218e-01
SecRow3CGFTpcRS_2009_pp200_g: 3.45247e-01 2.26954e-01
*/
const Double_t RowSigmaTrs[4] = {//
1.74024e-01, 0.,
1.21653e-01, 0.};
Float_t *b = &row.SecRowSigIW[0];
for (Int_t i = 0; i < 8; i++) {
b[i] = RowSigmaTrs[i%4];
}
row.PolyaInner = 1.38;
row.PolyaOuter = 1.38;
// 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.25; // 02/06/13 Xianglei Zhu from Run 9 pp 200 GeV embedding
row.T0offset = 0.25 - 0.56; // 03/08/14 Lukasz Fulek
tableSet->AddAt(&row);
// ----------------- end of code ---------------
return (TDataSet *)tableSet;
}
11 changes: 0 additions & 11 deletions StarDb/RunLog/onl/starClockOnl.devE.C
View file
Open in desktop

This file was deleted.

73 changes: 0 additions & 73 deletions StarDb/RunLog/onl/tpcRDOMasks.20071128.13225.C
View file
Open in desktop

This file was deleted.

73 changes: 0 additions & 73 deletions StarDb/RunLog/onl/tpcRDOMasks.20200625.121203.C
View file
Open in desktop

This file was deleted.

Loading
⚡