(r64) LogbookFourTopAnalysisPabloJacome < Sandbox

---++!! Pablo Jacome's Four Top Logbook

%TOC%

%TWISTY{<br /> mode="div"<br />  showimgleft="%ICONURLPATH{toggleopen}%"<br />  hideimgleft="%ICONURLPATH{toggleclose}%"<br />  }% 
%ENDTWISTY%

---+++ References

* <a target="_blank" href="http://arxiv.org/abs/1001.0221">http://arxiv.org/abs/1001.0221</a><br /> * &ldquo;Four-heavy-quark hadroproduction&rdquo;, V. Barger, A. L. Stange, and R. J. N. Phillips, Phys. Rev. D 44 (Oct, 1991) 1987&ndash;1996, doi:10.1103/PhysRevD.44.1987.

---+++ P) 22-Jul-2013 to 28-Jul-2013
   1 Playing with ROOT, trying to plot PrimaryVertices histogram: 
      1 <verbatim>TFile* file = new TFile("FourtopsPreSelectionCycle.DATA.DATA.root","READ");
TTree* mytree = (TTree*)(file -> Get("AnalysisTree"));
mytree->Draw("goodOfflinePrimaryVertices.m_nTracks");
</verbatim>
   1 Superimposed and adding histograms: 
      1 <verbatim>TCanvas *c1 = new TCanvas("c1","hists with different scales",600,400);

TFile* f1 = new TFile("FourtopsPreSelectionCycle.DATA.DATA.root","READ");
TTree* t1 = (TTree*)(f1 -> Get("AnalysisTree"));

TFile* f2 = new TFile("FourtopsPreSelectionCycle.MC.TTTT.root","READ");
TTree* t2 = (TTree*)(f2 -> Get("AnalysisTree"));


t1->Draw("goodOfflinePrimaryVertices.m_nTracks>>h1");
t2->Draw("goodOfflinePrimaryVertices.m_nTracks>>h2");
TH1F *h1 = (TH1F*)gDirectory->Get("h1");
TH1F *h2 = (TH1F*)gDirectory->Get("h2");
TH1F *h3 = (TH1F*)gDirectory->Get("h1");
h1->SetLineColor(kRed);
h3->Add(h2);
h3->Draw();

h1->Draw();
c1->Update();
h2->Draw("same");
</verbatim>

---+++ O) 15-Jul-2013 to 21-Jul-2013
   1 *QUESTIONS:* 
      1 
   1 Reading about PROOF: 
      1 http://indico.cern.ch/getFile.py/access?contribId=13&resId=0&materialId=slides&confId=71202 Write-up on SFrame: http://sourceforge.net/apps/mediawiki/sframe/ Configuring a small PROOF cluster: http://sourceforge.net/apps/mediawiki/sframe/index.php?title=SFrame-PROOF#Example_setup_of_a_small_PROOF_cluster
      1 Add the following in the config file, inside <cycle> tag definition: 
         1 <verbatim>  RunMode="PROOF" ProofServer="lite://"  ProofWorkDir="__CONDOR__"
RunMode="PROOF" ProofServer="lite"  ProofWorkDir=""
</verbatim>
      1 Files for sending CONDOR jobs. *sframe_condor_exe.csh*: 
         1 <verbatim>#! /bin/csh
#

# Modify SFrame config
cd ${_CONDOR_SCRATCH_DIR}
cp -r /uscms_data/d3/pjacome/_TEST2/CMSSW_5_3_3/src/SFrame/FITAnalysis/config .
cp config/FourtopsSelectionCycleMuons_config.xml ./temp.xml
sed -i s:__CONDOR__:${_CONDOR_SCRATCH_DIR}:g temp.xml
cp temp.xml config/.
cd config

# Source version of CMSSW associated to SFRAME
source /uscmst1/prod/sw/cms/cshrc uaf
cd /uscms_data/d3/pjacome/_TEST2/CMSSW_5_3_3/src
eval `scramv1 runtime -csh`

# Setting SFrame
cd /uscms_data/d3/pjacome/_TEST2/CMSSW_5_3_3/src/SFrame
source fullsetupSFrame.csh

# Run the sframe job
cd ${_CONDOR_SCRATCH_DIR}/config
sframe_main temp.xml
cd ${_CONDOR_SCRATCH_DIR}
rm -rf config
</verbatim>
      1 *sframe_condor_proof_cfg.txt* 
         1 <verbatim>universe = vanilla
Executable = sframe_condor_exe.csh
Requirements = Memory >= 199 &&OpSys == "LINUX"&& Machine == "cmswn1019.fnal.gov"
REQUEST_CPUS = 8
+BigMemoryJob = TRUE
Output = sframe_$(Cluster)_$(Process).stdout
Error = sframe_$(Cluster)_$(Process).stderr
Log = sframe_$(Cluster)_$(Process).log
notify_user = ${LOGNAME}@FNAL.GOV
Arguments =
Queue 1


