Commit b41ffd94 authored by Adrien Matta's avatar Adrien Matta
Browse files

* Progress on S034 analysis

parent c9371077
Pipeline #126510 passed with stages
in 8 minutes and 22 seconds
...@@ -28,7 +28,7 @@ void gamma(){ ...@@ -28,7 +28,7 @@ void gamma(){
chain->SetAlias("R","sqrt(Nebula.PosX*Nebula.PosX+Nebula.PosY*Nebula.PosY+(Nebula.PosZ+3774.7)*(Nebula.PosZ+3774.7))"); chain->SetAlias("R","sqrt(Nebula.PosX*Nebula.PosX+Nebula.PosY*Nebula.PosY+(Nebula.PosZ+3774.7)*(Nebula.PosZ+3774.7))");
new TCanvas(); new TCanvas();
unsigned int select =60; unsigned int select =91;
for(unsigned int i = 0 ; i < 150 ; i++){ for(unsigned int i = 0 ; i < 150 ; i++){
if(i!=select) if(i!=select)
process1bar(i); process1bar(i);
...@@ -49,16 +49,18 @@ void process1bar(int b){ ...@@ -49,16 +49,18 @@ void process1bar(int b){
double R = r1->GetBinCenter(r1->GetMaximumBin()); double R = r1->GetBinCenter(r1->GetMaximumBin());
auto h1 = GetV(b); auto h1 = GetV(b);
if(h1->GetEntries()<10) if(h1->GetEntries()<1)
return; return;
//h1->Rebin(8); //h1->Rebin(8);
double max = h1->GetBinCenter(h1->GetMaximumBin()); double max = h1->GetBinCenter(h1->GetMaximumBin());
//h1->Draw(); //h1->Draw();
auto f = new TF1("f","gaus(0)+pol0(3)",max-100,max+100); auto f = new TF1("f","crystalball(0)+pol0(5)",max-100,max+100);
f->SetParameter(0,h1->GetMaximum()); f->SetParameter(0,h1->GetMaximum());
f->SetParameter(1,max); f->SetParameter(1,max);
f->SetParameter(2,5); f->SetParameter(2,10);
f->SetParameter(3,5); f->SetParameter(3,1);
f->SetParameter(4,5);
f->SetParameter(5,5);
h1->Fit(f,"R"); h1->Fit(f,"R");
...@@ -68,10 +70,12 @@ void process1bar(int b){ ...@@ -68,10 +70,12 @@ void process1bar(int b){
// X=R*(1/c-1/Vbad) // X=R*(1/c-1/Vbad)
double offset=R*(1/c_light-1/f->GetParameter(1)) ; double offset=R*(1/c_light-1/f->GetParameter(1)) ;
cout << "hello " << max-f->GetParameter(1) << endl; cout << "hello " << max-f->GetParameter(1) << endl;
//double offset=R*(1/c_light-1/max) ; //double offset=R*(1/c_light-1/max) ;
cout <<f->GetParameter(1) << " " << offset << " " << R/(offset+R/f->GetParameter(1)) << endl; cout <<f->GetParameter(1) << " " << offset << " " << R/(offset+R/f->GetParameter(1)) << endl;
if(offset>0){ if(offset>0){
output << "NEBULA_T_ID" << b << " " << offset << endl; output << "NEBULA_T_ID" << b << " " << offset << endl;
graph->Set(graph->GetN()+1); graph->Set(graph->GetN()+1);
......
NEBULA_T_ID1 34.5922 NEBULA_T_ID1 34.9114
NEBULA_T_ID2 34.7694 NEBULA_T_ID2 34.8487
NEBULA_T_ID3 34.9723 NEBULA_T_ID3 35.0828
NEBULA_T_ID4 34.6604 NEBULA_T_ID4 34.7753
NEBULA_T_ID5 34.8851 NEBULA_T_ID5 35.1731
NEBULA_T_ID6 34.6936 NEBULA_T_ID6 34.8492
NEBULA_T_ID7 34.7044 NEBULA_T_ID7 34.8129
NEBULA_T_ID8 34.7642 NEBULA_T_ID8 35.1259
NEBULA_T_ID9 34.4816 NEBULA_T_ID9 34.957
NEBULA_T_ID10 34.5612 NEBULA_T_ID10 34.9146
NEBULA_T_ID11 34.6275 NEBULA_T_ID11 34.9044
NEBULA_T_ID12 34.619 NEBULA_T_ID12 34.9218
NEBULA_T_ID13 34.3853 NEBULA_T_ID13 34.6699
NEBULA_T_ID14 34.3388 NEBULA_T_ID14 35.1466
NEBULA_T_ID15 35.0631 NEBULA_T_ID15 35.0836
NEBULA_T_ID16 34.5556 NEBULA_T_ID16 34.6919
NEBULA_T_ID17 33.9783 NEBULA_T_ID17 34.796
NEBULA_T_ID18 34.2425 NEBULA_T_ID18 34.3545
NEBULA_T_ID19 34.5992 NEBULA_T_ID19 34.6873
NEBULA_T_ID20 34.2835 NEBULA_T_ID20 34.5355
NEBULA_T_ID21 34.4126 NEBULA_T_ID21 34.5093
NEBULA_T_ID22 34.4358 NEBULA_T_ID22 34.5631
NEBULA_T_ID23 34.0272 NEBULA_T_ID23 34.6092
NEBULA_T_ID24 34.1103 NEBULA_T_ID24 34.6205
NEBULA_T_ID25 34.5739 NEBULA_T_ID25 34.6199
NEBULA_T_ID26 35.0359 NEBULA_T_ID26 35.161
NEBULA_T_ID27 34.9156 NEBULA_T_ID27 35.1816
NEBULA_T_ID28 34.2105 NEBULA_T_ID28 35.0658
NEBULA_T_ID29 34.8138 NEBULA_T_ID29 35.0873
NEBULA_T_ID30 34.9965 NEBULA_T_ID30 35.1039
NEBULA_T_ID31 34.1339 NEBULA_T_ID31 34.6067
NEBULA_T_ID32 34.4452 NEBULA_T_ID32 34.6076
NEBULA_T_ID33 34.0956 NEBULA_T_ID33 34.862
NEBULA_T_ID34 34.3118 NEBULA_T_ID34 34.4877
NEBULA_T_ID35 34.3546 NEBULA_T_ID35 34.4975
NEBULA_T_ID36 34.3213 NEBULA_T_ID36 34.4971
NEBULA_T_ID37 34.7201 NEBULA_T_ID37 34.7477
NEBULA_T_ID38 34.2632 NEBULA_T_ID38 34.2818
NEBULA_T_ID39 34.4773 NEBULA_T_ID39 34.6409
NEBULA_T_ID40 33.719 NEBULA_T_ID40 34.2638
NEBULA_T_ID41 34.5495 NEBULA_T_ID41 34.5886
NEBULA_T_ID42 34.3425 NEBULA_T_ID42 34.4999
NEBULA_T_ID43 34.2744 NEBULA_T_ID43 35.0054
NEBULA_T_ID44 34.6795 NEBULA_T_ID44 34.7799
NEBULA_T_ID45 34.1254 NEBULA_T_ID45 34.7292
NEBULA_T_ID46 34.5549 NEBULA_T_ID46 34.5549
NEBULA_T_ID47 34.4596 NEBULA_T_ID47 34.6511
NEBULA_T_ID48 34.5931 NEBULA_T_ID48 34.9767
NEBULA_T_ID49 34.2254 NEBULA_T_ID49 34.4121
NEBULA_T_ID50 34.2657 NEBULA_T_ID50 34.7922
NEBULA_T_ID51 34.5111 NEBULA_T_ID51 34.5915
NEBULA_T_ID52 34.581 NEBULA_T_ID52 34.6975
NEBULA_T_ID53 34.099 NEBULA_T_ID53 35.0689
NEBULA_T_ID54 34.0757 NEBULA_T_ID54 34.7948
NEBULA_T_ID55 34.6373 NEBULA_T_ID55 34.6392
NEBULA_T_ID56 34.6463 NEBULA_T_ID56 34.7534
NEBULA_T_ID57 34.8271 NEBULA_T_ID57 35.1284
NEBULA_T_ID58 34.5796 NEBULA_T_ID58 35.4314
NEBULA_T_ID59 35.2348 NEBULA_T_ID59 35.3455
NEBULA_T_ID61 34.1897 NEBULA_T_ID60 35.3207
NEBULA_T_ID62 34.2115 NEBULA_T_ID61 34.2079
NEBULA_T_ID63 35.6809 NEBULA_T_ID62 34.3236
NEBULA_T_ID64 34.6327 NEBULA_T_ID63 35.7209
NEBULA_T_ID65 35.0425 NEBULA_T_ID64 34.7389
NEBULA_T_ID66 35.4113 NEBULA_T_ID65 35.0687
NEBULA_T_ID67 34.9083 NEBULA_T_ID66 35.4752
NEBULA_T_ID68 34.1289 NEBULA_T_ID67 35.0603
NEBULA_T_ID69 33.5752 NEBULA_T_ID68 34.2005
NEBULA_T_ID70 34.4984 NEBULA_T_ID69 33.6966
NEBULA_T_ID71 34.2975 NEBULA_T_ID70 34.8085
NEBULA_T_ID72 33.7079 NEBULA_T_ID71 34.4027
NEBULA_T_ID73 34.2989 NEBULA_T_ID72 34.1083
NEBULA_T_ID74 34.3563 NEBULA_T_ID73 34.7039
NEBULA_T_ID75 33.5771 NEBULA_T_ID74 34.4356
NEBULA_T_ID76 34.7965 NEBULA_T_ID75 34.3117
NEBULA_T_ID77 34.7209 NEBULA_T_ID76 35.085
NEBULA_T_ID78 33.7591 NEBULA_T_ID77 35.0264
NEBULA_T_ID79 34.059 NEBULA_T_ID78 33.7459
NEBULA_T_ID80 34.8741 NEBULA_T_ID79 34.3353
NEBULA_T_ID82 33.0554 NEBULA_T_ID80 35.0204
NEBULA_T_ID83 34.4457 NEBULA_T_ID81 34.8225
NEBULA_T_ID84 34.4195 NEBULA_T_ID82 33.1864
NEBULA_T_ID85 33.5292 NEBULA_T_ID83 34.8137
NEBULA_T_ID86 34.8708 NEBULA_T_ID84 34.404
NEBULA_T_ID87 33.9424 NEBULA_T_ID85 33.5315
NEBULA_T_ID88 33.1291 NEBULA_T_ID86 34.8705
NEBULA_T_ID89 34.668 NEBULA_T_ID87 34.3997
NEBULA_T_ID90 33.7628 NEBULA_T_ID88 33.2526
NEBULA_T_ID91 34.5667 NEBULA_T_ID89 34.7098
NEBULA_T_ID92 32.1052 NEBULA_T_ID90 34.3216
NEBULA_T_ID93 34.6439 NEBULA_T_ID92 31.9321
NEBULA_T_ID94 34.174 NEBULA_T_ID93 34.6501
NEBULA_T_ID95 34.4052 NEBULA_T_ID94 34.3427
NEBULA_T_ID96 34.777 NEBULA_T_ID95 34.6861
NEBULA_T_ID97 34.8321 NEBULA_T_ID96 34.8366
NEBULA_T_ID98 34.5154 NEBULA_T_ID97 34.832
NEBULA_T_ID99 34.8771 NEBULA_T_ID98 35.1193
NEBULA_T_ID100 34.82 NEBULA_T_ID99 35.172
NEBULA_T_ID101 34.3954 NEBULA_T_ID100 35.275
NEBULA_T_ID102 34.5721 NEBULA_T_ID101 34.459
NEBULA_T_ID103 34.6479 NEBULA_T_ID102 34.5771
NEBULA_T_ID104 34.5427 NEBULA_T_ID103 34.926
NEBULA_T_ID105 35.1569 NEBULA_T_ID104 34.548
NEBULA_T_ID106 34.8748 NEBULA_T_ID105 35.1097
NEBULA_T_ID107 35.0694 NEBULA_T_ID106 35.1887
NEBULA_T_ID108 33.8616 NEBULA_T_ID107 35.0866
NEBULA_T_ID109 34.5262 NEBULA_T_ID108 34.6665
NEBULA_T_ID110 33.4922 NEBULA_T_ID109 34.5351
NEBULA_T_ID111 34.367 NEBULA_T_ID110 33.65
NEBULA_T_ID112 34.8241 NEBULA_T_ID111 34.8446
NEBULA_T_ID113 33.6154 NEBULA_T_ID112 34.9961
NEBULA_T_ID114 35.2176 NEBULA_T_ID113 34.6419
NEBULA_T_ID115 33.723 NEBULA_T_ID114 35.2164
NEBULA_T_ID116 40.2125 NEBULA_T_ID115 33.7334
NEBULA_T_ID117 35.4705 NEBULA_T_ID116 34.3578
NEBULA_T_ID118 34.2625 NEBULA_T_ID117 35.7166
NEBULA_T_ID118 34.5191
NEBULA_T_ID119 35.471 NEBULA_T_ID119 35.471
NEBULA_T_ID120 34.8941 NEBULA_T_ID120 35.0177
NEBULA_T_ID121 33.6064 NEBULA_T_ID121 33.6312
NEBULA_T_ID122 34.3692 NEBULA_T_ID122 34.3798
NEBULA_T_ID123 32.8653 NEBULA_T_ID123 34.6003
NEBULA_T_ID124 34.6004 NEBULA_T_ID124 34.3961
NEBULA_T_ID125 35.4875 NEBULA_T_ID125 35.5018
NEBULA_T_ID126 23.867 NEBULA_T_ID126 23.4914
NEBULA_T_ID127 33.0042 NEBULA_T_ID127 34.5968
NEBULA_T_ID128 32.6078 NEBULA_T_ID128 32.6766
NEBULA_T_ID129 34.4132 NEBULA_T_ID129 34.45
NEBULA_T_ID130 34.1434 NEBULA_T_ID130 34.1432
NEBULA_T_ID131 23.8805 NEBULA_T_ID131 24.0027
NEBULA_T_ID132 32.6514 NEBULA_T_ID132 32.7579
NEBULA_T_ID133 35.5018 NEBULA_T_ID133 35.3537
NEBULA_T_ID134 34.633 NEBULA_T_ID134 34.55
NEBULA_T_ID135 36.0812 NEBULA_T_ID135 36.1924
NEBULA_T_ID136 35.0311 NEBULA_T_ID136 34.9019
NEBULA_T_ID137 32.5883 NEBULA_T_ID137 32.7547
NEBULA_T_ID138 25.9258 NEBULA_T_ID138 25.6774
NEBULA_T_ID139 28.8419 NEBULA_T_ID139 28.3729
NEBULA_T_ID140 25.3851 NEBULA_T_ID140 19.0246
NEBULA_T_ID141 32.6386 NEBULA_T_ID141 32.7578
NEBULA_T_ID142 28.5499 NEBULA_T_ID142 28.4096
NEBULA_T_ID143 29.1687 NEBULA_T_ID143 28.9315
NEBULA_T_ID144 60.234 NEBULA_T_ID144 31.2423
NEBULA_T_ID60 34.7372 NEBULA_T_ID91 34.6
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment