ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/event_count.C
Revision: 1.3
Committed: Thu Mar 5 17:30:21 2009 UTC (16 years, 2 months ago) by kukartse
Content type: text/plain
Branch: MAIN
CVS Tags: V00-03-01, ZMorph_BASE_20100408, gak040610_morphing, V00-02-02, gak011410, gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609, V00-01-15, V00-01-14, V00-01-13, V00-01-12, V00-01-11, V00-01-10, gak031009, gak030509, HEAD
Branch point for: ZMorph-V00-03-01, CMSSW_22X_branch
Changes since 1.2: +9 -3 lines
Error occurred while calculating annotation data.
Log Message: 
*** empty log message ***

File Contents

# Content
1 {
2 TChain chain("ttljets");
3 //chain . Add("/uscms_data/d1/lpcljm/MVA/Summer08/PAT_Layer1/InclusiveMu5Pt50_Summer08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/InclusiveMu5Pt50_Summer08_IDEAL_V9_v1_GEN-SIM-RECO-secondary-22dec2008.root");
4 chain . Add("/uscms/home/kukarzev/work/disk_ljmet/MVA/Summer08/Jan_23_2009/PAT_Layer1/InclusiveMuPt15_Summer08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_InclusiveMuPt15-12feb2009.root");
5 //chain . Add("/uscms_data/d1/lpcljm/MVA/Summer08/Jan_23_2009/PAT_Layer1/TTJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_ttjets-06feb2009.root");
6 //chain . Add("/uscms_data/d1/lpcljm/MVA/Summer08/Jan_23_2009/PAT_Layer1/WJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_wjets-06feb2009.root");
7 //chain . Add("/uscms_data/d1/lpcljm/MVA/Summer08/Jan_23_2009/PAT_Layer1/ZJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_zjets-06feb2009.root");
8
9 //gROOT->Macro("./cuts.C");
10
11 // selection
12 #include "cuts.C"
13 TCut nom_cut = semilep_cut && reversible_cut;
14 TCut rev_cut = semilep_cut && !reversible_cut;
15 the_cut = nom_cut;
16 the_cut = semilep_cut && lepton_jet_dr05_cut && !reversible_cut;;
17
18 TH1F h("h","Event weight",1000000,0, 1000);
19 TH1F h_weight("h_weight","Event weight",1000000,0, 1000);
20
21 //TCut the_cut = "lepton_calo_iso<7.0 && muon_isGood_AllGlobalMuons==1";
22 //TCut the_cut = "(1.0/(1.0+(lepton_calo_iso+lepton_track_iso)/lepton_pt))>0.98 && muon_isGood_AllGlobalMuons==1";
23 //TCut the_cut = "n_tagged_jets_trackCounting_loose>1";
24 //TCut the_cut = "";
25 //TCut the_cut = nominal_cut;
26
27 // weighted count
28 //chain . Draw("event_weight>>h_weight", weight*the_cut);
29 //unweighted count
30 chain . Draw("event_weight>>h", the_cut);
31
32 //cout << "Weighted events: " << h_weight . Integral() << endl;
33 cout << "Not weighted events: " << h . Integral() << endl;
34 }