Steps for forward jet production and analysis:

1. Inclusive forward jet analysis (generator level):

MC sample production / Input MC files:

CTEQ5L_QCD_Pt_15_20.cfg: Pythia config file (example) to produce a QCD sample (pT-hat bin 15-20 GeV/c) with CTEQ5L parton distribution function from PDFLIB

CTEQ6M_QCD_Pt_15_20.cfg: Pythia config file (example) to produce a QCD sample (pT-hat bin 15-20 GeV/c) with CTEQ6M parton distribution function from PDFLIB

MRST03_QCD_Pt_15_20.cfg: Pythia config file (example) to produce a QCD sample (pT-hat bin 15-20 GeV/c) with MRST03 parton distribution function from PDFLIB

The resulting root files are (in /castor/cern.ch/user/s/salim/pythia_data/CTEQ5L/ or MRST03/):

 QCD_Pt_15_20_CTEQ5L_AllEta.root, QCD_Pt_15_20_MRST03_AllEta.root, ..

Data analysis code:
ForwardParton.cc (CMSSW_1_6_6): Analyses the input MC files (see above) and (i) looks for the 2 partons just AFTER the hard scattering which are in HF & Castor acceptance, (ii) produces outgoing parton distributions: ET, Pt, Eta and Phi for HF&CASTOR.

ForwardParton.h : Header file for ForwardParton.cc

PythiaData_CTEQ5L_Pt_15_20.cfg : EDAnalyzer code to run: ForwardParton.cc

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/ForwardPartonData/):

Forward_QCD_CTEQ5L_Pt_15_20_AllEta_1.root, Forward_QCD_MRST03_Pt_15_20_AllEta_1.root...

Add_rootfile_and_norm_sigma_ForwardParton.C: This code converts the root files (see above) that were generated for each pT-hat bins into one root file for each model (CTEQ5L,CTEQ6M & MRST03) by using sigma normalization. It read sigma values from crossection_CTEQ5L_AllEta.dat, crossection_CTEQ6M_AllEta.dat & crossection_MRST03_AllEta.dat

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/ForwardPartonData/):

Forward_QCD_CTEQ5L_Pt_AllEta.root, Forward_QCD_MRST03_AllEta.root...

ForwardJet.cc: Analysis of GenJet BjoX1,X2 ranges HF & Castor (CMSSW_1_6_12). It produces HF&CASTOR GenJet BjoX1,X2 and highx parton PID.

The input files CSA07 QCD sample pT-hat bins: 15<pT<230 GeV/c [see here a typical GEN-SIM configuration file]

Code runing time for ~ 1. M : ~ 5 hour on Grid Computing ( DBS ,Ex:cmssrm.fnal.gov, gridka-dCache.fzk.de ...)

ForwardJet.h : Header file for ForwardJet.cc.

ForwardJet_QCD_pt_20_30.cfg : EDAnalyzer code to run: ForwardJet.cc . Code (example) for the Iterative Cone R=0.5 jet algorithm with CMSSW_1_6_12.

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/ForwardJet_Bjo_x1_x2/):

ForwardJet_QCD_Pt_20_30_1.root...

Add_rootfile_and_sigma_norm_HF_CASTOR_Jet_BjoXrange.C: This code converts the root files (see above) that were generated for each pT-hat bins into one root file by using sigma normalization. It read sigma values for pT-hat bins 15-230 GeV/c from CSA07_QCD_event_sigma_values.dat .
To see extra information about CSA07 QCD event sigma values of all pT-hat bins: CSA07_QCD_event_sigma_values_AllpT-hat_bins.pdf

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/ForwardJet_Bjo_x1_x2/)

HF_CASTOR_Jet_BjoXrange.root

Running time (for ~ 1.M evts):
~5 hour on Grid Computing ( DBS ,Ex:cmssrm.fnal.gov, gridka-dCache.fzk.de ...)

NLO Calculations:

