<<O>>  Difference Topic AnalisisDeLaDensidadDeHitsEnLosFTD (r1.5 - 03 Mar 2009 - Main.iglesias)

META TOPICPARENT WebHome

Analisis de la densidad de hits en los FTD

Line: 25 to 25

Estos plots muestran bastante bien de forma visual lo que queremos hacer, se aprecia el diferente tamaño de los discos y de sus correspondientes secciones.… sin embargo, no puedo obtener ningún valor numérico porque hay variables de la ntuple no definidas. Para conseguir tener histogramas del numero de hits por evento necesitaría tener definida la variable EVENT en la ntuple.

Deleted:
<
<
Segmentacion mas adecuada para los discos
Tratamos de elegir la segmentacion mas adecuada para los discos del FTD, teniendo en cuenta que los sensores tienen que caber en una oblea de 13.5x13.5 cm^2.

  • El disco 1, tiene Rint=39mm, y Rext=164mm, lo que le confiere un area de 797.17 cm2. Un cuarto de disco (1/4), lo que tu llamas un "cuadrante" tendría un area de 199.29 cm2 y con una oblea de 13.5x13.5 cm^2 podemos abarcar como mucho: 182.25cm2, por lo que en el disco 1, nos harian falta no 4, sino 5 petalos como minimo, y como queremos que sea un numero par, mejor pondremos 6, o incluso subir hasta 8 petalos.
  • Para el disco 2, tenemos una situación similar al disco anterior, es decir, conviene poner entre 6 y 8 petalos.
  • Para el disco 3, la situación se complica, pues al tener un Rext mucho mayor (308mm) nos da un area mucho mas grande: 644.906cm2, que para abarcarla se llegan a necesitar hasta 18.2 porciones de 13.5x13.5 cm^2, por lo que se requerirían minimo 20 petalos.
  • Para el disco 4, el area de 2683,58cm2, puede ser abarcada con 14,7 veces el area de 13.5x13.5 cm^2, por lo que queda bien poner 16 petalos.
  • Para los discos 5 y 6, se puede seguir usando 16 petalos para ser homogéneos con el disco 4, o bien reducirlo a 14 petalos.
  • Para el disco 7, de nuevo, podemos usar 16 petalos para mantener la homogeneidad, aunque bastarían 12 petalos para poder abarcar todo el area del disco con sensores.

En resumen, y dejándolo por simetría en múltiplos de 4, yo pondría:

  • Discos 1 y 2: 8 petalos
  • Disco 3: 20 petalos
  • Discos 4, 5, 6 y 7: 16 petalos.

Muestras utilizadas
  • Samples of single particles of muons at 1, 10 and 100 GeV?, reconstructed by Marcel with release v01-06. Las muestras estan en: /mnt/data2/ilcsoft/ilcinstall_prueba/ilcinstall/Samples/reconstructed_tracks/muones/Marcel_reco_v01-06.
Line: 58 to 41

encuentras una particula MC o una traza reconstruida. El metodo fillTransientHitsCollectionFTD simplemente vuelca el contenedor de hits en una n-tupla. Si te limitas a analizar la ntupla hits el codigo - tal y como esta ahora - te vale.
Changed:
<
<
Numero de hits en porcentaje
>
>
DISCOS COMPLETOS: Numero de hits en porcentaje


Changed:
<
<
Due to many problems with code, for now only the total number of hits is obtained. With this data only the percentage of hits in each section can be seen:
percentage_hits.jpg
>
>
The percentage of hits in each disk can be seen:

procentaje_hits.jpg


  • the most of the hits in 1 and 3 disks, more than 55% in the pixels sensors
  • the number of hist decrease for the farmost 4 rings

Changed:
<
<
Numero de Hits por evento
>
>
DISCOS COMPLETOS: Numero de Hits por evento


Changed:
<
<
Para cada muestra de muones y de ttbar, hago el plots del numero de eventos en cada disco (layer=0, 1, 2, 3, 4, 5 y 6).
  • nhits_evt_muones_100GeV.jpg:
    nhits_evt_muones_100GeV.jpg
