Pablo Jacome's Four Top Logbook

More... Close

References

* http://arxiv.org/abs/1001.0221
* “Four-heavy-quark hadroproduction”, V. Barger, A. L. Stange, and R. J. N. Phillips, Phys. Rev. D 44 (Oct, 1991) 1987–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. TFile* file = new TFile("FourtopsPreSelectionCycle.DATA.DATA.root","READ");
TTree* mytree = (TTree*)(file -> Get("AnalysisTree"));
mytree->Draw("goodOfflinePrimaryVertices.m_nTracks");
  2. Superimposed and adding histograms: 
    1. 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");

O) 15-Jul-2013 to 21-Jul-2013

  1. QUESTIONS:
  2. 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
    2. Add the following in the config file, inside tag definition: 
      1.   RunMode="PROOF" ProofServer="lite://"  ProofWorkDir="__CONDOR__"
RunMode="PROOF" ProofServer="lite"  ProofWorkDir=""
    3. Files for sending CONDOR jobs. sframe_condor_exe.csh: 
      1. #! /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
    4. sframe_condor_proof_cfg.txt 
      1. 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
    5. Send job: 
      condor_submit  sframe_condor_proof_cfg.txt
    6. State of the job: 
      condor_q -bet <JOB_ID>
    7. Killing 'idle' jobs that have been more than 24 hours: 
      condor_rm <JOB_ID>
    8. Changing settings of condor job: 
      1. Requirements = Memory >= 199 && OpSys == "LINUX"&& (Arch != "DUMMY" )&& Disk > 1000000
  3. Playing with ROOT: 
    1. AnalysisTree->Scan("event")
  4. Learning about Histos:
    1. http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/algomez/python/plothistos.py?view=log 
      1. python plotHistos.py -i /eos/uscms/store/user/yumiceva/fourtops/sframePreselection/FourtopsPreSelectionCycle.MC.TTTT.root  -o output.root -p 6
      2. Didn't work. It didn't throw output.root
    2. http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/Yumiceva/TreeAnalyzer/test/Analyzer.C?view=log
    3. 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
    2. What is BaseCycleContainer?
    3. Why Selection doesn't take input data drom PreSelection?
  2. Exploring src/FourtopsPreSelectionCycle.cxx and src/FourtopsSelectionCycle.cxx
    1. To compare: /uscms/home/yumiceva/work/sframe/CMSSW_5_3_3/src/SFrame/FITAnalysis/src
  3. Preselection FourtopsPreSelectionCycle.cxx :
    1. "At least one good muon"
      1. Changed in the function: FourtopsPreSelectionCycle::BeginInputData
      2. preselection->addSelectionModule(new NElectronSelection(0,0));//no electron preselection->addSelectionModule(new NMuonSelection(1,int_infinity()));
    2. "min pt, max eta , min reliso, and muon ID " for muon, electrons and jets.
      1. Changed in the function: FourtopsPreSelectionCycle::ExecuteEvent
  4. Trying Selection code: 
    1. 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
    2. No problems building. But when executing sframe_main, It got problems loading libFITAnalysis
    3. SFrame/lib/libFITAnalysis.so: undefined symbol: _ZN22FourtopsSelectionCycle11ShowMembersER16TMemberInspector
    4. Tried: c++filt _ZN22FourtopsSelectionCycle11ShowMembersER16TMemberInspector
    5. Message: FourtopsSelectionCycle::ShowMembers(TMemberInspector&) 
    6. nano include/FITAnalysis_LinkDef.h
                     #### ADD
                     #pragma link C++ class FourtopsSelectionCycle+;