</verbatim>
      1 Send job: <verbatim>condor_submit  sframe_condor_proof_cfg.txt</verbatim>
      1 State of the job: <verbatim>condor_q -bet <JOB_ID> </verbatim>
      1 Killing 'idle' jobs that have been more than 24 hours: <verbatim>condor_rm <JOB_ID> </verbatim>
      1 Changing settings of condor job: 
         1 <verbatim>Requirements = Memory >= 199 && OpSys == "LINUX"&& (Arch != "DUMMY" )&& Disk > 1000000</verbatim>
   1 Playing with ROOT: 
      1 <verbatim>AnalysisTree->Scan("event")</verbatim>
   1 Learning about Histos: 
      1 http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/algomez/python/plothistos.py?view=log 
         1 <verbatim>python plotHistos.py -i /eos/uscms/store/user/yumiceva/fourtops/sframePreselection/FourtopsPreSelectionCycle.MC.TTTT.root  -o output.root -p 6</verbatim>
         1 Didn't work. It didn't throw output.root
      1 http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/Yumiceva/TreeAnalyzer/test/Analyzer.C?view=log
      1 *https://twiki.cern.ch/twiki/bin/view/Main/FYumicevaTopFwk#Analysis_of_ntuples*

---+++ N) 08-Jul-2013 to 14-Jul-2013
   1 *QUESTIONS:* 
      1 Do I have to run the PreSelection for every type of filter? (... one for electrons, then one for muons and finally one for jets) 
         1 ... muons and electrons seems to be exclusive
      1 What is BaseCycleContainer?
      1 Why Selection doesn't take input data drom PreSelection?
   1 Exploring *src/FourtopsPreSelectionCycle.cxx* and *src/FourtopsSelectionCycle.cxx* 
      1 To compare: /uscms/home/yumiceva/work/sframe/CMSSW_5_3_3/src/SFrame/FITAnalysis/src
   1 Preselection *FourtopsPreSelectionCycle.cxx* : 
      1 "At least one good muon" 
         1 Changed in the function: FourtopsPreSelectionCycle::BeginInputData
         1 preselection-&gt;addSelectionModule(new NElectronSelection(0,0));//no electron preselection-&gt;addSelectionModule(new NMuonSelection(1,int_infinity()));
      1 "min pt, max eta , min reliso, and muon ID " for muon, electrons and jets. 
         1 Changed in the function: FourtopsPreSelectionCycle::ExecuteEvent
   1 Trying *Selection* code: 
      1 <verbatim>cp /uscms/home/yumiceva/work/sframe/CMSSW_5_3_3/src/SFrame/FITAnalysis/src/FourtopsSelectionCycle.cxx* ./src
cp /uscms/home/yumiceva/work/sframe/CMSSW_5_3_3/src/SFrame/FITAnalysis/include/FourtopsSelectionCycle.h* ./include
cp /uscms/home/yumiceva/work/sframe/CMSSW_5_3_3/src/SFrame/FITAnalysis/config/FourtopsSelectionCycleMuons_config.xml* ./config
make
</verbatim>
      1 No problems building. But when executing sframe_main, It got problems loading libFITAnalysis
      1 SFrame/lib/libFITAnalysis.so: undefined symbol: _ZN22FourtopsSelectionCycle11ShowMembersER16TMemberInspector
      1 Tried: c++filt _ZN22FourtopsSelectionCycle11ShowMembersER16TMemberInspector
      1 Message: FourtopsSelectionCycle::ShowMembers(TMemberInspector&)
      1 <verbatim>nano include/FITAnalysis_LinkDef.h
                     #### ADD
                     #pragma link C++ class FourtopsSelectionCycle+;
