KoreaCmsWiki> OldWebs/OldWebs/DAQ Web>JeewonSeoLog (2014-07-22, GunmoRyu) EditAttach

Jeewon's Log

Work Summary

ttbar channels

There are three channels, 'Full-hadronic', 'semi-leptonic' and 'di-leptonic', in ttbar production. At first, I choose the full-hadronic channel.

full-hadronic channel

The fully hadronic final state is characterized by the six-jets topology ttbar->ww+bbbar->6jets. However, this channel is affected by a large background from QCD. The ttbar production can be main background of SUSY. Before analyzing the ttbar production as a background, I try to analyze the ttbar as a signal.

Simple example

I did three example for training. Each work is performed different source. Main training is the making plot for analysis.

Pt and Eta distribution check

This work is for learning what SusyAnalyzer consists of and how it works. To learn how to use the module in SusyAnalyzer, I did default simulation at first. And next, I added eta-cut.
  • Performed job
    • Analysis tool : SusyAnalyzer
    • CMSSW version : CMSSW_1_5_2
    • Cut : |η| < 2.4
    • Source : I used sample source in Susyanalyzer.
  • Result
    RecomuonPt.gif RecomuonEta.gif
    Pt and Eta distribution of Reconstructed muon
    MatchedmuonPt.gif MatchedmuonEta.gif
    Pt and Eta distribution of Matched muon
  • Conclusion
    For this simulation, I used not pt-cut but only eta-cut. So, Pt-distribution has to be from 0 GeV/c to 100 GeV/c and Eta-distribution has to be smaller than 2.4. But, there are no plot from 0 GeV/c to 10 GeV/c. Because default pt-cut of SusyAnalyzer is more than 10 GeV/c. If I want to change this pt-cut for another event, I can do that. And there are not only pt-cut but also many variable for analysis.

Measurement efficiency & purity by CRAB

There are two purpose of this work. The first one is for learning to use CRAB and the second one is for learning to calculate the efficiency and purity.
  • Performed job
    • Analysis tool : SusyAnalyzer
    • CMSSW version : CMSSW_1_5_2
    • CRAB version : CRAB_1_5_4
    • Cuts :
      • |η| < 2.4
      • Pt > 10 GeV
    • Event information :
        electron measurement muon measurement
      event source RelVal152Z-EE DrellYan_mumu40
      event number 5500 15000
  • Result :
    • electron
      Efficiency2.gif
      efficiency histogram of pt and eta
      purityelec1.gif
      purity histogram of pt and eta
    • muon
      MuonEfficiency.gif
      efficiency histogram of pt and eta
      Purity1.gif
      purity histogram of pt and eta
  • Conclusion :
    There are some information In PTDR. The efficiency has to be about 98% for muon and 75% for electron. The purity have to be about 100% for muon and electron. And my result is clear.

Z boson mass measurement

The purpose of this work is measurement of Z boson mass by calculating the invariant mass of muon.
  • Performed job :
    • Analysis tool : SusyAnalyzer
    • CMSSW version : CMSSW_1_5_2
    • CRAB version : CRAB_1_5_4
    • Cuts :
      • |η| < 2.4
      • Pt > 10 GeV
    • Data set : DrellYan_mumu40
    • Event number : 15000
  • Result :
    RecoInvMass.gif
  • Conclusion :
    Invariant mass of muon is about 87.4 . This is similar Z boson mass.

Monthly Summary

March 2008

  • Schedule
    • Korea CMS group
      • 7. ~ 8. Presentation of result about experiment

February 2008

  • Schedule
    • Meeting with Filip & Luc
      • 21. 11:00 discussion about result
    • SUSY Group Meeting
      • 6. 15:00~19:00
      • 20. 15:00~17:30
        1. PAT update
        2. QCD Etmiss tail from data
        3. top background - update
        4. W+jets for Z->nunu+jets - update
        5. Study of Z+jets background to SUSY
    • KOREA CMS Meeting
      • 11. 11:00~14:00 DAQ meeting (G.N.Kim, J.H.Chung, S.U.Ahn, J.W.Seo)
      • 13. 5:00~13:00 2nd workshop(1st day)
      • 14. 3:00~13:00 2nd workshop(2nd day)

January 2008

  • Schedule
    • Meeting with Filip & Luc
      • 15. 11:00 hemi-analysis for jet pairing
    • SUSY Group Meeting
      • 30. 15:00~18:00
        1. Introduction and 2008 plan of work
        2. H->bb in SUSY cascades
        3. Z0 candle calibration
        4. W+jets for candle calibration of Z->nu-nu
        5. Z0+jets as background for SUSY
        6. t-tbar background for susy
        7. Trigger study for jets+MET
    • Korea CMS Meeting
      • 23. 11:00 DAQ meeting ( S.K.Oh, G.N.Kim, J.H.Chung, G.S.Kim, S.U.Ahn, S.H.Song, J.W.Seo )
        1. Have to meet and discuss month by month.
        2. Send the report about work to professor monthly.
        3. When I am faced with difficulties, doing something or some physics..., ask Dr. Ro and Dr. chung or... someone.
    • PHYSICS DAY
      • 24. 11:00~17:40
        1. WW, WZ, ZZ
        2. t-tbar
        3. Single top
        4. SUSY: what is needed from the Standard Model
        5. Exotica: what is needed from the Standard Model
        6. Higgs: what is needed from the Standard Model
        7. Analysis issues in searching for things

