1 / 18

Status of the Software and MC production

Status of the Software and MC production. Antoine Cazes Laboratoire de L’accélérateur Linéaire. Analysis. The Software chain. Beam Files. Available for the summer. Used for the present MC production. OpRoot. OpSym: OpGeom, OpVMC. OpTreeConvert. Simulation Files. OpDigit. OpNtuple.

havard
Download Presentation

Status of the Software and MC production

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Status of the Software and MC production Antoine Cazes Laboratoirede L’accélérateurLinéaire

  2. Analysis The Software chain Beam Files Available for the summer Used for the present MC production OpRoot OpSym: OpGeom, OpVMC OpTreeConvert Simulation Files OpDigit OpNtuple OpRec MyAna Reconstructed Files Antoine Cazes

  3. Software Status • OpROOT(D.Duchesneau) • Test have been performed • Problem for neutral charmed particle fixed • Few improvements • OpDigit • No Changes • OpRec(AC, JEC) • Pattern : No Changes • Tracking: • Tentative of including the tracking in the spectrometer. • But it is not working • Old version is used : use MC information for initialization of the track ending in a Spectrometer Antoine Cazes

  4. Tracking in the spectrometer: • Tracking in the Spectrometer: • Kalman filtering • parabolic • include matter effect. • include magnetic field • It has been coded. • X-Y reconstruction • Slope reconstruction • Momentum reconstruction OK OK Diverge!!! Antoine Cazes

  5. Temporary Solution • Use of the old Kalman filtering • Tracking in the Target Tracker. • Kalman filtering including matter effect • initialization with fit on the first points • Track going through the spectrometer • initialization using MC information • Tracking done only in the first hit TT • Already more realistic, since pattern do not use any MC information. Antoine Cazes

  6. Result : X-Y resolution (nmCC event) New Tracking Previous Tracking * Xtrue - Xfit Xtrue - Xfit Ytrue - Yfit Ytrue - Yfit * OPERA-NOTE 36 RMS = s = 12mm 8mm RMS = s = 16mm 9mm Antoine Cazes

  7. RMS = s = 37mrad 14mrad RMS = s = 47mrad 17mrad Results : Slope resolution (nmCC event) New Tracking Previous Tracking * Atrue - Afit Atrue - Afit Btrue - Bfit Btrue - Bfit * OPERA-NOTE 36 Antoine Cazes

  8. OpNtuple (A.Marotta) • OpNtuple is a converter from the official output of the software chain to Hbook files • It is available on cvs (v2) • Derived from MyAna • For compatibility reason, the information stored are (almost) the same as in the previous Ntuple production (may 2001) Antoine Cazes

  9. The Hbook file containt • Monte Carlo information * /MCINFO/ * NUMEVTevent id number* /MCINFO/ * ENEUTRIneutrino energy[0,3] * /MCINFO/ * VERTnumber of vertexes* /MCINFO/ * InPart(VERT)incoming particle* /MCINFO/ * XVTX(VERT)x,y,z of the interaction vertexes* /MCINFO/ * YVTX(VERT)* /MCINFO/ * ZVTX(VERT) [1,132] * /MCINFO/ * IPLBEST1-62 wall number if vertex in brick 1-31 sm=1; 32-62 sm=2) 64 interaction before the first sm 66 interaction in the 1st spec zone 68 interaction in the 2nd spec zone 71-132 wall number if vertex in TT 71-101 sm=1; 72-132 sm=2 Antoine Cazes

  10. Muon information [0,3] * /MCMUON/ * NMUONSnumber of (anti)muons generator level* /MCMUON/ * GeantCode(NMUONS)geant code of the (anti)muon* /MCMUON/ * energy(NMUONS)* /MCMUON/ * PXMUON(NMUONS)* /MCMUON/ * PYMUON(NMUONS)* /MCMUON/ * PZMUON(NMUONS)* /MCMUON/ * MUONVTX(NMUONS)link with vtx (1,2,3)* /MCMUON/ * XMUT(NMUONS)last hit coordinate in the TT* /MCMUON/ * YMUT(NMUONS)* /MCMUON/ * ZMUT(NMUONS)* /MCMUON/ * PXMUT(NMUONS) * /MCMUON/ * PYMUT(NMUONS) * /MCMUON/ * PZMUT(NMUONS)* /MCMUON/ * XMUS(NMUONS)last hit coordinate in the spect* /MCMUON/ * YMUS(NMUONS) * /MCMUON/ * ZMUS(NMUONS) * /MCMUON/ * PXMUS(NMUONS) * /MCMUON/ * PYMUS(NMUONS) * /MCMUON/ * PZMUS(NMUONS) [0,1000] * /MCOTH/ * NOTHnumber of particle generator level (no mu)* /MCOTH/ * CODOTH(NOTH)Geant Code of the particle * /MCOTH/ * PXOTH(NOTH) * /MCOTH/ * PYOTH(NOTH) * /MCOTH/ * PZOTH(NOTH) * /MCOTH/ * OTHVTX(NOTH)link with vtx (1,2,3)* /MCOTH/ * ENOTH(NOTH) Antoine Cazes

  11. Last emulsion sheet information [0,150] * /LASTES/ * NLASTES number of hits in the last sheet of the interaction brick * /LASTES/ * IndexLastes(NLASTES) index * /LASTES/ * IDLASTES(NLASTES) Geant Code of the particle * /LASTES/ * IDNUMBLASTES(NLASTES) Brick Id 10000*WALL+100*COLUMN+RPW (WALL:1-62) * /LASTES/ * XLASTES(NLASTES) * /LASTES/ * YLASTES(NLASTES) * /LASTES/ * ZLASTES(NLASTES) * /LASTES/ * PalphaXLASTES(NLASTES) alphaX * /LASTES/ * PalphaYLASTES(NLASTES) alphaY * /LASTES/ * MODPLASTES(NLASTES) momentum [0,1000] * /STRIPS/ * NSTRIPS number of strips * /STRIPS/ * IndexTT(NSTRIPS) index * /STRIPS/ * XSTR(NSTRIPS) * /STRIPS/ * YSTR(NSTRIPS) * /STRIPS/ * ZSTR(NSTRIPS) * /STRIPS/ * NUPLSTR(NSTRIPS) Id SM*100000+WALL*1000+BOX*100+BAR * /STRIPS/ * ESTRL(NSTRIPS) number of photons left side * /STRIPS/ * ESTRR(NSTRIPS) number of photons right side [0,1000] * /STRIPS/ * LKMUSTR(NSTRIPS) link with primary muons 0 if no muon linked +1 if a hit in the dgt by muon from vtx 1 +10 if a hit in the dgt by muon from vtx 2 +100 if a hit in the dgt by muon from vtx 3 Antoine Cazes

  12. Detector information [0,800] * /TUBES/ * NTUBES * /TUBES/ * IndexDT(NTUBES) index * /TUBES/ * IDTUBE(NTUBES) 100*SM+PLANE * /TUBES/ * XTUBE(NTUBES) * /TUBES/ * ZTUBE(NTUBES) [0,1000] * /TUBES/ * LKMUTUBE(NTUBES) link with primary muons [0,1000] * /SRPC/ * NHSRPC number of hit * /SRPC/ * IDHSRPC(NHSRPC) geant code of particle-hit * /SRPC/ * EHSRPC(NHSRPC) * /SRPC/ * PLHSRPC(NHSRPC) 1000*SM+100*halfSRPC+plane * /SRPC/ * XHSRPC(NHSRPC) * /SRPC/ * YHSRPC(NHSRPC) * /SRPC/ * ZHSRPC(NHSRPC) [0,1000] * /SRPC/ * NRSRPC number of reconstructed hit * /SRPC/ * IndexSRPC_H(NRSRPC) index of H dgt of the rec hit * /SRPC/ * IndexSRPC_V(NRSRPC) index of V dgt of the rec hit * /SRPC/ * XSRPC(NRSRPC) * /SRPC/ * YSRPC(NRSRPC) * /SRPC/ * ZSRPC(NRSRPC) * /SRPC/ * IDPLSRPC(NRSRPC) 1000*SM+100*halfSRPC+plane [0,800] * /XPC/ * NHXPC number of hit * /XPC/ * IDHXPC(NHXPC) geant code of particle-hit * /XPC/ * EHXPC(NHXPC) * /XPC/ * PLHXPC(NHXPC) 100*SM+plane * /XPC/ * XHXPC(NHXPC) * /XPC/ * YHXPC(NHXPC) * /XPC/ * ZHXPC(NHXPC) Antoine Cazes

  13. Detector information and tracking [0,800] * /XPC/ * NRXPC number of reconstructed hit * /XPC/ * IndexXPC_U(NRXPC) index of U dgt of the rec hit * /XPC/ * IndexXPC_V(NRXPC) index of V dgt of the rec hit * /XPC/ * XXPC(NRXPC) * /XPC/ * YXPC(NRXPC) * /XPC/ * ZXPC(NRXPC) * /XPC/ * IDPLXPC(NRXPC) 100*SM+plane * /TRACK/ * NTRACK Number of Track * /TRACK/ * NDIM(10)Dimention of the fit of track i (2 or 5) * /TRACK/ * NTT(10)Number of TT hit in the track * /TRACK/ * IDTT(150,10)IDTT(i,j)Index of the TT hit i of track j * /TRACK/ * NSRPC(10)idem for SRPC * /TRACK/ * IDSRPC(150,10) * /TRACK/ * NSDT(10) idem for DT * /TRACK/ * IDSDT(150,10) * /TRACK/ * PARA(5,10) PARA(i,j) result of the fit of track j * /TRACK/ * VARI(15,10) Variance Matrix of the fit of track j PARA(k,i): kth parameter fit at First TT Point of Track i 5Dim Model: (x,y,a=px/pz,b=py/pz,d=q/p) 2Dim Model: (x,a) or (y,b) VARI(.,it): element of the Variance Matrix. In case of 5 Parameter fit the storing in the following Vxx Vxy Vxa Vxb Vxd 1 2 3 4 5 Vyy Vya Vyb Vyd 6 7 8 9 Vaa Vab Vad = 10 11 12 Vbb Vbd 13 14 Vdd 15 Antoine Cazes

  14. Monte Carlo Production

  15. Production Status • December 2003 production was not perfect • We started a new production (may 2004) • Input • Beam Files from may2001 • Vertices only in lead. • Tracking is included. • Already 10,000 or 5,000 events of type: • nmCC • nmNC • Tests are under going to check the quality. • tm ntuple is investigated • Full production as soon as checks are OK • te • te QE • tm • tm QE Antoine Cazes

  16. File location: • Files are stored on the CASTOR area: • /castor/cern.ch/opera/mc/prodmay2004/ • simul/ : ouput of the simulation • digit/ : output of the digitization • reco/ : output of the Reconstruction: • ntuple/ : output of the Hbook converter • File name structure : numucc1_reco_1.root Simulation file number (1-6 for one Beam File) Type of event Beam file number (1-3) Level in the chain (simul, digit, reco) Antoine Cazes

  17. Crash problem • OpRec and OpNtuple experience crashes for big files when running in batch. • Apparently, there is a memory problem, in the OpIO package, maybe memory leak • Test using dedicated debuging software are foreseen. • Purify • Insure • Does not stop the production (just slow down ) • Temporary solution: number of event/file limited to less than 3000. Antoine Cazes

  18. Documentation • An html documentation is being written. • For each package, it will contain information about: • Installation • Contents • Action • Input and Output • Documentation about FEDRA already available. Antoine Cazes

More Related