Including the NN in h2gglobe

How to load the NN

Sanity-check : rerunning official results

(Instructions from Nicolas!)
  • The official code is in UserCode/HiggsAnalysis/HiggsTo2photons/h2gglobe
  • The tag suggested is baseline_workspace_08Dec2011_nw which should run the standard 4-categories version of the baseline analysis

  1. run the fitter.py python fitter.py -i datafiles.dat -n NUMBER -j JOB (so eg : for i in `seq 0 99`; do SUBMITTER( python fitter.py -i datafiles.dat -n 100 -j $i ) ; done (Note, you may have to set the variable DISPLAY to empty to avoid crashes: export DISPLAY
  2. run combiner.py: python combiner.py with the correct files in filestocombine.dat

Some notes :

  • The submission script to run on cern batch is : /afs/cern.ch/user/o/obondu/public/forMorgan/submission.sh
  • The output of python fitter.py -i datafiles.dat -n NUMBER -j JOB is stored on castor /castor/cern.ch/user/o/obondu/Higgs/7TeV/5fb/test/
  • The output of python combiner.py is /castor/cern.ch/user/o/obondu/Higgs/7TeV/5fb/test/CMS-HGG_4763pb.root

Now for the limit results :

  • The official datacard can be found here /afs/cern.ch/user/c/chanon/public/forMauro/cms-hgg-datacard.txt
  • The output of combine cms-hgg-datacard.txt -M Asymptotic -m 120 --generateBinnedWorkaround -S 1 ("-S 1" means that systematic uncertainties are incuded.) gives:

Sanity checks on the model: OK
Computing limit starting from observation
Will compute both limit(s) using minimizer Minuit2 with strategy 0 and tolerance 0.1
Median for expected limits: 1.54476
Sigma  for expected limits: 0.788156
Restricting r to positive values.

Make global fit of real data
NLL at global minimum of data: -84995.8 (r = 1.99846e-06)

Make global fit of asimov data
NLL at global minimum of asimov: -84994.9 (r = 3.56712e-05)
At r = 2.235139:   q_mu = 6.79976   q_A  = 7.43612   CLsb = 0.00456   CLb  = 0.54748   CLs  = 0.00833
At r = 1.117571:   q_mu = 2.29889   q_A  = 2.10419   CLsb = 0.06455   CLb  = 0.47325   CLs  = 0.13639
At r = 1.676355:   q_mu = 4.41133   q_A  = 4.46089   CLsb = 0.01785   CLb  = 0.50469   CLs  = 0.03537
At r = 1.396963:   q_mu = 3.30703   q_A  = 3.19116   CLsb = 0.03447   CLb  = 0.48706   CLs  = 0.07077
At r = 1.536659:   q_mu = 3.84982   q_A  = 3.80513   CLsb = 0.02487   CLb  = 0.49543   CLs  = 0.05021
At r = 1.606507:   q_mu = 4.12774   q_A  = 4.12767   CLsb = 0.02109   CLb  = 0.49999   CLs  = 0.04219
At r = 1.571583:   q_mu = 3.98801   q_A  = 3.96488   CLsb = 0.02291   CLb  = 0.49768   CLs  = 0.04604
At r = 1.554121:   q_mu = 3.91849   q_A  = 3.88447   CLsb = 0.02388   CLb  = 0.49656   CLs  = 0.04809
At r = 1.545390:   q_mu = 3.88387   q_A  = 3.84455   CLsb = 0.02437   CLb  = 0.49600   CLs  = 0.04914
At r = 1.541025:   q_mu = 3.86652   q_A  = 3.82456   CLsb = 0.02463   CLb  = 0.49572   CLs  = 0.04968

 -- Asymptotic -- 
Observed Limit: r < 1.5410
Expected  2.5%: r < 0.8381
Expected 16.0%: r < 1.1148
Expected 50.0%: r < 1.5448
Expected 84.0%: r < 2.1457
Expected 97.5%: r < 2.8508

Done in 0.62 min (cpu), 0.62 min (real)

Limit results from J. Tao

  • Steps: 1) prepared the RooWorkspace for mass point Mass between 110 to 150 with CMS-HGG_4763pb.root -- " root -b -q InterpolateMass.C\(${Mass}\)" with Normalization.C for SM Higgs; 2) "combine cms-hgg-datacard.txt -M Asymptotic -m ${Mass} --generateBinnedWorkaround -S 1", "-t Ntoys" can works; 3) "python limit-plotter-complete.py Asymptotic sm -r -s" to get the results; option -r mean the ratio and -s means smoothing, you can choose none(abs limit)/one(ratio or smoothing)/both (ratio and smoothing) of these 2 options.
  • Limit plots: limit_sm_Asymptotic.gif and ratio limit_sm_Asymptotic_ratio.gif, 3 strange points at mass ~128, ~137 and ~147GeV for the expected result, being checked. The Workspace and data for these point are there. Any other suggestions? (no effect from "-t Ntoys" I have tried (-t 10). : ( ) After smooth: limit_sm_Asymptotic_smooth.gif and ratio limit_sm_Asymptotic_ratio_smooth.gif

Modifications for including NN

(From Olivier, 30/01/12)

  • the HiggsAnalysis/HiggsTo2 photons compiles with the following packages:
Test Release based on: CMSSW_4_2_8_patch2
Base Release in: /afs/cern.ch/cms/slc5_amd64_gcc434/cms/cmssw-patch/CMSSW_4_2_8_patch2
Your Test release in: /afs/cern.ch/user/o/obondu/scratch0/Higgs/GlobeTest/CMSSW_4_2_8_patch2__baseline_workspace_08Dec2011_nw
--- Tag ---    -------- Package --------                        
V06-22-00      CondFormats/DataRecord                           
V00-03-01      CondFormats/EgammaObjects                        
baseline_workspace_08Dec2011_nw HiggsAnalysis/HiggsTo2photons                    
regression_Dec3d RecoEgamma/EgammaTools                           
---------------------------------------
total packages: 4 (4 displayed)

  • The NN function has to be included in GeneralFunctions_cc.h and GeneralFunctions_h.h files (see files in /afs/cern.ch/user/o/obondu/public/forHgg)
  • The NN weights (latest version from Hugues) is to be included in weights/TMVAClassification_MLP.class.C (see files in /afs/cern.ch/user/o/obondu/public/forHgg)
  • The input NN variables have to be recognized by globe (note that they ARE in V11_04_05 version of the non-reduced globetuples located in /castor/cern.ch/user/c/cmshgg/processed/V11_04_05): several files in the branchdef/ folder have to be modified (see files in /afs/cern.ch/user/o/obondu/public/forHgg)
  • This way both the make in the h2gglobe folder should do it and the scramv1 build in the CMSSW folder

  • Nicolas' core modifications are available here : /afs/cern.ch/user/c/chanon/public/forOlivier

Add NN variables to the reduction step

  • /sps/cms/chanon/save/ARCHIVE_OldCMSSWreleasesLxplus/CMSSW_4_2_8_patch7/src/HiggsAnalysis/HiggsTo2photons/h2gglobe/ contains an old installation of globe (LP?), with fixed pt cuts, and no energy regression
  • Example: go to /sps/cms/chanon/save/ARCHIVE_OldCMSSWreleasesLxplus/CMSSW_4_2_8_patch7/src/HiggsAnalysis/HiggsTo2photons/h2gglobe/PhotonAnalysis_scripts/ReduceGJet/ and create the jobs with createH2GGlobeReducer.bash then go in the Jobs directory, you can launch them with bsub -q 1nh RunCastorH2GGlobeReducer_0,1,2...

* (Tao) The root files locate in /castor/cern.ch/user/j/jtao/Hgg2011/h2gglobe_V11_04_05_reduction_jtao/ with the NN variables (pho_r19, pho_maxoraw, pho_cep, pho_lambdaratio, pho_lambdadivcov, pho_etawidth, pho_brem, pho_smaj, pho_e2x2 and pho_e5x5):

1) Data: PhotonRun2011A_Clean and PhotonRun2011B_Clean, 50 jobs each, with root files in /castor/cern.ch/user/j/jtao/Hgg2011/h2gglobe_V11_04_05_reduction_jtao/Data/ and the event statistic information data_stat.log. All the 100 jobs (50jobsx2samples) were copied to the castor successfully with the log file data_scp.log. (There are other two data samples, PhotonRun2011A and PhotonRun2011B in /castor/cern.ch/user/c/cmshgg/processed/V11_04_05/Data/, what's the difference? What's the meaning of "Clean"?)

