-- MinkyooChoi - 11 Feb 2008

Minkyoo's Log

JulyXbyMinkyoo

Elliptic flow

Date 07/15/2010

NegT_cen01.png NegT_cen12.png NegT_cen23.png
NegT_cen01_pt.png NegT_cen12_pt.png NegT_cen23_pt.png


Date 07/08/2010


hadd merge_1000708.root $(for i in $(cat file_list.txt); do echo -n "rfio:$CASTOR_HOME/crab/20100629/$i "; done)

Date 07/07/2010

LFN >> PFN 

>> ./fm_cli.py --lfn=/store/mc/Spring10/Hydjet_Quenched_MinBias_2760GeV/GEN-SIM-RECO/MC_36Y_V7A-v1/0005/E82E5989-BE70-DF11-9952-001EC9AA9FCC.root

Date 06/30/2010

>>crab -create
>>crab -submit
>>crab -status

Date 06/29/2010

[lxplus241] /afs/cern.ch/user/m/minkyoo > nsls
crab
[lxplus241] /afs/cern.ch/user/m/minkyoo > nsls crab
100208
100302
[lxplus241] /afs/cern.ch/user/m/minkyoo > nsmkdir crab/20100629
[lxplus241] /afs/cern.ch/user/m/minkyoo > nsls crab
100208
100302
20100629
[lxplus241] /afs/cern.ch/user/m/minkyoo > nsls -l
drwxrwxrwx   3 minkyoo  zh                        0 Jun 29 20:54 crab
[lxplus241] /afs/cern.ch/user/m/minkyoo > nsls -l crab
drwxr-xr-x   0 minkyoo  zh                        0 Feb 10 20:02 100208
drwxrwxr-x  97 cms003   zh                        0 Mar 03 02:48 100302
drwxr-xr-x   0 minkyoo  zh                        0 Jun 29 20:54 20100629
[lxplus241] /afs/cern.ch/user/m/minkyoo > nschmod 775 crab/20100629
[lxplus241] /afs/cern.ch/user/m/minkyoo > nsls -l
drwxrwxrwx   3 minkyoo  zh                        0 Jun 29 20:54 crab
[lxplus241] /afs/cern.ch/user/m/minkyoo > nsls -l crab
drwxr-xr-x   0 minkyoo  zh                        0 Feb 10 20:02 100208
drwxrwxr-x  97 cms003   zh                        0 Mar 03 02:48 100302
drwxrwxr-x   0 minkyoo  zh                        0 Jun 29 20:54 20100629
[lxplus241] /afs/cern.ch/user/m/minkyoo > 

Date 06/16/2010

**How to install root-sys for Mac
sudo port install Xft2 root

** auto start .bashrc 
>>vi .bash_login
write below it in .bash_login
------------------------
#!/bin/bash

[ -f ~/.bashrc ] && . ~/.bashrc
--------------------------


New dataset 

*** Discussion title: Heavy Ion Physics

Dear all,

Our 200K Hydjet Minbias production request has been completed:

/Hydjet_Quenched_MinBias_2760GeV/Spring10-MC_36Y_V7A-v1/GEN-SIM-RAW
/Hydjet_Quenched_MinBias_2760GeV/Spring10-MC_36Y_V7A-v1/GEN-SIM-RECO

They are currently stored at IN2P3 and being transferred to MIT.
Actually, the production has been done several days ago (took some time
to move to T1). This time, our request was processed all over the US for
only 3-4 days, not just at MIT. According to dataops, it went very well,
no memory issues even at sites with average memory of 1.5GB/core. Very
encouraging!

Cheers,
Wei

You have to search "like *" instead of "like VALID*". The transfers are
not done yet from IN2P3 to MIT so it's not valid yet. It's quite slow. I
am checking what's going on.

Date 05/25/2010


- 중이온 MC 현황 (이중 MB 사건을 이용)
   현재 제너레이션하는 데이터는 없는 상황

- 전체 사건을 다 읽고, centrality 를 플롯
   데이터 : summer 09 데이터(/store/mc/Summer09/Hydjet_MinBias_Unquen_4TeV/GEN-SIM-RAW/MC_31X_V3-GaussianVtx_311_ver1/0003/FE36E526-E17A-DE11-988F-0019B9CADC3D.root)

   
- centrality bin의 정의. (1-18인지, 1-10인지 등등, 각각의 빈이 몇%centrality인지 등등)
   Centrality of an event or impact parameter is defined as the distance between colliding two  nulei.
   Measure impact parameter or centrality of an event is extremely important for event characterisation.

    
    centrality table.

   Centrality.png

    Which one is a centrality?

  Centrality_reco.png
  
   
- centrality빈 별로 dn/deta 플롯, dn/dphi플롯
- event plane의 정의 (어제 미팅보니 7-8가지 sub-detector별로 정보를 종합하여 만들고 있음)
   총 34개 
   Track(10), Ecal(6),Hcal(6),Ecal+Hcal(24)
