1 / 15

SCT Simulation for CTB2004 - Status and HowTos

SCT Simulation for CTB2004 - Status and HowTos. Zdenka Broklova, Carlos Escobar, Peter Kodys. Our status , results and HowTos are summarized on http://ific.uv.es/~cescobar/simulation.html.

dewei
Télécharger la présentation

SCT Simulation for CTB2004 - Status and HowTos

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. SCT Simulation for CTB2004- Status and HowTos Zdenka Broklova, Carlos Escobar, Peter Kodys Our status, results and HowTos are summarized onhttp://ific.uv.es/~cescobar/simulation.html Special thanks to Thijs Cornelissen, Grant Gorfine, Pavel Nevski, Pavel Reznicek, Sasha Rozanov, Manuel Gallas, Vakho Tsulaia and many others for their help. Zdeňka Broklová

  2. History • December 03: - familiarize with Athena, simulations and NOVA database • January 04: - creation of endcap modules and testbeam geometric setup • February 04: - simulation and digitization • March 04: - reconstruction geometric setup G4  hits  digits  spacing  tracking Zdeňka Broklová

  3. SCT_GeoModel • Contents: • geometric parameters of all SCT components • set materials • properties of sensitive wafers • whole SCT description • prepare new geometric description of endcap modules • more realistic hybrid shape • two sensors on each side • adding subspines • improved precision of description • set right materials now it is integrated in SCT_GeoModel (made by Grant Gorfine) Zdeňka Broklová

  4. use components from SCT_GeoModel description of module and additional blocks position for CTB(all number come from NOVA database) geometric setup of various testbeams (May2003, Sep2003, CTB2004) SCT_TestBeamDetDescr Zdeňka Broklová

  5. Installation # setup cmt:# from csh, tcshsource /afs/cern.ch/sw/contrib/CMT/v1r14p20031120/mgr/setup.csh# or from zshsource /afs/cern.ch/sw/contrib/CMT/v1r14p20031120/mgr/setup.sh## use requirements file (if you do not have your favorite):cp -p /afs/cern.ch/user/k/kodys/public/NOVActb_GeoModel/requirements7.5.0 requirements# accomodate requirements to right release and good your paths# now do the cmt staffcmt configsource setup.csh# orsource setup.shcmt co InnerDetector/InDetTestBeam/InDetTestBeamDetDescr/SCT_TestBeamDetDescrcmt co TestReleasecd TestRelease/TestRelease-00-00-16/cmt# edit requirements # Now compilation cmt broadcast cmt configsource setup.csh# orsource setup.shcmt broadcast gmake Start in a new directory Only if you want to make some changes in this package Add all other packages which you want to work with In a new session: cd working_dir source setup.sh cd TestRelease/*/cmt source setup.sh cd ../run/# copy txt, mac and xml files (will be in CVS soon)cp /afs/cern.ch/user/c/cescobar/public/athena/last_release/share/* .cp /afs/cern.ch/user/c/cescobar/public/athena/last_release/data/* .cp /afs/cern.ch/user/c/cescobar/public/athena/last_release/G4macros/* . Zdeňka Broklová

  6. First Example // Load NOVA objects into the Detector store #include "$NOVACNVSVCROOT/share/NovaCnvSvc_jobOptions.txt" #include "$NOVACONVERTERSROOT/share/NovaConverters_jobOptions.txt" //load relevant libraries ApplicationMgr.DLLs += { "GeoModelSvc", "SCT_TestBeamDetDescr", "GeoModelGraphics"}; // Setup Geometry Model Service ApplicationMgr.ExtSvc += { "GeoModelSvc"}; // Setup Geometry Model Tools to be instantiated GeoModelSvc.Detectors = { "SCT_TBDetectorTool" }; // Choose a setup (1 = May 2003, 2 = Sep 2003, 3 = 2004) GeoModelSvc.SCT_TBDetectorTool.setup=3; // This is to include windows and additional blocks or not GeoModelSvc.SCT_TBDetectorTool.Block=true; // This variable must be set correctly according to // value of GeoModelSvc.SCT_TBDetectorTool.setup // for setup=1 and setup=2 //DetDescrCnvSvc.IdDictName="SCT_TB_IDS.xml"; // for setup=3 DetDescrCnvSvc.IdDictName="SCT_TB2004_IDS.xml"; // Setup example Detector Description Algorithm ApplicationMgr.TopAlg = { "GeoModelDisplay" }; TestBeamGraphics_jobOptions.txt Choose Geometry Tool (if you want create whole ATLAS SCT load SCT_GeoModel dll and put SCT_DetectorTool) Variables Setting Algorithm for displaying Zdeňka Broklová

  7. geometry in the GeoModel has to be translated it to a Geant4 geometry Geo2G4 conversion • CTB envelopes - coming from G4 macro ctb_envelopes.mac orctb_envelope_SCT.mac #include "G4Svc/jobOptions.G4Svc.txt" ApplicationMgr.DLLs += { "Geo2G4"}; ApplicationMgr.ExtSvc += { "Geo2G4Svc" }; SCT_GEO2G4_jobOptions.txt Graphic output to a vrml file Zdeňka Broklová

  8. From simulation … • ParticleGenerator • G4 simulation • the standard package G4SimAlg • connection between athena and Geant4 • G4 macros ApplicationMgr.DLLs +={"ParticleGenerator","TruthExamples"};ApplicationMgr.TopAlg = { "ParticleGenerator" , "DumpMC" , "G4SimAlg" };… and many parameters (see documentation) • Magnetic Field |__ ctb_Atlas.mac |__ ctb_envelopes.mac |__ ctb_SCT.mac | |__ ctb_SCTgeometry.mac | |__ ctb_SCTsensitivity.mac |__ ctb_visualization.mac | |__ ctb_SCTVis.mac |__ ctb_MagneticField.mac | |__ (ctb_magfield.xml) | |__ (ctb_magfield.dtd) |__ ctb_eventgraph.mac |__ ctb_physicslist.mac SCT_CTB_simulation_jobOptions.txt Zdeňka Broklová

  9. ****************************************************************************************************************************************************************************************************************** * G4Track Information: Particle = pi+, Track ID = 1, Parent ID = 0 ********************************************************************************************************* Step# X(mm) Y(mm) Z(mm) KinE(MeV) dE(MeV) StepLeng TrackLeng NextVolume ProcName 0 0 -9.98 24.9 1.8e+05 0 0 0 MBPS1::MBPS1 initStep 1 110 -9.98 24.9 1.8e+05 0.0225 110 110 MBPS1::MBPS1 hIoni :----- List of 2ndaries - #SpawnInStep= 1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 1 --------------- : 110 -9.98 24.9 0.00103 e- :----------------------------------------------------------------- EndOf2ndaries Info --------------- 2 264 -9.98 24.9 1.8e+05 0.0308 154 264 SCT::SCT Transportation 3 265 -9.98 24.9 1.8e+05 5.88e-05 0.5 265 Additional block Transportation 4 275 -9.98 24.9 1.8e+05 0.147 10 275 SCT::SCT Transportation 5 278 -9.98 24.9 1.8e+05 0.00108 3.68 278 Additional block Transportation 6 286 -9.98 24.9 1.8e+05 0.303 8 286 SCT::SCT Transportation 7 293 -9.98 24.9 1.8e+05 0.00136 6.3 293 SCT::SCT hIoni :----- List of 2ndaries - #SpawnInStep= 1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 2 --------------- : 293 -9.98 24.9 0.00378 e- :----------------------------------------------------------------- EndOf2ndaries Info --------------- 8 296 -9.98 24.9 1.8e+05 0.000406 3.07 296 SCT_FwdModule_00 Transportation 9 297 -9.98 24.9 1.8e+05 0.000186 1.25 297 SCT::SCT Transportation 10 301 -9.98 24.9 1.8e+05 0.000752 4.27 301 Sensor_Side#0 Transportation 11 302 -9.98 24.9 1.8e+05 0.0929 0.285 302 SCT_FwdModule_01 Transportation 12 302 -9.98 24.9 1.8e+05 6.85e-05 0.68 302 Sensor_Side#1 Transportation 13 303 -9.98 24.9 1.8e+05 0.0806 0.285 303 SCT::SCT Transportation 14 328 -9.98 24.9 1.8e+05 0.00411 25.8 328 SCT::SCT hIoni :----- List of 2ndaries - #SpawnInStep= 1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 3 --------------- : 328 -9.98 24.9 0.00104 e- :----------------------------------------------------------------- EndOf2ndaries Info --------------- 15 350 -9.98 24.9 1.8e+05 0.0052 21.6 350 Additional block Transportation 16 358 -9.98 24.9 1.8e+05 0.239 8 358 SCT::SCT Transportation 17 367 -9.98 24.9 1.8e+05 0.00288 9.38 367 SCT_FwdModule_10 Transportation 18 369 -9.98 24.9 1.8e+05 0.00131 1.25 369 SCT::SCT Transportation 19 373 -9.98 24.9 1.8e+05 0.00103 4.27 373 Sensor_Side#0 Transportation 20 373 -9.98 24.9 1.8e+05 0.0858 0.285 373 SCT_FwdModule_11 Transportation 21 374 -9.98 24.9 1.8e+05 5.95e-05 0.68 374 Sensor_Side#1 Transportation 22 374 -9.98 24.9 1.8e+05 0.0842 0.285 374 SCT::SCT Transportation 23 391 -9.98 24.9 1.8e+05 0.0046 17.3 391 SCT::SCT hIoni :----- List of 2ndaries - #SpawnInStep= 1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 4 --------------- : 391 -9.98 24.9 0.0017 e- :----------------------------------------------------------------- EndOf2ndaries Info --------------- 24 407 -9.98 24.9 1.8e+05 0.00397 15.4 407 SCT::SCT hIoni :----- List of 2ndaries - #SpawnInStep= 1(Rest= 0,Along= 0,Post= 1), #SpawnTotal= 5 --------------- : 407 -9.98 24.9 0.00141 e- :----------------------------------------------------------------- EndOf2ndaries Info --------------- 25 421 -9.98 24.9 1.8e+05 0.00496 14.7 421 Additional block Transportation 26 429 -9.98 24.9 1.8e+05 0.249 8 429 SCT::SCT Transportation 27 439 -9.98 24.9 1.8e+05 0.00158 9.88 439 SCT_FwdModule_20 Transportation 28 441 -9.98 24.9 1.8e+05 0.00116 1.25 441 SCT::SCT Transportation 29 445 -9.98 24.9 1.8e+05 0.0012 4.27 445 Sensor_Side#0 Transportation magnet input window support frame module support frame module Zdeňka Broklová

  10. … to reconstruction • digitization & clusterization • space points • tracking // Digitization and clusterization stuff // ---------------------------------------------- #include "SCT_Digitization/SCT_Digitization_jobOptions.txt" #include "SiClusterizationTool/SCT_ClustersOptions.txt" // Space points // ----------------------------------------------- #include "SiSpacePointFormation/SiTrackerTBSpacePointFinderOptions.txt" NewSiTrackerSpacePointFinder.ProcessPixels = false; NewSiTrackerSpacePointFinder.ProcessSCTs = true; NewSiTrackerSpacePointFinder.ProcessOverlaps = false; // If false, only process opposite wafer NewSiTrackerSpacePointFinder.AllClusters = true; // Tracking // ----------------------------------------------- ApplicationMgr.DLLs+= {"SCT_TestBeamAlgs"}; CBNT_Athena.NtupleLocID="/FILE1/MUTB/tree"; CBNT_Athena.Members += {"TestBeam::SCT_TrackFitter/SCT_TrackFitter"}; SCT_TrackFitter.Chi2cut=true; ApplicationMgr.TopAlg += { "CBNT_Athena" }; Zdeňka Broklová

  11. Some pictures Zdeňka Broklová

  12. Simple diagram NOVADatabase GeoModel Geo2G4 Reconstruction SpacepointMaker Geant4 Digitisation Clustering Tracking TransientEventStore Hits Digits Clusters Space points Tracks Read Algorithm ntuples, histograms, …. Zdeňka Broklová

  13. Outlook • verification • simulation of older testbeams and result comparisons • CTB 2004 tuning in future: • follow reality  reaction to real CTB conditions - comparison with real data Zdeňka Broklová

  14. Conclusion • whole sequenceseems to be working (in athena 7.5.0,7.7.0 and 8.0.0) • standard packages are being used • we know how to run it and make changes(e.g. how to put results into text file, histograms or ntuples) • up to now it hasn’t been much stable (due to massive development) hopefully it will be better Our status and results are summarized onhttp://ific.uv.es/~cescobar/simulation.html Zdeňka Broklová

  15. Conclusion II. • We are about to start completing physical tasks and answering physical questions like standard testbeam analysis (noise occupancy, efficiency, fake tracks, etc) • Practicle guides and examples for newcomers are on our web • We are able and willing to share our experience with others (feel free to contact us) Our status and results are summarized onhttp://ific.uv.es/~cescobar/simulation.html Zdeňka Broklová

More Related