rm src/FITAnalysis_Dict.*
rm obj/*
rm obj/dep/*</verbatim>
      1 Other issues (LumiFiles, PileUpHistos and JECFiles) where solved doing the following inside ./FITAnalysis directory: 
         1 <verbatim>cp /uscms/home/yumiceva/work/sframe/CMSSW_5_3_3/src/SFrame/ZprimeAnalysis2012/PileUpHistos/PileUpHistoCycle.DATA.ElectronHad_208686.root ./PileUpHistos/
cp /uscms/home/yumiceva/work/sframe/CMSSW_5_3_3/src/SFrame/ZprimeAnalysis2012/LumiFiles/HLT_Ele30_CaloIdVT_TrkIdT_PFNoPUJet100_PFNoPUJet25_208686.root ./LumiFiles/
cp -r /uscms/home/yumiceva/work/sframe/CMSSW_5_3_3/src/SFrame/ZprimeAnalysis2012/JECFiles/ ./
</verbatim>

---+++ M) 01-Jul-2013 to 07-Jul-2013
   1 *QUESTIONS:* 
      1 
   1 *<Cycle Name=* , in config file , stores the name of the class
   1 Trying to install SFrame from the sources:
      1 Following: *https://twiki.cern.ch/twiki/bin/viewauth/CMS/B2GZprimeSFrameInstall#Installation_of_sframe_Current_t*
         1 %TWISTY{<br /> mode="div"<br />  showimgleft="%ICONURLPATH{toggleopen}%"<br />  hideimgleft="%ICONURLPATH{toggleclose}%"<br />  }% \x{fffd}verbatim3\x{fffd}%ENDTWISTY%
   1 guide USING GITHUB:
      1 %TWISTY{<br /> mode="div"<br />  showimgleft="%ICONURLPATH{toggleopen}%"<br />  hideimgleft="%ICONURLPATH{toggleclose}%"<br />  }%<verbatim>setenv SCRAM_ARCH slc5_amd64_gcc462
scram project CMSSW CMSSW_5_3_3
cd CMSSW_5_3_3/src
cmsenv

git clone https://github.com/UHHAnalysis/NtupleWriter UHHAnalysis/NtupleWriter
cvs co -d EGamma/EGammaAnalysisTools -r V00-00-31 UserCode/EGamma/EGammaAnalysisTools
scram b -j 8

svn co https://sframe.svn.sourceforge.net/svnroot/sframe/SFrame/tags/SFrame-03-06-11  SFrame
cd SFrame
source setup.csh
cp ~/fullsetup_SFRAME_from_B2GZ.csh ./
source fullsetup_SFRAME_from_B2GZ.csh
sed -i 's/-lpcre//g' core/Makefile
make

git clone https://github.com/UHHAnalysis/NtupleWriter.git NtupleWriter
git clone https://github.com/UHHAnalysis/SFrameTools.git SFrameTools
git clone https://github.com/UHHAnalysis/SFrameAnalysis.git SFrameAnalysis

cd NtupleWriter
/bin/sh ./configure.sh
make
cd ..


cp NtupleWriter/include/*.h SFrameTools/include/
cd SFrameAnalysis/include
wget --no-check-certificate https://raw.github.com/UHHAnalysis/SFrameAnalysis/v1-00/include/ExampleHists.h
cd ..
cd src
wget --no-check-certificate https://raw.github.com/UHHAnalysis/SFrameAnalysis/v1-00/src/ExampleHists.cxx
cd ../..

cd SFrameTools/include
wget --no-check-certificate https://raw.github.com/UHHAnalysis/SFrameTools/v1-00/include/ObjectHandler.h
cd ..
cd src
wget --no-check-certificate https://raw.github.com/UHHAnalysis/SFrameTools/v1-00/src/ObjectHandler.cxx
wget --no-check-certificate https://raw.github.com/UHHAnalysis/SFrameTools/v1-00/src/QjetsPlugin.cxx

cd ..
make 
cd JetMETObjects
make
cd ../../SFrameAnalysis
make

git clone git://github.com/yumiceva/FITAnalysis.git
cd FITAnalysis
        nano  src/FourtopsPreSelectionCycle.cxx
                   comment:  if(bcc->electrons) cleaner.ElectronCleaner(35,2.5,0.1,m_reversed_electron_selection,m_reversed_electron_selection); // min pt, eta, iso, and ID)
                   THE PROBLEM IS BECAUSE : m_reversed_electron_selection
make</verbatim>%ENDTWISTY%
   1 Messages:
      1 libSFrameTools.so: undefined symbol: _ZN17LuminosityHandler11ShowMembersER16TMemberInspector
   1  *FINAL GUIDE!!*  Using CMSSW
      1 <verbatim>kserver_init
svn co https://sframe.svn.sourceforge.net/svnroot/sframe/SFrame/tags/SFrame-03-06-11  SFrame
cd SFrame
echo "#\!/bin/csh" > fullsetupSFrame.csh
echo "setenv FASTJETDIR /uscms_data/d3/pjacome/fastjet/3.0.2-install/lib" >> fullsetupSFrame.csh
echo "source setup.csh" >> fullsetupSFrame.csh
echo setenv LD_LIBRARY_PATH  \$\{FASTJETDIR\}:\$\{SFRAME_DIR\}/SFrameTools/JetMETObjects/lib:\$\{LD_LIBRARY_PATH\} >> fullsetupSFrame.csh
source fullsetupSFrame.csh
sed -i 's/-lpcre//g' core/Makefile
make

cvs co -d SFrameAnalysis -r Feb-15-2013-v1  UserCode/UHHAnalysis/SFrameAnalysis
cvs co -d SFrameTools -r Feb-15-2013-v1 UserCode/UHHAnalysis/SFrameTools
cd SFrameTools/include
cvs co -d Objects -r Feb-15-2013-v1  UserCode/UHHAnalysis/NtupleWriter/Objects
cvs co -d HEAD UserCode/UHHAnalysis/SFrameTools/include/
cp ./HEAD/BaseHists.h ./
cd ..
make
cd JetMETObjects 
make 
cd ../../SFrameAnalysis
make
cd ..

git clone git://github.com/yumiceva/FITAnalysis.git
cd FITAnalysis
make</verbatim>


---+++ L) 24-Jun-2013 to 30-Jun-2013
   1 *QUESTIONS:* 
      1 LHE_NtupleWriter.C in Newer versions of NtupleWriter
      1 In B2GZ they use other versions of software, ex: *git clone https://github.com/UHHAnalysis/NtupleWriter UHHAnalysis/NtupleWriter*
      1 What is PROOF run mode in config.xml files
      1 ERROR when executing SFrame_main: 
         1 \x{fffd}verbatim6\x{fffd}
      1 I think I know how to do the cuts with NTupleWriter: 
         1 What is minreliso ?
         1 At least 6 jets?
      1 Don't know how to do cuts with FITAnalysis
   1 Finally I solved all the problems: 
      1 Just copy the whole SFrame directory: 
         1 \x{fffd}verbatim7\x{fffd}
      1 With the previous step all works ok. But I also try to change some references without success so, don't do the following: 
         1 Go to SFrame directory: 
            1 \x{fffd}verbatim8\x{fffd}
         1 Change SFrame references: 
            1 \x{fffd}verbatim9\x{fffd}
         1 Changing FASTJET references: 
            1 \x{fffd}verbatim10\x{fffd}
         1 Now check again for other references: 
            1 \x{fffd}verbatim11\x{fffd}
         1 To keep PileUp information, Lumi files, etc. Copy again config files: 
            1 \x{fffd}verbatim12\x{fffd}
   1 Finish with: 
      1 \x{fffd}verbatim13\x{fffd}
   1 Creating Ntuples from PATuples using Francisco's recipe: 
      1 Using my PATuples: I got errors: 
         1 %TWISTY{<br /> mode="div"<br />  showimgleft="%ICONURLPATH{toggleopen}%"<br />  hideimgleft="%ICONURLPATH{toggleclose}%"<br />  }% \x{fffd}verbatim14\x{fffd}%ENDTWISTY%
      1 Using PATuples from Titas Roy, it succeded!: 
         1 \x{fffd}verbatim15\x{fffd}
         1 Note: I just changed the source file. I didn't change any other code line. 
            1 \x{fffd}verbatim16\x{fffd}
   1 *Skimmed Analysis ntuples (Preselection)* 
      1 Edit: config/FourtopsPreSelectionCycleMuons_config.xml 
         1 Leave *<!ENTITY TTTT* line. Delete the other ENTITIES definitions
         1 Edit *config/53xNTuples-v2/TTTT.xml* and put the location of the NTuple (Ntuple made using NTupleWriter)

---+++ K) 17-Jun-2013 to 23-Jun-2013
   1 Following new prescription: https://twiki.cern.ch/twiki/bin/view/Main/FYumicevaSFrameFwk 
      1 SFrame installation: 
         1 ROOT *(OK)*
         1 CMSSW 5 3 3 *(OK)*
         1 NTuple writer (in CMSSW) *(OK)*
         1 FastJet *(OK)* -> Next time could be installed in other directory, outside /CMSSW_5_3_3/src
      1 
         1 
            1 fastjet-3.0.4.tar.gz
            1 %TWISTY{<br /> mode="div"<br />  showimgleft="%ICONURLPATH{toggleopen}%"<br />  hideimgleft="%ICONURLPATH{toggleclose}%"<br />  }% \x{fffd}verbatim17\x{fffd} %ENDTWISTY%
         1 SFrame<strong> (OK)
<strong> </strong>

<strong>
---+++ J) 12-Jun-2013 to 16-Jun-2013

   1 Stat I - Hands-on Tutorial Session (HATS) at LPC: https://indico.cern.ch/conferenceDisplay.py?confId=253180 
      1 Is it mandatory to create a parameter with a defined range?
      1 Page 7, https://indico.cern.ch/getFile.py/access?sessionId=1&resId=1&materialId=0&confId=253180
      1 Created a signal Gaussian using parameters without ranges and works perfectly.
   1 Reading "Training and classification of BDT" for thesis plan.
   1 Reading "Hypothesis testing" from CLASHEP Perú slides: http://indico.cern.ch/getFile.py/access?contribId=20&resId=0&materialId=slides&confId=208901
   1 Analysis of Ntuples: 
      1 Source: https://twiki.cern.ch/twiki/bin/view/Main/FYumicevaTopFwk#Analysis_of_ntuples
      1 After executing "cvs co ..." code, do (this avoids dependencies problems): 
         1 cvs co -r V00-03-03 -d Yumiceva/Top7TeV UserCode/Yumiceva/Top7TeV
         1 scram b
      1 Another Error messages: 
         1 <verbatim>20:03:15  8637 Wrk-0.17 | Error in <TFile::TFile>: file /query-result.root does not exist
20:03:15  8637 Wrk-0.17 | Info in <TProofServLite::HandleArchive>: file cannot be open (/query-result.root)
20:03:15  8572 Wrk-0.10 | Error in <TFile::TFile>: file /query-result.root does not exist
20:03:15  8572 Wrk-0.10 | Info in <TProofServLite::HandleArchive>: file cannot be open (/query-result.root)
</verbatim>
      1 To plot root histograms: 
         1 root -l results_data.root 
            1 <verbatim>_file0->GetListOfKeys()->Print();
TH1F *h = (TH1F*)_file0->Get("h1test");
h->Draw();
</verbatim>
         1 Default example just plot h1test muon PT 
            1 <img width="350" alt="example_ntuples_analysis.png" src="https://twiki.cern.ch/twiki/pub/Sandbox/LogbookFourTopAnalysisPabloJacome/example_ntuples_analysis.png" style="vertical-align: middle;" title="example_ntuples_analysis.png" height="237" />
         1 Plotting leading muon histogram. Change *SimpleAnalyzer.C* : 
            1 <verbatim>  TopMuonEvent max_pt_muon;
  max_pt_muon.pt = 0.;
  for ( size_t imu=0; imu < total_muons; ++imu) {
    TopMuonEvent muon = muons[imu];
    if ( muon.pt > max_pt_muon.pt){
      max_pt_muon.pt = muon.pt;
    }
  }
  h1test->Fill( max_pt_muon.pt );
</verbatim>
            1 <img width="350" alt="leading_muon.png" src="https://twiki.cern.ch/twiki/pub/Sandbox/LogbookFourTopAnalysisPabloJacome/leading_muon.png" style="vertical-align: middle;" title="leading_muon.png" height="237" />

---+++ I) 27-May-2013 to 02-Jun-2013
   1 Due to problems in FourTop model for MadGraph (low Xsec) 
      1 
</strong>
<table align="center" border="1"> <tbody> <tr> <td> *Source* </td> <td> *X Sec. SM,<br />14TeV<br /> p p &gt; t t- t t-* </td> <td> *X Sec. GH 700,<br />14TeV,<br /> pp-&gt;GH GH-&gt; t t- t t-* </td> <td> *X Sec. GH 700,<br />14TeV,<br /> pp-&gt;GH GH* </td> <td> *X Sec. GH 700,<br />14TeV,<br /> pp-&gt;GH GH -&gt; b b- b b-* </td> <td> *Reference* </td> </tr> <tr> <td>MadGraph</td> <td>0.00779 pb</td> <td>0.00337 pb</td> <td>0.413 pb</td> <td>problems</td> </tr> <tr> <td> *Paper:* <br />Massive color-octet ..</td> <td>-</td> <td>-</td> <td>0.4 pb</td> <td>0.02 pb</td> <td>http://arxiv.org/pdf/0709.2378v3.pdf</td> </tr> <tr> <td> *Paper:* <br />Top compositeness ..</td> <td>0.0036 pb</td> <td>-</td> <td>-</td> <td>-</td> <td>http://arxiv.org/pdf/0712.3057v1.pdf</td> </tr> </tbody> </table>

<strong>

---+++ H) 13-May-2013 to 19-May-2013
   1 Reading about analysis script: https://twiki.cern.ch/twiki/bin/view/Main/FYumicevaTopFwk#Analysis_of_ntuples
   1 

---+++ G) 06-May-2013 to 12-May-2013
   1 Requesting severals letters (to EPN and IECE) for applying a bursary for SENESCYT
   1 Making the application
   1 Writing paper for EPN magazine
   1 Calculating number of events of Four Tops in SM and heavy colorons 
      1 pp collisions Xsec
</strong>
<table align="center" border="1"> <tbody> <tr> <td> </td> <td><latex> \sigma</latex>Tot</td> <td><latex> \sigma</latex> elastic</td> <td><latex> \sigma</latex> Inelastic</td> <td>Source</td> </tr> <tr> <td>8 Tev</td> <td>100 mb</td> <td>28 mb</td> <td>72 mb</td> <td>http://arxiv.org/pdf/1105.4916v3.pdf</td> </tr> <tr> <td>8 Tev</td> <td>101.7 mb</td> <td>27.1 mb</td> <td>74.7 mb</td> <td>http://arxiv.org/pdf/1303.2927v2.pdf</td> </tr> </tbody> </table> <table align="center" border="1"> <tbody> <tr> <td>Xsec/Model using MadGraph5<br /></td> <td>Standard Model</td> <td>Heavy scalar octects<br /> mGH=700GeV</td> </tr> <tr> <td><latex> \sigma pp \rightarrow t \bar{t}t \bar{t}</latex></td> <td>7.133E-4 [pb]</td> <td>2.32E-4 [pb]</td> </tr> </tbody> </table> <table align="center" border="1"> <tbody> <tr> <td>CMS Total<br /></td> <td>21.79 fb-1 * 101.7 mb</td> <td>= 2.21E15 events</td> </tr> <tr> <td>SM expected</td> <td>21.79 fb-1 * 0.7133 fb</td> <td>= 15.54 events</td> </tr> <tr> <td>Heavy scalar bosons expected</td> <td>21.79 fb-1 * 2.32E-4 pb</td> <td>= 5.05 events</td> </tr> </tbody> </table>

<strong>

---+++ F) 29-April-2013 to 05-May-2013
   1 Questions: 
      * Difference between : GEN,SIM and RECO vs GEN,FASTSIM, HLT
   1 https://twiki.cern.ch/twiki/bin/view/Main/FranciscoYumicevaSim
   1 Running code for plot histograms of GEN-FASTSIM-HLT root file : 
      1 <verbatim>cvs co UserCode/Yumiceva/Scripts
python GenParticleAnalyzer.py -i Hadronizer_TuneZ2star_8TeV_generic_LHE_pythia_tauola_cff_py_GEN_FASTSIM_HLT_PU.root</verbatim>
   1 Reading: https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookPATWorkflow 
      1 %RED%"The Analysis Tool (AT) group of CMS strongly recommends to use pat::Tuples instead of flat ntuples due to the following reasons:"%ENDCOLOR%
   1 Reading log output when created "tlbsm_53x_v2_mc.root" (PatTuple step), found an error: %TWISTY{<br /> mode="div"<br />  showimgleft="%ICONURLPATH{toggleopen}%"<br />  hideimgleft="%ICONURLPATH{toggleclose}%"<br />  }% <verbatim>Begin processing the 1st record. Run 1, Event 1, LumiSection 1 at 18-Apr-2013 09:36:03.058 CDT
   ----- Begin Fatal Exception 18-Apr-2013 09:36:03 CDT-----------------------
   An exception of category 'ProductNotFound' occurred while
      [0] Processing run: 1 lumi: 1 event: 1
      [1] Running path 'p0'
      [2] Calling event method for module HBHENoiseFilter/'HBHENoiseFilter'
   Exception Message:
    could not find HcalNoiseSummary.
      Additional Info:
         [a] If you wish to continue processing events after a ProductNotFound exception,
   add "SkipEvent = cms.untracked.vstring('ProductNotFound')" to the "options" PSet in the configuration.
   
   ----- End Fatal Exception -------------------------------------------------</verbatim> %ENDTWISTY% 
      1 Find the strings: <br /># process all the events<br />process.maxEvents.input = 100
      1 and add: *THIS JUST KEEP ITERATING OVER ALL EVENTS EVEN WHEN AN ERROR HAPPENED IN THE FIRST EVENT.* 
         1 <verbatim>process.options = cms.untracked.PSet(
    SkipEvent = cms.untracked.vstring('ProductNotFound')
)</verbatim>
      1 Add CMSSW package which contains HcalNoiseSummary.cc. *BUT THIS DIDN'T FIX THE PROBLEM*
      1 <blockquote> <verbatim>kserver_init    
addpkg DataFormats/METReco
scram b</verbatim> </blockquote>
      1 Removing all the lines containing: "HBHENoiseFilter" , "" 
         1 <verbatim>## The iso-based HBHE noise filter ___________________________________________||
#process.load('CommonTools.RecoAlgos.HBHENoiseFilter_cfi')

## The CSC beam halo tight filter ____________________________________________||
#process.load('RecoMET.METAnalyzers.CSCHaloFilter_cfi')

################################################################################################

process.filtersSeq = cms.Sequence(
   process.primaryVertexFilter *
   process.noscraping *
#   process.HBHENoiseFilter *
#   process.CSCTightHaloFilter *
   process.hcalLaserEventFilter *
   process.EcalDeadCellTriggerPrimitiveFilter *
   process.goodVertices * process.trackingFailureFilter *
   process.eeBadScFilter
)

##################################################################################################

#if options.runOnFastSim:
#    process.patseq.remove( process.HBHENoiseFilter )
#    process.patseq.remove( process.CSCTightHaloFilter )

</verbatim>
   1 Doing PatTuple Tutorial with my GEN_FASTSIM_HLT.root : 
      1 https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookPATTupleCreationExercise 
         1 PhysicsTools/PatAlgos/test/patTuple_standard_cfg.py
         1 edmDumpEventContent patTuple.root
      1 Exploring my GEN_FASTSIM_HLT.root file with ROOT/TBrowser I realized that there is just xxxx_HLT objects. ( No xxxx_RECO objects)
      1 Success creating *patTuple_standard.root*
   1 Now creating PatTuples following the tutorial: 
      1 https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookPATExampleTopQuarks 
         1 PhysicsTools/PatExamples/test/patTuple_topPreproduction_cfg.py
         1 edmDumpEventContent patTuple.root
   1 Doing NTuples: 
      1 https://twiki.cern.ch/twiki/bin/view/CMS/ExoticaWptb#Ntuple_Prescription
      1 Modifications: 
         1 $ nano Yumiceva/Top7TeV/test/TuplesFromPAT.py 
            1 <verbatim>if inputType=="MC":
    process.source = cms.Source("PoolSource",
                                fileNames = cms.untracked.vstring(
        'file:/uscms_data/d3/pjacome/exercises/ntuple_exoticaWptb/CMSSW_4_2_4/src/tlbsm_53x_v2_mc.root'
        ))</verbatim>
      1 $ cmsRun Yumiceva/Top7TeV/test/TuplesFromPAT.py useData=0 channel=muon events=100 &gt;& ntuple_100.log
      1 Error message: %TWISTY{<br /> mode="div"<br />  showimgleft="%ICONURLPATH{toggleopen}%"<br />  hideimgleft="%ICONURLPATH{toggleclose}%"<br />  }% 
         1 <verbatim>02-May-2013 11:33:19 CDT  Initiating request to open file file:/uscms_data/d3/pjacome/exercises/ntuple_exoticaWptb/CMSSW_4_2_4/src/tlbsm_53x_v2_mc.root
02-May-2013 11:33:20 CDT  Successfully opened file file:/uscms_data/d3/pjacome/exercises/ntuple_exoticaWptb/CMSSW_4_2_4/src/tlbsm_53x_v2_mc.root
02-May-2013 11:33:20 CDT  Closed file file:/uscms_data/d3/pjacome/exercises/ntuple_exoticaWptb/CMSSW_4_2_4/src/tlbsm_53x_v2_mc.root
%MSG-s CMSException:  AfterSourceConstruction 02-May-2013 11:33:20 CDT  pre-events
cms::Exception caught in cmsRun
---- FatalRootError BEGIN
Error occured while constructing primary input source.
Source is of type "PoolSource"
---- FatalRootError BEGIN
Fatal Root Error: @SUB=TStreamerInfo::Compile
Counter fNClusterRange should not be skipped from class TTree

Error occurred while creating source PoolSource
---- FatalRootError END
---- FatalRootError END


%MSG</verbatim>
         * %ENDTWISTY%
      * Now following new recipe from Francisco Yumiceva: <a target="_blank" href="Main.FYumicevaTopFwk">https://twiki.cern.ch/twiki/bin/view/Main/FYumicevaTopFwk</a> *IT WORKS*
   1 SM PatTuples: 
      1 /eos/uscms/store/user/troy2012/PaTTuples_v1/
   1 Questions: 
      1 DO I have to tunning HCAL signal response? 
         1 https://twiki.cern.ch/twiki/bin/view/CMSPublic/FastSimTuningHCALSignal
      1 Will I have to run EDAnalyzer for patTuple.root to get analyzePatMuons.root?
      1 I Did PatTuples using PhysicsTools/PatAlgos/test/patTuple_standard_cfg.py
      1 Event Selection? analyzerTopSelection_cfg.py
      1 What if I change CMSSW version of the prescription?
      1 Learning how to do NTuples: 
         1 https://twiki.cern.ch/twiki/bin/view/CMS/ExoticaWptb#Ntuple_Prescription
      1 Do I have to change/install pythia 8 ?
      1 If I began using an especific version of CMSSW, Do I have to keep working with this for all the processes (MG, FASTSIM, PatTuples, NTuples) ?
---+++ E) 22-April-2013 to 26-April-2013
   * Stuying QFT
   * Writing thesis plan.
   * Reading Dobrescu and 4 top papers: 
      * "Weak-triplet, color-octet scalars, and the CDF dijet excess", Dobrescu, Bogdan A. and Krnjaic, Gordan Z.
      * "Top compositeness at the Tevatron and LHC", Ben Lillie and Jing Shu and Tim M.P. Tait
      * "Massive color-octet bosons and pairs of resonances at hadron colliders", Dobresc
      * "Heavy octets and Tevatron signals with three or four b jets"
      * "Strong Dynamics and Electroweak Symmetry Breaking", Christopher T. Hill and Elizabeth H. Simmon
   * Questions: 
      1 WHat is current state of heavy scalar bosons since the higgs discovery? ( Heavy scalar bosons triy to explain electroweak symmetry breaking)
      1 Finding Dobrescu predictions for two pairs of t t~ , not just heavy scalar bosons x-sec.
      1 Found 3.6 fb for tt~tt~ in SM. What's energy? It is not visible over background?
      1 Dind't find yet the value for tt~tt~ x-sec in color-octet model.
      1 FIg 4 of "Massive color-octet bosons and pairs of resonances at hadron colliders"
      1 Fig 2. in "Weak-triplet, color-octet scalars, and the CDF dijet excess", Dobrescu, Bogdan A. and Krnjaic, Gordan Z., <br />doesn't seems to be for Heavy escalar bosons because it says " charged octo-triplet"
      1 It is b b~ decays greater than t t~ decay for GH bosons? in "Massive color-octet bosons and pairs of resonances at hadron colliders", Dobrescu, 
         1 <verbatim>(Section 5)
...
If the coefficient cs is flavor independent, GH decays predominantly 
into bb~ or tt~ for masses above 350 GeV. However, it is possible that
cs is nonzero only for down-type quarks, so that even above the
 tt~ threshold the dominant decay is into b b~
...</verbatim>
         1 <verbatim> 
 (Section 6)
...
and a 100% branching fraction
of GH to b b~(next-to-leading order corrections, which
may be sizable, are not included). The production
cross section is huge, more than 20 times that for
a quark of equal mass, allowing for early discovery.
...</verbatim>
---+++ D) 18-April-2013
   * PatTuple creation 
      * $ cp TopQuarkAnalysis/TopPairBSM/test/ttbsm_cfg.py TopQuarkAnalysis/TopPairBSM/test/4topbsm_cfg.py
      * $ nano TopQuarkAnalysis/TopPairBSM/test/4topbsm_cfg.py 
         * process.source.fileNames = [<br /># '/store/mc/Summer12_DR53X/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/AODSIM/PU_S10_START53_V7A-v1/0000/0076C8E3-9AE1-E111-917C-003048D439AA.root'<br /> 'file:/eos/uscms/store/user/pjacome/FOURTOPS_GH750_1000_1/Hadronizer_TuneZ2star_8TeV_generic_LHE_pythia_tauola_cff_py_GEN_FASTSIM_HLT_PU.root'<br /> ]
      * Edit input filename
      * $ cmsRun TopQuarkAnalysis/TopPairBSM/test/4topbsm_cfg.py useData=0 &gt;& output.log
---+++ C) 17-April-2013
   * Reading about Proca model
   * Reading about Local gauge invariance, massless bosons, chiral symmetries, helicity, etc. (... homework for a master's subject)
   * Learning how create PATtuple following: 
      * https://twiki.cern.ch/twiki/bin/viewauth/CMS/B2GTopLikeBSM53X#Version_2_53x_pre_Moriond_versio
---+++ B) 16-April-2013
   * Finish Chapter 2, QFT David Tong
   * Testing the scripts. 1 job of 1000 events lasts 2:30 hours.
---+++ A) 15-April-2013

Creating scripts to send several simulation jobs automatically.
   * 4top_mg5.csh , for madgraph
   * 4top_cmssw.csh , for GEN, FASTSIM and HLT 
      * 
</strong>
Attachments
Topic attachments
I	Attachment	History	Action	Size	Date	Who
png	SM-four-top_fig1.png	r1	manage	355.8 K	2013-05-14 - 17:32	PabloJacome
png	SM-four-top_fig2.png	r1	manage	454.3 K	2013-05-14 - 17:32	PabloJacome
png	example_ntuples_analysis.png	r1	manage	9.4 K	2013-06-14 - 16:58	PabloJacome
png	lead_muon_ntuples_analysis.png	r1	manage	9.4 K	2013-06-14 - 17:46	PabloJacome
png	leading_muon.png	r1	manage	9.1 K	2013-06-14 - 17:50	PabloJacome
Topic revision: r64 - 2013-09-13 - PabloJacome
Account
- Log In
- Cern Search
- TWiki Search
- Google Search
Sandbox All webs
~~Edit~~
~~Attach~~
Copyright &© 2008-2023 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