- event plane의 각분포 (플랫하게 나와야 겠지, 아니면 늘 일정한 방향으로 충돌한다는 이야기이니)
- 재구성된 event plane - true event plane 각의 분포 (가우시안 같은 분포가 나오겠지..)
- centrality bin별로의 event plane 분해능 분포

Date 03/12/2008

  • CMSSW_2_1_11 version 사용
  • CMSSW/GeneratorInterface/HydjetInterface 가져옴 2111
  • b=7 이벤트 50개 생성
081210-hydjettest-sourcecode.jpg

hydjettest.jpg

  • module 창출 eta/dn 작성중인데... 모듈은 mkedanlzr는 대략 만들었지만 현재 단계는 GEN 단계라 RECO단계까지 가야 만들수있는것인지 의문
  • GEN 단계만 거친 hydjettest1.root안에는 어떤게 들어있는가???
  • > edmDumpEventContent /afs/cern.ch/user/m/minkyoo/cms2111/src/GeneratorInterface/HydjetInterface/test/hydjet1.root | less 로 확인결과
  • GenHIEvent?HepMCProduct?만 존재..(흠..이게 멀까??)
  • 현재 작성한 모듈 소스
081210-minana-sourcecode1.jpg 081210-minana-sourcecode2.jpg

  • d eta /d n 플랏 만들기.. 하려면 일단 GEN단계에서 가능한지 판별..
  • 아직 전반적인 모듈 및 기타 등등의 잡다한 지식을 조금더.. 아직도 많이 헷갈림..

Date 16/02/2009

  • Hydjet gen단계 이벤트
  • gate.sscc.uos.ac.kr 접속
  • CERN과 연동
  • kinit minkyoo
  • CMSSW_2_1_4 버젼 GenerationInterface? 부분 다운로드
  • path>> //home/minkyoo/cmssw214/src/GeneratorInterface/HydjetInterface/test아래로 들어가 Hydjet_analysis.cfg
  • Hydjet_analysis.cfg 소스에서 untracked int32 input 넘버는 1로 해야됨
  • untracked uint32 sourceSeed 는 RANDOM
  • sourceSeed 초기값은 98765(의미없는값임)
  • replace HydjetSource?.bFixed = 0-13까지 옵션을 줄 수있음
  • 수정 소스는 아래와 같음

090217-Hydjet_analysis-sourcecode.jpg

  • >>cmsRun hydjet_analysis.cfg
  • gen단계가 돌아가면 hydjet1.root 와 myhistograms.root 두개의 파일 생성
  • >>root -l myhistograms.root
  • root 에서 myhistogram.root 에 들어가면 eventana;1 과 subana;1 두개의 디렉토리 보임
  • eventana;1 아래에 phdNdEta;1 , phdNdPhi;1 , phdNdPt;1 , phdNdY;1 있음
  • phdNdEta;1 플롯을 보면

090217-myhistogram-phdNdEta.jpg

  • phdNdPhi;1 플롯을 보면

090217-myhistogram-phdNdPhi.jpg

  • PhdNdPt? 플롯

090217-myhistogram-phdNdPt.jpg

  • phdNdY;1 플롯

090217-myhistogram-phdNdY.jpg

*>> vi HydjetAnalyzer.cc 에서 eta , e , pt , phi , pz , y 정의를 보면

  • part_eta = (*p)->momentum().eta();
  • part_e = (*p)->momentum().e();
  • part_pt = (*p)->momentum().perp();
  • part_phi = (*p)->momentum().phi();
  • part_pz = (*p)->momentum().z();
  • part_y = 0.5*log((part_e+part_pz)/(part_e-part_pz));
  • phdNdEta->Fill(part_eta);
  • phdNdY->Fill(part_y);
  • phdNdPt->Fill(part_pt);
  • phdNdPhi->Fill(part_phi);

Date 18/02/2009

  • condor Gen단계
  • >>./condorSubmitMaker.pl [돌리고싶은.cfg] [Event 갯수]
  • 생성된 디렉토리 (현재 condor-cfgname-생성날짜)
  • 정한 Event 갯수 많큼 cndor-cfgname-생성날짜 디렉토리안에 hydhet_analysis.cfg가 생성됨
  • 각 디렉토리를 콘도르로 돌리는 방법은 >>condor_submit hydjet_analysis.condor

Date 23/02/2009

  • 현재 impact parameter = 0,2,4,6,7,8,10,12,14까지 각각 Event 1000개씩 생성
  • 2009Feb12_230800 -> b=7
  • 2009Feb18_234108 -> b=0
  • 2009Feb20_230703 -> b=2
  • 2009Feb20_232207 -> b=4
  • 2009Feb20_234352 -> b=6
  • 2009Feb21_001222 -> b=8
  • 2009Feb21_002749 -> b=10
  • 2009Feb21_004513 -> b=12
  • 2009Feb21_005624 -> b=14

  • merge 하는 법 ::::>>hadd [outputfilename] [originalfile1] [originalfile2] . . . . .[originalfile N] (N = unlimited)

  • 예상대로 impact parameter value가 작을수록 데이터 사이즈와 Entree event Number는 증가함

  • flot은 아직...