>
>
Para cada muestra de muones y de ttbar, hago el plots del numero de eventos en cada disco (layer=0, 1, 2, 3, 4, 5 y 6) y en un lado del detector (side=0).
  • nhits_evt_muones_1GeV_all_disk.jpg:
    nhits_evt_muones_1GeV_all_disk.jpg

Changed:
<
<
  • nhits_evt_tt_bar_bbcssc.jpg:
    nhits_evt_tt_bar_bbcssc.jpg
>
>
  • nhits_evt_ttbar_all_disk.jpg:
    nhits_evt_ttbar_all_disk.jpg

Changed:
<
<
Para estimar el numero de hits por evento, como cojo el valor medio, me basta con dividir el numero total de entradas (hits) en cada disco por el numero total de eventos en ese disco.
>
>
Para estimar el numero de hits por evento, como cojo el valor medio, me basta con dividir el numero total de entradas (hits) en cada disco por el numero total de eventos en ese disco. También es interesante ver el max nº de hits que hay por event en cada disco, lo que hemos llamado "peak".

Deleted:
<
<
numero_hits_per_event.jpg

Deleted:
<
<
Los resultados para las muestras con ttbar, como esperábamos, al ser unas muestras con gran cantidad de trazas, el número de hits es mucho mas que para muones. De este modo, mientras en para muones en el primer disco tenemos de media solo 1/3 hits por cada evento, en los ttbar tenemos hasta 17 hits por evento.

Changed:
<
<
Densidad de Hits
>
>
Los resultados para las muestras con ttbar, como esperábamos, al ser unas muestras con gran cantidad de trazas, el número de hits es mucho mas que para muones. De este modo, mientras en para muones en el primer disco tenemos de media solo 0.16 hits por cada evento, en los ttbar tenemos hasta 8.6 hits por evento.

Respecto al maximo numero de hits por disco, para los 3 primeros discos en las muestras de muones tenemos entre 3-5 hits simultaneos por disco, mientras que para las muestras ttbar esta cantidad sube hasta 37 hits.

DISCOS COMPLETOS: Densidad de Hits



Los resultados de la densidad de hits, calculada como nhits/cm2xBX (1 BX=1 event), para las muestras de muones a 1, 10 y 100GeV, son similares a los obtenidos por Marcel http://ific.uv.es/~vos/ilc/ilcFastForward/MachineBackground.html, pero ahora aplicando la nueva geometría del detector (ILD_00), y el nuevo campo magnetico (3.5T), además de usar una versión mas actual del software.

Changed:
<
<
Density_hits_color.jpg
>
>
hit_density_all_disk.jpg

Por otro lado, los valores obtenidos para la densidad de hits para las muestras de ttbar, son 2 órdenes de magnitud mayores, que para los muones.

Changed:
<
<
Influencia del fondo de pares
>
>
DISCOS DIVIDIDOS EN 8 SECTIONES O PETALOS
  • Para poder acceder al FTD se hace necesaria una estructura en petalos.
  • Interesa que el numero de petalos sea par, y asi poder separar la parte superior de la inferior.
  • Para poder comparar unos discos con otros en las mismas condiciones, para empezar, divido cada disco en 8 petalos (ya veremos en el apartado siguiente, que puesto que los sensores estan determinados por el tamaño de la oblea, cada disco en funcion de su area, debera tener un nº diferente de petalos.)

DISCOS EN 8 PETALOS: Numero de Hits por evento


Para cada muestra de muones y de ttbar, hago el plots del nº de hits en cada uno de los 8 petalos (hit.Phi()>0&&hit.Phi()<1.57/2) en los que divido cada disco (layer=0, 1, 2, 3, 4, 5 y 6), a un lado del detector (side=0).
  • plots_nhits_layers_muons_1GeV_8_petalos.jpg:
plots_nhits_layers_muons_1GeV_8_petalos.jpg

  • plots_nhits_layers_ttbar_8_petalos.jpg:
plots_nhits_layers_ttbar_8_petalos.jpg

Para muones se obtienen muchos 0 y algun 1, mientras que para ttbar, hay mayor cantidad de hits depositados, e incluso de varios hits por evento.

En este caso obtengo el valor medio del nº de hits por evento en cada petalo de disco directamente del mean value del histograma.(*) He comprobado que el valor medio del nº de hits coindice con el que obtengo calculandolo a partir de nº total de entradas (hits) en cada petalo de cada disco dividio por el nº total de eventos en ese disco.

