470 likes | 589 Views
Muon Reconstruction/Trigger Performance with TGC 1 Bunch Readout Mode. Koji Sato MTSG meeting July 21 , 2014. Motivation. TGC ROD busy caused a lot of downtime in 2012. Significant amount (0.35 fb -1 ) of collisions were lost due to this problem. ↓
E N D
MuonReconstruction/Trigger Performance with TGC 1 Bunch Readout Mode Koji Sato MTSG meeting July 21, 2014
Motivation • TGC ROD busy caused a lot of downtime in 2012. • Significant amount (0.35 fb-1) of collisions were lost due to this problem. ↓ • Track down the source of the problem (→electrical disturbance). • Implement a scheme to reject events with burst of noise to the trigger system. • Minimize data size, and allow more events in the readout buffer. • TGC read outs hits from 3 bunches: • Previous, current and next bunch • Does the muon quality degradate, if we only read out the current bunch? (1 bunch readout) • I show a study of offline muon properties this time. • How the change affects the trigger is under study.
RunNumber 215473 EventNumber39074754 CAL Burst Event ② TGC + CAL RPC + TGC + CAL MDT + RPC + TGC + CAL From Masato’s slides: TDAQ Week, 8 July 2013
Analyzed Datasets Today, I will show plots showing all the combined muons (EF/L2 and Offline) found in the D3PDs. Run 215433: Nov 30, 2012 Peak lum=6.7e33 IntLum=123 pb-1
EF muon distributions • Plotting distributions of EF muons with mu6 bit and associated with CB tracks. Blue: 3 bunch readout Red: 1 bunch readout
EF muon distributions • Plotting distributions of EF muons with mu6 bit and associated with CB tracks. Blue: 3 bunch readout Red: 1 bunch readout
Difference in Simulated EF Trigger Muons, Z->mumu MC • Plotting (1bunch – 3bunch mode) EF muons with mu6 bit • Similar results from J/psi sample • No difference in # EF muons • Very slight change in kinematics for few events
L2muon distributions • Plotting distributions of L2 comb muons with mu6 bit. Blue: 3 bunch readout Red: 1 bunch readout
L2 muon distributions Blue: 3 bunch readout Red: 1 bunch readout
Difference in Simulated L2 Trigger Muons, Z->mumu MC • Plotting (1bunch – 3bunch mode) muons with mu6 bit • Similar results from J/psi sample • No difference in # L2 muons • No change in kinematics
L1muon distributions • Plotting L1 muons without requiring trigger bit. • D3PD only contains L1 bit for event. Blue: 3 bunch readout Red: 1 bunch readout
L1 muon distributions Blue: 3 bunch readout Red: 1 bunch readout
Difference in Simulated L1 Trigger Muons, Z->mumu MC • Plotting (1bunch – 3bunch mode) • We expect no difference for L1, since L1 muon trigger only considers Current bunch hits. • Similar results from J/psi sample • No difference in kinematics
Summary of Trigger Muon Properties • I compared simulated trigger muons between 1 and 3 bunch readout scheme for TGC. • Using and MC samples. • Only slight difference in kinematic properties of EF muons were seen for a few events in sample.
Comparison of offline reconstruction Previous presentation at MTSG on Jan 27, 2014, where reconstructed muons were looked at in data: https://indico.cern.ch/event/296921/ Updates are: Muon eta regioning is updated. Also looking at MC samples.
Offline Reconstructed Muons (Full Stat.) Blue: 3 bunch readout Red: 1 bunch readout • CB muons of MUONS collection, Pt>7 GeV
Offline Reconstructed Muons(small stat.) Blue: 3 bunch readout Red: 1 bunch readout • CB muons of MUONS collection in 1.1<|η|<1.6. MC MC
Offline Reconstructed Muons(small stat.) Blue: 3 bunch readout Red: 1 bunch readout • CB muons of MUONS collection in 1.1<|η|<1.6. MC MC
Offline Reconstructed Muons(small stat.) Blue: 3 bunch readout Red: 1 bunch readout • CB muons of MUONS collection in 1.1<|η|<1.6. MC MC MC MC
Offline Reconstructed Muons(small stat.) Blue: 3 bunch readout Red: 1 bunch readout • CB muons of MUONS collection in 1.1<|η|<1.6. MC MC
Dimuon Mass (Full stat.) Blue: 3 bunch readout Red: 1 bunch readout MC MC
Difference in Offl. Reco Variables, Z->mumu MC • Plotting (1bunch – 3bunch mode)
Difference in Offl. Reco Variables, J/psi->mumu MC • Plotting (1bunch – 3bunch mode)
Difference in Offl. Reco Variables, Data • Plotting (1bunch – 3bunch mode)
Muons of 3 BC scheme matched to 1 BC muon within dR<0.01 Data Pt eta phi author nTGCEtaHits nTGCPhiHits
Muons of 3 BC scheme matched to 1 BC muon within 0.01<dR<0.1 Data Pt eta phi author nTGCEtaHits nTGCPhiHits
Muons of 3 BC scheme having no matched 1 BC muonwithin dR <0.1 Data Pt eta phi author nTGCEtaHits nTGCPhiHits
Summary of Offline Muon Properties • Compared properties of offline reconstructed CB MUONS with Pt>7 GeV. • Using and MC samples. • Only found one muon reconstructed with ~0.001 change in eta and phi. • Using Muons.Merge Data of Run 215433. • 99.86% of muons of 3 BC scheme find a matched 1 BC muon within dR<0.01. • 3,139,971 (3,139,888) muons were found with 3 (1) bunch scheme, in 3512939 events. • 3,135,629 muons of 3 BC scheme have a matched muon of 1 BC scheme within dR<0.01. • 1,281 muons have a matched 1BC muon within 0.01<dR<0.1. • 3,061 muons don’t have a matched 1 BC muon within dR<0.1.
Summary • Almost no difference is seen between 1 and 3 bunch readout schemes in terms of CB muon quantities of trigger and offline muons in MC studies. • ~0.1% of 3 BC muons in data don’t get reconstructed in 1 BC scheme. • The results of offline studies of data and MC don’t look very consistent. • Is there anything MC samples doesn’t correctly simulate? Is its simulation of prev. and next bunch collisions tuned well? • Implementation of 1 bunch readout to Athena code is different between data and MC. • Plan to reprocess data but with MC Athena code, and see if the results change.
Motivation • TGC ROD busy caused a lot of downtime in 2012. • Significant amount (0.35 fb-1) of collisions were lost due to this problem. • Summary of 2012 Downtime: • Total: 1.6 fb-1 • Start/end of run: 0.6 fb-1 • Deadtime: 1.0 fb-1 • Luminosity in 2012: • Delivered: 23.3 fb-1 • Atlas Recoded: 21.7 fb-1
RunNumber 215473 EventNumber39074754 CAL Burst Event ② TGC + CAL RPC + TGC + CAL MDT + RPC + TGC + CAL From Masato’s slides: TDAQ Week, 8 July 2013
TGC problem in 2012 • TGC ROD busy caused a lot of downtime • 0.35fb-1 of online data was lost Frequency of ROD busy later 2012 A-side 2012 pp runs 2013 pPb runs C-side From Masato’s slides: TDAQ Week, 8 July 2013
Analysis Outline Modified code to kill hits from adjacent bunches. For MC, run the trigger simulation. Thanks to Susumu Oda for instruction and help. • All the processing done by Reco_trf.py. RAW data 1 bunch readout 3 bunch readout ESD ESD AOD AOD Compare NTUP_SMWZ NTUP_SMWZ NTUP_TRIGMU NTUP_TRIGMU Compare MuonSpectrometer/MuonCnv/MuonTGC_CnvTools/src/TgcRdoToPrepDataTool.h and .cxx were modified.
Change in Athena Code To Only Use Current Bunch Hits in TGC • Data analysis: • MuonSpectrometer/MuonCnv/MuonTGC_CnvTools/src/TgcRdoToPrepDataTool.cxx • Changed voidMuon::TgcRdoToPrepDataTool::selectDecoder(constTgcRdo::const_iterator& itD, constTgcRdo * rdoColl) to simply return if (*itD->bcTag()=!TgcDigit::BC_CURRENT). • This change only erases prior and next bunch hits in reconstruction, and the trigger information in the event record is unchanged. • epp.phys.Kyushu-u.ac.jp/~oda/pukiwiki/index.php?DecodeOnlyCurrentBC • MC • MuonSpectrometer/MuonCnv/MuonTGC_CnvTools/src/TGC_RodDecoderReadout.cxx: • Changed code not to produce TgcRawData object for prior and next bunches. • Also changed MuonRdoToMuonDigitTool::decodeTgc( … ) in MuonSpectrometer/MuonCnv/MuonByteStreamCnvTest/src/MuonRdoToMuonDigitTool.cxx • So that it executes collection->push_back(newDigit) only for newDigit->bcTag()==TgcDigit::BC_CURRENT. • https://www.tsukuba.jp.hep.net/twiki/bin/view/Main/KojiSatoTrigSimKillingAdjacentBunchHits
Offline Muon Selection • Following MCP Guidelines. • Analyzed combination of MUID/STACO/MUONS and CB/CB+ST(/Muons for High Pt Searches). • Results are similar across different muon selection. • Today, I show results on MUONS, CB selection.
MUONS Selection • medium+ selection && !(isStandAloneMuon) • has good ID hits • nPixHits>0 • nSCTHits>4 • nSCTHoles<3 • n>5&&nTRToutliers<n for 0.1<|eta|<1.9 • with n=nTRToutliers+nTRThits • |d0_exPV|<1. • |z0_exPV|<10. • Pt>7 GeV MCP guidelines from analysis web page
TGC Hits in Adjacent Bunches in Data • sl_nvariables in NTUP_TRIGMU are zero for prior and next bunches. tgc_currentBC_sl_n Prior and next bunch hits are safely killed with 1 bunch readout scheme. tgc_nextBC_sl_n tgc_priorBC_sl_n
TGC Hits in Adjacent Bunches in MC • TGC_prd_bunch in NTUP_L1TGC. • 1:prior, 2:current, 3:next bunch. Prior and next bunch hits are safely killed with 1 bunch readout scheme.
Details of trigger muon counting • I simply count the number of EF/L2/L1 or offline CB muons that are found in the D3PD. • Rejecting multiple entries: • For trigger muons, I only count one muon when I find multiple muons with . • I typically find multiple L2 muons pointing to the same direction even after requiring a specific trigger bit. • I don’t find multiple EF muons pointing to the same direction when I require a specific trigger bit. • EF • Only counting EF muons associated with CB track. Kinematic plots are with CB track quantities. • L2 • Only counting combfeature L2 muons. • L1 • For counting, I simply count all L1 muons regardless of associated trigger bit – D3PD doesn’t have trigger bits for each L1 muon. Only have event bit.
Trigger Muon Yields ( MC) • Only fails to find one EF muon for mu4T trigger with one bunch readout.
Trigger Muon Yields (J/psi MC) • Only fails to find one EF muon for mu4T trigger.
MuonYields (Data) • #MUONS found in 3512939 Muon stream data events. 3bunch reco, Pt(mu)> 7 GeV/c 1bunch reco , Pt(mu)> 7 GeV/c Difference for Endcap CB Difference=(1bunch-3bunch)/3bunch *100%: 3bunch reco , Pt(mu)> 30 GeV/c 1bunch reco , Pt(mu)> 30 GeV/c Difference for Endcap CB Difference=(1bunch-3bunch)/3bunch *100%:
Muon Yields (Z->mumu MC) • Analyzed 249500 MC events. • No change observed in MUONS CB yield (Pt cut >7 or >30 GeV).
Muon Yields (J/psi->mumu MC) • Analyzed 499000 MC events. • Only <~0.002% change is observed, corresponding to one muon for Pt>7 GeVMUONS CB selection.
Time of Flight at TGC location • TGCinnerlayer at R~13 m