Date 24/02/2009

  • Merge 스크립트
    #!/bin/bash
    
    step=0
    filename=""
    srcfiles=" "
    
    
    while [ $step -le 999 ]
    do
     filename="./${step}/myhistograms.root "
     srcfiles="${srcfiles}${filename}"
    # echo "$srcfiles"
     let step=$step+1
    done
    
    hadd Merge_test.root $srcfiles
    ~
    ~
    ~
    

  • impact parameter value에 따른 histogram source
    {
    gROOT -> Reset();
    TFile* b0 = new TFile("./condor-hydjet_analysis-2009Feb18_234108/merge_b0_5500GeV_EN1000.root");
    TFile* b2 = new TFile("./condor-hydjet_analysis-2009Feb20_230703/merge_b2_5500GeV_EN1000.root");
    TFile* b4 = new TFile("./condor-hydjet_analysis-2009Feb20_232207/merge_b4_5500GeV_EN1000.root");
    TFile* b6 = new TFile("./condor-hydjet_analysis-2009Feb20_234352/merge_b6_5500GeV_EN1000.root");
    TFile* b7 = new TFile("./condor-hydjet_analysis-2009Feb12_230800/merge_b7_5500GeV_EN1000.root");
    TFile* b8 = new TFile("./condor-hydjet_analysis-2009Feb21_001222/merge_b8_5500GeV_EN1000.root");
    TFile* b10 = new TFile("./condor-hydjet_analysis-2009Feb21_002749/merge_b10_5500GeV_EN1000.root");
    TFile* b12 = new TFile("./condor-hydjet_analysis-2009Feb21_004153/merge_b12_5500GeV_EN1000.root");
    TFile* b14 = new TFile("./condor-hydjet_analysis-2009Feb21_005624/merge_b14_5500GeV_EN1000.root");
    
    TH1F *b0_5500GeV_EN1000 = (TH1F*)b0->Get("eventana/phdNdEta");
    TH1F *b2_5500GeV_EN1000 = (TH1F*)b2->Get("eventana/phdNdEta");
    TH1F *b4_5500GeV_EN1000 = (TH1F*)b4->Get("eventana/phdNdEta");
    TH1F *b6_5500GeV_EN1000 = (TH1F*)b6->Get("eventana/phdNdEta");
    TH1F *b7_5500GeV_EN1000 = (TH1F*)b7->Get("eventana/phdNdEta");
    TH1F *b8_5500GeV_EN1000 = (TH1F*)b8->Get("eventana/phdNdEta");
    TH1F *b10_5500GeV_EN1000 = (TH1F*)b10->Get("eventana/phdNdEta");
    TH1F *b12_5500GeV_EN1000 = (TH1F*)b12->Get("eventana/phdNdEta");
    TH1F *b14_5500GeV_EN1000 = (TH1F*)b14->Get("eventana/phdNdEta");
    TCanvas* h1 = new TCanvas("h1","First canvas" , 800 , 600);
    b0_5500GeV_EN1000 -> Draw();
    b2_5500GeV_EN1000 -> Draw("same");
    b4_5500GeV_EN1000 -> Draw("same");
    b6_5500GeV_EN1000 -> Draw("same");
    b7_5500GeV_EN1000 -> Draw("same");
    b8_5500GeV_EN1000 -> Draw("same");
    b10_5500GeV_EN1000 -> Draw("same");
    b12_5500GeV_EN1000 -> Draw("same");
    b14_5500GeV_EN1000 -> Draw("same");
    }
    
    
    

  • dNdEta plot

0090224-mergedNdEta.jpg

Date 25/02/2009

  • root Draw file-dNdPhi
TFile* b6 = new TFile("./condor-hydjet_analysis-2009Feb20_234352/merge_b6_5500GeV_EN1000.root");
TFile* b7 = new TFile("./condor-hydjet_analysis-2009Feb12_230800/merge_b7_5500GeV_EN1000.root");
TFile* b8 = new TFile("./condor-hydjet_analysis-2009Feb21_001222/merge_b8_5500GeV_EN1000.root");
TFile* b10 = new TFile("./condor-hydjet_analysis-2009Feb21_002749/merge_b10_5500GeV_EN1000.root");
TFile* b12 = new TFile("./condor-hydjet_analysis-2009Feb21_004153/merge_b12_5500GeV_EN1000.root");
TFile* b14 = new TFile("./condor-hydjet_analysis-2009Feb21_005624/merge_b14_5500GeV_EN1000.root");

