antares software & km3net

Bologna 26-27 Aprile 2012

ANTARES software &

KM3NeT

R. Coniglione

Laboratori Nazionali del Sud -INFN


Ease of simulation of different detector lay-outs

Ease of implementation of PMTs and water parameters

Ease of adding modules and change part of the code

Output/input format friendly

Simulation in KM3NeT requirements


The KM3NeT optical sensors

la

ANTARES like DOM

1 large PMT (10” or 8“) inside a benthos sphere (17”)

31 little PMT (3”) inside a large benthos sphere (17”)


The KM3NeT Detection UNits

la

String like Tower -like


Large number of DUs (up to about 300)

Tipically 20 floors/DU with different spacing

Very large number of PMTs (up to about 200000) and hits (500000)

Simulation for nm from 100 GeV to about 108 GeV

Some numbers


The actual simulation scheme:

MupageAtmosheric m

generator

GenhenNeutrino generator

KM3Light & hits

GEASIMLight & hits

RECOReconstruction

code

MODK4040K Background hits

I/O format -> ASCII files

GENDETgeometry

Running on Lyon clusters

+

Only for low energy events <100GeV


Software organization

Different packages and different environments for different PMTs:

Antares_km3: for large PMTs

Antares_seawiet: for DOM


The versions

Antares_KM3 Antares_sewiet

GENDET v1r3 v1r3

MUPAGE v3r4 V3r4seawiet

GENHEN v6r3km3 v6r7seawiet

KM3 v3r6km3 up to v4r3seawiet Many versions

MODK40 v4r10km3 v4r10seawiet In geasim

ANTCC v3r7 v3r7seawiet

GEASIM v4r10km3 V4r10seawiet

IO v1r12 V1r12seawiet


IO librarymain modifications

Only DIMENSION in f_77.inc were changed

◦ PARAMETER (MAXTAG=150) it was 100◦ PARAMETER (MAXDAT=400000) it was 200000


GEASIM librarymain modifications (1)

• In the ANTARES like geometries a cluster consists of PMTs that can be repeated in space. In a tower based geometry with DOM a cluster consists of two perpendicular bars with 4 DOM. In total 124 PMTs per cluster. In a tower with 20 floors will be 10 cluster. 1540 cluster in a detector with 154 Towers.

hits.inc PARAMETER (MAXHIT=1000000) !it was 19500 - maximum number of hits

phomul.inc PARAMETER (NPMMAX=300000) maximal number of PMTs allowed PARAMETER (MAXCLU=125) maximal number of PMT cluster PARAMETER(MAXNCL=7000) maximal number of cluster


GEASIM librarymain modifications (2)

• gugepm.f is the subroutine that reads the geometry file. In a geometry based of towers the lines in the .det le with tag OM_cluster are too long. The section in this subroutine that reads the lines and calculates the number of PMT in the cluster has been modified

• guout.f • FORMAT has been modified

• lambda.f• PMT angular acceptance (function wang(ct)) and PMT Quantum

Efficiency changed. Not Called by KM3 but by the reconstruction code (wang

function) and by GEASIM code


GENHENmain modifications (1)

Added a new subroutine to take into account generation from extended source: extsource.f

NuPrecut.f: bugs found and a line added◦ HISTCAN was neither calculated nor transferred by a COMMON CALL DISTTOCAN(VX,VY,VZ,-VX,-VY,-VZ,DISTCAN,IFLAGCAN) !

added by R. Coniglione 26-1-2006 RangeEffective_music.f: : bugs found

◦ cc RANGEEFFECTIVE = 0. RANGEEFFECTIVE_music = 0. !!corretto da R. Coniglione 24-4-2007 cc RANGEEFFECTIVE = (A * EMU) + B RANGEEFFECTIVE_music = (A * EMU) + B !!corretto da R. Coniglione

24-4-2007 defcut.f

◦ cc MSTU(22) = 1000 MSTU(22) = 10000! number of errors allowed !! changed by R.

Coniglione 12-10-200


GENHENmain modification (2)

evpointdir.f For point-like generation back to antares version v6r3

◦ In the new v6r7 antares version +x is east and +y is north, soc EVCX = SIN(THETA)*SIN(PHI)c EVCY = SIN(THETA)*COS(PHI)c EVCZ = COS(THETA)cC R. Coniglione 18-7-2011 to go back to the old genhen v6r3c old ANTARES orientationc EVCX = -SIN(THETA)*COS(PHI) EVCY = -SIN(THETA)*SIN(-PHI) EVCZ = -COS(THETA)


KM3-GEN-HITmain modifications (1)

Water properties read from a file: gen-water-optical-partic-nemo-0.0075.dat

km3.inc DIMENSION changed◦ PARAMETER (MAX NUM OMS DET = 200000) ! Max number of OMs in

whole detector◦ PARAMETER (MAX NUM OMS CLS = 125) ! Max number of OMs in any

cluster◦ PARAMETER (MAX NUM STR = 512) ! Max number of strings in whole

detector hit-eff_area_pmt.f

◦ the PMT photocatode area changed PARAMETER ( pm diam = 2.99 ) 3” PMT◦ Changed IF statement on the photon energy IF (energy.gt.4.14334 .or. energy.lt.1.173847) THEN… the limits depend on the wave length reported in the water description

file


KM3-GEN-HITmain modifications (2)

hit-effang.f modified to include the new PMT angular acceptance (function effang(ct)): A function is required.

hit-ini_optic.f modified to take into account the new PMT quantum efficiency◦ DATA Q_EFF0/0.01/◦ DATA Q_EFF/0.8,…..

FORMAT changed in the write statement of the tag “hit:” inkm3-elec direct hits.fkm3-muon direct hits.fkm3-muon scattered hits.fkm3-phot direct hits.f:km3-phot scattered hits.fkm3-process k40.fkm3-write antares event.f110 FORMAT('hit:',1X,I7,1X,I7,1X,F6.1,1X,F10.2,3(1X,I6))


MODK40 in GEASIM• main-modk40.f

• a COMMON was introduced for the FFKEY cern subroutines (from 32bit to 64 bit)

COMMON /FFR/ ISEED,NEVE• addK40coinc.f new subroutine included to take into account the L1

coincidence due to 40K in DOM• The addraw.f (called also by Km3)

• FORMATs in the write statements have been changed


Input/output formatASCII files that can be converted in root

files with

Antviewer (read root files) & Anteater (create ROOT files) largely modified.

Some leaves deleted to reduce file size


Running time & storageDetector ref154_3inch31pm180_204006

5 108 events with 100<E<108 GeV spectral index 2

code Time (minutes)

Size

GENHEN 3 40M

KM3 20 65M

MODK40-L1 >120 2 G

GENHEN 8 21M

KM3 30 350M

MODK40-L1 >120 2.1G

• 2 108 events with 100<E<108 GeV spectral index 1.5


Ease of simulation of different detector lay-outs

Ease of implementation of PMTs and water parameters

Ease of adding modules and change part of the software

Output/input format friendly

Simulation in KM3NeT requirements


PMT angular acceptanceAngular acceptance large PMTs

Angular acceptance 3” PMTs


PMT quantum efficiency


absorption lenght

Toulon*Capo PasseroCapo Passero + 10%Smith&Baker

antares software & km3net

Documents