The NLO calculations for HF jets with NLOJET++/fastNLO done by Klaus Rabbertz. A short usage description of how to evaluate the inclusive jets calculations and the code in Fortran/HBOOK or C++/ROOT can be found iCSA08 jets Twiki. Klaus Rabbertz produced tables for HF jets with different jet algorithm (SiSCone R=0.5 & FaskT D=0.4 ) at 14 TeV and NLOJET++/fastNLO . We produced FastNLO single inclusive HF jets dsigma/dpTdEta spectrum root file (/castor/cern.ch/user/s/salim/FastNLO_HF_Jets_root_file/) by using these tables and runing code description on the iCSA08 jets Twiki.

Plots macros:

inp ut root files: (see above) in /castor/cern.ch/user/s/salim/ForwardPartonData/
, /castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt_Analysis_Data/
for produce_spectrum_plots_Phytia_cteq5l_fNLO_sisCone05_HFjet.C (see below)

/castor/cern.ch/user/s/salim/ForwardJet_Bjo_x1_x2/
for produce_plots_x_ranges.C (see below)

produce_spectrum_plots_Phytia_cteq5l_fNLO_sisCone05_HFjet.C : This macro produces plots with different PDFs model predictions FastNLO(CTEQ5L,CTEQ6M,MRST03) and Phytia-parton level(CTEQ5L).
Also, it compares Phytia-parton level(CTEQ5L) with three theoretical models (FastNLO ,cteq5l, cteq6m, mrst03).
(1) Invariant pT-differential cross-section: (1/2pi pT)dsigma/dpT
(2) pT-differential yields (for a given integrated luminosity): dN/dpT )
(3) Comparisons plots: [dsigma(PDF1) -dsigma(PDF2)]/dsigma(PDF2)

produce_plots_x_ranges.C: This macro produces the Bjorken-x ranges for CASTOR & HF GenJet. Also, it produce parton (jet) id histograms for high-x range.

2. Inclusive forward jet analysis (full simulation)

Input MC files:
The input files CSA07 QCD sample pT-hat bins: 15<pT<230 GeV/c [see here a typical GEN-SIM configuration file] and a special production of 400k events by Kostas Kousouris and Robert Haris with SiSCone05 and FastKt04 algorithms not available in the CSA07 production.
Data analysis code:
SingleIncForwardJet.cc : Analysis of the Single Forward Jet spectrum in HF (CMSSW_1_6_12). It produces Single ForwardJet pT, ET, Eta & phi spectrum for Corrected CaloJet and GenJet in forward CMS area (HF, CASTOR in the future) Code runing time for ~ 1.2 M : ~ 6 hour on Grid Computing ( DBS ,Ex:cmssrm.fnal.gov, gridka-dCache.fzk.de ...)

SingleIncForwardJet.h : Header file for SingleIncForwardJet.cc

SingleForwardJet_L2L3CorFastJet04_Pt_20_30.cfg : EDAnalyzer Code (example) to run which applies combined L2 & L3 Jet Energy corrections for the FastKt D=0.4 jet algorithm with CMSSW_1_6_12.

SingleForwardJet_L2L3CorSiSCone05_Pt_20_30.cfg : EDAnalyzer Code (example) which applies combined L2 & L3 Jet Energy corrections for the SISCone R=0.5 jet algorithm with CMSSW_1_6_12.

SingleForwardJet_L2L3CorJetsIC_Pt_20_30.cfg : EDAnalyzer Code (example) which applies combined L2 & L3 Jet Energy corrections for the case of the SISCone R=0.5 Iterative Cone R=0.5 jet algorithm with CMSSW_1_6_12.

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt.tar/)FastkT04,ICone05 & SiSCone05/):

SingleIncForwardJet_QCD_Pt_20_30_1.root, SingleIncForwardJet_QCD_Pt_30_40_1.root ....
Add_rootfile_and_sigma_norm_singleforwardjet_IC_SiSC_FastKt.C : This code converts the root files (see above) that were generated for each pT-hat bins into one root file by using sigma normalization. It read sigma values for pT-hat bins 15-230 GeV/c from CSA07_QCD_event_sigma_values.dat .

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt_Analysis_Data/):
Single_Inc_ForwardJet_IC.root,Single_Inc_ForwardJet_SisCone05.root,Single_Inc_ForwardJet_FastKt04.root