TH1F *b0_5500GeV_EN1000 = (TH1F*)b0->Get("eventana/phdNdPhi");
TH1F *b2_5500GeV_EN1000 = (TH1F*)b2->Get("eventana/phdNdPhi");
TH1F *b4_5500GeV_EN1000 = (TH1F*)b4->Get("eventana/phdNdPhi");
TH1F *b6_5500GeV_EN1000 = (TH1F*)b6->Get("eventana/phdNdPhi");
TH1F *b7_5500GeV_EN1000 = (TH1F*)b7->Get("eventana/phdNdPhi");
TH1F *b8_5500GeV_EN1000 = (TH1F*)b8->Get("eventana/phdNdPhi");
TH1F *b10_5500GeV_EN1000 = (TH1F*)b10->Get("eventana/phdNdPhi");
TH1F *b12_5500GeV_EN1000 = (TH1F*)b12->Get("eventana/phdNdPhi");
TH1F *b14_5500GeV_EN1000 = (TH1F*)b14->Get("eventana/phdNdPhi");
TCanvas* h1 = new TCanvas("h1","dNdPhi canvas" , 800 , 600);
gStyle->SetOptStat(111111);
b0_5500GeV_EN1000->SetLineColor(2);
b2_5500GeV_EN1000->SetLineColor(4);
b4_5500GeV_EN1000->SetLineColor(6);
b6_5500GeV_EN1000->SetLineColor(8);
b7_5500GeV_EN1000->SetLineColor(10);
b8_5500GeV_EN1000->SetLineColor(12);
b10_5500GeV_EN1000->SetLineColor(14);
b12_5500GeV_EN1000->SetLineColor(16);
b14_5500GeV_EN1000->SetLineColor(18);

b0_5500GeV_EN1000 -> Draw();
b2_5500GeV_EN1000 -> Draw("same");
b4_5500GeV_EN1000 -> Draw("same");
b6_5500GeV_EN1000 -> Draw("same");
b7_5500GeV_EN1000 -> Draw("same");
b8_5500GeV_EN1000 -> Draw("same");
b10_5500GeV_EN1000 -> Draw("same");
b12_5500GeV_EN1000 -> Draw("same");
b14_5500GeV_EN1000 -> Draw("same");
leg = new TLegend(0.5,0.6,0.89,0.89);
leg->AddEntry(b0_5500GeV_EN1000,"B=0","l");
leg->AddEntry(b2_5500GeV_EN1000,"B=2","1");
leg->AddEntry(b4_5500GeV_EN1000,"B=4","l");
leg->AddEntry(b6_5500GeV_EN1000,"B=6","l");
leg->AddEntry(b7_5500GeV_EN1000,"B=7","l");
leg->AddEntry(b8_5500GeV_EN1000,"B=8","l");
leg->AddEntry(b10_5500GeV_EN1000,"B=10","l");
leg->AddEntry(b12_5500GeV_EN1000,"B=12","l");
leg->AddEntry(b14_5500GeV_EN1000,"B=14","l");
leg->Draw();

histo->SetyAxis(0,100000);

090225-mergedNdPhi.jpg

Date 26/02/2009

  • dNdPhi

090226-mergedNdPhi.jpg

root [14] f->lsA()
Error: Can't call TFile::lsA() in current scope (tmpfile):1:
Possible candidates are...
(in TFile)
(in TDirectoryFile)
*** Interpreter error recovered ***
root [15] f->ls()
TFile**         myhistograms.root
 TFile*         myhistograms.root
  TDirectoryFile*               eventana        eventana (HydjetAnalyzer) folder
   KEY: TH1D    phdNdEta;1
   KEY: TH1D    phdNdY;1
   KEY: TH1D    phdNdPt;1
   KEY: TH1D    phdNdPhi;1
  TDirectoryFile*               subana  subana (SubEventAnalyzer) folder
   KEY: TNtuple vertices;1
   KEY: TNtuple subs;1
   KEY: TNtuple events;1
  KEY: TDirectoryFile   eventana;1      eventana (HydjetAnalyzer) folder
  KEY: TDirectoryFile   subana;1        subana (SubEventAnalyzer) folder
root [16] TBrowser b
root [17] : created default TCanvas with name c1
: created default TCanvas with name c1

root [17]
root [17]
root [17]
root [17] f->cd("vertices")
Error in : Unknown directory vertices
(Bool_t)0
root [18] f->cd("subana/vertices")
Error in : Object vertices is not a directory
(Bool_t)0
root [19] f->cd("subana/event")
Error in : Unknown directory event
(Bool_t)0
root [20]




