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. 
      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. 

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

1. Due to problems in FourTop model for MadGraph (low Xsec)

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

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
   • Now following new recipe from Francisco Yumiceva: https://twiki.cern.ch/twiki/bin/view/Main/FYumicevaTopFwk IT WORKS
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

• 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