Plots macros:
produce_spectrum_plots_HFjet.C : This macro produces plots with different PDFs model predictions and jets algo. (ICone05, SiSCone05&FastkT04) for CaloJet smeared spectrum and FastNLO Eeach one:
(1) Invariant pT-differential cross-section: (1/2pi pT)dsigma/dpT
(2) pT-differential yields (for a given integrated luminosity): dN/dpT
Also,
this macro produces " SmearedSpectrum _Single_Inc_ForwardJet_IC_SISCone_FastkT.root " for Unsmearing procedure in /castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt_Analysis_Data/
Unsmearing:
The reconstructed jet spectrum needs an unsmearing correction to account for energy resoulution of the HF calorimeter in a steeply falling spectrum with wide bins. To obtain the energy resolution correction factor (for SiSCone0 5) as a function of jet pT, we used the ansatz method describe in the CSA08 jet studies.
1. The ansatz function is smeared with the resolution (assumed Gaussian).
2. The convoluted function is used to fit the measured pT spectrum.
3. The unsmearing correction is obtained as the ratio of the spectrum predicted by the unsmeared function over that predicted by the smeared one.

Unsmearing.cc : This macro runs Unsmearing.h (below)

Unsmearing.h :
This macro fits the SmearedSpectrum with Ansatz function and creates Fitted Spectrum root file

Input File: "SmearedSpectrum_Single_Inc_ForwardJet_IC_SISCone_FastkT.root" in /castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt_Analysis_Data/
Output File: " FittedSpectrum_dNpT_HF_CaloJet_SISCone05.root , FittedSpectrum_dNpT_HF_CaloJet_FastkT04.root, FittedSpectrum_dNpT_HF_CaloJet_IC05.root " in /castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt_Analysis_Data/

SpectrumCorrection.cc : This macro calculates correction factors and corrects SmearedSpectrum by using correction factors. Also, creates UnsmearedSpectrum root file.

Input File: " FittedSpectrum_dNpT_HF_CaloJet_SISCone05.root , FittedSpectrum_dNpT_HF_CaloJet_FastkT04.root, FittedSpectrum_dNpT_HF_CaloJet_IC05.root " in /castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt_Analysis_Data/
Output
File: " UnsmearedSpectrum_dNdpT_HF_CaloJet_SISCone05.root , UnsmearedSpectrum_dNdpT_HF_CaloJet_FastkT04.root , UnsmearedSpectrum_dNdpT_HF_CaloJet_IC05.root " in/castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt_Analysis_Data/

produce_Unsmeared_spectrum_plots_HF.C :This macro produces plots with " Unsmeared " CaloJet (SiSCone R=0.5) and FastNLO (CTEQ6M & MRST03) Eeach one:
(1) Invariant pT-differential cross-section: (1/2pi pT)dsigma/dpT
(2) pT-differential yields (for a given integrated luminosity): dN/dpT

Jet energy scale (JES) error:

The most important source of systematic uncertainties of the single forward jet spectrum measurement is due to JES.

produce_plots_Comparison_PDF_and_JES_Uncertainty.C : This macro produces change of yield single Unsmeared forward jet spectrum (CaloJet SiSCone R=0.5) due to Jet Energy scale, also it compares data (Unsmeared CaloJet SiSCone R=0.5 ) with two theoretical models (FastNLO cteq6m & mrst03).
Input files : "UnsmearedSpectrum_dNdpT_HF_CaloJet_SISCone05.root" (above)
" fnlo_HF_eta_14TeV_SISCone_05_MRST2003cnlo.root, fnlo_HF_eta_14TeV_SISCone_05_cteq5l.root "
in castor/cern.ch/user/s/salim/SingleForwardJet_IC_SiSC_FastKt_Analysis_Data/

Inclusive forward jet HLT-L1Trigger:
Input MC files:

The input files are the CSA07/QCD_Pt_30_50/CMSSW_1_6_0-PreCSA07-HLT-A1/GEN-SIM-DIGI-RECO and /QCD_Pt_50_80/CMSSW_1_6_0-PreCSA07-HLT-B3/GEN-SIM-DIGI-RAW

Data analysis code:

SingleIncForwardJet_HLT.cc: This code determines which HLT triggers are fired by our QCD forward-jet events.We fill out histograms for the potentially relevant HLT streams (marked with "<-------" below): MB, QCD-jets: 0,13,14,15,16,87,88,89) as a function of the jet pT and eta

Trigger Name[ 0]: HLT1jet <--------- Threshold=200 GeV, Prescale = 1 Trigger Name[ 13]: HLT1jetPE1 <--------- Threshold=150 GeV, Prescale = 10 Trigger Name[ 14]: HLT1jetPE3 <--------- Threshold=100 GeV, Prescale = 100 Trigger Name[ 15]: HLT1jetPE5 <--------- Threshold= 60 GeV, Prescale = 10e4 Trigger Name[ 16]: CandHLT1jetPE7 <--------- Threshold=30 GeV, Prescale = 10e5

Trigger Name[ 87]: HLTMinBiasPixel <--------- Threshold= 0 GeV, Prescale = 3e6 (at L1) Trigger Name[ 88]: HLTMinBias <--------- Threshold= 0 GeV, Prescale = 3e6 (at L1) Trigger Name[ 89]: HLTZeroBias <--------- Threshold= 0 GeV, Prescale = 3e6 (at L1)

SingleIncForwardJet_HLT.h : Header file of SingleIncForwardJet _HLT.cc

ForwardJet_HLT_info_Pt_30_50.cfg : EDAnalyzer Code (example) to run which applies combined L2 & L3 Jet Energy corrections for the IterativeCone R=0.5 jet algorithm with CMSSW_1_6_12.:

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/ForwardJet_HLT_info_2/):

ForwardJet_HLT_info_QCD_Pt_50_80_1.root,ForwardJet_HLT_info_QCD_Pt_30_50_1.root....

Add_rootfile_and_sigma_norm_HLT.C : This code converts the root files in (castor/cern.ch/user/s/salim/ForwardJet_HLT_info/ForwardJet_HLT_info_QCD_Pt_50_80_1.root,.....) that were generated for each pT-hat bins into one root file by using sigma normalization. It read sigma values from CSA07_QCD_event_sigma_values.dat http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/SCerci/macros_note/jet_reco_performances/plots_macros/CSA07_QCD_event_sigma_values.dat

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/ForwardJet_HLT_info_1(pT_50_80) & ForwardJet _HLT_info_2(pT_30_50)/):

Single_ForwardJet_HLT_pT_30_80.root

SingleIncForwardJet _L1Trigger.cc: This code determines which L1Triggers are fired by our QCD forward-jet events.

We fill out histograms for the potentially relevant L1Triggers streams (marked with "<-------" below):

MB, QCD-jets: 21,23,25,26,118,119) as a function of the jet pT and eta

Trig #21 name L1_SingleJet15 <--------- Threshold=30 GeV, Prescale = 1e5 Trig #23 name L1_SingleJet30 <--------- Threshold=60 GeV, Prescale = 1e4 Trig #25 name L1_SingleJet70 <--------- Threshold=110 GeV, Prescale = 100 Trig #26 name L1_SingleJet100 <--------- Threshold=150 GeV, Prescale =1 Trig #27 name L1_SingleJet150 <--------- Threshold= ?GeV, Prescale = 1 Trig #118 name L1_CMS.MinBias_HTT10 <--------- Threshold=0 GeV, Prescale = 3e6 Trig #119 name L1_ZeroBias <--------- Threshold=0 GeV, Prescale = 3e6

SingleIncForwardJet_L1Trigger.h: Header file of SingleIncForwardJet _L1Trigger.cc

ForwardJet L1Trigger_info_Pt_30_50.cfg : EDAnalyzer Code (example) to run which applies combined L2 & L3 Jet Energy corrections for the IterativeCone _R=0.5 jet algorithm with CMSSW_1_6_12.:

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/ForwardJet_L1Trigger_info_1):

ForwardJet_L1Trigger_info_QCD_Pt_50_80_1.root ..,ForwardJet_L1Trigger_info_QCD_Pt_30_50_1.root ..