rm src/FITAnalysis_Dict.*
rm obj/*
rm obj/dep/*
    7. Other issues (LumiFiles, PileUpHistos and JECFiles) where solved doing the following inside ./FITAnalysis directory: 
      1. 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/ ./

M) 01-Jul-2013 to 07-Jul-2013

  1. QUESTIONS:
  2. <Cycle Name= , in config file , stores the name of the class
  3. Trying to install SFrame from the sources:
    1. Following: https://twiki.cern.ch/twiki/bin/viewauth/CMS/B2GZprimeSFrameInstall#Installation_of_sframe_Current_t
      1. More... Close \x{fffd}verbatim3\x{fffd}
  4. guide USING GITHUB:
    1. More... Close 
      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
  5. Messages:
    1. libSFrameTools.so: undefined symbol: _ZN17LuminosityHandler11ShowMembersER16TMemberInspector
  6. FINAL GUIDE!! Using CMSSW 
    1. 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

L) 24-Jun-2013 to 30-Jun-2013

  1. QUESTIONS:
    1. LHE_NtupleWriter.C in Newer versions of NtupleWriter
    2. In B2GZ they use other versions of software, ex: git clone https://github.com/UHHAnalysis/NtupleWriter UHHAnalysis/NtupleWriter
    3. What is PROOF run mode in config.xml files
    4. ERROR when executing SFrame_main:
      1. \x{fffd}verbatim6\x{fffd}
    5. I think I know how to do the cuts with NTupleWriter:
      1. What is minreliso ?
      2. At least 6 jets?
    6. Don't know how to do cuts with FITAnalysis
  2. Finally I solved all the problems:
    1. Just copy the whole SFrame directory:
      1. \x{fffd}verbatim7\x{fffd}
    2. 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}
      2. Change SFrame references:
        1. \x{fffd}verbatim9\x{fffd}
      3. Changing FASTJET references:
        1. \x{fffd}verbatim10\x{fffd}
      4. Now check again for other references:
        1. \x{fffd}verbatim11\x{fffd}
      5. To keep PileUp information, Lumi files, etc. Copy again config files:
        1. \x{fffd}verbatim12\x{fffd}
  3. Finish with:
    1. \x{fffd}verbatim13\x{fffd}
  4. Creating Ntuples from PATuples using Francisco's recipe:
    1. Using my PATuples: I got errors:
      1. More... Close \x{fffd}verbatim14\x{fffd}
    2. Using PATuples from Titas Roy, it succeded!:
      1. \x{fffd}verbatim15\x{fffd}
      2. Note: I just changed the source file. I didn't change any other code line.
        1. \x{fffd}verbatim16\x{fffd}
  5. Skimmed Analysis ntuples (Preselection)
    1. Edit: config/FourtopsPreSelectionCycleMuons_config.xml
      1. Leave <!ENTITY TTTT line. Delete the other ENTITIES definitions
      2. 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)
      2. CMSSW 5 3 3 (OK)
      3. NTuple writer (in CMSSW) (OK)
      4. FastJet (OK) -> Next time could be installed in other directory, outside /CMSSW_5_3_3/src
        1. fastjet-3.0.4.tar.gz
        2. More... Close \x{fffd}verbatim17\x{fffd}
      2. SFrame (OK)

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?
    2. Page 7, https://indico.cern.ch/getFile.py/access?sessionId=1&resId=1&materialId=0&confId=253180
    3. Created a signal Gaussian using parameters without ranges and works perfectly.
  2. Reading "Training and classification of BDT" for thesis plan.
  3. Reading "Hypothesis testing" from CLASHEP Perú slides: http://indico.cern.ch/getFile.py/access?contribId=20&resId=0&materialId=slides&confId=208901
  4. Analysis of Ntuples:
    1. Source: https://twiki.cern.ch/twiki/bin/view/Main/FYumicevaTopFwk#Analysis_of_ntuples
    2. After executing "cvs co ..." code, do (this avoids dependencies problems):
      1. cvs co -r V00-03-03 -d Yumiceva/Top7TeV UserCode/Yumiceva/Top7TeV
      2. scram b
    3. Another Error messages: 
      1. 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)
    4. To plot root histograms:
      1. root -l results_data.root 
        1. _file0->GetListOfKeys()->Print();
TH1F *h = (TH1F*)_file0->Get("h1test");
h->Draw();
      2. Default example just plot h1test muon PT
        1. example_ntuples_analysis.png
      3. Plotting leading muon histogram. Change SimpleAnalyzer.C : 
        1.   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 );
        2. leading_muon.png

I) 27-May-2013 to 02-Jun-2013

  1. Due to problems in FourTop model for MadGraph (low Xsec)
Source X Sec. SM,
14TeV
p p > t t- t t- 		X Sec. GH 700,
14TeV,
pp->GH GH-> t t- t t- 		X Sec. GH 700,
14TeV,
pp->GH GH 		X Sec. GH 700,
14TeV,
pp->GH GH -> b b- b b- 		Reference
MadGraph 0.00779 pb 0.00337 pb 0.413 pb problems
Paper:
Massive color-octet ..		 - - 0.4 pb 0.02 pb http://arxiv.org/pdf/0709.2378v3.pdf
Paper:
Top compositeness ..		 0.0036 pb - - - http://arxiv.org/pdf/0712.3057v1.pdf

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

G) 06-May-2013 to 12-May-2013

  1. Requesting severals letters (to EPN and IECE) for applying a bursary for SENESCYT
  2. Making the application
  3. Writing paper for EPN magazine
  4. Calculating number of events of Four Tops in SM and heavy colorons
    1. pp collisions Xsec
Tot elastic Inelastic Source
8 Tev 100 mb 28 mb 72 mb http://arxiv.org/pdf/1105.4916v3.pdf
8 Tev 101.7 mb 27.1 mb 74.7 mb http://arxiv.org/pdf/1303.2927v2.pdf
Xsec/Model using MadGraph5
 Standard Model Heavy scalar octects
mGH=700GeV
7.133E-4 [pb] 2.32E-4 [pb]
CMS Total
 21.79 fb-1 * 101.7 mb = 2.21E15 events
SM expected 21.79 fb-1 * 0.7133 fb = 15.54 events
Heavy scalar bosons expected 21.79 fb-1 * 2.32E-4 pb = 5.05 events

F) 29-April-2013 to 05-May-2013

  1. Questions:
    • Difference between : GEN,SIM and RECO vs GEN,FASTSIM, HLT
  2. https://twiki.cern.ch/twiki/bin/view/Main/FranciscoYumicevaSim
  3. Running code for plot histograms of GEN-FASTSIM-HLT root file : 
    1. cvs co UserCode/Yumiceva/Scripts
python GenParticleAnalyzer.py -i Hadronizer_TuneZ2star_8TeV_generic_LHE_pythia_tauola_cff_py_GEN_FASTSIM_HLT_PU.root
  4. Reading: https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookPATWorkflow
    1. "The Analysis Tool (AT) group of CMS strongly recommends to use pat::Tuples instead of flat ntuples due to the following reasons:"
  5. Reading log output when created "tlbsm_53x_v2_mc.root" (PatTuple step), found an error: More... Close 
    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 -------------------------------------------------
    1. Find the strings:
      # process all the events
      process.maxEvents.input = 100
    2. and add: THIS JUST KEEP ITERATING OVER ALL EVENTS EVEN WHEN AN ERROR HAPPENED IN THE FIRST EVENT. 
      1. process.options = cms.untracked.PSet(
    SkipEvent = cms.untracked.vstring('ProductNotFound')
)
    3. Add CMSSW package which contains HcalNoiseSummary.cc. BUT THIS DIDN'T FIX THE PROBLEM 
    4. kserver_init    
addpkg DataFormats/METReco
scram b
    5. Removing all the lines containing: "HBHENoiseFilter" , "" 
      1. ## 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 )
  6. 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
      2. edmDumpEventContent patTuple.root
    2. Exploring my GEN_FASTSIM_HLT.root file with ROOT/TBrowser I realized that there is just xxxx_HLT objects. ( No xxxx_RECO objects)
    3. Success creating patTuple_standard.root
  7. Now creating PatTuples following the tutorial:
    1. https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookPATExampleTopQuarks
      1. PhysicsTools/PatExamples/test/patTuple_topPreproduction_cfg.py
      2. edmDumpEventContent patTuple.root
  8. Doing NTuples:
    1. https://twiki.cern.ch/twiki/bin/view/CMS/ExoticaWptb#Ntuple_Prescription
    2. Modifications:
      1. $ nano Yumiceva/Top7TeV/test/TuplesFromPAT.py 
        1. 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'
        ))
    3. $ cmsRun Yumiceva/Top7TeV/test/TuplesFromPAT.py useData=0 channel=muon events=100 >& ntuple_100.log
    4. Error message: More... Close 
      1. 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
  9. SM PatTuples:
    1. /eos/uscms/store/user/troy2012/PaTTuples_v1/
  10. Questions:
    1. DO I have to tunning HCAL signal response?
      1. https://twiki.cern.ch/twiki/bin/view/CMSPublic/FastSimTuningHCALSignal
    2. Will I have to run EDAnalyzer for patTuple.root to get analyzePatMuons.root?
    3. I Did PatTuples using PhysicsTools/PatAlgos/test/patTuple_standard_cfg.py
    4. Event Selection? analyzerTopSelection_cfg.py
    5. What if I change CMSSW version of the prescription?
    6. Learning how to do NTuples:
      1. https://twiki.cern.ch/twiki/bin/view/CMS/ExoticaWptb#Ntuple_Prescription
    7. Do I have to change/install pythia 8 ?
    8. 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)
    2. Finding Dobrescu predictions for two pairs of t t~ , not just heavy scalar bosons x-sec.
    3. Found 3.6 fb for tt~tt~ in SM. What's energy? It is not visible over background?
    4. Dind't find yet the value for tt~tt~ x-sec in color-octet model.
    5. FIg 4 of "Massive color-octet bosons and pairs of resonances at hadron colliders"
    6. Fig 2. in "Weak-triplet, color-octet scalars, and the CDF dijet excess", Dobrescu, Bogdan A. and Krnjaic, Gordan Z.,
      doesn't seems to be for Heavy escalar bosons because it says " charged octo-triplet"
    7. 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. (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~
...
      2.  
 (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.
...

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 = [
        # '/store/mc/Summer12_DR53X/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/AODSIM/PU_S10_START53_V7A-v1/0000/0076C8E3-9AE1-E111-917C-003048D439AA.root'
        'file:/eos/uscms/store/user/pjacome/FOURTOPS_GH750_1000_1/Hadronizer_TuneZ2star_8TeV_generic_LHE_pythia_tauola_cff_py_GEN_FASTSIM_HLT_PU.root'
        ]
    • Edit input filename
    • $ cmsRun TopQuarkAnalysis/TopPairBSM/test/4topbsm_cfg.py useData=0 >& output.log

C) 17-April-2013

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
Attachments Attachments
Topic attachments
I AttachmentSorted ascending History Action Size Date Who Comment
PNGpng example_ntuples_analysis.png r1 manage 9.4 K 2013-06-14 - 16:58 PabloJacome  
PNGpng lead_muon_ntuples_analysis.png r1 manage 9.4 K 2013-06-14 - 17:46 PabloJacome  
PNGpng leading_muon.png r1 manage 9.1 K 2013-06-14 - 17:50 PabloJacome  
PNGpng SM-four-top_fig1.png r1 manage 355.8 K 2013-05-14 - 17:32 PabloJacome  
PNGpng SM-four-top_fig2.png r1 manage 454.3 K 2013-05-14 - 17:32 PabloJacome  
Edit | Attach | Watch | Print version | History: r66 < r65 < r64 < r63 < r62 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r64 - 2013-09-13 - PabloJacome
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

  • Edit
  • Attach
This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 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