2) MC samples incluing bkg and Hgg signal samples: 20 jobs for each bkg sample (BoxPt10to25/BoxPt25to250/BoxPt250/DiPhotonJets/DYJetsToLL_M50/GJet_Pt-20_pp /GJet_Pt-20_pf/QCDPt30to40_pp/QCDPt30to40_pf/QCDPt30to40_ff/QCDPt40_pp/QCDPt40_pf/QCDPt40_ff) and 1 job for each signal samples, GluGluToHToGG/VBF_HToGG/TTH_HToGG/WH_ZH_HToGG with mass points (90 95 100 105 110 115 120 121 123 125 130 135 140 145 150 155 160)GeV. All the MC root files locate in /castor/cern.ch/user/j/jtao/Hgg2011/h2gglobe_V11_04_05_reduction_jtao/MC_Sig_Fall11_S6. mc_stat.log show the event statistic information. All the 328 jobs (20jobsx13BkgSamples+4jobs*17MassPoints) were copied to the castor with the log file mc_scp.log.

Create diphoton minitrees

  • Nicolas hacked the globe to produce diphoton minitrees after the preselection. Changes are in h2gglobe/PhotonAnalysis/src/statanalysis.cc to create a minitree, including many isolation variables and the NN output for the leading and trailing photon. The selection was set from l.phoSUPERTIGHT to l.phoNOCUTS.
  • To run the minitree, you have to go in :
/sps/cms/chanon/save/ARCHIVE_OldCMSSWreleasesLxplus/CMSSW_4_2_8_patch7/src/HiggsAnalysis/HiggsTo2photons/h2gglobe/PhotonAnalysis_scripts/ReduceGJet/ you create the jobs with createH2GGlobeLooper.bash then go in Jobs_Looper directory, launch the jobs with bsub -q 1nh RunCastorH2GGlobeLooper_0,1,2...

Optimization code

  • Nicolas is working on it

Getting the brazilian flag plot

  1. Run interpolation between mass points in h2gglobe/Macros/InterpolateMass.C
  2. In recent versions of globe: h2gglobe/Macros/limit-plotter-complete.py.

Limit result with re-optimized CiC4+NN (by J. Tao)

* I compared the results with the following cases:

(1) Rerunning with the official RooWorkspace: CMS-HGG_4763pb.root (official)

(2) Running from the reduction steps (Data & MC processed samples V11_04_05 with additional NN variables), the reproduce the RooWorkspace and same CiC selections. (Re CiC)

(3) Analysis with re-optimized CiC4 and photon NN selections, input without shift (CiC+NN )

(4) Same as (3), but MC S4ratio rescaling, 1.007 (EB) and 1.009(EE) from Fan (CiC+NNSF, SF means shift)

* 20120202_GlobalAnaWithNN.pdf with a bug for the data analysis,, forgot to use the Jason files for the data analysis.

* 20120205GlobalAnaWithNNDebug.pdf with the bug fixd.

* The selections with reoptimized CiC4+NN CANNOT give a better even equivalent result. How to do the re-optimized CiC4+NN in the global framework?

Limit result with MVA (by J. Tao)

* Categories based on MIT BDT output in Globe (5 cats with keeping lowest BDT cat, no VBF cat): limit_sm_Asymptotic_ratio_MVA.gif.gif.