Add_rootfile_and_sigma_norm_L1Trigger.C: This code converts the root files in (castor/cern.ch/user/s/salim/ForwardJet_L1Trigger_info_1/ ForwardJet_L1Trigger_info_QCD_Pt_50_80_1.root .. ,ForwardJet_L1Trigger_info_QCD_Pt_30_50_1.root ..,) It read sigma values from CSA07_QCD_event_sigma_values.dat http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/SCerci/macros_note/jet_reco_performances/plots_macros/CSA07_QCD_event_sigma_values.dat

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/ForwardJet_L1Trigger_info_1/):

Single_ForwardJet_L1Trigger_pT_30_80.root

Plots macros:

input file: (see above) Single_ForwardJet_HLT_pT_30_80.root in castor/cern.ch/user/s/salim/ForwardJet_HLT_info_2/
produce_plots_HLT_TriggerEff.C : This macro produces the HLT Triggers Efficiency vs pT & CMS HF Eta

input file: (see above) Single_ForwardJet_HLT_pT_30_80.root in castor/cern.ch/user/s/salim/ForwardJet_L1Trigger_info_1/
produce_plots_L1TriggerEff.C : This macro produces the L1 Triggers Efficiency vs pT & CMS HF Eta


3. Jet reconstruction performance (HF) analysis:

Input MC files:

The input files are the CSA07 QCD sample pT-hat bins: "15-20", "20-30", "30-50", "50-80", "80-120", "120-170", "170-230" GeV/c [see here a typical GEN-SIM configuration file] and a special production of 400k events by Kostas Kousouris and Robert Haris with SiSCone05 and FastKt04 algorithms not available in the CSA07 production.

Data analysis code:

ForwardJet.cc: Analysis of the performances of jet reconstruction in HF with CMSSW_1_6_12. It produces CaloJetpT/GenJetpT, |CaloJetEta-GenJetEta|, CaloJetphi-GenJetphi distributions , jet response vs pT and HF eta range for HF Minus, HF Plus & AllHF, also pT, phi and eta spectrum for GenJet and Correct CaloJet.

Code runing time for ~ 1.2 M : ~ 6 hour on Grid Computing (DBS, Ex:cmssrm.fnal.gov, gridka-dCache.fzk.de ...)

ForwardJet.h: Header file of ForwardJet.cc.

L2L3CorJetsFastKt04_Pt_15_20.cfg: Code (example) which applies combined L2 & L3 Jet Energy corrections for the FastKt D=0.4 jet algorithm with CMSSW_1_6_12.

L2L3CorJetsIC_Pt_15_20.cfg : Code (example) which applies combined L2 & L3 Jet Energy corrections for the Iterative Cone R=0.5 jet algorithm with CMSSW_1_6_12.

L2L3CorJetsSisCone05_Pt_15_20.cfg : Code (example) which applies combined L2 & L3 Jet Energy corrections for the case of the SISCone R=0.5 jet algorithm with CMSSW_1_6_12.

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/HF_Jet_Performance_Data.tar/ICone05,SiSCone05 or FastKt):

ForwardJet_QCD_Pt_15_20_1.root,ForwardJet_QCD_Pt_20_30_1.root...

Add_rootfile_and_norm_sigma_pT_eta_phi_dist_IC_SiSC_FastKt.C: This code converts the root files (see above) that were generated for each pT-hat bins into one root file for three different algo (IterativeCone05,SISCone05 & FastKt04) by using sigma normalization. It read sigma values for pT-hat bins 15-230 GeV/c from CSA07_QCD_event_sigma_values.dat .
To see extra information about CSA07 QCD event sigma values of all pT-hat bins: CSA07_QCD_event_sigma_values_AllpT-hat_bins.pdf

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/HF_Jet_Performance_Analysis_Data/)

IC_ET_Eta_Phi_dist_L2_L3_corr.root, FastKt04_ET_Eta_Phi_dist_L2_L3_corr.root, SisCone05_ET_Eta_Phi_dist_L2_L3_corr.root