root [6]
root [6]
root [6] gSystem->Load("libFWCoreFWLite.so")
Error in : libFWCoreFWLite.so does not exist in .:/usr/local/root5.22.0/lib::/usr/local/root5.22.0/lib:/usr/local/root/lib:/opt/intel/impi/3.1/lib64:/opt/intel/cmkl/9.1.023/lib/em64t:/opt/intel/fce/9.1.051/lib:/opt/intel/cce/9.1.051/lib:/usr/local/KL_E14/pro/e14/lib/so:/opt/condor/lib:/usr/local/root/lib:/opt/intel/impi/3.1/lib64:/opt/intel/cmkl/9.1.023/lib/em64t:/opt/intel/fce/9.1.051/lib:/opt/intel/cce/9.1.051/lib:/usr/local/KL_E14/pro/e14/lib/so:/opt/condor/lib::/PBS/lib:/PBS/lib:/usr/local/root5.22.0/cint/cint/stl
(int)(-1)
root [7] .q
[minkyoo@gate test]$ cmsenv
[minkyoo@gate test]$ root -l
root [0] gSystem->Load("libFWCoreFWLite.so")
dlopen error: /usr/local/root/lib/vector.dll: cannot open shared object file: No such file or directory
Load Error: Failed to load Dynamic link library /usr/local/root/lib/vector.dll
Error in : failure loading library vector.dll for class vector >
(int)(-1)
*** Interpreter error recovered ***
root [1] gSystem->Load("libFWCoreFWLite.so")
(int)0
root [2] AutoLibraryLoader::enable()
root [3] TFile* f = new TFile("hydjet1.root")
root [4]
root [4]
root [4]

TTree 전체 구조 볼수 있는 방법이 멀까???????

브렌치 밖에 안보이..

root [5] TTree* g1=(TTree*)f->Get("Events")
root [6] g1->Print()
******************************************************************************
*Tree    :Events    :                                                        *
*Entries :        1 : Total =         6125496 bytes  File  Size =    1855448 *
*        :          : Tree compression factor =   3.30                       *
******************************************************************************
*Br    0 :EventAuxiliary :                                                   *
*Entries :        1 : Total  Size=        829 bytes  One basket in memory    *
*Baskets :        0 : Basket Size=      16384 bytes  Compression=   1.00     *
*............................................................................*
*Branch  :edmGenHIEvent_source__PROD.                                        *
*Entries :        1 : BranchElement (see below)                              *
*............................................................................*
*Br    1 :edmGenHIEvent_source__PROD.edm::EDProduct :                        *
*Entries :        1 : Total  Size=        709 bytes  One basket in memory    *
*Baskets :        0 : Basket Size=      16384 bytes  Compression=   1.00     *
*............................................................................*
*Br    2 :edmGenHIEvent_source__PROD.present :                               *
*Entries :        1 : Total  Size=        809 bytes  One basket in memory    *
*Baskets :        0 : Basket Size=      16384 bytes  Compression=   1.00     *
*............................................................................*
*Br    3 :edmGenHIEvent_source__PROD.obj :                                   *
*Entries :        1 : Total  Size=       3729 bytes  One basket in memory    *
*Baskets :        0 : Basket Size=      16384 bytes  Compression=   1.00     *
*............................................................................*
*Br    4 :edmGenHIEvent_source__PROD.obj.evts_ :                             *
*         | edmGenHIEvent_source__PROD.obj.evts__                            *
*Entries :        1 : Total  Size=       2668 bytes  One basket in memory    *
*Baskets :        0 : Basket Size=      16384 bytes  Compression=   1.00     *
*............................................................................*
*Br    5 :edmGenHIEvent_source__PROD.obj.evts_.sub_id :                      *
*         | sub_id[edmGenHIEvent_source__PROD.obj.evts__]                    *
*Entries :        1 : Total  Size=       1352 bytes  One basket in memory    *
*Baskets :        0 : Basket Size=      16384 bytes  Compression=   1.00     *
*............................................................................*
*Br    6 :edmGenHIEvent_source__PROD.obj.sel_ :                              *
*Entries :        1 : Total  Size=        818 bytes  One basket in memory    *
*Baskets :        0 : Basket Size=      16384 bytes  Compression=   1.00     *
*............................................................................*
*Branch  :edmHepMCProduct_source__PROD.                                      *
*Entries :        1 : BranchElement (see below)                              *



* seed 값이 다 틀려야 되는가??? random 넘버를 써야되는가??? 아니면 포문으로 1-1000까지해야되는가??

Date 05/03/2009

CERN(유럽공동핵과학연구소)에서 LHC에서 pbpb collision energy는 5.5TeV
로 예상되며 그에 따른 elltiptic flow 현상을 분석하고자 한다.
Elliptic flow는 과거 SPS,RHIC 등에서부터 꾸준히 분석되어왔고,
LHC에서는 핵자당 5.5TeV의 에너지를 갖기때문에 과거 연구결과에서
약하게나마 보였던 QGP 상태가 더욱 확고히 입증될 것으로 예상되며
이벤트 방법으로 분석할 예정이기에 그에 따른 현황을 보고하고자한다.  

Date 09/03/2009

 
제 4회 연구 성과 발표회 마쳤음

자료는 첨부자료에

Date 13/03/2009


-study for the 4 vector
-difference rapidity with sudo rapidity
-study a paper of v2(rapidity function)

Date 26/03/2009

Try to Creating NTuple but not yet

try again...

#include 

#include "HydjetAnalyzer.h"
#include "SimDataFormats/HepMCProduct/interface/HepMCProduct.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

