using jupiter and satellites

UsingJupiter and SatellitesAkiya Miyamoto

KEK

Feb 2007

A.Miyamoto, LCPAC, 2004/02/20

Target of Jupiter/Satellites/UranusJupiter and Satellites are tools for detector optimization based on Geant4 Full detector simulation.Implement an ideal GLD geometry in Jupiter ( other geometry is also possible )Study PFA performances by an ultimate conditionImplement tower calorimeterDevelop analysis tools Study physics performance vs. detector choice.Not for simulation studies for beam test studies

A study by detector simulation is also crucial for the design of the IR design of ILC.


Jupiter/Satellites ConceptsJUPITERJLC UnifiedParticle InteractionandTracking EmulatoRIOInput/Outputmodule setURANUSLEDAMonte-Calro Exact hits ToIntermediate Simulated outputUnified Reconstruction and ANalysis Utility SetLibrary Extention for Data AnalysisMETISSatellitesGeant4 basedSimulatorJSF/ROOT basedFrameworkJSF: the analysis flow controller based on ROOT The release includes event generators, Quick Simulator, and simple event display MC truth generatorEvent ReconstructionTools for simulation ToolsFor real data


JupiterbindocincludeconfiglibmacrosscriptssourcestmpbdcalcdctpcctitvtxirsrcincludekernlcexpsolmainJupiter.ccSub-Detector filesLinux-g++JupiterReadmecommon.gmkSample G4 scriptsDetector components to use are switched by CPP flags in Jupiter.cc or parameters in JSFJ4Jupiter Directory Structure


Standalone Jupiter$ Jupiter [options] [g4macro file] options: -v N : N is a message level. > 10 ; Write message every time new geometry is created. 5-10 ; Write message of geometry whose rank is less than N-4. < 5 ; No message is written at geometry creation. -f FILE : Read parameters from the file FILE More than one -f option can be specified. In this case a value in the first file is used. -w : Write default parameters defined in the program to the file, jupiter.defaults -help : Print help message Useful to display events and debug


G4macro fileCommands to Geant4 and Jupiter are provided as a g4macro file. A list of commands can be found at When Jupiter is executed as a standalone application interactively, help command can be used to get help information. http://www-jlc.kek.jp/~miyamoto/Jupiter/html/_.html


Jupiter parameter list fileParameters for detector configuration can be modified by a parameter file given by f option.

$JUPITERROOT/doc/HowToRun for more details

Standard configuration file as of Jan 2006 is at http://ilcphys.kek.jp/soft/samples_dod/dec05/index.html


Satellites/Uranus


JSFJ4Interface to run Jupiter in JSF environment.Allows to writeJupiter output by ROOT or LCIOUse generators in JSF or StdHep format file.Satellites package reads ROOT file created by JSFJ4


Metis packageMetis is a collection of reconstruction tools for Jupiter data.Current aim is to prepare a minimum set of Metis modules for studies of Particle Flow Algorithm.Novice users will be able to do physics analysis using information of PFO classes.As a first step, a cheated track finder and a cluster maker, etc are in preparation in order to know ultimate performance.

Each module is independent, thus shall be easy to implement different reconstruction algorithm according to interests Development of a real clustering code begins recently


make smeared TPC hits from exact hitmake tracks from TPCmake hybrid tracks ( TPC+IT+VTX)make smeared/merged CAL hits from exact hitmake cluster from CAL hitsmake Particle Flow Objectsjet clusteringMetis Analysis Flow


.cdcvtxSatellitesbinsrclibincludetestiojsfj4mctruthkernexamplesRun Jupiter in JSF to create a ROOTfile calS4xxxExactHit class= J4xxxHit classmetisSatellites Directory StructureLedaj42lcioOutput LCIO dataexamplesmacro


calmetistpchitmakertrackmakerhybthybtmakerhitmaker(JSFs) Modules for MC data analysismake smeared TPC hits from exact hitmake tracks from TPCmake hybrid tracks ( TPC+IT+VTX)make smeared/merged CAL hits from exact hitclustermakermake cluster from CAL hitspfopfomakermake Particle Flow ObjectsMetis Directory Structurejetjetmakermake jet Objects


Uranus Directory StructureUranusdatasrclibincludekernhybtVTXPackages for real data analysisITTPCdetconfighitmakertrackmaker


S4XXXTrackHitU4VTrackS4XXXTrackdouble fChi2U4VTrackMaker* fU4VJSFHelicalTrack* fHTS4XXXHitdouble fErrorS4XXXExactHit *HitRelation among TrackHits


Relation among CAL HitsS4CALHitS4VHitTObjectTAttLockable+GetExactHitsPtr()+GetAddress()+GetHitPosition()+ GetCALType()+GetClusterType()S4CALExactHit+ GetPreHitPtr()S4VExactHitGetAddressPtr()GetSmearedHitPtr()


Typical results


Momentum resolutionExact hit points created by single m were fitted by Kalman filter packagespt/pt2 (GeV -1)By A.Yamaguchi(Tsukuba)


By A.L.C.Sanchez (Niigata U.)Linearity is good, but to get energy resolution similar to beam test results,Randge cut of O(1) mm is required. It is very small and we dont know why we need O(1) mm.EM Cal Performance


Cheated Particle Flow AnalysisKey point of PFA is use Cal. signals only for neutral particles Remove CAL signal if connected to a track

Use Simulation information to connect track and cal. signals Cheated PFAWhat affects jet energy resolutionSignal sampling fluctuation in Cal.Tracker resolutionTreatment of V0, decays,and interactions before Cal. more Understand factors which affect resolutionA sample detector signal


Cheated PFO analysisZH event at Ecm=500 GeVBy K.Fujii(KEK), S.Yamamoto(GUAS), A.Yamaguchi(Tsukuba)- Exact hit points of TPC and CAL are displayed.Hits belong to the same PFO are shown with the same color

A framework of event display in JSF is used.


X3DROOTsX3d view of the same event


X3D-JetSame event, after a forced 4-jet clustering on PFObjects


Jet energy resolutione+e- Z q qbar at sqrt(E)=91.1GeVBy Sumie Yamamoto


DE for Perfect CAL./TrackerJet energy resolution when resolution of CAL and tracker are perfect.GeV/c2Contribution of each detector

Total 2.7 GeVHD Cal 2.0 EM Cal 1.2Tracker 0.4Others 1.4 (under investigation)By Sumie Yamamoto


Satellites examplesExam01Exam02Exam03Exam04Exam05Exam06Exam07Exam08Exam09

Analysis code examples


using jupiter and satellites

Documents