ForwardJet_Parton_Matching.cc : Analysis of jet-parton matching for GenJet & corrected CaloJet for HF (CMSSW_1_6_12). It produces outgoing HF parton (PYTHIA status == 3) pT and Eta spectrum, also matched GenJet&CaloJet outgoing parton pT & Eta spectrum, jet-parton response vs pT and Eta , DeltaR matching parameter for different pT bins.

Code runing time for ~ 0.7.M : ~ 4 hour on Grid Computing (DBS, Ex:cmssrm.fnal.gov, gridka-dCache.fzk.de ...)

ForwardJet_Parton_Matching.h :Header file of ForwardJet _Parton_Matching.cc.

ForwardJet_Parton_Matching_QCD_pt_20_30.cfg : Code (example) which applies combined L2 & L3 Jet Energy corrections for the ICone R=0.5 jet algorithm with CMSSW_1_6_12.

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/jet_parton_matching_data):

ForwardJet_QCD_Pt_15_20_1.root,ForwardJet_QCD_Pt_20_30_1.root...

Add_rootfile_and_sigma_norm_HF_parton_jet_matching.C: This code converts the root files (see above) that were generated for each pT-hat bins into one root file for IterativeCone05 by using sigma normalization. It read sigma values for pT-hat bins 15-230 GeV/c from CSA07_QCD_event_sigma_values.dat .

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/HF_Jet_Performance_Analysis_Data/)

HF_Parton_Jet_matching.root

Running time (for ~1.2 M evts):
~6 hour on Grid Computing (DBS, Ex:cmssrm.fnal.gov, gridka-dCache.fzk.de ...)

Plots macros:
input root files: (see above) in /castor/cern.ch/user/s/salim/HF_Jet_Performance_Analysis_Data/
"IC_pT_Eta_Phi_dist_L2_L3_corr.root, FastKt04 _pT_Eta_Phi_dist_L2_L3_corr.root, SisCone05_pT_Eta_Phi_dist_L2_L3_corr.root" for
produce_plots_fwd_jets_performance_HF.C & produce_plots_JetEff.C (below)
"HF_Parton_jet_matching.root" for produce_plots_parton_matching_Eff.C (below)

produce_plots_fwd_jets_performance_HF.C : This macro produces the pT, phi, eta resolution and response plots for ALLHF, HFPlus and HFMinus.

produce_plots_JetEff.C : This macro produces the Jet Matching Efficiency vs pT & CMS HF Eta

produce_plots_parton_matching_Eff.C : This macro produces the Jet-Parton Matching Efficiency and response vs pT & CMS Forward Eta . Also, It produces GenJet-CaloJet Matching parameter DeltaR distributions(parton-CaloJet)

4. Mueller-Navelet dijets analysis:

Input MC files:

The input files are the CSA07 QCD sample pT-hat bins: "15-20", "20-30", "30-50", "50-80", "80-120", GeV/c [see here a typical GEN-SIM configuration file] (GenJets-CMS.CaloJets) and

PreCSA07-HLT-A1/GEN-SIM-DIGI-RECO samples pT-hat bins :"30-50" GeV/c (only CaloJets).

Data analysis code:

dijet.cc : Analysis of dijet in HF (CMSSW_1_6_12). It produces dN_MNdijets_dET_Deta, dsigma_dijet_dDphi_Deta, cosDphi_vs_Deta cos2Dphi_vs_Deta plots for GenJet , CaloJet & CaloTower. Code runing time for full event (~ 4.7 M) : ~ 1 day on Grid Computing (DBS , Ex:cmssrm.fnal.gov, gridka-dCache.fzk.de ...)

dijet.h : Header file to use dijet.cc

dijet_QCD_pt_30_50.cfg : EDAnalyzer code to run: dijet.cc by using the combined L2 & L3 Jet Energy corrections for Iterative Cone R=0.5 with CMSSW_1_6_12.

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/MN_dijets_Data.tar/pT_30_50, pT_50_80 ...):

dijet_QCD_Pt_15_20_1.root,dijet_QCD_Pt_20_30_1.root...