n_hits_per_event_8_petalos.jpg

Tambien es interesante ver el max nº de hits que hay por event en un petalo de cada disco. Dicho valor, se ve mejor a partir de los picos de los plots del nº de eventos en 1 petalo de los 8 que hay en cada disco (layer=0, 1, 2, 3, 4, 5 y 6), en un lado del detector (side=0).

  • nhits_evt_muones_1GeV_8_petalos.jpg:
    nhits_evt_muones_1GeV_8_petalos.jpg

  • nhits_evt_ttbar_8_petalos.jpg:
    nhits_evt_ttbar_8_petalos.jpg

DISCOS EN 8 PETALOS: Densidad de Hits


Divido el nº de hits por evento en cada petalo del disco por el area total del petalo , y asi obtengo la densidad de hits, calculada como nhits/cm2xBX (1 BX=1 event).

hit_density_8_petalos.jpg

Segmentacion mas adecuada para los discos: DISCOS POR ZONAS


Tratamos de elegir la segmentacion mas adecuada para los discos del FTD, teniendo en cuenta que los sensores tienen que caber en una oblea de 13.5x13.5 cm^2.

  • El disco 1, tiene Rint=39mm, y Rext=164mm, lo que le confiere un area de 797.17 cm2. Un cuarto de disco (1/4), lo que tu llamas un "cuadrante" tendría un area de 199.29 cm2 y con una oblea de 13.5x13.5 cm^2 podemos abarcar como mucho: 182.25cm2, por lo que en el disco 1, nos harian falta no 4, sino 5 petalos como minimo, y como queremos que sea un numero par, mejor pondremos 6, o incluso subir hasta 8 petalos.
  • Para el disco 2, tenemos una situación similar al disco anterior, es decir, conviene poner entre 6 y 8 petalos.
  • Para el disco 3, la situación se complica, pues al tener un Rext mucho mayor (308mm) nos da un area mucho mas grande: 644.906cm2, que para abarcarla se llegan a necesitar hasta 18.2 porciones de 13.5x13.5 cm^2, por lo que se requerirían minimo 20 petalos.
  • Para el disco 4, el area de 2683,58cm2, puede ser abarcada con 14,7 veces el area de 13.5x13.5 cm^2, por lo que queda bien poner 16 petalos.
  • Para los discos 5 y 6, se puede seguir usando 16 petalos para ser homogéneos con el disco 4, o bien reducirlo a 14 petalos.
  • Para el disco 7, de nuevo, podemos usar 16 petalos para mantener la homogeneidad, aunque bastarían 12 petalos para poder abarcar todo el area del disco con sensores.

En resumen, y dejándolo por simetría en múltiplos de 4, yo pondría:

  • Discos 1 y 2: 8 petalos
  • Disco 3: 20 petalos
  • Discos 4, 5, 6 y 7: 16 petalos.

DISCOS POR ZONAS: Numero de Hits por evento



Changed:
<
<
-- Main.iglesias - 10 Feb 2009
>
>
Para cada muestra de muones y de ttbar, hago el plots del nº de hits en cada uno de los 8, 20 o 16 petalos los que divido los discos (layer=0, 1, 2, 3, 4, 5 y 6), a un lado del detector (side=0) y obtengo el valor medio del nº de hits por evento en cada petalo de disco directamente del mean value del histograma.(*) He comprobado que el valor medio del nº de hits coindice con el que obtengo calculandolo a partir de nº total de entradas (hits) en cada petalo de cada disco dividio por el nº total de eventos en ese disco.

Added:
>
>
n_hits_per_event_8_20_16_petalos.jpg

Tambien es interesante ver el max nº de hits que hay por event en un petalo de cada disco.

DISCOS POR ZONAS: Densidad de Hits


hit_density_8_20_16petalos.jpg

Influencia del fondo de pares



Line: 107 to 158

