Skip to topic | Skip to bottom
Home
ILC
ILC.AnalisisDeLaDensidadDeHitsEnLosFTDr1.4 - 26 Feb 2009 - 09:47 - Main.iglesiastopic end

Start of topic | Skip to actions

Analisis de la densidad de hits en los FTD

Division de los discos por sectores (petalos)


The idea is to analyse the hits density in FTD disks dividing them in sections (petals). In principle, Disks 1-3 in 4 sections in phi angle (phi section of 90 degrees or atan(y/x) below 1.57 rad) and Disks 4-7 in 16 sections (phi sections of 22.5 degrees or atan(y/x) below 0.3925, 0.785, 1.1775 and 1.57 rad), similar to IFCA design, but…
Area_sections.jpg

….It’s seems better divide the 3 Disk in 16 sections instead of only 4. NOTA: Si, como geometria el tercer disco es mas parecido a los cuatro ultimos que a los dos primeros. En una nueva version del detector habra que añadir un disco pequeñito y seguir con 3+5.

Plots show the x-y positions of hits in Forward Disk, for the whole FTD and for Disk 1 and 4, and their sections:

  • En el primer plot, se muestran los hits en el cojunto de discos del FTD, pero ya se aprecia que la mayoría de los hits caen en la zona central correspondiente a los 3 primeros discos.
  • El 2ºplot, 3º y 4º corresponden al primer disco, primero entero y luego divido en las 4 secciones (2 a 2)
  • El 5º plot y 6º y el 7º y 8º plots corresponden al segundo y tercer disco respectivamente, con sus secciones (2 a 2)
  • El 9º al 12º plots corresponden al cuarto disco y sus 16 secciones (2 a 2) , donde ya cambia algo el tamaño del disco y también disminuye el numero de hits que caen.
  • Del 13º al 16º plots , corresponden a los discos 5, 6 y 7. El numero de hits ha disminuido bastante.

more_entries.jpg

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.

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.

  • Samples with a great ammount of track: ttbar: from http://ilcsoft.desy.de/loi/6f_bbqqqq_500_ILD_00_REC.list. Hay vienen muchas y en principio no necesitaría todas, sino una de cada tipo de canal por ejemplo: Bbcsdu, bbcssc, bbuddu, bbudsc, ssssbb, uuccbb, ddssbb, ddddbb,… I have put them in fpmac116 machine in the directory: /mnt/data2/ilcsoft/ilcinstall_prueba/ilcinstall/Samples/reconstructed_tracks/ttbar

NOTA: El codigo de MyTracking? esta escrito, hoy por hoy, solo para analizar muestas con "single particle", la ntuple "track" no nos valdria para este analisis. Afortunadamente, la ntuple "hit" si sirve para cualquier tipo de muestra, luego si el analisis se limita a usar esta ntuple podria hacerlo con ttbar y fondo de pares.

Explicacion: getTracksFromMC() es valido solo donde haya una unica particula de interes en el suceso. Para muestras tt-bar NO sirve. El codigo de MyTracking? NO te servira para determinar eficiencias en muestras que no sean "single particle". En un futuro tendremos que crear una coleccion con nuestras trazas en el fichero de salida LCIO. Afortunadamente, la ntupla de hits se llenara independientemente de si 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.

Numero de hits en porcentaje


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 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

Numero de Hits por evento


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

  • nhits_evt_tt_bar_bbcssc.jpg:
    nhits_evt_tt_bar_bbcssc.jpg

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.

numero_hits_per_event.jpg

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.

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.

Density_hits_color.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.

Influencia del fondo de pares


-- Main.iglesias - 10 Feb 2009


to top

I Attachment sort Action Size Date Who Comment
Area_sections.jpg manage 37.9 K 10 Feb 2009 - 10:16 Main.iglesias  
more_entries.jpg manage 225.8 K 10 Feb 2009 - 10:34 Main.iglesias  
percentage_hits.jpg manage 34.8 K 10 Feb 2009 - 10:44 Main.iglesias  
numero_hits_per_event.jpg manage 65.7 K 18 Feb 2009 - 10:11 Main.iglesias  
nhits_evt_muones_10GeV.jpg manage 175.4 K 18 Feb 2009 - 10:16 Main.iglesias  
Density_hits.jpg manage 65.8 K 18 Feb 2009 - 10:29 Main.iglesias  
nhits_evt_tt_bar_bbcssc.jpg manage 274.1 K 18 Feb 2009 - 10:41 Main.iglesias  
nhits_evt_muones_100GeV.jpg manage 201.4 K 18 Feb 2009 - 10:43 Main.iglesias  
Density_hits_color.jpg manage 43.5 K 18 Feb 2009 - 10:49 Main.iglesias  

You are here: ILC > AnalisisDeLaDensidadDeHitsEnLosFTD

to top

Copyright © 1999-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding this material Send feedback