The Production spreadhseet

The spreadsheet contains up-to-date status of the 4 steps: IpnTreeProduction, photon miniTrees, diphoton miniTrees, plots

Cut-based analysis

MVA analysis

*The macro that Lyon uses for training TMVAClassification.C: TMVAClassification.C

* From Tao: Phontons were selected with ET>20 and |eta|<2.5 excluding [1.4442,1.566] , PixelSeed veto and some lSO selections, hoe<0.05 && IsoEcal<4.2 && IsoHcal<2.2, in PhotonCandDis.cxx after the event analysis with my own analyzer. And then selections, TrkISO<2.0 && SigIetaIeta<0.01(EB)/0.03(EE) && MCtruth selections, were used in the training TMVAClassification.C

* MC Samples used for training and test by Tao: Based on my presentation on 19th May, the "weight files" from the MVA analysis were obtained with Fall10 samples after the Jinan FCPPL workshop. The prompt photon were selected from the following DiphotonBox & Born MC samples:

/DiPhotonBorn_Pt*_TrackingParticles_7TeV-pythia6/Fall10-START38_V12-v1/GEN-SIM-RECODEBUG   (*=10to25, 25to250, 250toinf)       
/DiPhotonBox_Pt*_TrackingParticles_7TeV-pythia6/Fall10-START38_V12-v1/GEN-SIM-RECODEBUG   (*=10to25, 25to250, 250toinf)
The non-prompt photon were selected from the following QCD MC samples:
/QCD_Pt-*_EMEnriched_TuneZ2_7TeV-pythia6/Fall10-START38_V12-v1/GEN-SIM-RECO (*=20to30,30to80,80to170)

For the Spring11 samples, the GJet and QCD with DoubleEMenriched filter are available and further 2 photons requirement with pT>20GeV and |eta|<2.6. The statistic are not enough. These samples are used to validate the shape of the NN output, as shown in my presentation on 19th May (TaoSlides.ppt) or 20110525 NN Higgs. The corresponding Spring11 MC samples as Fall10 described above are still in process. Unfortunately always CRAB errors happened (60317).

Zmumugamma

data mc plots

reco files

  • /castor/cern.ch/user/o/obondu/Zmumugamma/Run2011/42X/Run2011A-May10ReReco-v1.root
  • /castor/cern.ch/user/o/obondu/Zmumugamma/Run2011/42X/Run2011A_Prompt_v4_June17.root
  • /castor/cern.ch/user/o/obondu/Zmumugamma/Run2011/42X/Run2011A_Prompt_v4_June10.root

Codes

Recipes

Fake photons from data

Data and selections

With the 2010 (36pb-1) and 2011 (3.013fb-1, up to run 177515 with 07Oct. jason file) data samples, after the kinematic requirements ET>30GeV, |eta|< 2.5 excluding the transition region (1.4442, 1.566), we select the photons except the 2 most isolated photons with the combinediso method. With the total 3.139fb-1 data, we select 14149 photons, including 4474/5589 unconverted/converted in Barrel and 1851/2235 unconverted/converted in Endcap. If we ask the Sigma_IetaIeta cuts, there are 2953 photons left, 953/930 unconverted/converted in Barrel and 503/567 unconverted/converted in Endcap. For the comparisons of the NN inputs with MC, we select the fake photons after the matching with MC truth, from the GJet (/GJet_Pt-20_doubleEMEnriched_TuneZ2_7TeV-pythia6/Summer11-PU_S4_START42_V11-v1/AODSIM) and QCD (/QCD_Pt-*_doubleEMEnriched_TuneZ2_7TeV-pythia6/Summer11-PU_S3_START42_V11-v2/AODSIM with *=30to40, 40) samples

Fake photons from MC were selected with the basic selections (ET, Eta, spike removal and pixel seed veto), the ISO selections and without matching any real photons from MC truth, which means the same bkg selections as we used in MVA analysis.

Data/MC comparisons plots:

* No EGM-10-006 loose selections (basic selections+two most isolated photons veto on real data): FakePhotonDistributionsDataMC_EB_Converted_NoISO.pdf, FakePhotonDistributionsDataMC_EB_Unconverted_NoISO.pdf, FakePhotonDistributionsDataMC_EE_Converted_NoISO.pdf and FakePhotonDistributionsDataMC_EE_Unconverted_NoISO.pdf. The distributions are much different.

* With only Sigma_IetaIeta cut (basic selections+two most isolated photons veto + Sigma_IetaIeta on real data): FakePhotonDistributionsDataMC_EB_Converted_OnlySigIeIe.pdf, FakePhotonDistributionsDataMC_EB_Unconverted_OnlySigIeIe.pdf, FakePhotonDistributionsDataMC_EE_Converted_OnlySigIeIe.pdf and FakePhotonDistributionsDataMC_EE_Unconverted_OnlySigIeIe.pdf. Some distributions are different, esp the etaWidth.

* I also draw the pT, R9 and NN outputsin these plots. The NN outputs were obtained by my TMVA analysis with 4 categories and CFMlpANN classifier with 8 inputs.

* Basic selections + two most isolated photons veto + Sigma_IetaIeta + EcalISO<4.2 on real data: FakePhotonDistributionsDataMC_UB_BasicSigVetoEcalISO.pdf, FakePhotonDistributionsDataMC_CB_BasicSigVetoEcalISO.pdf, FakePhotonDistributionsDataMC_UE_BasicSigVetoEcalISO.pdf, FakePhotonDistributionsDataMC_CE_BasicSigVetoEcalISO.pdf, with~370/300 for the unconverted/converted photon in barrel and ~710/790 unconverted/converted photon in endcap for all 2011 data.

check with CiC4 selections

1) The photons were selected with the CiC4 selections but failed any of the CiC ISO selections (RelativeCombinedISO_SelectedVertex,RelativeCombinedISO_WorstVertex,RelativeTrackISO). The same selections were used for both data and MC, to check the signal contamination from such selections. For the EB case, the signal fraction is about ~11% , seeing from plot FakePhotonDistributionsDataMC_FailAnyCiCisoEB.pdf.