Add_rootfile_and_sigma_norm_dijet.C : This code converts the root files (see above) that were generated for each pT-hat bins into one root file for jet algorithm (IterativeCone05) by using sigma normalization. It read sigma values for pT-hat bins 15-120 GeV/c from CSA07_QCD_event_sigma_values.dat .
To see extra information about CSA07 QCD event sigma values of all pT-hat bins: CSA07_QCD_event_sigma_values_AllpT-hat_bins.pdf

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/Analysis_MNdijet_Data/):

MN_dijet_genjet.root, MN_dijet_calojet.root

Running time (~ 4.7 M events):
~1 day on Grid Computing (DBS , Ex:cmssrm.fnal.gov, gridka-dCache.fzk.de ...)

Plots macros:
input root files: (see below) in /castor/cern.ch/user/s/salim/Analysis_MNdijet_Data/

produce_plots_mueller_navelet_dijets.C : This macro produces dN_MNdijets_dET_Deta, 1/N* dN_dDphi_Deta, cosDphi_vs_Deta,

cos2Dphi_vs_Deta and ratio of cos2Dphi and cosDphi vs Deta for Gen&Calo Jets.

MN dijets with Herwig Generator (genjet-level):
MC sample production / Input MC files:

Herwig_Genereted_QCDJets_pT_30_50.cfg : Herwig config file (example) to produce a QCD sample (pT-hat bin 30_50 GeV/c) with CTEQ5L parton distribution function from PDFLIB

The resulting root files are (in /castor/cern.ch/user/s/salim/HerwigQCDData_3, HerwigQCDData _4, HerwigQCDData _5):

 herwigQCD_pT_30_50_1.root ..,herwigQCD_pT_50_80_1.root ..,herwigQCD_pT_80_120_1.root ..

Data analysis code:

dijet.cc : Analysis of dijet in HF (CMSSW_1_6_12). It produces dN_MNdijets_dET_Deta, dsigma_dijet_dDphi_Deta, cosDphi_vs_Deta cos2Dphi_vs_Deta plots for GenJet and HERWIG Generator. Code runing time for full event (~ 4.6 M) : ~ 1 day (with "bsub" at CERN Computing)

dijet.h : Header file to use dijet.cc

dijet_QCD_pt_30_50_1.cfg : EDAnalyzer code to run: dijet.cc for Iterative Cone R=0.5 with CMSSW_1_6_12.

This analysis step produces the following files (in /castor/cern.ch/user/s/salim/HerwigQCDDataAnalysis_3, HerwigQCDDataAnalysis _5,HerwigQCDDataAnalysis_6):

dijet_QCD_Pt_30_50_1.root,dijet_QCD_Pt_50_80_1.root,dijet_QCD_Pt_80_120_1.root...

Add_rootfile_and_sigma_norm_dijet.C : This code converts the root files (see above )in (/castor/cern.ch/user/s/salim/Analysis_MN_dijet)

that were generated for each pT-hat bins into one root file for jetalgorithm (IterativeCone05) by using sigma normalization.It read sigma values from Herwig_Generated_QCDJets_pT_30_50_output.dat,Herwig_Generated_QCDJets_pT_50_80_output.dat,Herwig_Generated_QCDJets_pT_80_120_output.dat

on http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/SCerci/macros_note/Muller_Navelet_jets_studies/Herwig_MN_dijet/MC_sample_production/


This analysis step produces the following files (in /castor/cern.ch/user/s/salim/dijetData):

MN_dijet_genjet_herwig.root

Running time (~ 4.6 M events):

~1 day on with "bsub" at CERN Computing

Plots macros:

input root files: "MN_dijet_genjet_pythia.root", "MN_dijet_genjet_herwig.root" in /castor/cern.ch/user/s/salim/dijetData/

produce_plots_Herwig_and_Pyhtia_MN_dijets_match_parton_and_High_pT_GenJets.C : This macro produces ,cosDphi_vs_Deta and ratio of cos2Dphi and cos2Dphi vs Deta plots for GenJets for Mathed Parton, HighestpT and PYHTIA and HERWIG Generator



-- DavidDEnterria & Salim.Cerci- 17 Apr 2008

Edit | Attach | Watch | Print version | History: r48 < r47 < r46 < r45 < r44 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r48 - 2010-11-09 - DenizSunarCerci
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    CMSPublic All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback