search

ladamczy / STAR-Analysis

Repository for AGH-STAR analysis code
2 stars 0 forks source link

Prezentacja do Warszawy #20

Open ladamczy opened 1 month ago

ladamczy commented 1 month ago

ploty z pracy Pana Fulka które można zrobić już:

selekcja: 1 werteks przynajmniej dwa ślady pT>0.2 |eta|<0.7 (n_sel>=2) ciecie anty-elastyczne

do Warszawy chyba wystrczy tlko tryger bez BBCL veta (aby nie komplikować spraw)

3.4 tylko dane 3.5 tylko dane 3.6 dane + pythia8 3.7 tylko pythia8 (i w funkcji xi1xi2) 3.18(g) dane + pythia8 (i w funkcji xi1xi2) tylko n_sel>=2 3.19(a) dane + pythia8 (i w funkcji xi1*xi2) tylko n_sel>=2 3.20(a) dane + pythia8 tylko n_sel>=2 3.21(a) dane + pythia8 tylko n_sel>=2 3.60 3.62 bez fitów

do tego rozkład masy niezmienniczej, pt, eta kandydatów na K0_S, Lambda, anty-Lambda (to Pan ma)

i podobne ploty dla kandydatów na phi, K0*

i prosze przygotwac slajdy z opisem fizyki i motywacji:

ladamczy commented 1 month ago

pytania Pana Adama

Mam kilka pytań doprecyzowujących:

ladamczy commented 1 month ago

I dalsze

lus niektóre wykresy mają cząstki z selekcją którą nie jestem w stanie powtórzyć (np. jedno z cięć wymaga aby stosunek N_fit do wszystkich możliwych dla danego śladu był wyższy od 0.52, czego nie jestem w stanie policzyć w tak krótkim terminie), czy powinienem je po prostu pominąć?

ladamczy commented 1 month ago

I dalej Czy ta wartosc N_fit do wszystkich mozliwych nie jest dostepna po prostu jako jeden z parametrow sladu? Odp: To jest, ale ilość możliwych dla danego śladu już nie.

ladamczy commented 1 month ago

I dalej Czy ta wartosc N_fit do wszystkich mozliwych nie jest dostepna po prostu jako jeden z parametrow sladu? Odp: To jest, ale ilość możliwych dla danego śladu już nie.

N_fit pewnie nie ma updDST. Niech Pan selekcje sladów zostawi tak jak Pan miał do tej pory. Referencja do pracy P. Fulka była aby Pan ogarniał całość tego co trzeba zrobić ale nie powtarzał literalnie tej samej selekcji. Analiza 510 jest oparta o upcDST i pewnyce rzeczy uległy zmianie,

ladamczy commented 1 month ago 
* tryger bez veta w starych danych brać z ograniczeniem  runów czy bez

nie rozumiem pytania. Czy chodzi o to że tryger bez weta na początku był bez skalowania a póżniej ze skalowaniem? Jesli tak to wystarczy ten pierwszy kores bez skalowania.

* czy werteks brać z funkcji getNPrimVerices() czy getNumberOfVertices() (pamiętam że była różnica między nimi, bo poerwsza dawała wcześniej obliczoną liczbę, a druga po prostu sprawdza ile jest obecnie werteksów, ale nie pamiętam szczegółów)(poza tym i tak będę musiał brać do V0findera pozycję werteksu z listy, więc chyba druga opcja lepsza)

getNumberOfVertices()

* do normowania 1/N brać liczbę eventów sprzed czy już po selekcji (zakładam że po, ale wolę się upewnić)

po selekcji

* czy "tylko n_sel" oznacza że mam pominąć absolutnie wszystkie wymagania w danych oprócz tych do 2 śladów, czy chodzi o zastosowanie go w symulacji PYTHIA?

"tylko n_sel" to znaczy że ze strony TPC-ToF żadamy tylko n_sel>=2 . Natomiast oprócz tego normalne ciecia trygerowe i zakres kinemtyczny protonów + weto anty-elastyczne

* jako że w nowszych danych jest flaga kV0 i możliwość przypadkowej duplikacji danych, zakładam że te wszystkie płoty mam robić na starszych danych

Przypadkowa duplikacja sladów jest mało istotna. Istotne jest ze w nowych danych ma Pan preselekcje kV0. Doszliśmy juz do ustaleń że uzywamy starych danych (tylko ploty z K0, Lambda sa z nowych danych)

* zakładam że do śladów mam używać tych samych ograniczeń co w pracy Pana Fulka (w stylu miał napisane że używa N_dEdx, N_fit itp.)

Teraz to jest obojetne. Niech Pan używa swoich obecnych cięć ale jak Pan dodał dEdx to niech zostanie.

AddahDraggon commented 1 month ago

Wykresy wyszczególnione pownienem mieć na wtorek, póki co zrobiłem na szybko wykresy żeby sprawdzić czy faktycznie widać Kstar i phi:

AnaOutput_Inclusive_analysis_Kstar_phi_old_data.zip

Normalnie nie widać Kstar, a na masie par K+K- odzywa się inny rezonans - dopiero po identyfikacji cząstek dE/dx da się je znaleźć (na wykresie K+K- aż 4 piki, w tym phi(1020), w związku z tym czy jest jakiegoś rodzaju repozytorium które pozwala wyszukiwać cząstki po kanałach rozpadu? PDG nie jest tu zbyt pomocne).

Czy wykresy 3.4-3.6 mają być w postaci n-1 czy po prostu po cięciach na protony?

ladamczy commented 1 month ago

Jak tylko się da to wykresy robimy w posatci n-1 , ale jak jest problem to prosze teraz nie tracic na to czasu. I tak wszsytkich wykresów pan nie pokaże. Nie wiem nic na temat narzędzia które wypisze Panu wszsytkie czatki rozpadajace is na dany stan końcowy.

Jesli chodzi o KK to tam mamy phi ( lub f0(980) ) problem jest podobny jak dla analizy ekskluywnej. Jesli Pan zrobi histogram o wiekszej ilosci binów i dofituje cos co pasuje do phi to bedzie to dowód. Oprócz tego spodziewamy się f2(1270) i f0(1500) . Dodatkowo może pan widzieć coś co rozpada sie nie na KK tylko np. pipi (K0). Robiąc plot KK sugeruje dać veto na K0S ( w sensie masa pary pipip poza oknem K0S).

Jk Pan będzie miał wykresy do trygera bez weta BBL to może Pan jeszcze dodać rany z trygerem BBCL.

AddahDraggon commented 1 month ago

Napisałem wstępnie kod do powyższych wykresów, w środę przed spotkaniem zrobię analogiczną część do symulacji i je połączę. Nie stosowałem cięć n-1, jak będzie czas to przerobię. Póki co wyniki:

presentation_pdfs.zip

Po zastosowaniu veto na K0S wykresy bez identyfikacji cząstek się nie poprawiły, pojawił się za to nowy pik na wykresie Kpi w okolicach 0.73, po identyfikacji przeniósł się na 0.70.

Obecnie wszystkie ślady centralne (po upewnieniu się że n_sel jest spełniony tam gdzie żądany) są brane bez żadnych ograniczeń oprócz sygnału w ToF.

AddahDraggon commented 1 month ago

Zrobiłem ok. 90% treści prezentacji, brakujące wykresy i niezidentyfikowane rezonanse uzupełnię dzisiaj lub najpóźniej jutro. Prosiłbym o zerknięcie na prezentację.

Mam też kilka pytań:

ladamczy commented 1 month ago

Zrobiłem ok. 90% treści prezentacji, brakujące wykresy i niezidentyfikowane rezonanse uzupełnię dzisiaj lub najpóźniej jutro. Prosiłbym o zerknięcie na prezentację.

Mam też kilka pytań:

* Czy mógłbym prosić o zerknięcie na n_sel w danych Pythii, bo nie jestem pewien czy tak powinien wyglądać wykres wydajności n_sel>=2 ([link do miejsca brania n_sel](https://github.com/ladamczy/STAR-Analysis/blob/e4b1ec578cf434e857671d1465fee4c075d39c46/Inclusive_Analysis/pythia/src/Central_diffraction_presentation.cxx#L209))

nie rozumem dlaczego jest ten warunek TPCParticlesForNsel==abs(chargeForNsel)

Odrzuca Pan rzypadki gdy wszystkie ślady sa tego samego znaku. Dlaczego?

* czy dobrym pomysłem byłoby umieszczenie informacji o innych zmiennych na wykresach, jak w pracy Łukasza Fulka?

generalnie tak, ale z drugiej strony nie może tam być zbyt dużo informacji.

* czy obecne cięcia na K0 i Lambdę jeśli chodzi o V0 (hypoDecayLength itp.) są nadal aktualne w przypadku Kstar i phi(1020)? Pytanie głównie w kwestii tych wykresów różniczkowych.

Nie. Kstar i phi rozpadają się natychmiast w wierzchołku pierwotnym . Do ich rekonstrukcji nie korzystamy z V0. Tworzymy wszystkie pary z wierzchołka pierwotnego i sprawdzamy masę niezmiennicza.

AddahDraggon commented 1 month ago

nie rozumem dlaczego jest ten warunek TPCParticlesForNsel==abs(chargeForNsel)

Odrzuca Pan rzypadki gdy wszystkie ślady sa tego samego znaku. Dlaczego?

To jes najprostsza metoda zagwarantowania że w śladach znajdzie się co najmniej jedna para o przeciwnych znakach - do poszukiwań Delta++ faktycznie może być nieprzydatne, ale póki co tego nie szukamy

ladamczy commented 1 month ago

nie rozumem dlaczego jest ten warunek TPCParticlesForNsel==abs(chargeForNsel) Odrzuca Pan rzypadki gdy wszystkie ślady sa tego samego znaku. Dlaczego?

To jes najprostsza metoda zagwarantowania że w śladach znajdzie się co najmniej jedna para o przeciwnych znakach - do poszukiwań Delta++ faktycznie może być nieprzydatne, ale póki co tego nie szukamy

Ale dlaczego chce Pan być pewny zę jest para o przeciwnych znakach? W inkluzywnej analizie nie ma to znaczenia. Patzrymy na wszystkie kombinacje znaków, ważne aby były przynajmniej dwie cząstki naładowane.