2) With the CiC4 selections except the ISO selections, the relations between the etaWidth and the isolation variables in EB were plotted here for the signal photons and fake photons respectively, EtaWidth_RelCombinedISO_EB_True.gif/ EtaWidth_RelCombinedISO_EB_Fake.gif for the RelativeCombinedISO_SelectedVertex case with 2D scattering points and the TProfile distribution in the same plot. For the fake case, more relations between etaWidth and isolation variables are list as follows, EtaWidth_RelTrackISO_EB_Fake.gif for the relative track isolation in CiC4, EtaWidth_TrackISODR04RhoCorr_EB_Fake.gif for the absolute track isolation in DR=0.4 after rho correction, EtaWidth_EcalISODR04RhoCorr_EB_Fake.gif for the absolute ECAL isolation in DR=0.4 after rho correction, and EtaWidth_HcalISODR04RhoCorr_EB_Fake.gif for the absolute HCAL isolation in DR=0.4 after rho correction;

3) The results should/will be check with the global analysis.

code for cross check CIC

Where to find it

/UserCode/hbrun/higgsAnalysis

How to use it ? ( ALERT! in SPS )

  1. create the CMSSW area
scramv1 project -n CMSSW_4_2_3_test CMSSW CMSSW_4_2_3

  • go in src and download the macro
 cd CMSSW_4_2_3_test/src/ 
 cvs co -r defautTag -d higgsAnalysis UserCode/hbrun/higgsAnalysis
 cd higgsAnalysis

  • set the environement and compile
 source prepareEnv.(c)sh 
 source compile_dir.csh 

  • run the code
 ./theEXE.exe 

what is the result ?

list of passing event are here passingEvents.log

we have 22299 event (should have 28500)

what are the input files :

sample name = /GluGluToHToGG_M-120_7TeV-powheg-pythia6/Summer11-PU_S3_START42_V11-v2/AODSIM

IPN tree upples : /pnfs/in2p3.fr/data/cms/t2data/store/user/obondu/GluGlu_120_v2/GluGlu_120_*.root

list of the CIC variables of the first 10000 events

the file is here : CICvarsValue.log

var0 combIso default vtx *50/Et
var1 worst comb Iso *50 /Etworst vtx
var2 iso tracker default vtx*50/Et
var3 sig Ieta Ieta
var4 HoE
var5 R9
var6 dR to track

Coef to fix the cluster shape variables

More details in ClusterShapeFix

BARREL ENDCAP
R9 value 1.0048 value 1.00492
E2x2/E5x5 value 1.00708 value 1.008074

Recipes to run everything

How to run IpnTreeProduction

How to run photon miniTrees

How to run diphoton miniTrees

How to run the plots

How to run discovery/exclusion plots

Presentations

* NNworkingslides.pdf: Fan Jiawei and Hugues slides in 19 of may meeting

* TaoSlides.ppt: Tao slides in 19 of may meeting

* 20110525 NN Higgs: Tao slides in 26 of may meeting. For the Hgg part, there is a bug for the normalization of Born MC samples. I forgot to multiply the matching efficiency ( 0.35, madgraph samples). I corrected this in my presentation for the IHEP weekly meeting 20110527 NN Higgs. With K-factor, data and MC agree well. The mass distribution of Mgg changed with the NN cut, please have a look at these 2 pdf files, TaoDiPhotonDataMC_mH115_NN.pdf and TaoDiPhotonDataMC_mH120_NN.pdf.

* NNshapeTests.pdf: Fan Jiawei and Hugues slides with tests about NN output

* 20110601NNadditinalVar.ppt: Some plots from J. Tao and H. Xiao for the additional variable. Jiawei and Hugues, please check the additional variable λ-/σηη2, including the data/mc agreement and its performance in MVA analysis.

* 201000615 NN Higgs.pdf: given by Tao on 15th June. I have checked the definition of Lambda+-. It's true that we should use the Cee, Cpp and Cep for the calculation, not the sigma ones. I am trying with the corrected variables and redoing the MVA analysis.

* NNshapeTests.pdf: Some slides from Fan and Hugues in order to understand the NN output shape, some DATA/MC comparisons and some performance checks

Cross check with training the NN with diphotonborn

* signal dataset /DiPhotonBorn_Pt10to25_TrackingParticles_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

/DiPhotonBorn_Pt25to250_TrackingParticles_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

/DiPhotonBorn_Pt250toinf_TrackingParticles_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

(dataset for comparison /GJet_Pt-20_doubleEMEnriched_TuneZ2_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM )

*background dataset /QCD_Pt-20to30_EMEnriched_TuneZ2_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

/QCD_Pt-30to80_EMEnriched_TuneZ2_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

/QCD_Pt-80to170_EMEnriched_TuneZ2_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

*statistic = DiPhotonBorn reduced to have after cut the same stat that QCD EM Enriched after cuts

*I put some weights for training to take in account the differents X section of each sample

*what are the cuts = pho_et>20 && (abs(pho_eta)<=2.5)&&(!((abs(pho_eta)>1.4442)&&(abs(pho_eta)<1.566))) && pho_hoe<0.05 && pho_hasPixelSeed==0 && pho_IsoHollowTrkCone<2 && pho_IsoEcalRechit<4.2 && pho_IsoHcalRechit<2.2 && ((pho_isEB==1&&pho_sigmaIetaIeta<0.01)||(pho_isEE==1&&pho_sigmaIetaIeta<0.03))

*variables : cEP, phi_width/eta_width, lambdaRatio, eMax/eSCraw, eta_width, e2x2/e5x5, lamdba-/cEE, r19

cross check for the pt categories test

* signal dataset /GJet_Pt-20_doubleEMEnriched_TuneZ2_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

*background dataset /QCD_Pt-20to30_EMEnriched_TuneZ2_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

