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

* Progress on Nebula Claibration

parent 952b282e
Pipeline #125107 passed with stages
in 3 minutes and 27 seconds
...@@ -31,11 +31,11 @@ void gamma(){ ...@@ -31,11 +31,11 @@ void gamma(){
h->Draw("colz"); h->Draw("colz");
new TCanvas(); new TCanvas();
for(unsigned int i = 0 ; i < 160 ; i++){ for(unsigned int i = 0 ; i < 160 ; i++){
if(i!=60) if(i!=61)
process1bar(i); process1bar(i);
} }
process1bar(60); process1bar(61);
output.close(); output.close();
new TCanvas(); new TCanvas();
graph->Draw("ap"); graph->Draw("ap");
...@@ -49,7 +49,7 @@ void process1bar(int b){ ...@@ -49,7 +49,7 @@ void process1bar(int b){
double R = r1->GetBinCenter(r1->GetMaximumBin()); double R = r1->GetBinCenter(r1->GetMaximumBin());
auto h1 = h->ProjectionY(Form("h%d",b),b,b+1); auto h1 = h->ProjectionY(Form("h%d",b),b,b+1);
//h1->Rebin(4); h1->Rebin(4);
double max = h1->GetBinCenter(h1->GetMaximumBin()); double max = h1->GetBinCenter(h1->GetMaximumBin());
//h1->Draw(); //h1->Draw();
auto f = new TF1("f","crystalball(0)+pol1(5)",max-50,max+50); auto f = new TF1("f","crystalball(0)+pol1(5)",max-50,max+50);
...@@ -66,8 +66,8 @@ void process1bar(int b){ ...@@ -66,8 +66,8 @@ void process1bar(int b){
// 1/Vbad +X/R = 1/c // 1/Vbad +X/R = 1/c
// 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)) ;
// double offset=R*(1/c_light-1/max) ; double offset=R*(1/c_light-1/h1->GetMean()) ;
cout <<f->GetParameter(1) << " " << offset << " " << R/(f->GetParameter(1)-offset) << endl; cout <<f->GetParameter(1) << " " << offset << " " << R/(f->GetParameter(1)-offset) << endl;
if(offset>0){ if(offset>0){
......
NEBULA_T_ID1 35.0208 NEBULA_T_ID1 33.4905
NEBULA_T_ID2 34.9058 NEBULA_T_ID2 33.4387
NEBULA_T_ID3 34.8484 NEBULA_T_ID3 33.2579
NEBULA_T_ID4 35.0285 NEBULA_T_ID4 33.1961
NEBULA_T_ID5 34.5798 NEBULA_T_ID5 33.3351
NEBULA_T_ID6 35.1383 NEBULA_T_ID6 33.4268
NEBULA_T_ID7 34.6951 NEBULA_T_ID7 33.3499
NEBULA_T_ID8 34.6655 NEBULA_T_ID8 33.2878
NEBULA_T_ID9 35.1679 NEBULA_T_ID9 33.2271
NEBULA_T_ID10 34.902 NEBULA_T_ID10 33.1101
NEBULA_T_ID11 34.929 NEBULA_T_ID11 33.1321
NEBULA_T_ID12 34.929 NEBULA_T_ID12 33.1597
NEBULA_T_ID13 34.8411 NEBULA_T_ID13 33.1065
NEBULA_T_ID14 34.9144 NEBULA_T_ID14 33.1347
NEBULA_T_ID15 35.1963 NEBULA_T_ID15 33.2256
NEBULA_T_ID16 34.8331 NEBULA_T_ID16 33.146
NEBULA_T_ID17 34.6084 NEBULA_T_ID17 33.0473
NEBULA_T_ID18 34.4628 NEBULA_T_ID18 32.9851
NEBULA_T_ID19 34.3097 NEBULA_T_ID19 32.9804
NEBULA_T_ID20 34.6923 NEBULA_T_ID20 33.0581
NEBULA_T_ID21 34.5083 NEBULA_T_ID21 33.0919
NEBULA_T_ID22 34.5421 NEBULA_T_ID22 33.0102
NEBULA_T_ID23 34.7158 NEBULA_T_ID23 32.9368
NEBULA_T_ID24 34.641 NEBULA_T_ID24 32.8589
NEBULA_T_ID25 34.2835 NEBULA_T_ID25 32.8442
NEBULA_T_ID26 34.8876 NEBULA_T_ID26 33.0195
NEBULA_T_ID27 35.1372 NEBULA_T_ID27 33.1863
NEBULA_T_ID28 35.1372 NEBULA_T_ID28 33.187
NEBULA_T_ID29 35.227 NEBULA_T_ID29 33.3282
NEBULA_T_ID30 35.0475 NEBULA_T_ID30 33.195
NEBULA_T_ID31 34.5307 NEBULA_T_ID31 32.978
NEBULA_T_ID32 34.4025 NEBULA_T_ID32 33.3618
NEBULA_T_ID33 34.6834 NEBULA_T_ID33 33.3447
NEBULA_T_ID34 34.6565 NEBULA_T_ID34 33.3058
NEBULA_T_ID35 34.4795 NEBULA_T_ID35 33.3803
NEBULA_T_ID36 34.478 NEBULA_T_ID36 33.4223
NEBULA_T_ID37 34.6938 NEBULA_T_ID37 33.5491
NEBULA_T_ID38 34.2743 NEBULA_T_ID38 33.2656
NEBULA_T_ID39 34.4684 NEBULA_T_ID39 33.1152
NEBULA_T_ID40 34.4441 NEBULA_T_ID40 33.1433
NEBULA_T_ID41 34.4249 NEBULA_T_ID41 33.192
NEBULA_T_ID42 34.4866 NEBULA_T_ID42 33.442
NEBULA_T_ID43 34.646 NEBULA_T_ID43 33.5688
NEBULA_T_ID44 34.9125 NEBULA_T_ID44 33.5181
NEBULA_T_ID45 34.5875 NEBULA_T_ID45 33.34
NEBULA_T_ID46 34.914 NEBULA_T_ID46 33.214
NEBULA_T_ID47 34.6781 NEBULA_T_ID47 33.3359
NEBULA_T_ID48 34.8383 NEBULA_T_ID48 33.4933
NEBULA_T_ID49 34.5358 NEBULA_T_ID49 33.2248
NEBULA_T_ID50 34.5403 NEBULA_T_ID50 33.0666
NEBULA_T_ID51 34.5465 NEBULA_T_ID51 33.0149
NEBULA_T_ID52 34.7831 NEBULA_T_ID52 33.111
NEBULA_T_ID53 34.8545 NEBULA_T_ID53 33.4055
NEBULA_T_ID54 34.6405 NEBULA_T_ID54 33.3783
NEBULA_T_ID55 34.8021 NEBULA_T_ID55 33.263
NEBULA_T_ID56 35.1102 NEBULA_T_ID56 33.0668
NEBULA_T_ID57 35.1386 NEBULA_T_ID57 33.1642
NEBULA_T_ID58 35.1951 NEBULA_T_ID58 33.4036
NEBULA_T_ID59 35.3143 NEBULA_T_ID59 33.3579
NEBULA_T_ID61 32.9368 NEBULA_T_ID60 33.2143
NEBULA_T_ID62 34.1741 NEBULA_T_ID62 34.5329
NEBULA_T_ID63 34.189 NEBULA_T_ID63 34.9431
NEBULA_T_ID64 35.8693 NEBULA_T_ID64 35.1015
NEBULA_T_ID65 35.0106 NEBULA_T_ID65 35.1265
NEBULA_T_ID66 35.3523 NEBULA_T_ID66 35.1442
NEBULA_T_ID67 35.3218 NEBULA_T_ID67 34.9116
NEBULA_T_ID68 34.805 NEBULA_T_ID68 34.6604
NEBULA_T_ID69 33.8813 NEBULA_T_ID69 34.4462
NEBULA_T_ID70 33.8717 NEBULA_T_ID70 34.5221
NEBULA_T_ID71 34.3202 NEBULA_T_ID71 35.0051
NEBULA_T_ID72 34.3202 NEBULA_T_ID72 35.1518
NEBULA_T_ID73 33.3849 NEBULA_T_ID73 34.9795
NEBULA_T_ID74 34.5386 NEBULA_T_ID74 34.941
NEBULA_T_ID75 34.3927 NEBULA_T_ID75 34.864
NEBULA_T_ID76 34.9837 NEBULA_T_ID76 35.2727
NEBULA_T_ID77 35.0867 NEBULA_T_ID77 35.4467
NEBULA_T_ID78 35.0796 NEBULA_T_ID78 35.1904
NEBULA_T_ID79 33.2896 NEBULA_T_ID79 35.2455
NEBULA_T_ID80 34.8329 NEBULA_T_ID80 35.4184
NEBULA_T_ID81 34.7558 NEBULA_T_ID81 35.4121
NEBULA_T_ID82 37.2023 NEBULA_T_ID82 35.1133
NEBULA_T_ID83 34.8499 NEBULA_T_ID83 35.1475
NEBULA_T_ID84 34.8686 NEBULA_T_ID84 35.1741
NEBULA_T_ID85 33.5995 NEBULA_T_ID85 34.7787
NEBULA_T_ID86 33.632 NEBULA_T_ID86 34.914
NEBULA_T_ID87 34.3618 NEBULA_T_ID87 34.9607
NEBULA_T_ID88 33.7309 NEBULA_T_ID88 34.5766
NEBULA_T_ID89 33.4773 NEBULA_T_ID89 34.5742
NEBULA_T_ID90 34.3177 NEBULA_T_ID90 34.7277
NEBULA_T_ID91 34.336 NEBULA_T_ID91 34.6957
NEBULA_T_ID92 34.2217 NEBULA_T_ID92 35.0481
NEBULA_T_ID93 27.3398 NEBULA_T_ID93 35.2177
NEBULA_T_ID94 34.6431 NEBULA_T_ID94 35.4574
NEBULA_T_ID95 34.8253 NEBULA_T_ID95 35.4506
NEBULA_T_ID96 34.7884 NEBULA_T_ID96 35.5051
NEBULA_T_ID97 34.9366 NEBULA_T_ID97 35.8906
NEBULA_T_ID98 34.9366 NEBULA_T_ID98 36.0723
NEBULA_T_ID99 35.0008 NEBULA_T_ID99 35.617
NEBULA_T_ID100 35.0008 NEBULA_T_ID100 35.4569
NEBULA_T_ID101 34.8941 NEBULA_T_ID101 35.5729
NEBULA_T_ID102 34.3539 NEBULA_T_ID102 35.3768
NEBULA_T_ID103 34.6346 NEBULA_T_ID103 35.482
NEBULA_T_ID104 34.6135 NEBULA_T_ID104 35.3516
NEBULA_T_ID105 34.7392 NEBULA_T_ID105 35.4782
NEBULA_T_ID106 34.8739 NEBULA_T_ID106 35.6773
NEBULA_T_ID107 35.2196 NEBULA_T_ID107 35.7758
NEBULA_T_ID108 34.8644 NEBULA_T_ID108 35.6255
NEBULA_T_ID109 34.6989 NEBULA_T_ID109 35.3739
NEBULA_T_ID110 33.8876 NEBULA_T_ID110 35.141
NEBULA_T_ID111 33.9751 NEBULA_T_ID111 35.2823
NEBULA_T_ID112 34.7086 NEBULA_T_ID112 35.8331
NEBULA_T_ID113 34.7238 NEBULA_T_ID113 35.3862
NEBULA_T_ID114 34.6405 NEBULA_T_ID114 34.9962
NEBULA_T_ID115 34.0107 NEBULA_T_ID115 35.2729
NEBULA_T_ID116 34.3824 NEBULA_T_ID116 35.3663
NEBULA_T_ID117 34.1178 NEBULA_T_ID117 35.2767
NEBULA_T_ID118 34.3166 NEBULA_T_ID118 35.1023
NEBULA_T_ID119 34.9824 NEBULA_T_ID119 35.1875
NEBULA_T_ID120 35.7801 NEBULA_T_ID120 35.46
NEBULA_T_ID121 33.2615 NEBULA_T_ID121 32.8085
NEBULA_T_ID122 33.6554 NEBULA_T_ID122 32.7149
NEBULA_T_ID123 34.4064 NEBULA_T_ID123 32.7373
NEBULA_T_ID124 34.6484 NEBULA_T_ID124 32.7586
NEBULA_T_ID125 34.67 NEBULA_T_ID125 32.7595
NEBULA_T_ID126 34.6457 NEBULA_T_ID126 32.5096
NEBULA_T_ID127 34.5139 NEBULA_T_ID127 32.5093
NEBULA_T_ID128 34.4616 NEBULA_T_ID128 32.7319
NEBULA_T_ID129 34.3498 NEBULA_T_ID129 32.7193
NEBULA_T_ID130 34.2412 NEBULA_T_ID130 32.6077
NEBULA_T_ID131 34.0883 NEBULA_T_ID131 32.5304
NEBULA_T_ID132 33.48 NEBULA_T_ID132 32.2799
NEBULA_T_ID133 33.0478 NEBULA_T_ID133 32.5294
NEBULA_T_ID134 34.7691 NEBULA_T_ID134 34.6905
NEBULA_T_ID135 34.4737 NEBULA_T_ID135 34.8101
NEBULA_T_ID136 34.7838 NEBULA_T_ID136 35.0837
NEBULA_T_ID137 35.9547 NEBULA_T_ID137 35.0707
NEBULA_T_ID138 35.7866 NEBULA_T_ID138 34.4088
NEBULA_T_ID139 35.6697 NEBULA_T_ID139 34.235
NEBULA_T_ID140 34.5767 NEBULA_T_ID140 34.2653
NEBULA_T_ID141 35.8165 NEBULA_T_ID141 34.3418
NEBULA_T_ID142 33.4558 NEBULA_T_ID142 34.6938
NEBULA_T_ID143 35.1309 NEBULA_T_ID143 34.4649
NEBULA_T_ID144 34.5652 NEBULA_T_ID144 34.6477
NEBULA_T_ID145 21.315 NEBULA_T_ID145 34.8407
NEBULA_T_ID60 35.3134 NEBULA_T_ID61 33.1146
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