// essentials !!!
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "HepMC/HeavyIon.h"

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "PhysicsTools/UtilAlgos/interface/TFileService.h"

#include "TFile.h"
#include "TH1.h"
#include "TH2.h"
#include "TNtuple.h" 
using namespace edm;
using namespace std;

/************************************************************************
 *
 *  This code is analyzing pyquen events produced with pyquen+analysis.cfg
 *
 *  Author: cmironov@lanl.gov
 *
 *************************************************************************/

 
HydjetAnalyzer::HydjetAnalyzer(const ParameterSet& pset)
  : phdNdEta(0), phdNdY(0), phdNdPt(0), phdNdPhi(0),
   modLabel_(pset.getUntrackedParameter("moduleLabel","source"))
{
  // constructor
}


//_______________________________________________________________________
void HydjetAnalyzer::beginJob( const EventSetup& )
{
  //runs at the begining of the job
std::cout << "======= this is begin job section ========" << std::endl;
  edm::Service fs;
  TH1::SetDefaultSumw2(true);

  phdNdEta      = fs->make("phdNdEta",";#eta;",100,-10.,10.);
  phdNdY        = fs->make("phdNdY",";y;",100,-10.,10.) ;
  phdNdPt       = fs->make("phdNdPt",";p_{T}(GeV/c);",100, 0.,10.) ;    
  phdNdPhi      = fs->make("phdNdPhi",";d#phi(rad);",100,-3.15,3.15);



  return ;
}

/* 
void HydjetAnalyzer::beginJob( const EventSetup& )
{
  //runs at the begining of the job

  edm::Service fs;
  TH2::SetDefaultSumw2(true);

  phdPtdPl      = fs->make("phdPtdPl",";PtPl;",100,0,10,100,0,10) ;
  return ;
}
*/                   
//______________________________________________________________________
void HydjetAnalyzer::analyze( const Event& e, const EventSetup& )
{
  //runs every event

   Handle EvtHandle;
   
   // find initial (unsmeared, unfiltered,...) HepMCProduct
   // by its label - HydjetSource, that is
// e.getByLabel( "source", EvtHandle ) ;
   e.getByLabel( modLabel_, EvtHandle ) ;
   
   //double part_eta, part_y, part_pt, part_phi, part_e, part_pz;
   const HepMC::GenEvent* myEvt = EvtHandle->GetEvent() ;
  
   	   	
  // 	TFile f("NTuple_Gen_b0.root","recreate");
   	int pid ;
   	double px,py,pz,part_eta,part_e,part_pt,part_phi,part_pz ;
   	
   	TNtuple *Ntu = new TNtuple("minple","Min_NTuple","pid:px:py:pz");
   	 
   	  for( HepMC::GenEvent::particle_const_iterator p = myEvt->particles_begin(); p != myEvt->particles_end(); p++ )
	 			 {
				  pid = (*p)->pdg_id() ;
	  
	 				px = (*p)->momentum().x();
	 				py = (*p)->momentum().y();
    			pz = (*p)->momentum().z();
    /*
    px       = (*p)->momentum().px();
    py       = (*p)->momentum().py();
    pz       = (*p)->momentum().pz();
    
    part_eta = (*p)->momentum().eta();
		part_e   = (*p)->momentum().e();
		part_pt  = (*p)->momentum().perp();
		part_phi = (*p)->momentum().phi();
		part_pz  = (*p)->momentum().z();
		*/
		
					Ntu->Fill(pid,px,py,pz);
    
	  		}
    
   // Ntu->Write();	 
   	                          
   	


   	
   	/*
   if(myEvt)
     {
       for( HepMC::GenEvent::particle_const_iterator p = myEvt->particles_begin();
	    p != myEvt->particles_end(); p++ )
	  {
	   int pid = (*p)->pdg_id() ; 
  //     std::cout << "PDG partle ID = "<end_vertex() && abs( (*p)->pdg_id() ) == 211)
//	   if( !(*p)->end_vertex() && abs(pid) == 211)
	     {
		part_eta = (*p)->momentum().eta();
		part_e   = (*p)->momentum().e();
		part_pt  = (*p)->momentum().perp();
		part_phi = (*p)->momentum().phi();
		part_pz  = (*p)->momentum().z();
		part_y = 0.5*log((part_e+part_pz)/(part_e-part_pz));

	       phdNdEta->Fill(part_eta);
	       phdNdY->Fill(part_y);
	       phdNdPt->Fill(part_pt);
	       phdNdPhi->Fill(part_phi);
           phdPtdPl->Fill(part_pt,part_pz);
  //         std::cout << pid << "\t" << part_pt << "\t" << part_pz << std::endl;
         }
	  }
     }

  HepMC::HeavyIon *hi = myEvt->heavy_ion();
  if ( hi ) {
     std::cout << "B = " << hi->impact_parameter() << std::endl;
  }
   return ;   
}

*/
}
//_____________________________________________________________
void HydjetAnalyzer::endJob()
{
  // executed at the end of the job 
  phdNdEta->Scale(phdNdEta->GetBinWidth(0));
  phdNdY->Scale(phdNdY->GetBinWidth(0));
  phdNdPt->Scale(phdNdPt->GetBinWidth(0));
  phdNdPhi->Scale(phdNdPhi->GetBinWidth(0));  
  //phdPtdPl->Scale(phdPtdPl->GetBinWidth(0));
  return ;
}

