KADD 2020 Laboratorium 12 - Revision history

Lgraczyk: /* Kroki postępowania */

2020-05-18T14:17:47Z

Kroki postępowania

← Older revision		Revision as of 14:17, 18 May 2020
Line 6:		Line 6:
	== Kroki postępowania ==		== Kroki postępowania ==

-	1. Ściągamy plik z danymi: [http://www.if.pw.edu.pl/~lgraczyk/KADD2020/lab12/example_data.root example_data.root] i wczytujemy go w funkcji głównej makra.	+	1. Ściągamy plik z danymi: [http://www.if.pw.edu.pl/~lgraczyk/KADD2020/lab12/example_data.root example_data.root] i wczytujemy go w funkcji głównej makra.<br>
		+	Wczytywanie pliku:
		+	TFile *ifile = new TFile("file.root");
		+	TH1D hist = (TH1D)ifile->Get("hist");

	2. Do naszych danych dopasujemy taką oto funkcję zdefiniowaną jako funkcja globalna (nasz model opisujący dane doświadczalne):		2. Do naszych danych dopasujemy taką oto funkcję zdefiniowaną jako funkcja globalna (nasz model opisujący dane doświadczalne):

Lgraczyk: /* Zadanie */

2020-05-18T09:25:01Z

Zadanie

← Older revision		Revision as of 09:25, 18 May 2020
Line 2:		Line 2:
	'''Dopasowanie funkcji za pomocą paietu <code>MINUIT</code>''' (5 pkt.)		'''Dopasowanie funkcji za pomocą paietu <code>MINUIT</code>''' (5 pkt.)

-	Chcemy wykonać dopasowanie do danych doświadczalnych (wczytywanych z pliku) '''nie''' używając funkci <code>Fit</code>, tylko za pomocą pakietu <code>MINUIT</code> ([https://root.cern.ch/root/html534/TMinuit.html dokumentacja]), który daje nam zdecydowanie większą kontrolę nad procesem dopasowania funkcji.	+	Chcemy wykonać dopasowanie do danych doświadczalnych (wczytywanych z pliku) '''nie''' używając funkci <code>Fit</code>, tylko za pomocą pakietu <code>MINUIT</code> ([https://root.cern.ch/root/html534/TMinuit.html dokumentacja]), który daje nam zdecydowanie większą kontrolę nad procesem dopasowania funkcji. Przede wszystkim pakiet <code>MINUIT</code> pozwala na zdefiniowanie przez nas minimalizowanej wielkości (może to być chi-kwadrat, funkcja wiarygodności, lub cokolwiek innego co wymyślimy), i na podstawie tej minimalizacji pakiet dostarcza nam finalne parametry dopasowania.

	== Kroki postępowania ==		== Kroki postępowania ==

Lgraczyk: /* Kroki postępowania */

2020-05-18T09:17:17Z

Kroki postępowania

← Older revision		Revision as of 09:17, 18 May 2020
Line 39:		Line 39:
	<code>MINOS</code> to algorytm używany po zminimalizowaniu <code>fcn</code>, który pozwala na dokładną estymację niepwności (w tym asymetrycznych).<br>		<code>MINOS</code> to algorytm używany po zminimalizowaniu <code>fcn</code>, który pozwala na dokładną estymację niepwności (w tym asymetrycznych).<br>
	'''Uwaga!''' Nie musimy znać szczegółów obu algorytmów. Dla zainteresowanych można poczytać [https://root.cern.ch/root/htmldoc/guides/minuit2/Minuit2Letter.pdf Minuit2Letter.pdf]		'''Uwaga!''' Nie musimy znać szczegółów obu algorytmów. Dla zainteresowanych można poczytać [https://root.cern.ch/root/htmldoc/guides/minuit2/Minuit2Letter.pdf Minuit2Letter.pdf]
		+
		+	8. Za pomocą:
		+	gMinuit->GetParameter(...)
		+	wyciągamy otrzymane parametry dopasowania, tworzymy funkcję <code>TF1</code> z tymi parametrami i rysujemy wszystko na wykresie.

	== Wynik ==		== Wynik ==

Lgraczyk: /* Kroki postępowania */

2020-05-18T09:15:51Z

Kroki postępowania

← Older revision		Revision as of 09:15, 18 May 2020
Line 27:		Line 27:
	'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat. To oznacza, że '''wewnątrz''' funkcji <code>fcn</code> musimy wykonać iterację bin po binie z wczytanego pliku i dla każdego binu policzyć kwadrat różnic ( (model-wartosc_binu)^2/niepwnosc_binu^2 ). Ostateczna wartość <code>f</code>, zgodnie z definicją chi-kwadrat, to po prostu suma kwadratów różnic dla wszystkich binów.		'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat. To oznacza, że '''wewnątrz''' funkcji <code>fcn</code> musimy wykonać iterację bin po binie z wczytanego pliku i dla każdego binu policzyć kwadrat różnic ( (model-wartosc_binu)^2/niepwnosc_binu^2 ). Ostateczna wartość <code>f</code>, zgodnie z definicją chi-kwadrat, to po prostu suma kwadratów różnic dla wszystkich binów.

-	5. W funkcji głównej, mówimy naszemu obiektowi <code>gMinuit>, że ma skorzystać z funkcji FCN zdefiniowanej powyżej:	+	5. W funkcji głównej, mówimy naszemu obiektowi <code>gMinuit</code>, że ma skorzystać z funkcji FCN zdefiniowanej powyżej:
	gMinuit->SetFCN(fcn);		gMinuit->SetFCN(fcn);

Line 33:		Line 33:
	gMinuit->DefineParameter(0, "p0", 100., 1., 0., 200.);		gMinuit->DefineParameter(0, "p0", 100., 1., 0., 200.);

-	7. Główna część zadania to wykonanie ~~trzech~~ komend:	+	7. Główna część, gdzie dokonuje się zadania to wykonanie dwóch komend:
-	gMinuit->Command("MIGRAD");	+	gMinuit->Command("MIGRAD");
-	~~gMinuit->Command("MIGRAD");~~	+
	gMinuit->Command("MINOS");		gMinuit->Command("MINOS");
		+	<code>MIGRAD</code> to algorytm szukający minimum funkcji <code>fcn</code> (czyli szuka najniższej wartości liczby <code>f</code>, którą definiujemy w określony przez nas sposób).<br>
		+	<code>MINOS</code> to algorytm używany po zminimalizowaniu <code>fcn</code>, który pozwala na dokładną estymację niepwności (w tym asymetrycznych).<br>
		+	'''Uwaga!''' Nie musimy znać szczegółów obu algorytmów. Dla zainteresowanych można poczytać [https://root.cern.ch/root/htmldoc/guides/minuit2/Minuit2Letter.pdf Minuit2Letter.pdf]

	== Wynik ==		== Wynik ==

Lgraczyk: /* Kroki postępowania */

2020-05-18T09:07:26Z

Kroki postępowania

← Older revision		Revision as of 09:07, 18 May 2020
Line 20:		Line 20:

	3. W funkcji głównej makra tworzymy obiekt TMinuit, przykładowo:		3. W funkcji głównej makra tworzymy obiekt TMinuit, przykładowo:
-	TMinuit *gMinuit = new TMinuit(liczba_parametrow_fitu);	+	TMinuit *gMinuit = new TMinuit(liczba_parametrow_fitu); //liczba parametrów wynika z definicji funkcji teoretycznej

	4. Musimy stworzyć funkcję globalną (definiujemy ją przed funkcją główną), którą <code>MINUIT</code> będzie starał się zminimalizować manipulując naszymi parametrami. Musi ona mieć postać typu:		4. Musimy stworzyć funkcję globalną (definiujemy ją przed funkcją główną), którą <code>MINUIT</code> będzie starał się zminimalizować manipulując naszymi parametrami. Musi ona mieć postać typu:
Line 26:		Line 26:
	'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>		'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>
	'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat. To oznacza, że '''wewnątrz''' funkcji <code>fcn</code> musimy wykonać iterację bin po binie z wczytanego pliku i dla każdego binu policzyć kwadrat różnic ( (model-wartosc_binu)^2/niepwnosc_binu^2 ). Ostateczna wartość <code>f</code>, zgodnie z definicją chi-kwadrat, to po prostu suma kwadratów różnic dla wszystkich binów.		'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat. To oznacza, że '''wewnątrz''' funkcji <code>fcn</code> musimy wykonać iterację bin po binie z wczytanego pliku i dla każdego binu policzyć kwadrat różnic ( (model-wartosc_binu)^2/niepwnosc_binu^2 ). Ostateczna wartość <code>f</code>, zgodnie z definicją chi-kwadrat, to po prostu suma kwadratów różnic dla wszystkich binów.
		+
		+	5. W funkcji głównej, mówimy naszemu obiektowi <code>gMinuit>, że ma skorzystać z funkcji FCN zdefiniowanej powyżej:
		+	gMinuit->SetFCN(fcn);
		+
		+	6. Następnie ustawiamy parametry funkcji, np.:
		+	gMinuit->DefineParameter(0, "p0", 100., 1., 0., 200.);
		+
		+	7. Główna część zadania to wykonanie trzech komend:
		+	gMinuit->Command("MIGRAD");
		+	gMinuit->Command("MIGRAD");
		+	gMinuit->Command("MINOS");

	== Wynik ==		== Wynik ==

Lgraczyk: /* Kroki postępowania */

2020-05-18T09:04:21Z

Kroki postępowania

← Older revision		Revision as of 09:04, 18 May 2020
Line 22:		Line 22:
	TMinuit *gMinuit = new TMinuit(liczba_parametrow_fitu);		TMinuit *gMinuit = new TMinuit(liczba_parametrow_fitu);

-	4. Musimy stworzyć funkcję globalną, którą <code>MINUIT</code> będzie starał się zminimalizować manipulując naszymi parametrami. Musi ona mieć postać typu:	+	4. Musimy stworzyć funkcję globalną (definiujemy ją przed funkcją główną), którą <code>MINUIT</code> będzie starał się zminimalizować manipulując naszymi parametrami. Musi ona mieć postać typu:
	void fcn(int &npar, double gin, double &f, double par, int iflag)		void fcn(int &npar, double gin, double &f, double par, int iflag)
	'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>		'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>

Lgraczyk: /* Kroki postępowania */

2020-05-18T09:03:42Z

Kroki postępowania

← Older revision		Revision as of 09:03, 18 May 2020
Line 6:		Line 6:
	== Kroki postępowania ==		== Kroki postępowania ==

-	1. Ściągamy plik z danymi: [http://www.if.pw.edu.pl/~lgraczyk/KADD2020/lab12/example_data.root example_data.root]	+	1. Ściągamy plik z danymi: [http://www.if.pw.edu.pl/~lgraczyk/KADD2020/lab12/example_data.root example_data.root] i wczytujemy go w funkcji głównej makra.

	2. Do naszych danych dopasujemy taką oto funkcję zdefiniowaną jako funkcja globalna (nasz model opisujący dane doświadczalne):		2. Do naszych danych dopasujemy taką oto funkcję zdefiniowaną jako funkcja globalna (nasz model opisujący dane doświadczalne):
Line 19:		Line 19:
	}		}

-	3. ~~Tworzymy~~ obiekt TMinuit, przykładowo:	+	3. W funkcji głównej makra tworzymy obiekt TMinuit, przykładowo:
	TMinuit *gMinuit = new TMinuit(liczba_parametrow_fitu);		TMinuit *gMinuit = new TMinuit(liczba_parametrow_fitu);

Lgraczyk: /* Kroki postępowania */

2020-05-18T09:02:56Z

Kroki postępowania

← Older revision		Revision as of 09:02, 18 May 2020
Line 25:		Line 25:
	void fcn(int &npar, double gin, double &f, double par, int iflag)		void fcn(int &npar, double gin, double &f, double par, int iflag)
	'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>		'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>
-	'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat. To oznacza, że '''wewnątrz''' funkcji <code>fcn</code> musimy wykonać iterację bin po binie z wczytanego pliku i dla każdego binu policzyć kwadrat różnic ( (model-wartosc_binu)^2/niepwnosc_binu^2 ). Ostateczna wartość <code>f</code>, zgodnie z definicją chi-kwadrat, to po prostu suma kwadratów różnic.	+	'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat. To oznacza, że '''wewnątrz''' funkcji <code>fcn</code> musimy wykonać iterację bin po binie z wczytanego pliku i dla każdego binu policzyć kwadrat różnic ( (model-wartosc_binu)^2/niepwnosc_binu^2 ). Ostateczna wartość <code>f</code>, zgodnie z definicją chi-kwadrat, to po prostu suma kwadratów różnic dla wszystkich binów.

	== Wynik ==		== Wynik ==

Lgraczyk: /* Kroki postępowania */

2020-05-18T09:02:22Z

Kroki postępowania

← Older revision		Revision as of 09:02, 18 May 2020
Line 25:		Line 25:
	void fcn(int &npar, double gin, double &f, double par, int iflag)		void fcn(int &npar, double gin, double &f, double par, int iflag)
	'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>		'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>
-	'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat.	+	'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat. To oznacza, że '''wewnątrz''' funkcji <code>fcn</code> musimy wykonać iterację bin po binie z wczytanego pliku i dla każdego binu policzyć kwadrat różnic ( (model-wartosc_binu)^2/niepwnosc_binu^2 ). Ostateczna wartość <code>f</code>, zgodnie z definicją chi-kwadrat, to po prostu suma kwadratów różnic.

	== Wynik ==		== Wynik ==

Lgraczyk: /* Kroki postępowania */

2020-05-18T08:58:42Z

Kroki postępowania

← Older revision		Revision as of 08:58, 18 May 2020
Line 24:		Line 24:
	4. Musimy stworzyć funkcję globalną, którą <code>MINUIT</code> będzie starał się zminimalizować manipulując naszymi parametrami. Musi ona mieć postać typu:		4. Musimy stworzyć funkcję globalną, którą <code>MINUIT</code> będzie starał się zminimalizować manipulując naszymi parametrami. Musi ona mieć postać typu:
	void fcn(int &npar, double gin, double &f, double par, int iflag)		void fcn(int &npar, double gin, double &f, double par, int iflag)
-	'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).</br>	+	'''Uwaga 1!''' Nie przerażamy się parametrami (to są automatyczne parametry, które muszą być i które potrzebuje <code>MINUIT</code>, my ich w ogóle nie użyjemy).<br>
	'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat.		'''Uwaga 2!''' Jedynym parametrem, który nas interesuje to <code>double &f</code> - to jest minimalizowana wielkość. W naszym przypadku chcemy, by <code>f</code> było po prostu tożsame z chi-kwadrat.