META FILEATTACHMENT nhits_evt_tt_bar_bbcssc.jpg attr="" comment="" date="1234953680" path="C:\Users\iglesias\Documents\My_Work_USC\My_Work\ILC_my_work\Density_Hits\Plots\tt_bar_bbcsdu\nhits_evt_tt_bar_bbcssc.jpg" size="280667" user="iglesias" version="1.1"
META FILEATTACHMENT nhits_evt_muones_100GeV.jpg attr="" comment="" date="1234953790" path="C:\Users\iglesias\Documents\My_Work_USC\My_Work\ILC_my_work\Density_Hits\Plots\muones\nhits_evt_muones_100GeV.jpg" size="206213" user="iglesias" version="1.1"
META FILEATTACHMENT Density_hits_color.jpg attr="" comment="" date="1234954144" path="C:\Users\iglesias\Documents\My_Work_USC\My_Work\ILC_my_work\Density_Hits\Density_hits_color.jpg" size="44559" user="iglesias" version="1.1"
Added:
>
>
META FILEATTACHMENT procentaje_hits.jpg attr="" comment="" date="1236077371" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\procentaje_hits.jpg" size="57784" user="iglesias" version="1.2"
META FILEATTACHMENT nhits_evt_muones_1GeV_all_disk.jpg attr="" comment="" date="1236074528" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\nhits_evt_muones_1GeV_all_disk.jpg" size="204649" user="iglesias" version="1.1"
META FILEATTACHMENT nhits_evt_ttbar_all_disk.jpg attr="" comment="" date="1236074548" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\nhits_evt_ttbar_all_disk.jpg" size="275782" user="iglesias" version="1.1"
META FILEATTACHMENT hit_density_all_disk.jpg attr="" comment="" date="1236075376" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\hit_density_all_disk.jpg" size="63243" user="iglesias" version="1.1"
META FILEATTACHMENT plots_nhits_layers_muons_1GeV_8_petalos.jpg attr="" comment="" date="1236076109" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\plots_nhits_layers_muons_1GeV_8_petalos.jpg" size="113990" user="iglesias" version="1.1"
META FILEATTACHMENT plots_nhits_layers_ttbar_8_petalos.jpg attr="" comment="" date="1236076126" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\plots_nhits_layers_ttbar_8_petalos.jpg" size="116133" user="iglesias" version="1.1"
META FILEATTACHMENT n_hits_per_event_8_petalos.jpg attr="" comment="" date="1236076482" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\n_hits_per_event_8_petalos.jpg" size="57276" user="iglesias" version="1.1"
META FILEATTACHMENT nhits_evt_muones_1GeV_8_petalos.jpg attr="" comment="" date="1236076850" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\nhits_evt_muones_1GeV_8_petalos.jpg" size="220859" user="iglesias" version="1.1"
META FILEATTACHMENT nhits_evt_ttbar_8_petalos.jpg attr="" comment="" date="1236076880" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\nhits_evt_ttbar_8_petalos.jpg" size="230090" user="iglesias" version="1.1"
META FILEATTACHMENT n_hits_per_event_all_disk.jpg attr="" comment="" date="1236078486" moveby="iglesias" movedto="ILC.AnalisisDeLaDensidadDeHitsEnLosFTD" movedwhen="1236078892" movefrom="ILC.AnalisisDeLaDensidadDeHitsEnLosFTD" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\n_hits_per_event_all_disk.jpg" size="80518" user="iglesias" version="1.2"
META FILEATTACHMENT hit_density_8_petalos.jpg attr="" comment="" date="1236078980" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\hit_density_8_petalos.jpg" size="86832" user="iglesias" version="1.2"
META FILEATTACHMENT hit_density_8_20_16petalos.jpg attr="" comment="" date="1236080080" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\hit_density_8_20_16petalos.jpg" size="90138" user="iglesias" version="1.2"
META FILEATTACHMENT n_hits_per_event_8_20_16_petalos.jpg attr="" comment="" date="1236080067" path="C:\Users\iglesias\Documents\My_Work_USC\Personal\Twiki_web_ILC_USC\n_hits_per_event_8_20_16_petalos.jpg" size="63250" user="iglesias" version="1.2"
View topic | Diffs | r1.6 | > | r1.5 | > | r1.4 | More
Revision r1.4 - 26 Feb 2009 - 09:47 - Main.iglesias
Revision r1.5 - 03 Mar 2009 - 11:06 - Main.iglesias