/QCD_Pt-30to80_EMEnriched_TuneZ2_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

/QCD_Pt-80to170_EMEnriched_TuneZ2_7TeV-pythia6/Spring11-PU_S1_START311_V1G1-v1/AODSIM

*statistic = same stat in the 2 samples

*I put some weights for training to take in account the differents X section of each sample

*what are the cuts = pho_et>20 && (abs(pho_eta)<=2.5)&&(!((abs(pho_eta)>1.4442)&&(abs(pho_eta)<1.566))) && pho_hoe<0.05 && pho_hasPixelSeed==0 && pho_IsoHollowTrkCone<2 && pho_IsoEcalRechit<4.2 && pho_IsoHcalRechit<2.2 && ((pho_isEB==1&&pho_sigmaIetaIeta<0.01)||(pho_isEE==1&&pho_sigmaIetaIeta<0.03))

*variables : cEP, phi_width/eta_width, lambdaRatio, eMax/eSCraw, eta_width, e2x2/e5x5, lamdba-/cEE, r19

*pt categories : 20to30, 30to40, 40to50, 50to06, pt>60 TIP one training with not cat, only cut aplied for effi-rej plot TIP

* 20110706_NN.ppt and the update based on the discussion 20110706_NN_update.ppt

* 20110817_NN_PT30.ppt with ET >30GeV. Bkg rejection is about 50%-60% with ET>30 GeV when keeping 90% signal efficiency

* 20110908HggIHEPIPNL.ppt: A simple presentation prepared in a short time.

* 20110918FakPhotonAndDataMCComparisonHgg.ppt: All photon ET>30GeV. With EGM-10-006 Loose selections, only 17 photons were selected except the 2 most isolated with combIso. With EGM-10-006 Loose selections except Sig_IeIe, 33 photons were selected except the 2 most isolated photons. The pileup reweighting and FastJet Rho correction were made for the Data/MC comparisons.

-- OlivierBondu - 20-May-2011

nw
Tag --- -------- Package -------- V06-22-00 CondFormats/DataRecord V00-03-01 CondFormats/EgammaObjects baseline20 and EtaWidth

