Closed ladamczy closed 9 months ago
w plikach ROOT przesyłam R dla sub-leading kaon i zmieniłam binowanie dla pT miss, dodałam również warunek na pT(pion). yniki zostały wykonane dla nowej selekcji masy - para o najmniejszej sumie kwadratów odległości od tablicowej masy kaonu.
Rozkłady vertexu, dca córek, pTmiss , różnych konfiguracji pędu poprzecznego, N_{tof cluter} są zrobione dla wąskiego okna masy.
Prosze dodac:
W zasadzie my znamy masy tych stanów, więc moglibyśmy zamiast eta juz rysowac rapidity (y)
czy te rozkłady wypełniać pod warunkiem że obydwa kaony są w wąskim oknie masy?
Tak,
A czemu tych przypadków tak dużo? To jest chyba bez ciecia na pT_mis i n_ToF (?) Rozkład masy prosze pociagna do 2 GeV ( a histogram mieć nawet 5 GeV aby być pewnym ze dalej juz nie ma przypadków)
tak były bez cięć, tylko na wąskie okno masy
A dlaczego rozklad masy K0K0 jest obciety na 1.35 GeV? Jak on wyglada dla wiekszych mas?
Rozkłady po nałożeniu ciecia na nTOF i pTmiss
czy rozkład PtK0K0 powinien tak wyglądać?
HistEtaK0K0.pdf HistEtaLeadingK0.pdf HistEtaLeadingSubLeadingDiff.pdf HistEtaSubLeadingK0.pdf HistMassK0K0.pdf HistPtK0K0.pdf HistPtLeadingK0.pdf HistPtSubLeadingK0.pdf
To jest rozkład pT na poziomie detektorowym po wszystkich naszych cięciach. Przypuszczam, że tutaj ograniczona akceptancja protonów nam psuje oczekiwany eksponencjalny spadek. Żeby to dokładnie zrozumieć to potrzebna byłaby symulacja MC. Proszę do tego zestawu rysunków dodać rozkład PhiK0K0, i jeszcze plot 2D eta vs phi dla par K0K0.
Wygląda dobrze, poza tym że nie widać żadnych rezonansów na roizkładzie masy. Pierwsz górka na rozkładzie pT to przypadki z protonamy o przeciwnych znajach py a ta dróga do przypadki o tych samych znajach pY. Dobrze jest dodać tutaj MC ale normalizując go do danych a nie luminosity.
Do MC można dodać opis DiMe * 0.01. tak aby było te rozkłady widac na plotach. Przyczym powinno sie raczej używąc jednego skalowania dla całej próbki ale my musimy zrobić dwa skalowania jedno dla przypakdów z pypy<0. a drugie dla przypadków pypy>0
Dodam MC, tak wygląda plot z fitem dla pT miss HistPtMissData.pdf
Gdzie to MC. Dodawnie MC na rozkładach pT,_miss px_miss, py_miss ma jedynie sens jak będ tru protonu zostanie rozmyty
px_reco = gRandom->Gauss(px_true,sigma) py_reco = gRandom->Gauss(py_true,sigma)
sigma trezba dobrac aby rozkłady sie zgadzały powinno być około 0.05 GeV
Rozkład masy prosze zrobić dla szerokości binu 0.1 od 0.9
Tak wyglądają rozkłady MC. Nie rozumiem czemu MC pokazuje że K0K0 ma tylko dodatnią etę HistEtaK0K0Data.pdf HistMassK0K0Data.pdf HistPtK0K0Data.pdf
Wydaje mi sie, ze gdyby na rozkladzie masy zrobic rebin(2) to znacznie lepiej byloby widac piki tam gdzie ich oczekujemy. Prosze sprobowac to zrobic.
rozkłady pTmiss, cluster, px, py,R są wykonane dla wąskiego okna masy : HistRLeadingKaonData.pdf HistPtProtonPtKaonsX1dOppYData.pdf HistPtProtonPtKaonsX1dSameYData.pdf HistPtProtonPtKaonsY1dOppYData.pdf HistPtProtonPtKaonsY1dSameYData.pdf HistNumOfClustersData.pdf HistPtMissData.pdf
Ten rozklad masy nie wyglada lepiej... Chcialem zeby Pani zrebinowala te biny, ktore sa widoczne na oryginalnym histogramie. Czyli 0+1=1, 2+4=6, 4+6=10, itd
rozkłady pTmiss, cluster, px, py,R są wykonane dla wąskiego okna masy : HistRLeadingKaonData.pdf HistPtProtonPtKaonsX1dOppYData.pdf HistPtProtonPtKaonsX1dSameYData.pdf HistPtProtonPtKaonsY1dOppYData.pdf HistPtProtonPtKaonsY1dSameYData.pdf HistNumOfClustersData.pdf HistPtMissData.pdf
Tutaj proponuje zrebinowac wszystkie rozklady HistPt po dwa biny, ale tak, zeby to bylo symetryczne wokol zera. Na tych rysunkach opis osi poziomej nie jest jasny - proponuje cos takiego p_{T,x}^{miss,same}
Symetryczne wokol zera w sensie, ze albo zero w srodku binu, albo zero oddzielajace dwa sasiednie biny. W przypadku tej wielkosci najlepsza by byla ta druga opcja.
HistPtProtonPtKaonsX1dOppYData.pdf HistPtProtonPtKaonsX1dSameYData.pdf HistPtProtonPtKaonsY1dOppYData.pdf HistPtProtonPtKaonsY1dSameYData.pdf
To wyglada znacznie lepiej. Generalnie, jesli nie chcemy wylapac jakiejs specyficznej cechy w danym rozkladzie, a mamy jedynie do czynienia z fluktuacjami statystycznymi, to lepiej zwiekszyc szerokosc binu, tak zeby zmniejszyc bledy i wyeliminowac flutuacje statystyczne, ktore nic nie wnosza.
Tak wyglądają rozkłady MC. Nie rozumiem czemu MC pokazuje że K0K0 ma tylko dodatnią etę HistEtaK0K0Data.pdf HistMassK0K0Data.pdf HistPtK0K0Data.pdf
A tutaj raczej nie powinnam pokazywać MC, jest bardzo mało przypadków
Jeżeli chodzi o rozkład masy to zmieniłam binowanie na poziomie wypełniania: HistMassK0K0Data.pdf
Jeżeli chodzi o rozkład masy to zmieniłam binowanie na poziomie wypełniania: HistMassK0K0Data.pdf
Ten rysunek wyglada zle... Chodzilo mi o przebinowanie tego histogramu https://github.com/ladamczy/STAR-Analysis/files/12891221/HistMassK0K0Data.pdf w takim zakresie w jakim jest pokazany, czyli w nowych binach powinno byc 0+1=1, 2+4=6, 4+6=10, itd
Czy moglaby Pani dodatkowo zrobic rozkłąd masy inwariantnej czterech pionow 2Pi+2Pi- (tych z ktorych Pani probuje rekonstruowac K0), ale bez żądania żeby pośrednio dało się z nich utworzyć dwa K0 w oknie masy - po prostu rozklad masy iwariantnej czterech pionow (wszystkie inne ciecia zostaja).
zakładam osobne issue dla plotu masy(KK) bo jest zbyt duze zamieszanie.
Rozkład rapidity(KK) wygląda faktycznie dziwnie dla MC . Czy moze Pani popatrzeć na rozkład true bez cieć ?
MC możemy zawsze wyrzucić z rysunków. Zdecydujemy dzień przed prezentacją. Generuje więcej MC ale nie wiem kiedy będzie.
Powinna Pani znowu zrobić taki zbiorczy plik ze wszystkimi plotami abyśmy mogli jeszcze zdecydowac które ploty pokazać w Kairze i na takich plotach powinnsmy sie koncentrować bo czasu jest juz mało.
Powinniśmy pokazac coś z cieć topologicznych (dca między cząstakmi, R - odległość werteksu(K) od beamlinu
Można by jeszcze zrobić plot kąt w płaszczyżnie poprzecznej między wektorem od werteksu(K) do beamline a wektorem pędu poprzecznego (wektor 2D) od werteksu(K) w kierunku pędu(K).
Miałem to dołozyć do StUPCV0 ale nie mam tam informacj o beamlinie. Albo mi Pani opisze jak dostać x_beam i y_beam dla danego z i to zrobie. Albo sama to samie policzy tak jak opisałem wyżej.
już dodałam do overleaf wstępne wykresy, to może na ich podstawie będziemy decydować co zostawić/zmienić/odrzucić
Zamykam to issue i przejdziemy do nowych issues dla dodatkowych analiz
nie widze rozkladu R dla sub-leading
Troche przesadziłem , prosze w rozkładzie pt_mis przyjąć szerokość binu 50 MeV