December 2007

  • Schedule
    • Work presentation 3. 11:30
    • 10~14 CMS week
      • 10. 11:00 OPENING PLENARY SESSION
      • 12. 11:00 PH-SUSY
    • Report
      • 14. Classical dynamics-canonical transformation
      • 17. Particle physics - the Lagrangian formulation of mechanics
      • 20. Quantum mechanics - time-independent perturbation theory

November 2007

  • Schedule
    • SUSY/BSM meeting
      • 2. 14:30 WG3/WG4 meeting
      • 9. 14:20 Dijet analysis Pre-Approval
    • 7. 16:30 CMS Meeting

DAQ.October 2007

  • Schedule
    • Meeting with Filip
      • 3. 11:00 simple example
      • 10. 10:00 ETHZ meeting - Start working with ETHZ group
      • 16. 11:00 ttbar+jet events
    • SUSY/BSM meeting
      • 23. RS Graviton to ZZ to 4 muons(J.H.Chung)
    • 13. First Korea CMS group physics meeting
    • 25. Physics Plenary meeting - measurements of top quark pair events

September 2007

  • Schedule
    • Meeting with Filip
      • 4. 11:00 Simple example - to draw histograms of muon and electron
      • 10. 10:00 Changed the event source
      • 12. 11:00 to draw efficiency and purity plots
      • 26. 11:00 muon efficiency
    • CMS week
      • 17. Opening session
      • 18. PH-muon
      • 19. PH-SUSY/BSM

Tip

CRAB

  • Set up 
       scramv1 b
   source /afs/cern.ch/cms/LCG/LCG-2/UI/cms_ui_env.csh
   eval `scramv1 runtime -csh`
   source /afs/cern.ch/cms/ccs/wm/scripts/Crab/CRAB_x_x_x/crab.csh
   voms-proxy-init -voms cms
  • Perform Job
    Job creation crab -create -cfg (configuation file name)

    Job Submission crab -submit all -c

    Job Status Check crab -status -c

    Job Output retrieval crab -getoutput all -c

    Kill jobs crab -kill (job number) -c

  • ce & se status check 
       lcg-infosites --vo cms ce 
   lcg-infosites --vo cms se
    or 
       lcg-infosites --vo cms ce | grep xxx
   lcg-infosites --vo cms se | grep xxx

Castor

  • copy from local directory to castor
    rfcp (file name) /castor/cern.ch/user/s/seoclara/(directory name)/(new file name)

  • copy from castor to local directory
    rfcp /castor/cern.ch/user/s/seoclara/(directory name)/(file name) (new file name)

  • delete a file(or directory)in castor 
       rfrm /castor/cern.ch/user/s/seoclara/(directory name)/(file name)
   rfrm -r /castor/cern.ch/user/s/seoclara/(directory name)
  • create a directory in castor
    rfmkdir /castor/cern.ch/user/s/seoclara/(new directory name)

  • Show directory(or file) list 
       rfdir /castor/cern.ch/user/s/seoclara/
   nsls -l /castor/cern.ch/user/s/seoclara/(directory name)

  • If use the .csh file, write this command, ./(filename).

ROOT macro

  • choose the .root file for analysis 
       TFile f1(" filename ");
   TFile f2(" filename ");
  • Show the current directory 
       gDirectory->pwd();
    or 
       .pwd
  • Show the list in current directory 
       .ls
  • Move directory f1 
     f1->cd();
    If i want to move to subdirectory of f1, type like this cd( "subdirectory name") instead of cd().

  • Go to the home directory 
     gROOT->Ld();
  • Sum of histogram 
       TFile f1("filename");
   TFile f2("filename");
   TCanvas *c1=new TCanvas("c1","name",200,10,600,400);
   TH1D *h1=(TH1D*)f1.Get("plotname");
   TH1D *h2=(TH1D*)f2.Get("plotname");
   c1->SetLogy();     // if you don't want make the histogram in log-scale, ignore this line.
   h1->SetLineColor(4);
   h1->Draw();
   h2->SetLineColor(3);
   h2->Draw("same");

  • normalize 
       Double_t scale1 = (number)/h1->Integral();
   h1->Scale(scale1);

Stage request and check

Type following commands in working area. 
   stager_get -M /castor/cern.ch/cms/(source address) -M /castor/cern.ch/cms/(source address) ... -U (name)
   stager_qry -U (name)
In this case, response can be expected :

EQUEST_RAILED : This message means that there isn't such file or directory in using server.

EQUEST_READY : There is such file or directory.

STAGEIN : The request is success.

STAGED : The source file is staged and can be used.

Batch system

  1. Make .csh file

    This is code in batch file (.csh) 
       #!/bin/csh
   cd /afs/cern.ch/user/(myid)/(mydirectory)/...
   eval `scramv1 runtime -csh`
   cmsRun (batch file name)
  2. In working area, type below commands 
       chmod 744 (batch file name)
   bsub -q 1nm -J job1 < (batch file name)
   bjobs
    In this command, 1nm means time of running job. m : minute, h : hour, d : day, w : week. ( ex. 1nm, 1nh, 8nh, 1nd, 1nw )

-- JeewonSeo - 14 Jan 2008

Edit | Attach | Watch | Print version | History: r30 < r29 < r28 < r27 < r26 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r30 - 2014-07-22 - GunmoRyu
 
KOREA-CMS
This site is powered by the TWiki collaboration platform Powered by Perl This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding KoreaCmsWiki? Send feedback