Topic attachments
I Attachment History Action Size Date Who Comment
Unknown file formatpptx 0529_PhotonIDComparisons.pptx r1 manage 1101.2 K 2012-05-29 - 22:44 JunquanTao  
Unknown file formatpptx 0529_PhotonIDComparisons_Outputs.pptx r1 manage 243.3 K 2012-06-05 - 23:56 JunquanTao  
PDFpdf 201000615_NN_Higgs.pdf r1 manage 2181.7 K 2011-06-17 - 01:54 JunquanTao  
PowerPointppt 20110525_NN_Higgs.ppt r1 manage 7831.0 K 2011-05-29 - 11:57 JunquanTao  
PowerPointppt 20110601NNadditinalVar.ppt r1 manage 2182.0 K 2011-06-01 - 11:23 JunquanTao  
PowerPointppt 20110706_NN.ppt r1 manage 5494.5 K 2011-07-06 - 14:14 JunquanTao  
PowerPointppt 20110706_NN_update.ppt r1 manage 9039.0 K 2011-07-11 - 08:17 JunquanTao  
PowerPointppt 20110817_NN_PT30.ppt r1 manage 2397.0 K 2011-08-18 - 09:59 JunquanTao  
PowerPointppt 20110908HggIHEPIPNL.ppt r1 manage 1566.0 K 2011-09-08 - 15:24 JunquanTao  
PDFpdf 20110916.pdf r6 r5 r4 r3 r2 manage 852.2 K 2011-09-22 - 12:40 FanJiawei  
PowerPointppt 20110918FakPhotonAndDataMCComparisonHgg.ppt r3 r2 r1 manage 587.0 K 2011-09-19 - 12:09 JunquanTao  
PowerPointppt 20111014_HggEMG006andCiC.ppt r1 manage 552.0 K 2011-10-14 - 04:16 JunquanTao  
Unknown file formatrar 2011DataMCpdf.rar r1 manage 352.6 K 2012-06-05 - 09:39 JunquanTao  
PDFpdf 2011_DataMC_Photon_MITPre_EB_AfterCorrections.pdf r1 manage 365.7 K 2012-05-31 - 20:43 JunquanTao  
PDFpdf 2011_DataMC_Photon_MITPre_EE_AfterCorrections.pdf r1 manage 385.7 K 2012-05-31 - 20:43 JunquanTao  
PDFpdf 20120202_GlobalAnaWithNN.pdf r2 r1 manage 718.7 K 2012-02-03 - 18:09 JunquanTao  
PDFpdf 20120205GlobalAnaWithNNDebug.pdf r3 r2 r1 manage 588.7 K 2012-02-08 - 19:43 JunquanTao  
PDFpdf 20120308_NNTry.pdf r1 manage 421.4 K 2012-03-10 - 14:28 JunquanTao  
Unknown file formatpptx 20120308_NNTry.pptx r1 manage 1263.8 K 2012-03-22 - 10:49 JunquanTao  
Unknown file formatpptx 20120323_PhotonIDwithNNShowerShape.pptx r1 manage 2408.6 K 2012-03-25 - 13:49 JunquanTao  
PDFpdf 20120329_MVAinGlobeAndNN.pdf r1 manage 290.8 K 2012-03-29 - 10:00 JunquanTao  
Unknown file formatpptx 20120405_HggGlobePhotonIDMVA.pptx r1 manage 4406.0 K 2012-04-06 - 12:15 JunquanTao  
Unknown file formatpptx 20120420HggGlobePhotonIDMVA.pptx r1 manage 3332.1 K 2012-04-19 - 14:29 JunquanTao  
PDFpdf 20120423HggGlobePhotonIDMVA.pdf r1 manage 2763.1 K 2012-04-23 - 16:18 JunquanTao  
Unknown file formatpptx 20120423HggGlobePhotonIDMVA.pptx r2 r1 manage 3424.0 K 2012-04-22 - 09:25 JunquanTao  
PDFpdf 20120615SigIeIeCheckElectronPhoton.pdf r2 r1 manage 309.1 K 2012-06-15 - 18:06 JunquanTao  
PDFpdf 20120616SigIeIeEB.pdf r1 manage 107.7 K 2012-06-16 - 17:13 JunquanTao  
PDFpdf 20120907_HggSpinCosThetaStarJHU_JTao.pdf r2 r1 manage 180.7 K 2012-09-11 - 15:46 JunquanTao  
PDFpdf 20120918_2DFitComparisons.pdf r1 manage 366.8 K 2012-10-02 - 21:28 JunquanTao  
PDFpdf 20120918_2DandMVAComparisons.pdf r1 manage 372.7 K 2012-09-20 - 09:55 JunquanTao  
PDFpdf 20121206_2DFitStatus.pdf r1 manage 1058.1 K 2012-12-06 - 10:55 JunquanTao  
PDFpdf 2012FCPPL_JTao_IHEP_IPNL_CMS_.pdf r1 manage 3129.1 K 2012-03-16 - 11:17 JunquanTao  
Unknown file formatpptx 2012FCPPL_JTao_IHEP_IPNL_CMS_.pptx r3 r2 r1 manage 3002.9 K 2012-03-21 - 12:22 JunquanTao  
Unknown file formatpptx 2012FCPPL_JTao_IHEP_IPNL_CMS_Suzanne.pptx r1 manage 2719.4 K 2012-03-21 - 10:13 SuzanneGascon Version partially commented by Suzanne
PDFpdf 2012MCFactors_Zmumugamma.pdf r1 manage 484.8 K 2012-05-31 - 10:01 JunquanTao  
Unknown file formatpptx 2012MCFactors_Zmumugamma.pptx r2 r1 manage 739.5 K 2012-05-30 - 23:44 JunquanTao  
PDFpdf 2012MCFactors_Zmumugamma_V2.pdf r1 manage 514.5 K 2012-05-31 - 18:58 JunquanTao  
PDFpdf 2012MCFactors_Zmumugamma_V3.pdf r6 r5 r4 r3 r2 manage 1278.5 K 2012-06-04 - 10:44 JunquanTao  
PDFpdf 2012MCFactors_Zmumugamma_V4.pdf r1 manage 1216.2 K 2012-06-05 - 16:05 JunquanTao  
PDFpdf 2012MCFactors_Zmumugamma_V5_Linear.pdf r4 r3 r2 r1 manage 828.3 K 2012-06-08 - 11:52 JunquanTao  
PDFpdf 2012MCFactors_Zmumugamma_V6_Linear.pdf r2 r1 manage 722.5 K 2012-06-14 - 11:47 JunquanTao  
Unknown file formattop 2gNNLO_WRT2010Data_All.top r1 manage 23.8 K 2012-12-03 - 14:40 JunquanTao  
PDFpdf AN2008_063_v1Neutral_Pion_rejection_for_isolated_and_unconverted_photon_candidates_using_CMS_ECAL_and_Preshower_detector.pdf r1 manage 3809.8 K 2012-05-13 - 23:52 JunquanTao  
PDFpdf BkgReweight_FirtPTThenEta_Both_GJet2BinsXS.pdf r1 manage 30.9 K 2012-05-25 - 00:01 JunquanTao  
PDFpdf Comparisons_2gNNLO_ResBos_2010Data.pdf r1 manage 166.0 K 2012-12-03 - 15:35 JunquanTao  
Texttxt CosThetaStarBin_All.txt r1 manage 7.5 K 2012-04-24 - 21:20 JunquanTao  
Texttxt CosThetaStarBin_Weights_All.txt r1 manage 7.2 K 2012-04-24 - 21:21 JunquanTao  
Texttxt CosThetaStar_All.txt r1 manage 8.6 K 2012-04-24 - 21:20 JunquanTao  
Texttxt CosThetaStar_Weights_All.txt r1 manage 8.4 K 2012-04-24 - 21:20 JunquanTao  
PDFpdf DiPhoton2011_IHEP.pdf r1 manage 1673.3 K 2012-12-10 - 10:45 JunquanTao  
Microsoft Word filedoc DidierRAPIPNL20082009_1607Higgs.doc r2 r1 manage 1488.5 K 2011-09-21 - 11:43 SuzanneGascon New version: Pour Patrice: RA 2008-2009 avec Higgs + photons de Suzanne
PDFpdf EB_Eff_PT_EtaorPTPlusEtaReweighting.pdf r2 r1 manage 14.9 K 2012-05-24 - 00:44 JunquanTao  
PDFpdf EE_Eff_PT_EtaorPTPlusEtaReweighting.pdf r2 r1 manage 15.1 K 2012-05-24 - 00:45 JunquanTao  
GIFgif EE_input1.gif r1 manage 28.2 K 2012-05-24 - 13:47 JunquanTao  
GIFgif Eb_input1.gif r1 manage 28.4 K 2012-05-24 - 13:46 JunquanTao  
GIFgif EtaWidth_EcalISODR04RhoCorr_EB_Fake.gif r1 manage 23.1 K 2012-01-29 - 13:22 JunquanTao  
GIFgif EtaWidth_HcalISODR04Rho_EB_Fake.gif r1 manage 24.5 K 2012-01-29 - 13:23 JunquanTao  
GIFgif EtaWidth_RelCombinedISO_EB_Fake.gif r1 manage 21.0 K 2012-01-29 - 13:19 JunquanTao  
GIFgif EtaWidth_RelCombinedISO_EB_True.gif r1 manage 20.4 K 2012-01-29 - 13:19 JunquanTao  
GIFgif EtaWidth_RelTrackISO_EB_Fake.gif r1 manage 22.6 K 2012-01-29 - 13:24 JunquanTao  
GIFgif EtaWidth_TrackISODR04RhoCorr_EB_Fake.gif r1 manage 20.9 K 2012-01-29 - 13:22 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_CB_BasicSigVetoEcalISO.pdf r1 manage 37.7 K 2012-01-06 - 22:50 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_CE_BasicSigVetoEcalISO.pdf r1 manage 41.0 K 2012-01-06 - 22:50 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_EB_Converted_NoISO.pdf r1 manage 124.8 K 2011-10-09 - 15:54 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_EB_Converted_OnlySigIeIe.pdf r1 manage 104.1 K 2011-10-09 - 15:55 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_EB_Unconverted_NoISO.pdf r1 manage 123.3 K 2011-10-09 - 15:56 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_EB_Unconverted_OnlySigIeIe.pdf r1 manage 102.4 K 2011-10-09 - 15:56 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_EE_Converted_NoISO.pdf r1 manage 116.8 K 2011-10-09 - 15:57 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_EE_Converted_OnlySigIeIe.pdf r1 manage 103.8 K 2011-10-09 - 15:57 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_EE_Unconverted_NoISO.pdf r1 manage 119.5 K 2011-10-09 - 15:57 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_EE_Unconverted_OnlySigIeIe.pdf r1 manage 105.0 K 2011-10-09 - 15:58 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_FailAnyCiCisoEB.pdf r1 manage 36.2 K 2012-01-29 - 12:53 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_UB_BasicSigVetoEcalISO.pdf r1 manage 38.3 K 2012-01-06 - 22:50 JunquanTao  
PDFpdf FakePhotonDistributionsDataMC_UE_BasicSigVetoEcalISO.pdf r1 manage 40.5 K 2012-01-06 - 22:50 JunquanTao  
PDFpdf HGammaSlidesDraft.pdf r4 r3 r2 r1 manage 5942.1 K 2011-09-26 - 13:30 HuguesBrun  
PDFpdf Hgg2Dfit_Fabrice.pdf r1 manage 3133.6 K 2012-10-31 - 02:52 JunquanTao  
PDFpdf HggInterferenceCosThetaStar.pdf r2 r1 manage 292.6 K 2012-04-24 - 21:27 JunquanTao  
PDFpdf HggInterferenceCosThetaStarOrThetaStar.pdf r2 r1 manage 1004.5 K 2012-04-26 - 16:18 JunquanTao  
PDFpdf HggInterferenceWithThetaStar.pdf r2 r1 manage 924.9 K 2012-04-26 - 22:13 JunquanTao  
PDFpdf HggInterferenceWithThetaStar_0428.pdf r3 r2 r1 manage 1187.3 K 2012-04-29 - 23:11 JunquanTao  
PDFpdf HggInterferenceWithThetaStar_Jan16.pdf r6 r5 r4 r3 r2 manage 1318.2 K 2012-05-02 - 22:44 JunquanTao  
PDFpdf Hugues6ofJuly.pdf r1 manage 849.7 K 2011-07-06 - 15:04 HuguesBrun  
PDFpdf HuguesNNhiggsAnalysis18juillet11.pdf r1 manage 13092.9 K 2011-08-04 - 00:33 SuzanneGascon Hugues' presentation to the IPNL CMS group
PDFpdf HuguescutOptiPhotonsNew9oct11.pdf r1 manage 8854.5 K 2011-10-20 - 00:00 SuzanneGascon Hugues' results,9oct11
PDFpdf HuguescutOptiPhotons_16oct11.pdf r1 manage 8745.0 K 2011-10-20 - 00:03 SuzanneGascon Hugues' results,16oct11
PDFpdf HugueslesSlidesEnPlus1aout11.pdf r1 manage 2614.9 K 2011-08-04 - 00:35 SuzanneGascon Extra slides from Hugues, August 1
Unknown file formatpptx MVA12_Eff_PTEtaNvtx.pptx r1 manage 263.5 K 2012-05-31 - 16:31 JunquanTao  
Texttxt Mass100160_Angle1090_1Degree.txt r1 manage 33.2 K 2012-04-28 - 23:03 JunquanTao  
PNGpng Mass_2012DataMC_MVA_photonID005_V7.png r1 manage 21.3 K 2012-06-21 - 19:24 JunquanTao  
PNGpng Mass_2012DataMC_MVA_photonIDm02_V7.png r1 manage 20.4 K 2012-06-21 - 19:29 JunquanTao  
PNGpng Mass_2012DataMC_OldCiC4_V7.png r1 manage 23.3 K 2012-06-21 - 19:22 JunquanTao  
PDFpdf NNshapeTests.pdf r3 r2 r1 manage 5701.4 K 2011-06-14 - 09:30 HuguesBrun  
PDFpdf NNworkingslides.pdf r1 manage 12833.3 K 2011-05-22 - 16:50 HuguesBrun  
Unknown file formatcxx PhotonCandDis.cxx r1 manage 62.4 K 2011-05-25 - 05:55 JunquanTao  
PDFpdf PhotonID52X_0511.pdf r1 manage 1103.7 K 2012-05-11 - 15:00 JunquanTao  
PDFpdf PhotonID52X_0514.pdf r1 manage 1218.2 K 2012-05-15 - 11:37 JunquanTao  
Unknown file formatpptx PhotonID52X_0521.pptx r7 r6 r5 r4 r3 manage 2825.0 K 2012-05-29 - 22:20 JunquanTao  
Unknown file formatpptx PhotonID52X_0523.pptx r1 manage 1011.5 K 2012-05-23 - 11:10 JunquanTao  
PDFpdf PhotonID52X_0525Egamma.pdf r2 r1 manage 1609.1 K 2012-05-25 - 12:14 JunquanTao  
Unknown file formatpptx PhotonID52X_0525Egamma.pptx r1 manage 2171.0 K 2012-05-29 - 22:39 JunquanTao  
PDFpdf PhotonIDMVAGlobe20110612.pdf r1 manage 489.9 K 2012-06-12 - 16:22 JunquanTao  
PowerPointppt SuzanneEtAlNNForHgg26Sept11.ppt r1 manage 16072.0 K 2011-09-26 - 10:13 SuzanneGascon Draft of Hgg NN presentation for 26 sep 11
PowerPointppt SuzanneEtAlNNForHgg7nov11.ppt r1 manage 18933.0 K 2011-11-07 - 08:31 SuzanneGascon Draft of Hgg NN presentation for 7 November 2011
PowerPointppt SuzanneHuguesFanOlivieretAlPhotonCommissioning2012EGM27avr12.ppt r2 r1 manage 1784.0 K 2012-04-27 - 14:14 HuguesBrun  
PowerPointppt SuzanneHuguesFanOlivieretAlPhotonCommissioning2012EGM27avr12Final_4mai27aprJSON.ppt r1 manage 2431.0 K 2012-05-05 - 23:43 SuzanneGascon 2012 Data-MC and 2012-2011 Data-data comparisons for cluster-shape variables
Unknown file formatpptx SuzanneIHEP-IPNL_CMS_Nanjing27Mars13New.pptx r2 r1 manage 12718.1 K 2013-03-27 - 04:23 SuzanneGascon Cleaned version of FCPPL talk
Unknown file formatpptx SuzanneTao20120323_PhotonIDwithNNShowerShape.pptx r1 manage 2370.8 K 2012-03-26 - 17:00 SuzanneGascon Draft of Hgg talk for today (Tao's slides with style etc modifications)
PowerPointppt SuzanneTaoFanClusterShapesROCs5may12.ppt r1 manage 336.5 K 2012-05-09 - 03:29 SuzanneGascon Roc plots in 4_2_x for cluster shape and detector Isolation variables
C source code filec TMVAClassification.C r1 manage 28.0 K 2011-05-24 - 16:43 HuguesBrun  
PDFpdf TaoDiPhotonDataMC_mH115_NN.pdf r1 manage 196.2 K 2011-05-29 - 12:03 JunquanTao  
PDFpdf TaoDiPhotonDataMC_mH120_NN.pdf r1 manage 193.7 K 2011-05-29 - 12:03 JunquanTao  
PowerPointppt TaoSlides.ppt r1 manage 3482.0 K 2011-05-22 - 16:55 HuguesBrun  
C source code filec TaoTMVAClassification.C r1 manage 26.4 K 2011-05-25 - 05:56 JunquanTao  
PDFpdf Tao_MCFactors_Zmumugamma_Linear.pdf r2 r1 manage 787.7 K 2013-02-08 - 08:55 JunquanTao  
PDFpdf Tao_MCFactors_Zmumugamma_Pol1Plo2NewZee.pdf r1 manage 1077.2 K 2013-02-21 - 11:13 JunquanTao  
PDFpdf Taojq_PhotonDiPhoton_ShowerShapeTemplate.pdf r1 manage 862.6 K 2012-06-01 - 13:38 JunquanTao  
Texttxt ThetaStar_125.txt r1 manage 4.6 K 2012-04-26 - 11:52 JunquanTao  
Texttxt ThetaStar_Delta_1DegreeBin.txt r4 r3 r2 r1 manage 25.7 K 2012-04-27 - 21:28 JunquanTao  
Texttxt ThetaStar_Weights_1DegreeBin.txt r4 r3 r2 r1 manage 25.4 K 2012-04-27 - 21:28 JunquanTao  
Unknown file formatlog data_scp.log r1 manage 11.4 K 2012-01-30 - 23:24 JunquanTao  
Unknown file formatlog data_stat.log r1 manage 239.5 K 2012-01-30 - 23:04 JunquanTao  
C source code filec do_NLLFit.C r1 manage 20.7 K 2012-05-28 - 03:39 FanJiawei Codes for correction factors ---FAN
PDFpdf lastNNslides.pdf r1 manage 14013.8 K 2011-06-14 - 09:27 HuguesBrun  
GIFgif limit_sm_Asymptotic.gif r1 manage 42.8 K 2012-01-20 - 09:41 JunquanTao  
GIFgif limit_sm_Asymptotic_ratio.gif r1 manage 40.3 K 2012-01-20 - 09:48 JunquanTao  
GIFgif limit_sm_Asymptotic_ratio_MVA.gif r1 manage 39.7 K 2012-02-23 - 17:54 JunquanTao  
GIFgif limit_sm_Asymptotic_ratio_smooth.gif r3 r2 r1 manage 35.1 K 2012-01-20 - 09:58 JunquanTao  
GIFgif limit_sm_Asymptotic_smooth.gif r1 manage 39.7 K 2012-01-20 - 09:42 JunquanTao  
Unknown file formatlog mc_scp.log r2 r1 manage 37.9 K 2012-01-31 - 08:29 JunquanTao  
Unknown file formatlog mc_stat.log r2 r1 manage 1193.2 K 2012-01-31 - 08:29 JunquanTao  
PDFpdf moreDATAMCcomp.pdf r1 manage 6160.7 K 2011-06-15 - 15:25 HuguesBrun  
Unknown file formatlog passingEvents.log r2 r1 manage 2678.1 K 2011-09-21 - 22:37 HuguesBrun  
PNGpng sieie_EB_Corr1_DavidTriedZeeCorr.png r1 manage 16.6 K 2012-06-15 - 00:01 JunquanTao  
PNGpng sieie_EB_Corr1_TaoSuggestZee.png r1 manage 16.3 K 2012-06-15 - 00:02 JunquanTao  
PNGpng sieie_EB_Corr1_ZeeCorrFromMatteo.png r1 manage 16.8 K 2012-06-15 - 00:24 JunquanTao  
PDFpdf slidesHugues15ofJune.pdf r1 manage 15817.3 K 2011-06-15 - 15:25 HuguesBrun  
Edit | Attach | Watch | Print version | History: r172 < r171 < r170 < r169 < r168 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r172 - 2013-11-20 - JunquanTao
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

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