export CVSROOT=:ext:pjacomeh@lxplus.cern.ch:/afs/cern.ch/user/c/cvscmssw/public/CMSSW cvs co -r V00-02-06 StatisticalTools/RooStatsRoutines cd StatisticalTools/RooStatsRoutines
gSystem -> SetIncludePath( "-I$ROOFITSYS/include" );
cd StatisticalTools/RooStatsRoutines/root root -l .L roostats_cl95.C+ LimitResult mylimit = GetClsLimit(17775.6, 711.0, 0.5, 0.1, 1887.1, 10.0, 1991.0); Double_t obs_limit = mylimit.GetObservedLimit(); cout << obs_limit << endl;
SAMPLE | PreSel | HLT_Mu40_ eta2p1_v TRIGGER | PV≥1 | Nmuon≥1 PT≥26GeV | Nmuon=1 PT≥26GeV | Nelect=0 | Njet≥6 PT≥40GeV | b-tags≥2 CSVM | HT≥400GeV |
---|---|---|---|---|---|---|---|---|---|
DATA |
9000.00 | 6610.00 | 6610.00 | 6610.00 | 6452.00 | 6452.00 | 6086.00 | 2478.00 | 2319.00 |
TTbar |
7717.84 | 5358.03 | 5358.03 | 5358.03 | 5314.84 | 5314.84 | 5043.89 | 2426.53 | 2305.05 |
WJets |
842.13 | 552.65 | 552.65 | 552.65 | 552.65 | 552.65 | 517.56 | 17.54 | 17.54 |
DYJets_50toInf |
229.15 | 180.05 | 180.05 | 180.05 | 128.90 | 128.90 | 112.53 | 8.18 | 8.18 |
SingleTop _T_sch |
1.82 | 1.04 | 1.04 | 1.04 | 1.04 | 1.04 | 1.04 | 0.52 | 0.52 |
SingleTop_T_tWch |
84.22 | 61.57 | 61.57 | 61.57 | 60.78 | 60.78 | 57.20 | 24.63 | 23.83 |
SingleTop _T_tch |
16.29 | 9.08 | 9.08 | 9.08 | 9.08 | 9.08 | 7.48 | 4.01 | 3.74 |
SingleTop _Tbar_tch |
6.49 | 3.11 | 3.11 | 3.11 | 3.11 | 3.11 | 2.82 | 1.41 | 1.41 |
SingleTop_Tbar_tWch |
77.71 | 54.88 | 54.88 | 54.88 | 54.88 | 54.88 | 52.48 | 24.44 | 22.03 |
SingleTop_Tbar_sch |
1.34 | 0.67 | 0.67 | 0.67 | 0.67 | 0.67 | 0.67 | 0.22 | 0.22 |
TOTAL MC | 8976.99 | 6221.07 | 6221.07 | 6221.07 | 6125.93 | 6125.93 | 5795.67 | 2507.48 | 2382.54 |
cd /uscms_data/d3/pjacome/sframe cd ZprimeAnalysis nano src/ZprimeSelectionCycle.cxx
cd ZprimeAnalysis rm obj/ZprimeSelectionCycle.o obj/dep/ZprimeSelectionCycle.d rm ../lib/ZprimeAnalysis.par ../lib/libZprimeAnalysis.so make
wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/multiConfigs_SELECTION_NO-PROOF.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/sframe_CONDOR_NO-PROOF.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/mySelectionCycle_config.xml wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/template.condor bash multiConfigs_NO-PROOF.sh
SFrame_new/SFramePlotter$ ./bin/Plots -f myExampleCyclePlots.steer
Sample | Events | Luminosity[pb-1] |
---|---|---|
MuDataA | 7956466 | 889 |
MuDataB | 43087517 | 4429 |
MuDataC | 63804279 | 7152 |
MuDataD | 66607606 | 7318 |
TOTAL DATA | 181455868 | 19788 |
Sample | TOTAL | SELECTED | Xsec[pb] | +Error[pb] | -Error[pb] | order | Luminosity[pb-1] | ||
---|---|---|---|---|---|---|---|---|---|
TTTT_GH400_1muo | 10000 | 10000 | 0.00001439 | 4.066E-008 | 4.066E-008 | 694927032.66157 | |||
TTTT_GH400 | 50000 | 50000 | 0.0001758 | 2.945E-007 | 2.945E-007 | 284414106.939704 | |||
TTTT_GH500_1muo | 10000 | 10000 | 0.00004151 | 1.185E-007 | 1.185E-007 | 240905805.82992 | |||
TTTT_GH500 | 50000 | 50000 | 0.0005019 | 6.669E-007 | 6.669E-007 | 99621438.5335724 | |||
TTTT_GH600_1muo | 10000 | 10000 | 0.000033057 | 302507789.575582 | |||||
TTTT_GH600 | 50000 | 50000 | 0.0003966 | 0.000000557 | 0.000000557 | 126071608.673727 | |||
TTTT_GH700_1muo | 10000 | 10000 | 0.00001928 | 5.244E-008 | 5.244E-008 | 518672199.170124 | |||
TTTT_GH700 | 50000 | 50000 | 0.0002323 | 3.568E-007 | 3.568E-007 | 215238915.195867 | |||
TTTT_GH800_1muo | 10000 | 10000 | 0.000010225 | 977995110.02445 | |||||
TTTT_GH800 | 50000 | 50000 | 0.0001227 | 1.958E-007 | 1.958E-007 | 407497962.510187 | |||
TTTT_GH900_1muo | 10000 | 10000 | 0.000005136 | 1.438E-008 | 1.438E-008 | 1947040498.44237 | |||
TTTT_GH900 | 50000 | 50000 | 0.00006222 | 0.000000086 | 0.000000086 | 803600128.576021 | |||
TTTT_GH1000_1muo | 10000 | 10000 | 0.000002558 | 7.388E-009 | 7.388E-009 | 3909304143.86239 | |||
TTTT_GH1000 | 50000 | 50000 | 0.00003076 | 3.925E-008 | 3.925E-008 | 1625487646.29389 | |||
TTTT_SM_1muo | 10000 | 10000 | 0.00003946 | 1.069E-007 | 1.069E-007 | 253421186.011151 | |||
TTTT_SM | 50000 | 50000 | 0.0007136 | 0.000002498 | 0.000002498 | 70067264.573991 | |||
TTbar | 21675970 | 21560109 | 245.8 | 0.025 | -0.034 | 2.6 | -2.6 | NNLO | 88185.3946297803 |
SingleTop_T_tch | 3758227 | 3748155 | 56.4 | 2.1 | 0.3 | 1.1 | -1.1 | NNLO | 66635.2304964539 |
SingleTop_T_tWch | 497658 | 495559 | 11.1 | 0.3 | 0.3 | 0.7 | -0.7 | NNLO | 44834.0540540541 |
SingleTop_T_sch | 259961 | 259176 | 3.79 | 0.07 | 0.07 | 0.13 | -0.13 | NNLO | 68591.2928759894 |
SingleTop_Tbar_tch | 1935072 | 1930185 | 30.7 | 0.7 | 0.7 | 0.9 | -1.1 | NNLO | 63031.661237785 |
SingleTop_Tbar_tWch | 493460 | 491463 | 11.1 | 0.3 | 0.3 | 0.7 | -0.7 | NNLO | 44455.8558558559 |
SingleTop_Tbar_sch | 139974 | 139604 | 1.76 | 0.01 | 0.01 | 0.08 | 0.08 | NNLO | 79530.6818181818 |
DYJetsToLL_M-50 | 30459503 | 30404232 | 1177.3 | 5.9 | -3.6 | 38.8 | 38.8 | NNLO | 25872.3375520258 |
DYJetsToLL_M-50 | 30459503 | 30404232 | 3531.9 | 17.7 | -10.8 | 116.4 | 116.4 | NNLO | 8624.112517342 |
DY_10to50 | 37835275 | 37803439 | |||||||
WJetsToLNu_v1 | 18393090 | 18375310 | 12234.4 | 79 | -39.7 | 414.7 | NNLO | 1503.3912574381 | |
WJetsToLNu_v2 | 57709905 | 57653686 | 12234.4 | 79 | -39.7 | 414.7 | NNLO | 4717.0196331655 | |
WJetsToLNu_v1 + WJetsToLNu _2 | 76102995 | 76028996 | 36703.2 | 237 | -119.1 | 1244.1 | NNLO | 2073.4702968679 | |
WW | 10000431 | 9982151 | 56 | 2.3 | -2.3 | 2 | -2 | NLO | 178579.125 |
WZ | 10000283 | 9981352 | 33.21 | NLO | 301122.643781993 | ||||
ZZ | 9799908 | 9780255 | 8.2 | 0.3 | -0.3 | 0.3 | -0.3 | NLO | 1195110.73170732 |
wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/multiConfigs_NO-PROOF.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/sframe_CONDOR_NO-PROOF.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/myPreSelection_Muo_NO-PROOF.xml wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/template.condor bash multiConfigs_NO-PROOF.sh
Proof.Sandbox: /uscms_data/d3/pjacome/proof
<Cycle Name="FourtopsSelectionCycle" RunMode="PROOF" ProofServer="lite://" ProofWorkDir="/uscms_data/d3/pjacome/proof/" OutputDirectory="/eos/uscms/store/user/pjacome/Selection/from_Zprime/Muons/" PostFix="" TargetLumi="1.0" >
cd ~ ln -s /eos/uscms/store/user/pjacome/rootrc .rootrc
cd /eos/uscms/store/user/pjacome/preselection/muon/run1 cp /uscms_data/d3/pjacome/sframe/ZprimeAnalysis/config/JobConfig.dtd . wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/myPreSelection_Muo_PROOF.xml wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/multiConfigs_PROOF.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/template.condor wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/allEOS_no-flag_condorPROOF_SFrame.sh bash multiConfigs_PROOF.sh
source fullsetup.sh cd ZprimeAnalysis rm obj/ZprimePreSelectionCycle.o obj/dep/ZprimePreSelectionCycle.d rm ../lib/ZprimeAnalysis.par ../lib/libZprimeAnalysis.so
Selection* preselection = new Selection("preselection"); if(m_Electron_Or_Muon_Selection=="Electrons" || m_Electron_Or_Muon_Selection=="Electron" || m_Electron_Or_Muon_Selection=="Ele" || m_Electron_Or_Muon_Selection=="ELE") { preselection->addSelectionModule(new NElectronSelection(1,int_infinity(),35,2.5));//at least one electron preselection->addSelectionModule(new NMuonSelection(0,0));//no muons } else if(m_Electron_Or_Muon_Selection=="Muon" || m_Electron_Or_Muon_Selection=="Muons" || m_Electron_Or_Muon_Selection=="Mu" || m_Electron_Or_Muon_Selection=="MU") { preselection->addSelectionModule(new NElectronSelection(0,0));//no electron preselection->addSelectionModule(new NMuonSelection(1,int_infinity(),26,2.1));//at least one muon } else { m_logger << ERROR << "Electron_Or_Muon_Selection is not defined in your xml config file --- should be either `ELE` or `MU`" << SLogger::endmsg; } preselection->addSelectionModule(new NJetSelection(6,int_infinity(),40,2.4));//at least two jets RegisterSelection(preselection); ## . . . //clean collections here if(bcc->muons) cleaner.MuonCleaner(26,2.1,0.12); // MuonCleaner(double ptmin, double etamax, double relisomax) if(bcc->electrons) cleaner.ElectronCleaner(35,2.5,0.1,m_reversed_electron_selection,false,true); // ElectronCleaner(double ptmin, double etamax, double relisomax, bool reverseID, bool reverseIso, bool egmId) if(bcc->jets) cleaner.JetLeptonSubtractor(m_corrector,false); if(!bcc->isRealData && bcc->jets) cleaner.JetEnergyResolutionShifter(); if(bcc->jets) cleaner.JetCleaner(40,2.4,true); // JetCleaner(double ptmin, double etamax, bool doPFID)
make
<Item Name="JECDataGlobalTag" Value="FT_53_V21_AN6" /> <Item Name="JECMCGlobalTag" Value="START53_V27" /> <Item Name="JECJetCollection" Value="AK5PFchs" />
cd /eos/uscms/store/user/pjacome/ntuples wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/D_Ntuples.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/template.condor wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/cmsrun_condor.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/myNtupleWriter_cfg.py.txt mv myNtupleWriter_cfg.py.txt myNtupleWriter_cfg.py nano D_Ntuples.sh ./D_Ntuples.sh
cd /eos/uscms/store/user/pjacome/count_ntuple grep -v "_x_OK" list_*.xml > bad_files
cd /uscms_data/d3/pjacome/sframe/SFrameAnalysis/config/Samples_TTBSM53_FNAL cp /eos/uscms/store/user/pjacome/count_ntuple/bad_files . wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/removing_bad_files.sh bash removing_bad_files.sh
cp /uscms_data/d3/pjacome/CMSSW_5_3_3/src/TopQuarkAnalysis/TopPairBSM/test/ttbsm_cfg.py my_pat.py ## USAGE: cmsRun my_pat.py useData=0 runOnFastSim=1 usePythia8=1 2>err 1> log
cd /eos/uscms/store/user/pjacome/pattuples wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/template.condor wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/C_PATtuples.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/cmsrun_condor.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/my_pat.py.txt mv my_pat.py.txt my_pat.py nano C_PATtuples.sh ./C_PATtuples.sh
./generate_GEN-FASTSIM-HLT_config_file.sh
cd /eos/uscms/store/user/pjacome/fastsim wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/template.condor wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/B2_GEN_FASTSIM_HLT.sh wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/cmsrun_condor.sh ./B2_GEN_FASTSIM_HLT.sh
edmDumpEventContent /eos/uscms/store/user/pjacome/fastsim/SM_tttt_8TeV_10000ev_r01.root
cd /eos/uscms/store/user/pjacome/mg5 wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/template.condor wget --no-check-certificate https://twiki.cern.ch/twiki/pub/Sandbox/FourTopAnalysisPabloJacomeV2/A2_GH_Madgraph.sh cond_f='gh_mg' cp template.condor ${cond_f}.condor sed -i 's|__script__|'${PWD}'/'A2_GH_Madgraph.sh'|g' ${cond_f}.condor sed -i 's|__name__|'${cond_f}'|g' ${cond_f}.condor sed -i 's|__config__| |g' ${cond_f}.condor condor_submit ${cond_f}.condor
SM ttbarttbar | 20000 | 7.132E-4 | |
GH600 ttbarttbar | 20000 | 3.959E-4 | 4.549E-7 |
GH800 ttbarttbar | 20000 | 1.227E-4 | 1.366E-7 |
GH1000 ttbarttbar | 20000 | 3.073E-5 | 3.732E-8 |
cd /uscms_data/d3/pjacome/MadGraph5_v1_5_9/madanalysis5 ./bin/ma5
cd ExRootAnalysis ExRootLHEFConverter my.LHE out.root
condor_q <usr> condor_q -analyze <usr> condor_rm _proc_id_
_file0->ls(); TDirectory* Event_Postsel = _file0->GetDirectory("Event_Postsel;1"); Event_Postsel->ls(); TH1F * HT = (TH1F*)Event_Postsel->Get("HT;1"); # HT->Scale(0.1); ## example HT->Draw(); TTree *myTree = (TTree*)_file0->Get("AnalysisTree;1");