//define as a plug-in
DEFINE_FWK_MODULE(HydjetAnalyzer);

DQM

Date 02/04/2010


How to install CMSHI

 scramv1 p CMSSW CMSSW_3_1_1
   cd CMSSW_3_1_1/src/
   cvs co UserCode/CmsHi
   ./UserCode/CmsHi/Utilities/scripts/setup31X.sh

Date 09/29/2009

minkyoo]$ cmsDriver.py step2_MC1_2_HI -s DQM -n 3 --filein file:dcache:///pnfs/cmsaf.mit.edu/hibat/cms/mc/hydj
et_b4_310/HYDJET_MC31XV3_pt25_RECO_9.root  --eventcontent RECOSIM --conditions FrontierConditions_GlobalTag,MC_31X_V2::All --mc
DQM,ENDJOB
DQM
ENDJOB
Starting  cmsRun step2_MC1_2_HI_DQM_MC.py
29-Sep-2009 10:16:38 EDT  Initiating request to open file file:dcache:///pnfs/cmsaf.mit.edu/hibat/cms/mc/hydjet_b4_310/HYDJET_MC3
1XV3_pt25_RECO_9.root
29-Sep-2009 10:16:57 EDT  Successfully opened file file:dcache:///pnfs/cmsaf.mit.edu/hibat/cms/mc/hydjet_b4_310/HYDJET_MC31XV3_pt
25_RECO_9.root
%MSG-w Missing Dictionary:  PoolSource:source{*ctor*} 29-Sep-2009 10:16:58 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed
 to be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
%MSG-w Missing Dictionary:  PoolSource:source{*ctor*} 29-Sep-2009 10:16:58 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed
 to be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
%MSG-w Missing Dictionary:  PoolSource:source{*ctor*} 29-Sep-2009 10:16:58 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed t
o be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
%MSG-w Missing Dictionary:  PoolSource:source{*ctor*} 29-Sep-2009 10:16:58 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed
 to be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
%MSG-w Missing Dictionary:  PoolSource:source{*ctor*} 29-Sep-2009 10:16:58 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed
 to be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
%MSG-w Missing Dictionary:  PoolSource:source{*ctor*} 29-Sep-2009 10:16:58 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed t
o be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
%MSG-w Missing Dictionary:  AfterModuleConstruction 29-Sep-2009 10:17:08 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed
 to be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
%MSG-w Missing Dictionary:  AfterModuleConstruction 29-Sep-2009 10:17:08 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed to be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
%MSG-w Missing Dictionary:  AfterModuleConstruction 29-Sep-2009 10:11:17 EDT  pre-events
Could not find a Reflex dictionary for class 'std::vector'.  This class was registered as one which is supposed to be held by an edm::Event, LuminosityBlock, or Run but will not be available. Please check
 1) was a Reflex dictionary created for the class,
 2) if so was the package with the dictionary linked with all plugins that use that class,
 3) the file is from an old release and this data type has been removed from the present release.
%MSG
Begin processing the 1st record. Run 1, Event 1, LumiSection 1 at 29-Sep-2009 10:12:08 EDT
29-Sep-2009 10:12:11 EDT  Closed file file:dcache:///pnfs/cmsaf.mit.edu/hibat/cms/mc/hydjet_b4_310/HYDJET_MC31XV3_pt25_RECO_9.root
%MSG-s CMSException:  PostModule 29-Sep-2009 10:12:17 EDT Run: 1 Event: 1
cms::Exception caught in cmsRun
---- EventProcessorFailure BEGIN
EventProcessingStopped
---- ScheduleExecutionFailure BEGIN
ProcessingStopped
---- ProductNotFound BEGIN
getByLabel: Found zero products matching all criteria
Looking for type: FEDRawDataCollection
Looking for module label: rawDataCollector
Looking for productInstanceName:

cms::Exception going through module SiStripFEDMonitorPlugin/siStripFEDMonitor run: 1 event: 1
---- ProductNotFound END
Exception going through path dqmoffline_step
---- ScheduleExecutionFailure END
an exception occurred during current event processing
cms::Exception caught in EventProcessor and rethrown
---- EventProcessorFailure END


%MSG

=============================================

