Four Tops Analysis - Logbook for Aug-Sep 2011

My undergrad Project.
  • Leyend: Instructions (A,B,C), Steps (1,2,3), Running code (monospace), Results (i. ii. iii.), Notes (bullets)

For further studies, please look here.

2011-09-30 - Pythia8 fix

  1. After a long session, it is a good idea to close terminal and start again. After a fresh start, everything works. smile
  2. In Pythia8, Event listing looks different. After a better sight, test file seems right. pythia8ex3_cfg.log
  3. Generate 10k 4top events in MG5.
    1. fourtop500_v4_proc_card_mg5.dat
    2. fourtop1000_v9_proc_card_mg5.dat

2011-09-29 - Do all again.

  1. rm -R /uscms/home/algomez/nobackup/pythia8/*
  2. wget http://home.thep.lu.se/~torbjorn/pythia8/pythia8153.tgz .
  3. tar xvf pythia8153.tgz
  4. cmsrel CMSSW_4_2_5
  5. cd CMSSW_4_2_5/src/ & cmsenv
  6. cd ../../pythia8153/
  7. export USRCXXFLAGS=-fPIC
  8. ./configure --enable-64bits --with-hepmc=/uscmst1/prod/sw/cms/slc5_amd64_gcc434/external/hepmc/2.05.01-cms2
  9. make
  10. cd ../CMSSW_4_2_5/src & cmsenv
  11. kserver_init
  12. addpkg GeneratorInterface/Pythia8Interface
  13. cd  GeneratorInterface/Pythia8Interface/doc/
  14. cp pythia8_107 pythia8_153 and change folders to /uscms/home/algomez/nobackup/pythia8/pythia8153/
  15. =scramv1 setup pythia8_153 1.0 file:./pythia8_153=
  16. Change pythia8 to pythia8_153 in ../BuildFile.xml
  17. scramv1 b -r doesn't compile. Error:
    algomez@cmslpc16:~/nobackup/pythia8/CMSSW_4_2_5/src/GeneratorInterface/Pythia8Interface$ scramv1 b -r
    Resetting caches
    >> Local Products Rules ..... started
    >> Local Products Rules ..... done
    >> Entering Package GeneratorInterface/Pythia8Interface
    >>Creating project symlinks
    Entering library rule at GeneratorInterface/Pythia8Interface
    >> Compiling /uscms_data/d3/algomez/pythia8/CMSSW_4_2_5/src/GeneratorInterface/Pythia8Interface/src/LHAupLesHouches.cc 
    >> Compiling /uscms_data/d3/algomez/pythia8/CMSSW_4_2_5/src/GeneratorInterface/Pythia8Interface/src/Pythia8Hadronizer.cc 
    >> Compiling /uscms_data/d3/algomez/pythia8/CMSSW_4_2_5/src/GeneratorInterface/Pythia8Interface/src/RandomP8.cc 
    >> Building shared library tmp/slc5_amd64_gcc434/src/GeneratorInterface/Pythia8Interface/src/GeneratorInterfacePythia8Interface/libGeneratorInterfacePythia8Interface.so
    @@@@ Checking shared library for missing symbols: libGeneratorInterfacePythia8Interface.so
    tmp/slc5_amd64_gcc434/src/GeneratorInterface/Pythia8Interface/src/GeneratorInterfacePythia8Interface/libGeneratorInterfacePythia8Interface.so: undefined reference to `Pythia8::UserHooks::biasSelectionBy(Pythia8::SigmaProcess const*, Pythia8::PhaseSpace const*, bool)'
    collect2: ld returned 1 exit status
    gmake: *** [tmp/slc5_amd64_gcc434/src/GeneratorInterface/Pythia8Interface/src/GeneratorInterfacePythia8Interface/libGeneratorInterfacePythia8Interface.so] Error 1
             

2011-09-28 - Change lhe file

  1. Francisco and Stephen fix the problem.
    What Pythia does not know is the ID number so Steve suggested a very
    easy fix by just changing the particle ID of the GH scalar:
    
    sed -i -e 's/    10025    2/  1000021    2/' ~/lpc1/fourtop1000.lhe
    
    I copied your file to my area, and I run it with
    
    /uscms/home/yumiceva/work/CMSSW_4_2_5/src/GeneratorInterface/Pythia8Interface/test/pythia8ex3_cfg.py
    
    you can run the Particle list drawer with the file dump.py located in
    the same area.
    
    This solve the problem with Pythia8. Now Steve suggested two options
    to decay the tops: either we have the full decay chain of tops in MG
    ie. top->Wb->Wdecay,b or we tell pythia to decay the tops. We cannot
    have the intermediate step where Pythia make W decays.
    
    The best is to have MG decay fully the top quarks because it will take
    into account spin correlations. You have a scalar so one top is LH and
    the topbar is RH, and this is propagated to the W decay.
    
    Please try to generate a LHE file withe MG fully decay the tops, fix
    the ID number on the LHE file, and pass it to new Pythia8.
             
  2. Change paramenters in former lhe file fourtop1000_n10000.lhe:
    1. In /uscms/home/algomez/nobackup/pythia8150/CMSSW_4_2_5/src/
    2. cp /uscms/home/yumiceva/work/CMSSW_4_2_5/src/GeneratorInterface/Pythia8Interface/test/pythia8ex3_cfg.py .
    3. This is how pythia8ex3_cfg.py looks like.
    4. ln -s ~/nobackup/files/fourtop/fourtop1000_n10000.lhe
    5. sed -i -e 's/    10025    2/  1000021    2/' fourtop1000_n10000.lhe
    6. cmsRun pythia8ex3_cfg.py > fourtop1000_n10000.log
    • It doesn't seem to be fine.
  3. Generate again a 4top lhe file:
    1. fourtop500_v4: proc_card_mg5.dat & fourtop500_v4.lhe
    2. ln -s ~/nobackup/files/fourtop/fourtop500_v4.lhe
    3. sed -i -e 's/    10025    2/  1000021    2/' fourtop500_v4.lhe
    4. cmsRun pythia8ex3_cfg.py > fourtop500_v4.log
    • Doesn't work either.

2011-09-25 - Build Pythia8_150 in CMSSW

  1. Download Pythia8_150, compiled and got the Pythia8 interface compiled with the new version.
    1. cd pythia8150
    2. wget http://home.thep.lu.se/~torbjorn/pythia8/pythia8150.tgz .
    3. tar xvf pythia8150.tgz
    4. cd pythia8150/
    5. export USRCXXFLAGS=-fPIC
    6. ./configure --enable-64bits --with-hepmc=/uscmst1/prod/sw/cms/slc5_amd64_gcc434/external/hepmc/2.05.01-cms2
    7. make
    8. cd ../ & cmsrel CMSSW_4_2_5
    9. cd CMSSW_4_2_5/src/
    10. addpkg GeneratorInterface/Pythia8Interface
    11. cd  GeneratorInterface/Pythia8Interface/doc/
    12. cp pythia8_107 pythia8_153
      <doc type="BuildSystem::ToolDoc" version="1.0"/>
      <Tool name="pythia8_150" version="1.0">
      
      <Environment name="LIBDIR" default="/uscms/home/algomez/nobackup/pythia8150/pythia8150/lib">
      </Environment>
      
      <Environment name="INCLUDE" default="/uscms/home/algomez/nobackup/pythia8150/pythia8150/include">
      </Environment>
      
      <lib name="pythia8"/>
      <lib name="hepmcinterface"/>
      
      <Runtime name="PYTHIA8DATA" value="/uscms/home/algomez/nobackup/pythia8150/pythia8150/xmldoc"/>
      
      </Tool>
               
    13. scram setup pythia8_150 1.0 file:./pythia8_150
    14. vi ../BuildFile.xml
      <use name="FWCore/Framework"/>
      <use name="FWCore/ParameterSet"/>
      <use name="FWCore/MessageLogger"/>
      <use name="SimDataFormats/GeneratorProducts"/>
      <use name="GeneratorInterface/Core"/>
      <use name="GeneratorInterface/ExternalDecays"/>
      <use name="boost"/>
      <use name="clhep"/>
      <use name="hepmc"/>
      <use name="lhapdf"/>
      <use name="pythia8_150"/>
      <flags EDM_PLUGIN="1"/>
      <export>
         <lib name="GeneratorInterfacePythia8Interface"/>
      </export>
               
    15. cd ../../ & scramv1 b
  2. After compiled:
    1. cd GeneratorInterface/Pythia8Interface/test
    2. Test run cmsRun pythia8ex3_cfg.py works fine!.
    3. cp pythia8ex3_cfg.py fourtop1000_n10000.py It seems like:
      import FWCore.ParameterSet.Config as cms
      
      process = cms.Process("PROD")
      
      process.load("Configuration.StandardSequences.SimulationRandomNumberGeneratorSeeds_cff")
      
      #process.source = cms.Source("EmptySource")
      process.source = cms.Source("LHESource",
          fileNames = cms.untracked.vstring('file:fourtop1000_n10000.lhe')
      )
      
      process.generator = cms.EDFilter("Pythia8HadronizerFilter",
          maxEventsToPrint = cms.untracked.int32(1),
          pythiaPylistVerbosity = cms.untracked.int32(1),
          filterEfficiency = cms.untracked.double(1.0),
          pythiaHepMCVerbosity = cms.untracked.bool(False),
          comEnergy = cms.double(7000.),
          LHEInputFileName = cms.untracked.string('fourtop1000_n10000.lhe'),
          PythiaParameters = cms.PSet(
              pythia8_example03 = cms.vstring(
              'Top:gg2ttbar = on',
              'Top:qqbar2ttbar = on'),
              parameterSets = cms.vstring('pythia8_example03')
          )
      )
      
      process.MessageLogger = cms.Service("MessageLogger",
          cout = cms.untracked.PSet(
              default = cms.untracked.PSet(
                  limit = cms.untracked.int32(0)
              )
          ),
          destinations = cms.untracked.vstring('cout')
      )
      
      process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
          moduleSeeds = cms.PSet(
              generator = cms.untracked.uint32(123456),
              g4SimHits = cms.untracked.uint32(123456788),
              VtxSmeared = cms.untracked.uint32(123456789)
          ),
      )
      
      process.maxEvents = cms.untracked.PSet(
          input = cms.untracked.int32(10)
      )
      
      process.GEN = cms.OutputModule("PoolOutputModule",
          fileName = cms.untracked.string('fourtop1000_n10000.root')
      )
      
      process.p = cms.Path(process.generator)
      process.outpath = cms.EndPath(process.GEN)
      
      process.schedule = cms.Schedule(process.p, process.outpath)
               
    • Pythiaparameters from a ttbar decay, just for test.
    1. ln -s ~/nobackup/files/fourtop/fourtop1000_n10000.lhe
    • Adding ttbar parameters, seems to works fine.

2011-09-23 - Add Pythia8 150 in CMSSW

  1. Changed my fnal shell to bash smile (Finally)
  2. Build package again
    1. rm -R GeneratorInterface
    2. scramv1 b clean
    3. export PYTHIASUBVERSION=150
    4. export CMSSWVERSION=CMSSW_4_2_6
    5. source /afs/cern.ch/cms/sw/cmsset_default.sh
    6. scramv1 project CMSSW $CMSSWVERSION
    7. cd $CMSSWVERSION/src
    8. addpkg GeneratorInterface/Pythia8Interface
    9. cd  GeneratorInterface/Pythia8Interface/doc/
    10. cp pythia8_107 pythia8_150 pythia8_150 as before.
    11. scram setup pythia8_150 1.0 file: pythia8_150
    12. scramv1 tool list to check that scram register my pythia8150 tool
    13. Any change in BuildFile.xml
    14. scramv1 b Error:
      >> Local Products Rules ..... started
      >> Local Products Rules ..... done
      >> Building CMSSW version CMSSW_4_2_6 ----
      >> Entering Package GeneratorInterface/Pythia8Interface
      >>Creating project symlinks
      >> Building shared library
      tmp/slc5_amd64_gcc434/src/GeneratorInterface/Pythia8Interface/src/GeneratorInterfacePythia8Interface/libGeneratorInterfacePythia8Interface.so
      @@@@ Checking shared library for missing symbols:
      libGeneratorInterfacePythia8Interface.so
      /uscms_data/d3/algomez/pythia8150/CMSSW_4_2_6/external/slc5_amd64_gcc434/lib/libhepmcinterface.so:
      undefined reference to `HepMC::GenEvent::GenEvent(int, int,
      HepMC::GenVertex*, HepMC::WeightContainer const&, std::vector<long,
      std::allocator<long> > const&)'
      collect2: ld returned 1 exit status
      gmake: ***
      [tmp/slc5_amd64_gcc434/src/GeneratorInterface/Pythia8Interface/src/GeneratorInterfacePythia8Interface/libGeneratorInterfacePythia8Interface.so]
      Error 1
               
  3. Ask about this in HN: https://hypernews.cern.ch/HyperNews/CMS/get/generators/1236/2/2.html

2011-09-22 - Add Pythia8 150 in CMSSW

  1. After found a error in generation, Francisco ask Stephen Mrenna about it. He suggested to use Pythia8. But then, they said: Steve and I tried to run the 4top LHE file you produced with pythia8 but it crashes. Steve looked at the error and found out that it works with pythia version pythia8150
  2. He suggested to build Pythia8150 in my area. I followed the steps shown here: https://twiki.cern.ch/twiki/bin/view/CMS/Pythia8Interface
    1. In my /uscms_data/d3/algomez/ area: mkdir pythia8150
    2. setenv PYTHIASUBVERSION 150
    3. setenv CMSSWVERSION CMSSW_4_2_6
    4. bash /afs/cern.ch/cms/sw/cmsset_default.sh
    5. scramv1 project CMSSW $CMSSWVERSION
    6. cd $CMSSWVERSION/src
    7. cvs co -r HEAD GeneratorInterface/Pythia8Interface
    8. cd  GeneratorInterface/Pythia8Interface/doc/
    9. cp pythia8_107 pythia8_150
    10. pythia8_150 looks like:
      <doc type="BuildSystem::ToolDoc" version="1.0"/>
      <Tool name="pythia8_150" version="1.0">
      
      <Environment name="LIBDIR" default="/afs/cern.ch/sw/lcg/external/MCGenerators/pythia8/150/x86_64-slc5-gcc43-opt/lib">
      </Environment>
      
      <Environment name="INCLUDE" default="/afs/cern.ch/sw/lcg/external/MCGenerators/pythia8/150/x86_64-slc5-gcc43-opt/include">
      </Environment>
      
      <lib name="pythia8"/>
      <lib name="hepmcinterface"/>
      
      <Runtime name="PYTHIA8DATA" value="/afs/cern.ch/sw/lcg/external/MCGenerators/pythia8/150/share/xmldoc"/>
      
      </Tool>
               
    11. mv pythia8_107 pythia8_107_backup
    12. scram setup pythia8_150 1.0 file: pythia8_150
    13. scramv1 tool list to check that scram register my pythia8150 tool
    14. change BuildFile.xml from GeneratorInterface/Pythia8Interface
      <use name="FWCore/Framework"/>
      <use name="FWCore/ParameterSet"/>
      <use name="FWCore/MessageLogger"/>
      <use name="SimDataFormats/GeneratorProducts"/>
      <use name="GeneratorInterface/Core"/>
      <use name="GeneratorInterface/ExternalDecays"/>
      <use name="boost"/>
      <use name="clhep"/>
      <use name="pythia8"/>
      <use name="pythia8_150"/>
      <flags EDM_PLUGIN="1"/>
               
    15. scramv1 b Error:
      >> Local Products Rules ..... started
      >> Local Products Rules ..... done
      >> Entering Package GeneratorInterface/Pythia8Interface
      >>Creating project symlinks
      >> Compiling /uscms_data/d3/algomez/pythia8150/CMSSW_4_2_6/src/GeneratorInterface/Pythia8Interface/src/Pythia8Hadronizer.cc 
      In file included from /uscms_data/d3/algomez/pythia8150/CMSSW_4_2_6/src/GeneratorInterface/Pythia8Interface/src/Pythia8Hadronizer.cc:15:
      /uscms_data/d3/algomez/pythia8150/CMSSW_4_2_6/src/GeneratorInterface/Pythia8Interface/interface/UserHooks.h: In member function 'virtual double PtHatReweightUserHook::biasSelectionBy(const Pythia8::SigmaProcess*, const Pythia8::PhaseSpace*, bool)':
      /uscms_data/d3/algomez/pythia8150/CMSSW_4_2_6/src/GeneratorInterface/Pythia8Interface/interface/UserHooks.h:17: error: 'selBias' was not declared in this scope
      /uscms_data/d3/algomez/pythia8150/CMSSW_4_2_6/src/GeneratorInterface/Pythia8Interface/interface/UserHooks.h:20: error: 'selBias' was not declared in this scope
      /uscmst1/prod/sw/cms/slc5_amd64_gcc434/external/boost/1.44.0-cms3/include/boost/system/error_code.hpp: At global scope:
      /uscmst1/prod/sw/cms/slc5_amd64_gcc434/external/boost/1.44.0-cms3/include/boost/system/error_code.hpp:214: warning: 'boost::system::posix_category' defined but not used [-Wunused-variable]
      /uscmst1/prod/sw/cms/slc5_amd64_gcc434/external/boost/1.44.0-cms3/include/boost/system/error_code.hpp:215: warning: 'boost::system::errno_ecat' defined but not used [-Wunused-variable]
      /uscmst1/prod/sw/cms/slc5_amd64_gcc434/external/boost/1.44.0-cms3/include/boost/system/error_code.hpp:216: warning: 'boost::system::native_ecat' defined but not used [-Wunused-variable]
      gmake: *** [tmp/slc5_amd64_gcc434/src/GeneratorInterface/Pythia8Interface/src/GeneratorInterfacePythia8Interface/Pythia8Hadronizer.o] Error 1
               

2011-09-21 - New 4top Generation in MG

  1. Because the previous generation in MG5 takes a huge amount of time, Francisco suggested a simpler generation in MG5. With this, we will leave W decay for Pythia.
    1. Generate again 10k for 4tops with 500 GeV and 1 TeV
      1. Now at: /uscms_data/d3/algomez/MadGraph5_v1_3_16
      2. Sadly, I generated this sample in many different ways but it doesn't seem to be as we expect. A bunch of proc_cards with the log files with tree decay after hadronization (with the former method) is below.
        No. Proc Card Log File
        1 fourtop1000_proc_card_mg5.dat fourtop1000.log
        2 fourtop1000_v2_proc_card_mg5.dat fourtop1000_v2.log
        3 fourtop1000_v4_proc_card_mg5.dat fourtop1000_v4.log
        4 fourtop1000_v5_proc_card_mg5.dat fourtop1000_v5.log
        5 fourtop1000_v6_proc_card_mg5.dat fourtop1000_v6.log
        6 fourtop1000_v7_proc_card_mg5.dat fourtop1000_v7.log
        7 fourtop1000_v8_proc_card_mg5.dat fourtop1000_v8.log
        8 fourtop500_v1_proc_card_mg5.dat fourtop500_v1.log

2011-09-20 - 4top generation problem and tree decay.

  1. After printing the tree of decays generated in Pythia, there are not difference between leptons with status 1 and 3:
    1. Tree with particle information - here seems like status 1 and 3 is the same for leptons.
    2. Complete log file
    3. Python script - fourtops with tree decay. Important changes are:
      process.generator = cms.EDFilter("Pythia6HadronizerFilter",
          pythiaHepMCVerbosity = cms.untracked.bool(True),
          maxEventsToPrint = cms.untracked.int32(5),
          pythiaPylistVerbosity = cms.untracked.int32(1),
          comEnergy = cms.double(7000.0),
          PythiaParameters = cms.PSet(
            
  2. BR with leptons in status 3:
    1. BR for fourtop with 500 GeV
    2. BR for fourtop with 1 TeV
  3. Generation problem: After looking the decays of GH in ttbar, I realized that W and b are NOT decaying from ttbar. frown I made a mistake in the generation. frown
    1. Generate again 10k for 4tops with 500 GeV and 1 TeV
      1. import model FourTop
      2. define l+ = l+ ta+
      3. define l- = l- ta-
      4. generate p p > GH GH, (GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl))
      5. add process p p > GH GH, (GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j)) @2
      6. add process p p > GH GH, (GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j)) @3
      7. add process p p > GH GH, (GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl)) @4
      8. output fourtop500 and output fourtop100
      9. launch (nevents = 10000 and without cuts)
      10. For 500 GeV proc_card_mg5.dat and ttbar_tau_run_card.dat
      11. For 1 TeV proc_card_mg5.dat and ttbar_tau_run_card.dat
    2. Hadronized as cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:fourtop500.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 10000
    3. Run process in condor.
    4. Result files: fourtop500_n10000.root and fourtop1000_n10000.root
    5. After a considerably large time, this process had been canceled.
  4. In the meantime, some improvements in the code.

2011-09-18 - Plots

  1. Modified raw code for ttbar acceptance and some pt and eta plots:
    #######################################
    def acceptance(infile,outfile):
    #######################################
            events = Events (infile)
            handleGen = Handle ("vector<reco::GenParticle>")
            label = "genParticles"
    
            hAE = TH1F('h1','Acceptance',12,0,12)
            hAE.SetYTitle('percentage')
            hAE.SetStats(0)
    
            hpt = TH1F('partons pt for all jets', 'partons pt for all jets',100,0,300)
            heta = TH1F('partons eta for all jets','partons eta for all jets',100,-3,3)
    
            llpt = TH1F('leptons pt for lepton + jets', 'leptons pt for lepton + jets',100,0,100)
            lleta = TH1F('leptons eta for lepton + jets','leptons eta for lepton + jets',100,-10,10)
            lppt = TH1F('partons pt for lepton + jets', 'partons pt for lepton + jets',100,0,300)
            lpeta = TH1F('partons eta for lepton + jets','partons eta for lepton + jets',100,-3,3)
    
            dlpt = TH1F('leptons pt for dilepton', 'leptons pt for dilepton',100,0,100)
            dleta = TH1F('leptons eta for dilepton','leptons eta for dilepton',100,-10,10)
            dppt = TH1F('partons pt for dilepton', 'partons pt for dilepton',100,0,300)
            dpeta = TH1F('partons eta for dilepton','partons eta for dilepton',100,-3,3)
    
            nevent = 0 
            hadronically = 0
            haccepJets = 0
            leptonJets = 0
            laccepMuons = 0
            laccepElectrons = 0
            laccepJets = 0
            ltotal = 0
            dilepton = 0
            daccepMuons = 0
            daccepElectrons = 0
            daccepJets = 0          
            dtotal = 0
            for event in events:
                    nevent += 1
                    leptonsTau = 0
                    leptons = 0
                    hjets = 0
                    lmuons = 0 
                    lelectrons = 0
                    ljets = 0
                    dmuons = 0 
                    delectrons = 0
                    djets = 0
    
    #               print event
                    event.getByLabel (label,handleGen)
                    gens = handleGen.product()
    
                    for p in gens:
                            if p.mother():
                                    #print e, p.status(), p.pdgId(), mother
                                    if p.status() == 1:
                                            if ((abs(p.pdgId()) == 11) or (abs(p.pdgId()) == 13) or (abs(p.pdgId()) == 15)) and ((abs(p.mother().pdgId()) == 11) or (abs(p.mother().pdgId()) == 13) or (abs(p.mother().pdgId()) == 15)):
                                                    leptonsTau +=1
                                                    if not (abs(p.pdgId()) == 15) : leptons += 1
    #               print nevent, leptons
                    if leptonsTau == 0 : 
                            hadronically += 1
                            for p in gens:
                                    if p.status() == 3: 
                                            if (abs(p.pdgId()) in range (1,7)) and (abs(p.eta()) < 2.4):
                                                    hjets += 1
                                                    hpt.Fill(float(p.pt()))
                                                    heta.Fill(float(p.eta()))
                            if hjets: haccepJets += 1
                            #print nevent, htotal
    
                    if leptons == 1:
                            leptonJets += 1
                            for p in gens:
                                    if p.status() == 1:
                                            if (abs(p.pdgId()) == 13) and (abs(p.eta()) < 2.1): lmuons += 1
                                            elif (abs(p.pdgId()) == 11) and (abs(p.eta()) < 2.5): lelectrons += 1
                                            elif (abs(p.pdgId()) in range(11,15)):
                                                    llpt.Fill(float(p.pt()))
                                                    lleta.Fill(float(p.eta()))
                                    if p.status() == 3:
                                            if (abs(p.pdgId()) in range (1,7)) and (abs(p.eta()) < 2.4):
                                                    ljets += 1
                                                    lppt.Fill(float(p.pt()))
                                                    lpeta.Fill(float(p.eta()))
                            if lmuons: laccepMuons += 1
                            if lelectrons: laccepElectrons += 1
                            if ljets: laccepJets += 1
                            if (lmuons and lelectrons and ljets): ltotal +=1
    #                       print nevent, ltotal
    
                    elif leptons == 2:
                            dilepton += 1
                            for p in gens:
                                    if p.status() == 1:
                                            if (abs(p.pdgId()) == 13) and (abs(p.eta()) < 2.1): dmuons += 1
                                            elif (abs(p.pdgId()) == 11) and (abs(p.eta()) < 2.5): delectrons += 1
                                            elif (abs(p.pdgId()) in range(11,15)):
                                                    dlpt.Fill(float(p.pt()))
                                                    dleta.Fill(float(p.eta()))
                                    if p.status() == 3:
                                            if (abs(p.pdgId()) in range (1,7)) and (abs(p.eta()) < 2.4):
                                                    djets += 1
                                                    dppt.Fill(float(p.pt()))
                                                    dpeta.Fill(float(p.eta()))
                            if dmuons: daccepMuons += 1
                            if delectrons: daccepElectrons += 1
                            if djets: daccepJets += 1
                            if (dmuons and delectrons and djets): dtotal +=1
            hAE.Fill('alljets |eta| jets < 2.4', (haccepJets*100/hadronically))
            hAE.Fill('alljets without cuts', (hadronically*100/hadronically))
            hAE.Fill('leptonJet |eta| muons < 2.1', ((laccepMuons*100)/leptonJets))
            hAE.Fill('leptonJet |eta| electrons < 2.5', (laccepElectrons*100/leptonJets))
            hAE.Fill('leptonJet |eta| jets < 2.4', (laccepJets*100/leptonJets))
            hAE.Fill('leptonJet combined', (ltotal*100/leptonJets)) 
            hAE.Fill('leptonJet without cuts', (leptonJets*100/leptonJets))
            hAE.Fill('dilepton |eta| muons < 2.1', ((daccepMuons*100)/dilepton))
            hAE.Fill('dilepton |eta| electrons < 2.5', (daccepElectrons*100/dilepton))
            hAE.Fill('dilepton |eta| jets < 2.4', (daccepJets*100/dilepton))
            hAE.Fill('diletpon combined', (dtotal*100/dilepton))
            hAE.Fill('dilepton without cuts', (dilepton*100/dilepton))
            can1 = TCanvas('can1', '', 800, 600)
            can1.cd()
            hAE.Draw()
            can1.SetGridx()
            can1.SetGridy()
            can1.SetBottomMargin(0.2)
            can1.SaveAs(outfile+'_acceptance.png')
            del can1
    
            c1 = TCanvas('c1', '', 800, 600)
            heta.Draw()
            c1.SaveAs(outfile+'_alljets_eta.png') 
            hpt.Draw()
            c1.SaveAs(outfile+'_alljets_pt.png') 
            llpt.Draw()
            c1.SaveAs(outfile+'_leptonjets_lepton_pt.png')
            lleta.Draw()
            c1.SaveAs(outfile+'_leptonjets_lepton_eta.png')
            lppt.Draw()
            c1.SaveAs(outfile+'_leptonjets_parton_pt.png')
            lpeta.Draw()
            c1.SaveAs(outfile+'_leptonjets_parton_eta.png')
            dlpt.Draw()
            c1.SaveAs(outfile+'_dilepton_lepton_pt.png')
            dleta.Draw()
            c1.SaveAs(outfile+'_dilepton_lepton_eta.png')
            dppt.Draw()
            c1.SaveAs(outfile+'_dilepton_parton_pt.png')
            dpeta.Draw()
          
    1. Complete raw code
    2. For ttbar:
      1. BR in ttbar_V12.
      2. Total Acceptance Efficiency in ttbar_V12.
      3. Acceptance for different cases in ttbar_V12. Percentages compared with total number of each case.
      4. eta for partons in all+jets in ttbar_V12.
      5. pt for partons in all+jets in ttbar_V12.
      6. eta for leptons in lepton+jets in ttbar_V12.
      7. pt for leptons in lepton+jets in ttbar_V12.
      8. eta for partons in lepton+jets in ttbar_V12.
      9. pt for partons in lepton+jets in ttbar_V12.
      10. deltaR in lepton+jets after acceptance cuts in ttbar_V12.
      11. eta for leptons in dilepton in ttbar_V12.
      12. pt for leptons in dilepton in ttbar_V12.
      13. eta for partons in dilepton in ttbar_V12.
      14. pt for partons in dilepton in ttbar_V12.
      15. deltaR in dilepton after acceptance cuts in ttbar_V12.
  2. For 10k fourtop with GH 500 GeV:
    1. Complete raw code
    2. With any cuts:
      1. eta for tops.
      2. pt for tops.
      3. mass for tops.
      4. number of tops.
    3. Acceptance and BR:
      1. Branching Ratio.
      2. Total Acceptance Efficiency.
      3. Acceptance for different cases. Percentages compared with total number of each case.
    4. All+jets
      1. eta for partons in all+jets.
      2. pt for partons in all+jets.
    5. Lepton+Jets
      1. eta for leptons in lepton+jets.
      2. pt for leptons in lepton+jets.
      3. eta for partons in lepton+jets.
      4. pt for partons in lepton+jets.
      5. deltaR in lepton+jets after acceptance cuts.
    6. Dilepton+Jets
      1. eta for leptons in dilepton.
      2. pt for leptons in dilepton.
      3. eta for partons in dilepton.
      4. pt for partons in dilepton.
      5. deltaR in dilepton after acceptance cuts.
    7. 3 lepton+Jets
      1. eta for leptons in 3 lepton.
      2. pt for leptons in 3 lepton.
      3. eta for partons in 3 lepton.
      4. pt for partons in 3 lepton.
      5. deltaR in 3 lepton after acceptance cuts.
    8. 4 lepton+Jets
      1. eta for leptons in 4 lepton.
      2. pt for leptons in 4 lepton.
      3. eta for partons in 4 lepton.
      4. pt for partons in 4 lepton.
      5. deltaR in 4 lepton after acceptance cuts.
  3. For 10k fourtop with GH 1 TeV:
    1. With any cuts:
      1. eta for tops.
      2. pt for tops.
      3. mass for tops.
      4. number of tops.
    2. Acceptance and BR:
      1. Branching Ratio.
      2. Total Acceptance Efficiency.
      3. Acceptance for different cases. Percentages compared with total number of each case.
    3. All+jets
      1. eta for partons in all+jets.
      2. pt for partons in all+jets.
    4. Lepton+Jets
      1. eta for leptons in lepton+jets.
      2. pt for leptons in lepton+jets.
      3. eta for partons in lepton+jets.
      4. pt for partons in lepton+jets.
      5. deltaR in lepton+jets after acceptance cuts.
    5. Dilepton+Jets
      1. eta for leptons in dilepton.
      2. pt for leptons in dilepton.
      3. eta for partons in dilepton.
      4. pt for partons in dilepton.
      5. deltaR in dilepton after acceptance cuts.
    6. 3 lepton+Jets
      1. eta for leptons in 3 lepton.
      2. pt for leptons in 3 lepton.
      3. eta for partons in 3 lepton.
      4. pt for partons in 3 lepton.
      5. deltaR in 3 lepton after acceptance cuts.
    7. 4 lepton+Jets
      1. eta for leptons in 4 lepton.
      2. pt for leptons in 4 lepton.
      3. eta for partons in 4 lepton.
      4. pt for partons in 4 lepton.
      5. deltaR in 4 lepton after acceptance cuts.

2011-09-16 - Generation problem solve

  1. After a couple of days and headaches, what Johan said solve the day!!. smile MG5 doesn't define taus by default!!!. Defining them, the BR problem is solve.
    1. In MG5:
      1. define l+ = l+ ta+
      2. define l- = l- ta-
      3. generate p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
      4. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
      5. add process p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j) @3
      6. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl) @4
      7. output ttbar_v12
      8. launch (nevents = 5000, ebeam1/2 = 3500, cut_decays = F)
      9. proc_card_mg5.dat and ttbar_tau_run_card.dat
    2. Hadronized as cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:ttbar_v12.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 5000
    3. Run process in condor.
    4. ttbarAnalyzer Code:
      def getBR(infile,outfile):
      #######################################
              events = Events (infile)
              handleGen = Handle ("vector<reco::GenParticle>")
              label = "genParticles"
      
              hBR = TH1F('h1','Branching Ratio',3,0,3)
              hBR.SetStats(0)
      
              nevent = 0
              hadronically = 0
              dilepton = 0
              leptonJets = 0
              for event in events:
      #               print event
                      nevent += 1
                      leptons = 0
                      leptonsTau = 0
                      event.getByLabel (label, handleGen)
                      gens = handleGen.product()
      
                      for p in gens:
                              if p.mother():
                                      #print e, p.status(), p.pdgId(), mother
                                      if p.status() == 1:
                                              if ((abs(p.pdgId()) == 11) or (abs(p.pdgId()) == 13) or (abs(p.pdgId()) == 15)) and ((abs(p.mother().pdgId()) == 11) or (abs(p.mother().pdgId()) == 13) or (abs(p.mother().pdgId()) == 15)):
                                                      leptonsTau +=1
                                                      if not (abs(p.pdgId()) == 15) : leptons += 1
      #               print e, leptons
                      if leptonsTau == 0 : hadronically += 1
                      elif leptons == 1: leptonJets += 1
                      elif leptons == 2: dilepton += 1
      
              #print e, hadronically, dilepton, leptonJets
              hBR.Fill('hadronically', ((hadronically*100)/nevent))
              hBR.Fill('lepton + Jets', (leptonJets*100/nevent))
              hBR.Fill('dilepton', (dilepton*100/nevent))
              hBR.Draw()
              c1.SetGridx()
              c1.SetGridy()
              c1.SaveAs(outfile+'_BR.png')
            
  2. Generate 10k Fourtop process with GH 500MeV:
    1. In MG5:
      1. import model FourTop
      2. define l+ = l+ ta+
      3. define l- = l- ta-
      4. generate p p > GH GH, GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
      5. add process p p > GH GH, GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
      6. add process p p > GH GH, GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j) @3
      7. add process p p > GH GH, GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl) @4
      8. output fourtop_v12
      9. launch (nevents = 10000)
      10. proc_card_mg5.dat and ttbar_tau_run_card.dat
    2. Hadronized as cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:fourtop500_n10000.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 10000
    3. Run process in condor.
    4. Result files: fourtop500_n10000.root
  3. Same process for GH 1 TeV:
    1. proc_card_mg5.dat and ttbar_tau_run_card.dat
    2. Result files: fourtop1000_n10000.root

2011-09-15 - Generation problems

  1. Tried with a different Hadronizer:
    1. cmsDriver.py Configuration/GenProduction/python/Z2onLHE_noQEDfsrElectrons_7TeV_pythia6_cff.py --filein=file:ttbar_v15.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 100
    2. Result files: Z2onLHE_noQEDfsrElectrons_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU py and root
    3. Z2onLHE_noQEDfsrElectrons_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU_nleptons.png
    • Same results.
  2. In BU Meeting, they suggested me to use status = 3 with my leptons, because there is some problems with pythia. I obtain this two plots:

2011-09-14 - MG5 problems

  1. Johan Alwall (MG Team) confirm that, at least, my generation process in MG5 is correct. He said: Using your process generation, I get results in perfect agreement with BR(t->b e/mu/ta nu)= 0.33, BR(t->bjj)=0.67. Note that you need to run without cuts on the decay products (e.g. by setting cut_decays=F) to get the expected cross section ratios.
  2. Generate a sample with cut_decays=F, as before, in MG5.
    1. Result files: ttbar_v14 (1000 events)
    2. Number of leptons in ttbar_V14
    • Nothing seems to change.
  3. Trying a different way to generate ttbar in MG5.
    1. generate p p > t~ t QED=0, (t~ > W- b~, W- > l- vl~), t > j j b @2, t~ > j j b~ @3, (t > W+ b, W+ > l+ vl) @4
    2. output ttbar_v15 (100 events)
    3. Number of leptons in ttbar_V15
    • Nothing seems to change.

2011-09-13 - MG5 problems

  1. After analysed 5k ttbar events with my python script, I realized that something weird is going on. BR of this sample (ttbar_v10.root) drew this plot, where it seems clearly that it is wrong. As Francisco said The BF of W to hadrons is ~0.67 and for lv is 0.11 Therefore in the case of all jets the BF(ttbar to hadrons) is 0.67*0.67= 45% and in the case of 4tops is 0.67^4 = 20%. In the case of lepton+jets (only and only one muon or electron), the BF(ttbar to muon+jets) = 0.67*0.11+0.11*0.67 = ~15% Similarly for the case of electron+jets BF(ttbar to e+jets)= 15%.
  2. Started to look up the sample in more detail:
    1. Number of leptons in ttbar_V10 and Percentage of leptons in ttbar_V10 seems that the generated process prefers one or two leptons instead of zero leptons.
    2. Just in case, I drew number of leptons in ttbar_V10 with state 3. But It gave me not a clue of what is going on.
  3. Generate 1k events for ttbar in MG5 as:
    1. Change environment to CMSSW_4_2_6 in /uscms_data/d3/algomez/work/
    2. In MG5
      1. generate p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
      2. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
      3. output ttbar_v11
      4. launch (nevents = 5000)
    3. Hadronized as cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:ttbar_v11.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 1000
    4. Result plots: Branching ratio in ttbar_V11, Number of leptons in ttbar_V11 and Percentage of leptons in ttbar_V11 doesn't seem to fix the problem.
  4. Ask about the generation steps in HN and here. There they suggested me to define tau and tau~ in MG5 before generated the process.
    1. In MG5
      1. define l+ = l+ tau~
      2. define l- = l- tau
      3. generate p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
      4. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
      5. add process p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j) @3
      6. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl) @4
      7. output ttbar_v12
      8. launch (nevents = 5000)
    2. Hadronized as cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:ttbar_v12.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 100
    3. Result plots: Branching ratio in ttbar_V12, Number of leptons in ttbar_V12 and Percentage of leptons in ttbar_V12 doesn't seem to fix the problem.
  5. Tried generating a more simple process:
    1. In MG%:
      1. generate p p > t t~
      2. output ttbar_13
      3. launch (nevents = 5000)
    2. Hadronized as cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:ttbar_v12.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 500
    3. Result plots: Branching ratio in ttbar_V13, Number of leptons in ttbar_V13 doesn't seem to fix the problem

2011-09-11 - Generate 5k ttbar events and RECO in Condor

  1. Generate 5k events for ttbar in MG5 with previous proceding
    • Results: ttbar_v10.lhe
  2. Send to Condor for 5k RECO events.
    • Results: ttbar_v10.root
  3. Move files to /uscms_data/d3/algomez/files/ttbar

2011-09-10 - Submit RECO process to Condor

  1. Prepare files to submit job to Condor because I need more events for analysis.
    1. Files to submit:
  2. Generate 5k events for GH 500 GeV and 1 TeV.
    1. Results: fourtop500.root and fourtop1000.root.
    2. Organize everything in /uscms_data/d3/algomez/files/fourtop/

2011-09-08 - Fixing some errors (deltaR, Branching Ratio, Acceptance Efficiency)

  1. Correct deltaR code:
    def deltaR(infile,outfile):
    ######################################
            events = Events (infile) 
            handleGen = Handle ("vector<reco::GenParticle>")
            label = "genParticles"
    
            hdeltaR = TH1F('deltaR','deltaR',100,0,2)
    
            for event in events:
    #               print event
                    count2 = 0 
                    count = 0
                    L1 = {}
                    Q1 = {}
                    dR = {}
                    Min = {}
                    event.getByLabel (label, handleGen)
                    gens = handleGen.product()
                     
                    for p in gens:
    #                       print 'particle', p.pdgId(), 'status', p.status()
                            if p.status() == 1:
                                    if (p.pdgId() == 11) or (p.pdgId() == 13):
    #                                       print 'name', p.pdgId(), p.status(), p.p4()
                                            count += 1
                                            L1[count] = p.p4()
                            if p.status() == 3:
                                    if (abs(p.pdgId()) in range(1,7)) or (p.pdgId() == 21) :
    #                                       print 'name', p.pdgId(), p.status()
                                            count2 += 1
                                            Q1[count2] = p.p4()
    #               print L1
    #               print Q1
                    for lepton, pLepton in L1.iteritems():
                            dR[lepton] = {}
                            for quark, pQuark in Q1.iteritems():
                                    deltaR = Math.VectorUtil.DeltaR(pLepton, pQuark)
                                    dR[lepton][quark] = deltaR
    #                               print lepton, quark, deltaR
    #               print dR
                    for l in dR:
                            myMin = min(dR[l].values())
                    for i in Min:
                            myMin2 = min(Min.values())
    #               print myMin2
                    hdeltaR.Fill(myMin2)
    #
            c1 = TCanvas('c1','c1',800,600)
            hdeltaR.Draw()
            c1.SetLogy()
            c1.SaveAs(outfile+'_deltaR.png')
          
  2. Correct branching ratio:
    def getBR(infile,outfile):
    #######################################
            events = Events (infile)
            handleGen = Handle ("vector<reco::GenParticle>")
            label = "genParticles"
    
            hBR = TH1F('h1','Branching Ratio',3,0,3)
    #       hBR.SetYTitle('BR')
            hBR.SetStats(0)
    
            e = 0
            hadronically = 0
            dilepton = 0
            leptonJets = 0
            for event in events:
    #               print event
                    e += 1
                    ntotal = 0
                    leptons = 0
                    quarks = 0
                    others = 0
                    event.getByLabel (label, handleGen)
                    gens = handleGen.product()
    
                    for p in gens:
                            if p.status() == 1:
    #                               print e, p.status(), p.pdgId()
                                    if (p.pdgId() == 11) or (p.pdgId() == 13):
                                            leptons +=1
                            if p.status() == 3:
                            #       print e, p.status(), p.pdgId()
                                    if abs(p.pdgId()) in range(1,7):
                                            quarks += 1
    #               print e, ntotal, quarks, leptons
                    if quarks > 16 :
                            hadronically += 1
                    elif leptons == 4:
                            dilepton += 1
                    elif leptons == 2:
                            leptonJets += 1
            hBR.Fill('hadronically', ((hadronically*100)/e))
            hBR.Fill('dilepton', (dilepton*100/e))
            hBR.Fill('lepton + Jets', (leptonJets*100/e))
            hBR.Draw()
            c1.SetGridx()
            c1.SetGridy()
            c1.SaveAs(outfile+'_ttbarBR.png')
       
  3. Acceptance Efficiency code:
    #######################################
    def accepEfficiency(infile,outfile):
    #######################################
            events = Events (infile)
            handleGen = Handle ("vector<reco::GenParticle>")
            label = "genParticles"
    
            hAE = TH1F('h1','Acceptance Efficiency',4,0,4)
            hAE.SetYTitle('percentage')
            hAE.SetStats(0)
    
            e = 0
            accepMuons = 0
            accepElectrons = 0
            accepJets = 0
            total = 0
            for event in events:
                    e += 1
                    muons = 0
                    electrons = 0
                    jets = 0
    #               print event
                    event.getByLabel (label,handleGen)
                    gens = handleGen.product()
    
                    for p in gens:
                            if p.status() == 1:
                                    if (abs(p.pdgId()) == 13) and (abs(p.eta()) < 2.1):
                                            muons += 1
                                    elif (abs(p.pdgId()) == 11) and (abs(p.eta()) < 2.5):
                                            electrons += 1
                            if p.status() == 3:
                                    if (abs(p.pdgId()) in range (1,7)) and (abs(p.eta()) < 2.4):
                                            jets += 1
                    if muons: accepMuons += 1
                    if electrons: accepElectrons += 1
                    if jets: accepJets += 1
                    if (muons and electrons and jets): total +=1
    
            hAE.Fill('|eta| muons < 2.1', ((accepMuons*100)/e))
            hAE.Fill('|eta| electrons < 2.5', (accepElectrons*100/e))
            hAE.Fill('|eta| jets < 2.4', (accepJets*100/e))
            hAE.Fill('combined', (total*100/e))
            hAE.Draw()
            c1.SetGridx()
            c1.SetGridy()
            c1.SaveAs(outfile+'_accepEfficiency.png')
       

2011-09-06 - Fixing some errors

  1. I was doing some errors that I have to fix and I need to reconstruct more events with cmsDriver.
  2. For fourtop_v10.lhe, ran the same process as before.
    • Number of events: 949/950 (finished for root cmslpc)
    • New files: Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU_v12
  3. For fourtop_v11.lhe, ran the same process as before.
    • Number of events: 767/767 (finished for root cmslpc)
    • New files: Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU_v13
  4. Trying to fix deltaR, but it doesn't work. frown

2011-09-01 - 4t DeltaR

  1. Finally it works!!. My mistake was not to match leptons and partons in every event instead of the whole process.
    1. Code plothistos.py
    2. fourtop_v10_deltaR.png for GH mass = 500 GeV and lepton status 3.
    3. fourtop_v11_deltaR.png for GH mass = 1000 GeV and lepton status 3.
    4. fourtop_v10_deltaR_l1.png for GH mass = 500 GeV and lepton status 1.
    5. fourtop_v11_deltaR_l1.png for GH mass = 1000 GeV and lepton status 1.

2011-08-31 - Generating different events in MG, Hadronizing and some analysis.

  1. After two days fighting with MG in cmslpc, finally I can run MG5 there. My mistake, not to setup the environment (cmsenv) BEFORE trying to run MG5. :P
    1. cd work
    2. wget http://launchpad.net/madgraph5/trunk/1.3.0/+download/MadGraph5_v1.3.11.tar.gz
    3. tar xvf MadGraph5_v1.3.11.tar.gz
    4. cd MadGraph5_v1_3_11/
    5. Copy FourTop model to folder: models/
    6. ./bin/mg5 (to run MG5)
  2. Generate 5k events for HG mass for 500 GeV.
    1. import model FourTop
    2. generate p p > GH GH, GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
    3. add process p p >  GH GH, GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @1
    4. add process p p >  GH GH, GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j) @2
    5. add process p p >  GH GH, GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl) @3
    6. output fourtop_v10
    7. launch (removed all cuts in run_card.dat)
      1. Madevent fourtop_v10.tar.gz
      2. run_card_fourtop_v10.dat
    8. Move lhe file to my fourtop work area: mv Events/run_01_unweighted_events.lhe.gz ../fourtop/CMSSW_4_2_5/src/. Unzip (gunzip run_01_unweighted_events.lhe.gz) and change name to fourtop_v10.lhe
    9. In /uscms/home/algomez/work/fourtop/CMSSW_4_2_5/src/: cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:fourtop_v10.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 100
      1. 100 over 100 events
      2. Results: Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU_v10 python and root files.
    10. Run my python script for some plots:
      1. fourtop_v10_top_eta.png
      2. fourtop_v10_top_mass.png
      3. fourtop_v10_top_nparticles.png
      4. fourtop_v10_top_pt.png
      5. fourtop_v10_BR.png
  3. Same previous process but changing HG mass to 1 TeV. Same results files.
    1. Results:
      1. fourtop_v11_top_eta.png
      2. fourtop_v11_top_mass.png
      3. fourtop_v11_top_nparticles.png
      4. fourtop_v11_top_pt.png
      5. fourtop_v11_BR.png

2011-08-28 - Removing Cuts and Increase mass in MG

  1. We don't need to jet matched our process. As Francisco and Bogdan suggested: In this case, we don't actually need matching. We are only considering the LO (leading order) processes. It is truth that with tops, the next LO process could be as high as 40%. That is why in the LHC we use MG to generate SM ttbar+jets, where matching of the additional jets is required in the simulation.
  2. Reorganize my files:
    1. Create folders ttbar and fourtop in my nobackup area (/uscms_data/d1/algomez)
    2. Move all lhe, py and root files to nobackuparea, organize in each folder.
  3. Next steps: Remove any default cut in MG, like pT or deltaRij in the partons. Or make sure that the default cuts are much smaller than the ones we are going to apply in the analysis. We don't know yet about this so better make zero the pT cut and relax to large values the deltaR cuts.
    1. Generate same process in MG:
      1. import model FourTop
      2. generate p p >  GH GH, GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
      3. add process p p >  GH GH, GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
      4. add process p p >  GH GH, GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j) @3
      5. add process p p >  GH GH, GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl) @4
      6. output ttbar_v6 (in my laptop)
      7. launch
    2. Hadronizer:
      1. cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:fourtop_v6.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 100
        • 100 over 100 events
        • Results: Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU_v6 python and root files.
  4. Next steps: The model he (Bodgan) sent is for a scalar mass of 500 GeV. We need to test various mass points. We can start with a couple of mass points: GH=500 GeV, and GH=1000 GeV.
    1. Generate same precious process in MG.
    2. Hadronizer:
      1. cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:fourtop_v7.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 100
        • 100 over 100 events
        • Results: Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU_v7 python and root files.

2011-08-26 - Hadronizer

  1. Same previous process but only with 100 events. cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:fourtop_v4.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 100
    • I got 100 over 100 events and the process takes like a 13 min long.
    1. Results: Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU python and root files.
  2. Tried to run hadronizer again. cmsDriver.py Configuration/GenProduction/python/Hadronizer_MgmMatchTuneZ2_7TeV_madgraph_tauola_cff.py --filein=file:ttbar_v5.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 100
    • Same error message as before. (The MGParamCMS header does not specify the jet matching parameter "etaclmax", but it is requested by the CMSSW configuration.)
    1. Log File

2011-08-25 - Matching Efficiency 4t

  1. Asked Elim Cheung about mgPostProc. She answer:
The mgPostProc.py is a script given by Aram. It changes the format of an input LHE file. For example, we found that in our LHE files, the neutrino mass, which is supposed to be zero according to the SM, appears to be a very small number (something times e^-33). This error occurs because we used an older version of MadGraph. But I am not sure if it would help...
    • Maybe, it is not the right way to solve the problem.
  1. Generate 5k events for 4t. Change energy per beam to 3.5 TeV.
    1. Generation in MG fourtop_v4.lhe (in my laptop)
    2. Check for number of events generated in MG: 5000 events.
    3. Upload fourtop_v4.lhe to my fnal area.
    4. Run cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:fourtop_v4.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 5000
      • It seems like it works. Maybe I had some disk space problems. But the process takes a lot of time and perhaps I should split it in shorter chunks.
      • The process stops after 956 events, because of some cmslpc policy about large computing process. But It seems as it works!!!

2011-08-24 - Matching Efficiency ttbar/4t

  1. Task: Generate 5k events ttbar events. Make tops decay in MG. Pass the LHE file to the FastSim using this hadronizer: Hadronizer_MgmMatchTuneZ2_7TeV_madgraph_tauola_cff.py instead of the Z2onLHE_7TeV_pythia6_cff.py. Then generate 5k events 4-tops events. Make the tops decay in MG. Use the same hadronizer and check the matching efficiency. (Francisco's instructions).

TTbar

    1. Generate 5k events:
      1. generate p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
      2. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
      3. add process p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j) @3
      4. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl) @4
      5. output ttbar_v5 (in my laptop)
      6. launch (nevents = 5000)
    2. Upload ttbar_v5.lhe to /uscms/home/algomez/work/ttbar/CMSSW_4_2_5/src
    3. cmsDriver.py Configuration/GenProduction/python/Hadronizer_MgmMatchTuneZ2_7TeV_madgraph_tauola_cff.py --filein=file:ttbar_v5.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO
      1. Error Message:
        GEN,FASTSIM,HLT:GRun
        We have determined that this is simulation (if not, rerun cmsDriver.py with --data)
        GEN
        FASTSIM
        HLT:GRun
        Starting  cmsRun Hadronizer_MgmMatchTuneZ2_7TeV_madgraph_tauola_cff_py_GEN_FASTSIM_HLT_PU.py
             MSTU(12)       changed from              0 to          12345
         FastSimulation/PileUpProducer -> distribution from configuration file 
         FastSimulation/PileUpProducer -> An histogram will be created with 25 bins in the range (0,25).
         FastSimulation/PileUpProducer -> MinBias events taken from MinBias7TeV_001.root et al.,  with an average number of events of 7.35814
        Number of Cells in file is 0 - probably bad file format
        
        Process and cross-section statistics
        ------------------------------------
        Process   events   tried   xsec [pb]      accepted [%]
        0   0   0   n/a         n/a
        2   0   0   n/a         n/a
        3   0   0   n/a         n/a
        4   0   0   n/a         n/a
        Total   0   0   0.000e+00 +/- 0.000e+00   nan
             MSTU(12)       changed from              0 to          12345
        1****************** PYINIT: initialization of PYTHIA routines *****************
         ==== PYTHIA WILL USE LHAPDF ====
         *************************************
         *       LHAPDF Version 5.6.0        *
         *************************************
        
         >>>>>> PDF description: <<<<<<
         CTEQ6L1 - LO with LO alpha_s                                    
         Reference:                                                      
         J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky,       
         W.K. Tung                                                       
         hep-ph/0201195                                                  
         >>>>>>                   <<<<<<
        
         Parametrization: CTEQ6           
        
         ==============================================
         PDFset name /uscmst1/prod/sw/cms/slc5_amd64_gcc434/external/lhapdf/5.6.0-cms4/share/lhapdf/P
         with          1 members
         ====  initialized. ===========================
         Strong coupling at Mz for PDF is:  0.12978
        
         ==============================================================================
         I                                                                            I
         I         PYTHIA will be initialized for p+ on p+ user configuration         I
         I            with   7000.000 GeV on   7000.000 GeV beam energies             I
         I                                                                            I
         I           corresponding to  14000.000 GeV center-of-mass energy            I
         I                                                                            I
         ==============================================================================
        
         ******** PYMAXI: summary of differential cross-section maximum search ********
        
                   ==========================================================
                   I                                      I                 I
                   I  ISUB  Subprocess name               I  Maximum value  I
                   I                                      I                 I
                   ==========================================================
                   I                                      I                 I
                   I    4   User process 0                I    1.9032D-08   I
                   I    6   User process 2                I    1.2168D-07   I
                   I    7   User process 3                I    4.8076D-08   I
                   I    9   User process 4                I    4.7397D-08   I
                   I   96   Semihard QCD 2 -> 2           I    8.9669D+03   I
                   I                                      I                 I
                   ==========================================================
        
         ****** PYMULT: initialization of multiple interactions for MSTP(82) = 4 ******
                pT0 = 3.22 GeV gives sigma(parton-parton) = 4.89D+02 mb: accepted
        
         ****** PYMIGN: initialization of multiple interactions for MSTP(82) = 4 ******
                pT0 = 3.22 GeV gives sigma(parton-parton) = 1.91D+02 mb: accepted
        
         ********************** PYINIT: initialization completed **********************
         ****************** PYSLHA v1.14: SUSY/BSM SPECTRUM INTERFACE *****************
         * (PYSLHA:) Last Change 10 Jun 2010 - P.Z. Skands
         * (PYSLHA:) Ignoring MASS  entry for KF =        5 (SLHA read-in not allowed)
         * (PYSLHA:) Reading  MASS  entry for KF =        6, pole mass =     173.000
         * (PYSLHA:) Ignoring MASS  entry for KF =       15 (SLHA read-in not allowed)
         * (PYSLHA:) Reading  MASS  entry for KF =       23, pole mass =      91.188
         * (PYSLHA:) Reading  MASS  entry for KF =       25, pole mass =     120.000
         * (PYSLHA:) Ignoring MASS  entry for KF =       12 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =       14 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =       16 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =        2 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =        4 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =        1 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =        3 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =       22 (SLHA read-in not allowed)
         * (PYSLHA:) Reading  MASS  entry for KF =       24, pole mass =      80.419
         * (PYSLHA:) Ignoring MASS  entry for KF =       21 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =       11 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring MASS  entry for KF =       13 (SLHA read-in not allowed)
         * (PYSLHA:) Allowing DECAY table for KF =        6 but this is NOT recommended.
         * (PYSLHA:) Allowing DECAY table for KF =       23 but this is NOT recommended.
         * (PYSLHA:) Allowing DECAY table for KF =       24 but this is NOT recommended.
         * (PYSLHA:) Ignoring DECAY table for KF =       12 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =       14 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =       16 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =        2 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =        4 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =        1 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =        3 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =        5 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =       22 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =       21 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =       11 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =       13 (SLHA read-in not allowed)
         * (PYSLHA:) Ignoring DECAY table for KF =       15 (SLHA read-in not allowed)
        1********* PYSTAT:  Statistics on Number of Events and Cross-sections *********
        
         ==============================================================================
         I                                  I                            I            I
         I            Subprocess            I      Number of points      I    Sigma   I
         I                                  I                            I            I
         I----------------------------------I----------------------------I    (mb)    I
         I                                  I                            I            I
         I N:o Type                         I    Generated         Tried I            I
         I                                  I                            I            I
         ==============================================================================
         I                                  I                            I            I
         I   0 All included subprocesses    I            0             0 I  0.000D+00 I
         I   4 User process 0               I            0             0 I  0.000D+00 I
         I   6 User process 2               I            0             0 I  0.000D+00 I
         I   7 User process 3               I            0             0 I  0.000D+00 I
         I   9 User process 4               I            0             0 I  0.000D+00 I
         I                                  I                            I            I
         ==============================================================================
        
         ********* Total number of errors, excluding junctions =        0 *************
         ********* Total number of errors, including junctions =        0 *************
         ********* Total number of warnings =                           0 *************
         ********* Fraction of events that fail fragmentation cuts =  1.00000 *********
        
        
        **********************************************************************
            PRETAUOLA_05 - summary: 
            TAU+ force to decay into:  all decay modes              
            TAU- force to decay into:  all decay modes              
            Nev =         0 N(taus) =         0 ADDITIONAL BR =        NaN
        
         ***************************************************************************
         *                         *****TAUOLA UNIVERSAL INTERFACE: ******         *
         *                         *****VERSION 1.21, September2005 ******         *
         *                         **AUTHORS: P. Golonka, B. Kersevan, ***         *
         *                         **T. Pierzchala, E. Richter-Was, ******         *
         *                         ****** Z. Was, M. Worek ***************         *
         *                         **USEFUL DISCUSSIONS, IN PARTICULAR ***         *
         *                         *WITH C. Biscarat and S. Slabospitsky**         *
         *                         ****** are warmly acknowledged ********         *
         *                         ****** END OF MODULE OPERATION ********         *
         ***************************************************************************
        
        
        Process and cross-section statistics
        ------------------------------------
        Process   events   tried   xsec [pb]      accepted [%]
        0   0   0   n/a         n/a
        2   0   0   n/a         n/a
        3   0   0   n/a         n/a
        4   0   0   n/a         n/a
        Total   0   0   0.000e+00 +/- 0.000e+00   nan
        %MSG-s CMSException:  PostEndRun 24-Aug-2011 23:23:05 CDT PostEndRun
        cms::Exception caught in cmsRun
        ---- Generator|PartonShowerVeto BEGIN
        The MGParamCMS header does not specify the jet matching parameter "qcut", but it is requested by the CMSSW configuration.
        cms::Exception going through module Pythia6HadronizerFilter/generator run: 1 luminosityBlock: 1
        ProcessingStopped
        Exception going through path generation_step
        EventProcessingStopped
        an exception occurred during current event processing
        cms::Exception caught in EventProcessor and rethrown
        ---- Generator|PartonShowerVeto END
        
        
        %MSG
        
        =============================================
        
        MessageLogger Summary
        
         type     category        sev    module        subroutine        count    total
         ---- -------------------- -- ---------------- ----------------  -----    -----
            1 Generator            -w AfterSource                            1        1
            2 LHEInterface         -w AfterSource                            1        1
            3 CMSException         -s PostEndRun                             1        1
        
         type    category    Examples: run/evt        run/evt          run/evt
         ---- -------------------- ---------------- ---------------- ----------------
            1 Generator            BeforeEvents                      
            2 LHEInterface         BeforeEvents                      
            3 CMSException         PostEndRun                        
        
        Severity    # Occurrences   Total Occurrences
        --------    -------------   -----------------
        Warning                 2                   2
        System                  1                   1
                 
      • Main Problem: The MGParamCMS header does not specify the jet matching parameter "qcut", but it is requested by the CMSSW configuration.
    4. Tried with different python files:
      1. =cmsDriver.py Configuration/GenProduction/python/Hadronizer_MgmMatchTuneZ2_7TeV_madgraph_tauola_cff_topmass175_5_py_GEN_FASTSIM_HLT_PU.py --filein=file:ttbar_v5.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO =
      2. cmsDriver.py Configuration/GenProduction/python/Hadronizer_MgmMatchTuneZ2_7TeV_madgraph_cff_py_GEN_FASTSIM_HLT_PU.py --filein=file:ttbar_v5.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO
      3. Same error message.
    5. Looking in internet for an answer, I found this: http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/elim/src/mgPostProc.py?revision=1.1&view=markup that seems like it solve the problem.
      1. python mgPostProc.py -o ttbar_v5.lhe -e 1 -q 1 run_01_unweighted_events.lhe (in my laptop, run_01_.... is the name of the process.)

4t

    1. Generate 5k events:
      1. generate p p > GH GH, GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
      2. add process p p > GH GH, GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
      3. add process p p > GH GH, GH > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j) @3
      4. add process p p > GH GH, GH > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl) @4
      5. output fourtop_v3 (in my laptop)
      6. launch (nevents = 5000)
    2. Upload fourtop_v3.lhe to /uscms/home/algomez/work/fourtop/CMSSW_4_2_5/src

DeltaR

  1. Francisco suggested: for partons you should use status==1. At the generator level, you can have only partons. Pythia takes these partons through the showering and hadronization models to simulate real life collisions.
  2. Changed deltaR function in plothistos.py:
    1. Change in if abs(p.pdgId()) in range(11,14): if p.status() == 1: shows a root error. Apparently leptons in status 1 don't have p4() (Lorentz Vector).
    2. There is no quark or gluon at status==1. Perhaps I should look for hadrons instead of quarks and gluons.
    3. Change in if (abs(p.pdgId()) in range(1,7)) or (p.pdgId() = 21) :  if p.status() = 3 shows a root error. Apparently gluons don't have p4().
    4. It seems like the code runs only with, for both leptons and quarks, status=3.

2011-08-23 - DeltaR, Branching Ratio and Fist 4t generation

  1. After some problems with the code, this are the results:
    1. Code: plothistos.py.
    2. Branching Ratio: ttbarBR_v1.png
      • It seems like always are 16 particles before parton showering.
      • With 100 events: 84 decay hadronically, 16 dilepton and 0 lepton+jet.
      • To prove this results I generate again the lhe file and run the same process again. I got the same results.
        1. Branching Ratio: ttbarBR_v1_1.png
  2. To avoid any statistical problem, I repeated the process again for 1000 events. The plots are:
    1. topmass_v4.png
    2. toppt_v4.png
    3. topeta_v4.png
    4. topnparticles_v4.png
    5. With lepton before parton showering deltaRTotal_v4.png
    6. Branching Ratio: ttbarBR_v4.png
    • Same results. frown
  3. Generate again ttbar process as:
    1. generate p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
    2. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
    3. add process p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > j j) @3
    4. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > l+ vl) @4
    5. output ttbar_v4 (in my laptop)
    6. launch (nevents = 2000)
    7. Same process as before
    8. Results:
      1. topmass_v5.png
      2. toppt_v5.png
      3. topeta_v5.png
      4. topnparticles_v5.png
      5. With lepton before parton showering deltaRTotal_v5.png
      6. Branching Ratio: ttbarBR_v5.png

Four Top

  1. First Instructions from Bodgan Dobrescu and Francisco about four top:
    1. Bodgan create the model: FourTop.tgz
    2. MG5 in my laptop:
      1. import model FourTop
      2. generate p p > GH GH , GH > t t~
      3. output fourtop_v1
      4. launch
    3. Inside my fnal area:
      1. mkdir work/fourtop
      2. scramv1 project CMSSW CMSSW_4_2_5
      3. = cd CMSSW_4_2_5/src= ; cmsenv
    4. Upload the lhe file: fourtop.lhe to ~/work/fourtop/CMSSW_4_2_5/src/
      1. cvs co Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py
      2. scramv1 b -j4
      3. cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:fourtops_v1.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 10000
        • log file:
          GEN,FASTSIM,HLT:GRun
          We have determined that this is simulation (if not, rerun cmsDriver.py with --data)
          GEN
          FASTSIM
          HLT:GRun
          Starting  cmsRun Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU.py
               MSTU(12)       changed from              0 to          12345
           FastSimulation/PileUpProducer -> distribution from configuration file 
           FastSimulation/PileUpProducer -> An histogram will be created with 25 bins in the range (0,25).
           FastSimulation/PileUpProducer -> MinBias events taken from MinBias7TeV_001.root et al.,  with an average number of events of 7.35814
          %MSG-w Generator:  AfterSource 23-Aug-2011 23:04:18 CDT BeforeEvents
          Event weight not set to one for abs(IDWTUP) == 3
          %MSG
          %MSG-w LHEInterface:  AfterSource 23-Aug-2011 23:04:18 CDT BeforeEvents
          Event weight not set to one for abs(IDWTUP) == 3
          %MSG
          Number of Cells in file is 0 - probably bad file format
          1****************** PYINIT: initialization of PYTHIA routines *****************
           ==== PYTHIA WILL USE LHAPDF ====
           *************************************
           *       LHAPDF Version 5.6.0        *
           *************************************
          
           >>>>>> PDF description: <<<<<<
           CTEQ6L1 - LO with LO alpha_s                                    
           Reference:                                                      
           J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky,       
           W.K. Tung                                                       
           hep-ph/0201195                                                  
           >>>>>>                   <<<<<<
          
           Parametrization: CTEQ6           
          
           ==============================================
           PDFset name /uscmst1/prod/sw/cms/slc5_amd64_gcc434/external/lhapdf/5.6.0-cms4/share/lhapdf/P
           with          1 members
           ====  initialized. ===========================
           Strong coupling at Mz for PDF is:  0.12978
          
           ==============================================================================
           I                                                                            I
           I         PYTHIA will be initialized for p+ on p+ user configuration         I
           I            with   7000.000 GeV on   7000.000 GeV beam energies             I
           I                                                                            I
           I           corresponding to  14000.000 GeV center-of-mass energy            I
           I                                                                            I
           ==============================================================================
          
           ******** PYMAXI: summary of differential cross-section maximum search ********
          
                     ==========================================================
                     I                                      I                 I
                     I  ISUB  Subprocess name               I  Maximum value  I
                     I                                      I                 I
                     ==========================================================
                     I                                      I                 I
                     I    4   User process 0                I    4.5597D-12   I
                     I   96   Semihard QCD 2 -> 2           I    8.9668D+03   I
                     I                                      I                 I
                     ==========================================================
          
           ****** PYMULT: initialization of multiple interactions for MSTP(82) = 4 ******
                  pT0 = 3.22 GeV gives sigma(parton-parton) = 4.89D+02 mb: accepted
          
           ****** PYMIGN: initialization of multiple interactions for MSTP(82) = 4 ******
                  pT0 = 3.22 GeV gives sigma(parton-parton) = 1.92D+02 mb: accepted
          
           ********************** PYINIT: initialization completed **********************
           ****************** PYSLHA v1.14: SUSY/BSM SPECTRUM INTERFACE *****************
           * (PYSLHA:) Last Change 10 Jun 2010 - P.Z. Skands
           * (PYSLHA:) Reading  QNUMBERS    for KF =    10025
           * (PYSLHA:) Ignoring MASS  entry for KF =        1 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =        2 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =        3 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =        4 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =        5 (SLHA read-in not allowed)
           * (PYSLHA:) Reading  MASS  entry for KF =        6, pole mass =     172.000
           * (PYSLHA:) Ignoring MASS  entry for KF =       11 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =       13 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =       15 (SLHA read-in not allowed)
           * (PYSLHA:) Reading  MASS  entry for KF =       23, pole mass =      91.188
           * (PYSLHA:) Reading  MASS  entry for KF =       25, pole mass =     120.000
           * (PYSLHA:) Ignoring MASS  entry for KF =    10025 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =       12 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =       14 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =       16 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring MASS  entry for KF =       22 (SLHA read-in not allowed)
           * (PYSLHA:) Reading  MASS  entry for KF =       24, pole mass =      79.824
           * (PYSLHA:) Ignoring MASS  entry for KF =       21 (SLHA read-in not allowed)
           * (PYSLHA:) Allowing DECAY table for KF =        6 but this is NOT recommended.
           * (PYSLHA:) Allowing DECAY table for KF =       23 but this is NOT recommended.
           * (PYSLHA:) Allowing DECAY table for KF =       24 but this is NOT recommended.
           * (PYSLHA:) Ignoring DECAY table for KF =    10025 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =       12 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =       14 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =       16 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =       11 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =       13 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =       15 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =        2 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =        4 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =        1 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =        3 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =        5 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =       22 (SLHA read-in not allowed)
           * (PYSLHA:) Ignoring DECAY table for KF =       21 (SLHA read-in not allowed)
          Begin processing the 1st record. Run 1, Event 1, LumiSection 1 at 23-Aug-2011 23:04:35.381 CDT
          %MSG-w SiStripQuality:  TrackProducer:iterativeFirstTracksWithTriplets  23-Aug-2011 23:04:44 CDT Run: 1 Event: 1
          The cabling should always include the active feds in runInfo and possibly have some morethere are instead 59 feds only active in runInfo
          %MSG
          %MSG-w Generator:  source 23-Aug-2011 23:04:53 CDT PostProcessEvent
          Event weight not set to one for abs(IDWTUP) == 3
          %MSG
          %MSG-w LHEInterface:  source 23-Aug-2011 23:04:53 CDT PostProcessEvent
          Event weight not set to one for abs(IDWTUP) == 3
          %MSG
          Begin processing the 2nd record. Run 1, Event 2, LumiSection 1 at 23-Aug-2011 23:04:53.245 CDT
          .........................
          ........................
          Begin processing the 168th record. Run 1, Event 168, LumiSection 1 at 23-Aug-2011 23:20:38.845 CDT
          %MSG-w Generator:  source 23-Aug-2011 23:20:46 CDT PostProcessEvent
          Event weight not set to one for abs(IDWTUP) == 3
          %MSG
          %MSG-w LHEInterface:  source 23-Aug-2011 23:20:46 CDT PostProcessEvent
          Event weight not set to one for abs(IDWTUP) == 3
          %MSG
          1********* PYSTAT:  Statistics on Number of Events and Cross-sections *********
           ==============================================================================
           I                                  I                            I            I
           I            Subprocess            I      Number of points      I    Sigma   I
           I                                  I                            I            I
           I----------------------------------I----------------------------I    (mb)    I
           I                                  I                            I            I
           I N:o Type                         I    Generated         Tried I            I
           I                                  I                            I            I
           ==============================================================================
           I                                  I                            I            I
           I   0 All included subprocesses    I          168           168 I  4.560D-12 I
           I   4 User process 0               I          168           168 I  4.560D-12 I
           I                                  I                            I            I
           ==============================================================================
          
           ********* Total number of errors, excluding junctions =        0 *************
           ********* Total number of errors, including junctions =        0 *************
           ********* Total number of warnings =                           1 *************
           ********* Fraction of events that fail fragmentation cuts =  0.00000 *********
          
          
          Process and cross-section statistics
          ------------------------------------
          Process   events   tried   xsec [pb]      accepted [%]
          0   168   168   4.560e-03 +/- 7.159e-06   100.0
          Total   168   168   4.560e-03 +/- 7.159e-06   100.0
          
          =============================================
          
          MessageLogger Summary
          
           type     category        sev    module        subroutine        count    total
           ---- -------------------- -- ---------------- ----------------  -----    -----
              1 AdaptiveVertexFitter -w PFDisplacedVerte                       1        1
              2 Generator            -w AfterSource                            1        1
              3 Generator            -w source                               168      168
              4 HcalHitMaker         -w FamosProducer:fa                       1        1
              5 LHEInterface         -w AfterSource                            1        1
              6 LHEInterface         -w source                               168      168
              7 SiStripQuality       -w TrackProducer:it                       1        1
              8 TwoTrackMinimumDista -w ConversionProduc                       3        3
          
           type    category    Examples: run/evt        run/evt          run/evt
           ---- -------------------- ---------------- ---------------- ----------------
              1 AdaptiveVertexFitter 1/10                              
              2 Generator            BeforeEvents                      
              3 Generator            PostProcessEvent PostProcessEvent PostProcessEvent
              4 HcalHitMaker         1/111                             
              5 LHEInterface         BeforeEvents                      
              6 LHEInterface         PostProcessEvent PostProcessEvent PostProcessEvent
              7 SiStripQuality       1/1                               
              8 TwoTrackMinimumDistanceHelixHelix 1/26             1/124            1/136
          
          Severity    # Occurrences   Total Occurrences
          --------    -------------   -----------------
          Warning               344                 344
                   

2011-08-21 - DeltaR

  1. After a lot of problems with the code, I think I got what we want for deltaR. As I am not sure if leptons should be from final state (1) or before parton showering (3).
    1. Code
      events = Events (infile)
              handleGen = Handle ("vector<reco::GenParticle>")
              label = "genParticles"
      
              leptons = [11, 13, 15, -11, -13, -15]
              quarks = [1, 2, 3, 4, 5, 6, -1, -2, -3, -4, -5, -6]
      
              hdeltaR = TH1F('deltaR','deltaR',100,0,2)
      
              L = {}
              Q = {}
              dR = {}
              count = 0
              count2 = 0
              for event in events:
      #               print event
                      event.getByLabel (label, handleGen)
                      gens = handleGen.product()
      
                      for p in gens:
                              if p.pdgId() in leptons:
                                      if p.status() == 1:                     # for final particles
                                              count += 1
                                              L[count] = {}
                                              L[count][int(p.pdgId())] = p.p4()
                              if p.pdgId() in quarks:         # for particles before parton showering
                                      if p.status() == 3:
                                              count2 += 1
                                              Q[count2] = {}
                                              Q[count2][int(p.pdgId())] = p.p4()
                      for i in L.keys():
                              dR[i] = {}
                              for lepton, ivalue in L[i].iteritems():
                                      for j in Q.keys():
                                              dR[i][j] = {}
                                              for quark, jvalue in Q[j].iteritems():
                                                      deltaR = Math.VectorUtil.DeltaR(ivalue, jvalue)
                                                      if (0 < deltaR < 50):
                                                              dR[i][j] = deltaR
                                                      elif deltaR > 50:
                                                              dR[i][j] = 50
                                                      else:
                                                              dR[i][j] = 1000
                      for k in dR.keys():
                              for x in dR[k].keys():
                                      myMin = min(dR[k].values())
                              hdeltaR.Fill(myMin)
      
              c1 = TCanvas('c1','c1',800,600)
              hdeltaR.Draw()
              c1.SetLogy()
              c1.SaveAs('deltaRTotal.png')
               
    2. With lepton before parton showering deltaRTotal_v2.png
    3. With lepton as final state deltaRTotal_v3.png
    4. Preliminar code: plothistos.py.

2011-08-20 - DeltaR

  1. Finally the number of particles plot is fine. It seems like I am initializing the nparticles in a wrong place.
    1. Code:
                      for event in events:
                              #entry += 1
                              nparticles = 0
                              event.getByLabel (label, handleGen)
                              gens = handleGen.product()
      
                              for p in gens:
                                      if abs( int(p.pdgId()) ) == int(PDG[particle]):
                                              if p.status() == 3: 
                                                      nparticles += 1
                                                      hpt.Fill(float(p.pt()))
                                                      heta.Fill(float(p.eta()))
                                                      hmass1.Fill(float(p.mass()))
                              hnparticles.Fill(nparticles)
               
    2. top_nparticles_v3.png
  2. deltaR calculation was wrong. I am experiencing some problems with this.
  3. Branching Ratio is easier than expected. I have to do it again.

2011-08-18 - DeltaR and Branching Ratios

  1. Fix the number of particles plot. (stupid mistake smile )
    1. hnparticles_v2_2.png, but seems to me that it is wrong or perhaps I don't understand what I am doing.
  2. Check Pythia Particle Status:
    status meaning
    1 stable particle (handed over to Geant4 for detector simulation)
    2 particle after parton showering and ISR/FSR
    3 particle before parton showering and ISR/FSR (i.e. as it comes from matrix element calculation)
  3. Trying to plot the deltaR(lepton,parton) between lepton and the closest parton.
    1. Add a function in plothistos.py to draw deltaR.
      1. Code:
                events = Events (infile)
                handleGen = Handle ("vector<reco::GenParticle>")
                label = "genParticles"
        
                leptons = [11, 13, 15, -11, -13, -15]
                quarks = [1, 2, 3, 4, 5, 6]
        
                hdeltaR = TH1F('deltaR','deltaR',100,0,10)
        
                for event in events:
                        event.getByLabel (label, handleGen)
                        gens = handleGen.product()
                        for p in gens:
                                if p.status() == 1:               # for final particles
                                        if p.pdgId() in leptons:
                                                lepton = p.p4()   #LorentzVector
                                if p.status() == 3:  
                                        if p.pdgId() in quarks: # for particles before parton showering
                                                quark = p.p4()
                        deltaR = Math.VectorUtil.DeltaR(lepton,quark)
                        hdeltaR.Fill(deltaR)
                hdeltaR.Draw()
                c1.SaveAs('deltaR.png')
                 
        • Something weird is going on, in my opinion.
      2. deltaR.png
    2. I didn't iterate over all posible particles!!!.
      1. Code:
        #######################################
        def deltaR(infile):
        ######################################
           events = Events (infile)
           handleGen = Handle ("vector<reco::GenParticle>")
           label = "genParticles"
        
           leptons = [11, 13, 15, -11, -13, -15]
           quarks = [1, 2, 3, 4, 5, 6, -1, -2, -3, -4, -5, -6]
        
           hdeltaR = TH1F('deltaR','deltaR',100,0,10)
           hdeltaR2 = TH1F('deltaR < 0.5','deltaR < 0.5',100,0,0.5)
        
           L = {}
           Q = {}
           for event in events:
              event.getByLabel (label, handleGen)
              gens = handleGen.product()
        
              for p in gens:
                 if p.status() == 1:         # for final particles
                    if p.pdgId() in leptons:
                       L[p.pdgId()] = p.p4()
                 if p.status() == 3: 
                    if p.pdgId() in quarks:    # for particles before parton showering
                       Q[p.pdgId()] = p.p4()
        
              for i in L:
                 for j in Q:
                    deltaR = Math.VectorUtil.DeltaR(L[i],Q[j])
                    if deltaR != 0:
                       hdeltaR.Fill(deltaR)
                       if deltaR < 0.5:
                          hdeltaR2.Fill(deltaR)
        
           hdeltaR.Draw()
           c1.SaveAs('deltaRTotal.png')
           hdeltaR2.Draw()
           c1.SaveAs('deltaR05.png')
                 
      2. deltaRtotal.png
      3. deltaR05.png
  4. Looking for the branching ratio of the event.
    1. Based on http://hep.uchicago.edu/cdf/howto/grapher2.cpp
    2. Create another function on plothistos.py
    3. I didn't finish it today.
  5. Read paper:
    1. Cacciapaglia, G. et al., 2011. Four tops on the real projective plane at LHC. Arxiv preprint arXiv:1107.4616, p.25. Available at: http://arxiv.org/abs/1107.4616 [Accessed August 4, 2011].

2011-08-17 - FastSim and some plots

  1. I should check twice every instruction that I've read!!! frown Z2onLHE_7TeV_pythia6_cff.py is on Configuration/GenProduction/python
  2. Moved all previous files to /uscms/home/algomez/work/ttbar/CMSSW_4_2_5/src/Tauola
  3. Working with Z2onLHE_7TeV_pythia6_cff.py
    1. cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:ttbar.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM  --datatier GEN-SIM-DIGI-RECO -n 100 --no_exec
    2. cmsRun Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU.py
      1. Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU.root
    3. edmDumpEventContent Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU.root > Z2onLHE_7TeV_pythia6_cff_py_GEN_FASTSIM_HLT_PU_edmDumpEventContent.log
  4. Modify gen.py to get the following histograms:
    1. topmass.png
    2. toppt.png
    3. topeta.png
    4. topnparticles.png
    5. ttbarmass.png
  5. As I would have to plot this histos again and again, I transform gen.py from Francisco to plothistos.py.

2011-08-16 - First FastSim

  1. Following instructions of Francisco.
  2. Following the Fast SImulation SW Guide inside my fnal area.
    1. Working directory: /uscms/home/algomez/work/ttbar/
    2. Working release CMSSW_4_2_5 as suggested in the Guide.
    3. cvs co Configuration/Generator/python
    4. Francisco suggeted to run: cmsDriver.py Configuration/GenProduction/python/Z2onLHE_7TeV_pythia6_cff.py --filein=file:myLHE.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI-RECO -n 100 --no_exec
      • But apparently there is no Z2onLHE_7TeV_pythia6_cff.py file in Configuration/Generator/python.
    5. Work with the TTbar_Tauola_7TeV_cfi.py as suggested in the Guide.
    6. cmsDriver.py Configuration/Generator/python/ZTT_Tauola_7TeV_cfi.py --filein=file:ttbar.lhe -s GEN,FASTSIM,HLT:GRun --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI=RECO -n 100 --no_exec
      1. ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_PU python script
      2. cmsDriver.py parameters:
        parameter What it mean
        ZTT_Tauola_7TeV_cfi.py Remove all unnecessary EmptySource and ProductionFilterSequence declarations for new 7TeV RelVals.
        --filein=file for infile
        -s for sequence
        --pileup=FlatDist10_2011EarlyData_50ns for flat distribution up to 10, then decreasing smoothly (since 4_2_X)
        --conditions START42_V11::All As the 'MC' GT (including bad channels), but with realistic startup alignment and calibration constants. More Info here
        --beamspot Realistic7TeVCollision ---
        --eventcontent=RECOSIM ---
        --datatier GEN-SIM-DIGI=RECO Data Tiers
        -n 100 number of events
        --no_exec argument tells cmsDriver.py to write out the configuration file. If you do not specify this argument, cmsDriver.py will proceed to executing cmsRun
      3. More info here
    7. cmsRun ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_PU.py > ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_PU_v1.log - change the output root file to ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_PU_v1.root
        1. log file
        2. ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_PU_v1.root
        3. NuclearInteractionOutputFile.txt
        4. PileUpOutputFile.txt
    8. To produce Release Validation plots:
      1. cmsDriver.py Configuration/Generator/python/ZTT_Tauola_7TeV_cfi.py --filein=file:ttbar.lhe -s GEN,FASTSIM,HLT:GRun,VALIDATION --pileup=FlatDist10_2011EarlyData_50ns --conditions START42_V11::All --beamspot Realistic7TeVCollision --eventcontent=RECOSIM --datatier GEN-SIM-DIGI=RECO -n 100 --no_exec
      2. cmsRun ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_VALIDATION_PU.py > ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_VALIDATION_PU_v2.log
        1. ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_VALIDATION_PU.root
        2. NuclearInteractionOutputFile.txt
        3. PileUpOutputFile.txt
      3. cmsDriver.py step3 -s HARVESTING:validationHarvestingFS --harvesting AtJobEnd --conditions=auto:startup --filein_file:ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_VALIDATION_PU.root
      4. cmsRun step3_HARVESTING.py
        1. DQM_V0001_R000000001__Global__CMSSW_X_Y_Z__RECO.root
        • I didn't found any difference between "normal" and "validation" files. - I will continue to work with the "normal" file.

  1. Review Francisco's python script gen.py
    1. Change root input file to ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_PU_v1.root
    2. python gen.py
      1. Histogram without restrictions (i.e. all particles) topmass.png
      2. Trying to plot the mass for any other particle gives me an empty plot. Weird!! frown

  1. Just for fun smile Check my root file in Fireworks:
    1. cmsShow ZTT_Tauola_7TeV_cfi_py_GEN_FASTSIM_HLT_PU_v1.root
      1. ttbar-1_1_1_RhoPhi.png
      2. ttbar-1_1_1_Lego.png
      3. ttbar-1_1_1_3DTower.png

  1. Quick review of the files:
    1. /uscms_data/d2/vasurang/EJTerm3Photon/MakeGenFile.py
    2. /uscms_data/d2/vasurang/EJTerm3Photon/FastSimExample.py
  2. Reading again the Event Gen & Sim Chapter in Workbook for a better understanding of what I am doing.

smile Seems like a very productive day!!! smile

2011-08-15 - Generation MG5

  1. Cannot run Madgraph inside my fnal area. frown
    1. Error Message:
      • [algomez@cmslpc17 MadGraph5_v1_3_9]$ ./bin/mg5
      • MadGraph 5 works only with python 2.6 or later (but not python 3.X).
      • Please upgrate your version of python.
  2. Run Madgraph inside my laptop. smile
    1. In Madgraph:
      1. generate p p > t t~, (t~ > W- b~, W- > l- vl~), (t > W+ b, W+ > l+ vl)
      2. add process p p > t t~, (t~ > W- b~, W- > j j), (t > W+ b, W+ > j j) @2
      3. output ttbar_v2
      4. launch - changing nevents = 2000 in the run_card.dat
    2. Results:

2011-08-14 - First generation process in MG5

  1. Starts the real work smile
  2. Read some papers:
    1. DOBRESCU, B., KONG, K. & MAHBUBANI, R., 2008. Massive color-octet bosons and pairs of resonances at hadron colliders. Physics Letters B, 670(2), pp.119-123. Available at: http://arxiv.org/abs/0709.2378 [Accessed August 10, 2011].
    2. Gregoire, T., Katz, E. & Sanz, V., 2011. Toptet. Arxiv preprint arXiv:1101.1294. Available at: http://arxiv.org/abs/1101.1294 [Accessed August 4, 2011].
    3. Gerbush, M. et al., 2008. Color-octet scalars at the CERN LHC. Physical Review D, 77(9), p.15. Available at: http://arxiv.org/abs/0710.3133 [Accessed August 10, 2011].
  3. Generate the process in my laptop and my fnal area:
    1. Follow the Madgraph tutorial, generate a ttbar and get the lhe file.
    2. Generate the process: p p > t t~ , (t > w+ b , w+ > l+ vl) , (t~ > w- b~ , w- > l- vl~)
      1. Feynman Diagrams: gg> t tbar and uubar > t tbar
      2. LHE file
  4. Start to read about hadronize the HLE file with Pythia in CMSSW.

Previous Work - Madgraph

  1. Download Madgraph to my laptop.
  2. Read Madgrapdh documentation.
  3. Produce the first lhe file for ttbar in SM: run_01_unweighted_events.lhe (temporarily available).

-- AlejandroGomez - 24-Aug-2011

Edit | Attach | Watch | Print version | History: r61 < r60 < r59 < r58 < r57 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r61 - 2012-10-26 - AlejandroGomez
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2021 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