Istruzioni per l'analisi in fortran
Le ntuple necessarie (.hbook generati da CMSJET) si possono trovare su pccmsfarm01, in /cmsraid2/giammanc/top. Quelle attualmente utilizzate sono:
- ewtop_s.hbook (20k single top s-channel -> 1l)
- ewtop_t.hbook (20k single top t-channel -> 1l)
- wt.hbook (20k single top Wt-channel -> 1l)
- wt_2l.hbook (20k single top Wt-channel -> 2l)
- ttbar.hbook (200k ttbar ufficiali -> 1l)
- tt_2l.hbook (20k ttbar -> 2l)
- wbb.hbook (20k Wbb -> 1l)
- wcc.hbook (10k Wcc -> 1l)
- ww_incl (20k WW inclusivi)
- zw_incl (20k ZW inclusivi)
- zz_incl (20k ZZ inclusivi)
Gli altri fondi ho visto che sembrano trascurabili gia' a livello di preselezione e per ora non li ho inclusi nell'analisi completa...
Per eseguire il programma hai bisogno di:
Per girare l'analisi su tutti i sample e avere una tabella degli eventi che passano trigger, preselezione e selezione completa, riscalati a 30 fb-1:
Questo script fa un loop in cui modifica il file nomefile.inc che e' incluso in readcwn.inc e chiama readcwn.csh per ogni sample; alla fine del loop esegue sn.f che legge gli output di readcwn.f e calcola Signal/Background (per le tre diverse definizioni di "signal": single top, s-channel, t-channel).
Parametri in readcwn.inc:
- selectionmode=1,2,3: serve a decidere quale analisi compiere (1: single top, 2: s-channel, 3: t-channel). In pratica regola quale dei cuts_*.inc sara' incluso.
- mrecoswitch=1,2,3: per la ricostruzione del top nell'evento (massa e vettore impulso) occorre ricostruire il neutrino. Per fare cio' uso l'impulso trasverso mancante per le componenti trasverse dell'impulso del neutrino, ed estraggo la componente longitudinale dell'impulso sfruttando l'impulso del leptone e imponendo il vincolo della massa del W. Questo pero' da' due soluzioni. 1: scelgo la prima soluzione (che ci azzecca il 50% delle volte; utile per confronti e debug), 2: scelgo la soluzione che da' il minimo Delta R tra neutrino e leptone, 3: scelgo la soluzione che da' la massa ricostruita del top piu' vicina a 175 GeV.
- mrecoswitch2=0,1: escludendo il caso del single top in t-channel, in cui il b associato al top tende a essere collineare al gluone padre e quindi e' quasi sempre troppo in avanti per essere rivelato, i b nell'evento sono due. In preselezione quindi richiediamo due jet. Quale dei due utilizzare nella ricostruzione del top? 0: il jet piu' b-taggato (questo criterio va benissimo per il t-channel ma non per l's-channel), 1: il jet che da' la massa ricostruita del top piu' vicina a 175 GeV.
- mrecoswitch3=0,1: 0: utilizza mrecoswitch e mrecoswitch2, 1: ignorali, e delle 4 possibili combinazioni di soluzioni neutriniche e jet utilizza quella che da' la massa ricostruita del top piu' vicina a 175 GeV.
- bestworst=1,2: 1: chiama "jet 1" il jet piu' b-taggato, 2: il contrario
- preselswitch=0,1,2: 1: preselezione normale, 0: nessuna preselezione, 2: preselezione alternativa di test
- dotrigger=.false.,.true.: applica o no la selezione di trigger
- lumi=.false.,.true.: decide quali tagli usare per il trigger: .false.: bassa luminosita', .true.: alta luminosita'
Tagli in cuts_*.inc:
- btagcut: taglio (inferiore) sulla probabilita' di b-tag per il jet 1 (vedi bestworst=1,2 per la numerazione dei jet)
- btagcut2: taglio (inferiore) sulla probabilita' di b-tag per il jet 2; serve per arricchire il campione in s-channel
- antibtagcut: taglio (superiore) sulla probabilita' di b-tag per il jet 2; serve per arricchire il campione in t-channel
- ptjetveto: sotto i 20 GeV di pT i jet non dovrebbero essere ricostruibili, ma CMSJET li ricostruisce lo stesso; questo taglio e' solo un trucco per non considerarli nel conteggio (e non va cambiato)
- ptjetcut: taglio (inferiore) sul pT dei jet
- ptjetcut_b: taglio (inferiore) sul pT dei jet b-taggati
- ptelecut: taglio (inferiore) sul pT degli elettroni
- ptmucut: taglio (inferiore) sul pT dei muoni
- mreco_low,mreco_high: finestra di massa ricostruita del top
- deltamreco_high: taglio (superiore) su M(lnj1)-M(lnj2); serve a rigettare gli eventi t-channel quando si cerca di arricchire in s-channel (il taglio complementare per arricchire in t-channel non occorre, e' gia' abbastanza facile anche senza questa variabile)
- ht_low,ht_high: finestra di Ht (somma degli impulsi trasversi di tutti gli oggetti: leptoni, jet, mancante)
- ptmiss_low,ptmiss_high: finestra di energia mancante
- pj1_low,pj1_high: finestra di impulso del jet 1
- pj2_low,pj2_high: finestra di impulso del jet 2
- etanonb_high: taglio (superiore) di pseudorapidita' del jet 2
- rbnonb_low,rbnonb_high: finestra di Delta R tra i due jet (per discriminazione s-channel/t-channel)
- mj1j2_low: taglio (inferiore) sulla massa invariante della coppia di jet. Utile per rigettare fondi in cui due b-jet vengono da una Z.
- ptlep_in_jet: taglio (superiore) sul pT di un leptone contenuto in un jet (i leptoni da decadimento del b superano raramente i 5 GeV; leptoni con impulsi di molto superiori probabilmente provengono da W o Z)
Per vedere "in faccia" queste variabili, si puo' usare top2.kumac. Per istruzioni sul suo uso, exe top2 help.
Per l'analisi di polarizzazione (dopo aver selezionato t-channel): polar.kumac.
Per la risoluzione in impulso ed energia del top: ptopresol.kumac.
Per il confronto tra l'uso del jet piu' e meno b-taggato (per il Delta R tra lui e il leptone, o per la massa ricostruita del top): bnonb.kumac.