MessageLogger Summary

 type     category        sev    module        subroutine        count    total
 ---- -------------------- -- ---------------- ----------------  -----    -----
    1 Missing Dictionary   -w AfterModuleConst                       3        3
    2 Missing Dictionary   -w PoolSource:sourc                       6        6
    3 HLTConfigProvider    -e                                        4*       4
    4 HLTConfigProvider    -e AfterModuleConst                       1*       1
    5 HLTConfigProvider    -e OfflineDQMMuonTr                       1*       1
    6 CMSException         -s PostModule                             1        1
    7 fileAction           -s PoolSource:sourc                       2        2
    8 fileAction           -s PostModule                             1        1

* Some occurrences of this message were suppressed in all logs, due to limits.

 type    category    Examples: run/evt        run/evt          run/evt
 ---- -------------------- ---------------- ---------------- ----------------
    1 Missing Dictionary   pre-events       pre-events       pre-events
    2 Missing Dictionary   pre-events       pre-events       pre-events
    3 HLTConfigProvider    BeforeEvents     BeforeEvents     BeforeEvents
    4 HLTConfigProvider    pre-events
    5 HLTConfigProvider    pre-events
    6 CMSException         1/1
    7 fileAction           pre-events       pre-events
    8 fileAction           1/1

Severity    # Occurrences   Total Occurrences
--------    -------------   -----------------
Warning                 9                   9
Error                   6                   6
System                  4                   4

Minkyoo's Logbook

test

Topic attachments
I Attachment History Action Size Date Who Comment
JPEGjpg 081210-hydjettest-sourcecode.jpg r1 manage 89.2 K 2008-12-10 - 00:01 UnknownUser  
JPEGjpg 081210-hydjettest.jpg r1 manage 32.7 K 2008-12-10 - 00:01 UnknownUser  
JPEGjpg 081210-minana-sourcecode1.jpg r1 manage 65.0 K 2008-12-10 - 00:15 UnknownUser  
JPEGjpg 081210-minana-sourcecode2.jpg r1 manage 56.8 K 2008-12-10 - 00:16 UnknownUser  
JPEGjpg 090217-Hydjet_analysis-sourcecode.jpg r1 manage 177.9 K 2009-02-17 - 15:50 UnknownUser  
JPEGjpg 090217-myhistogram-phdNdEta.jpg r1 manage 55.4 K 2009-02-17 - 16:19 UnknownUser  
JPEGjpg 090217-myhistogram-phdNdPhi.jpg r1 manage 65.8 K 2009-02-17 - 16:25 UnknownUser  
JPEGjpg 090217-myhistogram-phdNdPt.jpg r1 manage 47.6 K 2009-02-17 - 16:57 UnknownUser  
JPEGjpg 090217-myhistogram-phdNdY.jpg r1 manage 54.3 K 2009-02-17 - 17:01 UnknownUser  
JPEGjpg 090224-mergedNdEta.jpg r1 manage 94.5 K 2009-02-24 - 00:13 UnknownUser  
JPEGjpg 090225-mergedNdPhi.jpg r1 manage 69.0 K 2009-02-25 - 08:51 UnknownUser  
JPEGjpg 090226-mergedNdPhi.jpg r1 manage 80.3 K 2009-02-25 - 21:30 UnknownUser  
JPEGjpg 090819-쌍제니-교수님과다같이식사.JPG r1 manage 8812.3 K 2009-08-21 - 20:53 UnknownUser  
JPEGjpg 090820-CAF-현용건모.JPG r1 manage 6554.4 K 2009-08-21 - 20:56 UnknownUser  
JPEGjpg 090820-CAF-현용건모1.JPG r1 manage 7866.0 K 2009-08-21 - 20:59 UnknownUser  
JPEGjpg 090820-빔파이프모두.JPG r1 manage 9825.5 K 2009-08-21 - 21:03 UnknownUser  
PNGpng Centrality.png r1 manage 128.0 K 2010-05-25 - 07:33 UnknownUser centrality-text
PNGpng Centrality_reco.png r1 manage 54.7 K 2010-05-25 - 07:34 UnknownUser centrality-reco
PowerPointppt Measurement_of_Elliptic_flow.ppt r1 manage 1070.0 K 2009-03-09 - 07:52 UnknownUser  
PNGpng NegT_cen01.png r1 manage 36.0 K 2010-07-16 - 03:24 UnknownUser  
PNGpng NegT_cen01_pt.png r1 manage 32.8 K 2010-07-16 - 03:25 UnknownUser  
PNGpng NegT_cen12.png r1 manage 31.9 K 2010-07-16 - 03:28 UnknownUser  
PNGpng NegT_cen12_pt.png r1 manage 31.9 K 2010-07-16 - 03:28 UnknownUser  
PNGpng NegT_cen23.png r1 manage 36.1 K 2010-07-16 - 03:29 UnknownUser  
PNGpng NegT_cen23_pt.png r1 manage 36.2 K 2010-07-16 - 03:29 UnknownUser  
Compressed Zip archivezip wi.zip r1 manage 11238.1 K 2010-04-21 - 05:33 UnknownUser  
Edit | Attach | Watch | Print version | History: r28 < r27 < r26 < r25 < r24 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r28 - 2010-07-16 - MinkyooChoi
 
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