1 / 8

Access to SICB data

Access to SICB data. Modifications in Gaudi v5 release: Pile up. Other new features: Pre-staging of several tapes Access to new data base. Pile Up Implementation. Pile up infra-structure was already in Gaudi v4 But it uses internally SICB and it was not working.

mari
Download Presentation

Access to SICB data

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. Access to SICB data Modifications in Gaudi v5 release: Pile up. Other new features: Pre-staging of several tapes Access to new data base Gonzalo Gracia Abril

  2. Pile Up Implementation • Pile up infra-structure was already in Gaudi v4 • But it uses internally SICB and it was not working. • Some more debugging was required. • Pile up is implemented as a separate algorithm which has to run before any other algorithm which access data: JobOptions: ApplicationMgr.TopAlg = { "PileUpAlg",”myAlg1",”myAlg2" }; • Number of pile up events is calculated as a function of the luminosity. • RELUMI re-implemented in C++: Reads beam.cdf file to get input parameters. Uses new random generator service (See later) • Two files read in, one for the main event and the second for pile up events.

  3. SicbEventSelector “ EventSetector ” ApplicationMgr taSvc “ EventDataSvc ” icnbSicbPileupAlg “ PileupAlg ” SicbPileupSvc EventDataSvc “ PileupSvc ” “ PileupDataSvc ” SicbEventSelector “ PileupSelector ” IProperty IProperty IEventSelector SicbEventSelector IEventSelector ApplicationMgr IDataManagerService IDataManagerSvc ZEBRA EventDataSvc SicbPileUpAlg IDataProviderSvc IDataProviderSvc ReMerge SICB specific IDataProviderSvc SicbPileUpSvc EventDataSvcvc ZEBRA IDataManagerSvc IDataManagerService IEventSelector SicbEventSelector The Algorithm and Service are merged into one Algorithm for the Sicb version Not needed for the SICB version The Algorithm and IEventSelector Not needed for Service are merged the Sicb version into one Algorithm for IProperty IProperty the Sibc version

  4. Random Numbers • Reminder : Two random numbers are generated for every event to get the number of pile up events: • First number generated following a flat distribution (0,1) to calculate the current luminosity (from the fill duration and beam decay time, the current time in the fill is the random number). • Second (number of interactions): Follows a Poisson distribution whose mean is equal to the average number of interactions for the current luminosity (extracted from the cross section and the bunch crossing rate) • If the number of iterations is 0 we repeat the process until is different than 0 • Random number service is now available in Gaudi and we have chosen to use it.

  5. Histogram: Gaudi Points: SICBDST 500 Events Number of interactions.

  6. Histogram: Gaudi Points: SICBDST 500 Events Average number of interactions per bunch crossing.

  7. Pile Up: Some comments • The actual merging of events is done by REMERGE • ATMC, AVMC and raw hits banks are merged. (CALO raw banks are not converted) • When you retrieve, from your Gaudi Algorithm, data converted from one of these banks you will get the merged data. • In the conversion Gaudi follows every reference, we find always the missing or wrong references!!. (It runs in 500 events, several data samples tried… but watch the DST2 problems!!) • Two instances of the event selector which are configured separately

  8. Other new features • Pre-staging of several tapes: • You wait for the first tape. The others are staged while you read the first one. If you finish reading the first file and the second tape is not staged you will get a message and will have to wait. • Access to tapes is not possible from NT • Access to new data base • Gaudi sends a request to the ORACLE web server • Executes the same procedure used by mcstagein to get the volume number, file sequence number and size.

More Related