17 Aug 2010
I am trying to calculate the Higgs Cross Section using some fortran programs found at: http://wwweth.cern.ch/HiggsCrossSections/ The problem is their calculations only go to 600 GeV and I need to go higher.

higlu To get this to work I had to install lhapdf-5.8.4 which can be found at http://projects.hepforge.org/lhapdf/ then I had to change the LIBS line to

LIBS = /usr/local/lib/libLHAPDF.so
trying ldd run gave
linux-vdso.so.1 =>  (0x00007fff8b9ff000)
   libLHAPDF.so.0 => not found
   libgfortran.so.3 => /usr/lib/libgfortran.so.3 (0x00007f707c886000)
   libm.so.6 => /lib/libm.so.6 (0x00007f707c602000)
   libgcc_s.so.1 => /lib/libgcc_s.so.1 (0x00007f707c3eb000)
   libc.so.6 => /lib/libc.so.6 (0x00007f707c068000)
   /lib64/ld-linux-x86-64.so.2 (0x00007f707cb9f000)

so then I ran sudo ldconfig -v

then when running it I get

./run: error while loading shared libraries: libLHAPDF.so.0: cannot map zero-fill pages: Cannot allocate memory

Change in the FFLAGS line -march=pentiumpro to -march=native also delete the -fno-fixed-form flag

Then

 make 
and it installed.

10 Aug 2010
Here are the cross sections for HtoZZto2L2Nu for the various higgs masses that I calculated
105.00 0.00339363
110.00 0.00639475
115.00 0.0117958
120.00 0.0201309
125.00 0.0312436
130.00 0.0440731
135.00 0.0567509
140.00 0.0670555
145.00 0.0727671
150.00 0.0716391
155.00 0.0604
160.00 0.0324218
165.00 0.0158218
170.00 0.0157241
175.00 0.0201839
180.00 0.0352601
185.00 0.0826275
190.00 0.108895
195.00 0.117728
200.00 0.119991
205.00 0.11946
210.00 0.117407
215.00 0.114678
220.00 0.111584
225.00 0.108351
230.00 0.105194
235.00 0.102028
240.00 0.0989891
245.00 0.0960891
250.00 0.0933463
255.00 0.0907284
260.00 0.0882344
265.00 0.0859175
270.00 0.0837383
275.00 0.0816947
280.00 0.0797791
285.00 0.0780472
290.00 0.076418
295.00 0.0749407
300.00 0.073617
305.00 0.0724285
310.00 0.0714083
315.00 0.0705673
320.00 0.0698997
325.00 0.0694504
330.00 0.069276
335.00 0.0694732
340.00 0.0702583
345.00 0.0720113
350.00 0.0724766
355.00 0.0721054
360.00 0.0711599
365.00 0.0697845
370.00 0.0681678
375.00 0.0663271
380.00 0.0643836
385.00 0.0623637
390.00 0.0602918
395.00 0.0582335
400.00 0.0561582
405.00 0.0541391
410.00 0.0521285
415.00 0.0501823
420.00 0.0482859
425.00 0.0464449
430.00 0.0446559
435.00 0.0429323
440.00 0.0412644
445.00 0.0396485
450.00 0.0380957
455.00 0.0366011
460.00 0.035166
465.00 0.0337777
470.00 0.0324585
475.00 0.0311832
480.00 0.0299533
485.00 0.0287793
490.00 0.027659
495.00 0.0265789
500.00 0.0255403
505.00 0.0245404
510.00 0.023589
515.00 0.022673
520.00 0.0218036
525.00 0.0209589
530.00 0.0201569
535.00 0.0193856
540.00 0.0186449
545.00 0.0179353
550.00 0.0172526
555.00 0.0166042
560.00 0.0159765
565.00 0.0153809
570.00 0.0148028
575.00 0.014254
580.00 0.0137265
585.00 0.0132207
590.00 0.0127317
595.00 0.0122655
600.00 0.0118195

10 Aug 2010
I am looking on what cuts are done for soft muons. I need to talk to Zoltan and ask what they are doing with J/psi since they work with soft muons. This will be looked at as perhaps one of the cuts for the 2l2nu analysis. Also I will look at the 3D impact parameter. Nicola seems to think that this will remove the need to do the B tagging.

This is the website for the b-tagging from the cmssw workbook https://twiki.cern.ch/twiki/bin/view/CMS/WorkBookBTagging

Also at yesterdays meeting I am supposed to find which z+jet sample is the best one to use.

8 July 2010
Sorry I haven't been keeping a good log of what is going on. Anyways, here is the summary part of one of the 500 parts of the 150 GeV Higgs to ZZ to 2L2NU sample. I think this is the cross section and I am going to add these up for all the files and get a total amount. This is think is before the filters.

 
 ==============================================================================
 I                                                                            I
 I              PYTHIA will be initialized for a p on p collider              I
 I                  at   7000.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   96   Semihard QCD 2 -> 2           I    9.6420D+03   I
           I  102   g + g -> h0                   I    6.2908D-11   I
           I  123   f + f' -> f + f' + h0         I    4.7414D-11   I
           I  124   f + f' -> f" + f"' + h0       I    1.2203D-10   I
           I                                      I                 I
           ==========================================================

 ****** PYMULT: initialization of multiple interactions for MSTP(82) = 4 ******
        pT0 = 2.25 GeV gives sigma(parton-parton) = 5.52D+02 mb: accepted

 ****** PYMIGN: initialization of multiple interactions for MSTP(82) = 4 ******
        pT0 = 2.25 GeV gives sigma(parton-parton) = 2.50D+02 mb: accepted

 ********************** PYINIT: initialization completed **********************

  PRETAUOLA: additional BR after          1 tau+- decays =   1.000000


  PRETAUOLA: additional BR after       1001 tau+- decays =   1.000000


  PRETAUOLA: additional BR after       2001 tau+- decays =   1.000000


  PRETAUOLA: additional BR after       3001 tau+- decays =   1.000000


  PRETAUOLA: additional BR after       4001 tau+- decays =   1.000000


  PRETAUOLA: additional BR after       5001 tau+- decays =   1.000000

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

MessageLogger Summary

 type     category        sev    module        subroutine        count    total
 ---- -------------------- -- ---------------- ----------------  -----    -----
    1 CSCTFPtMethods::Pt2S -w CSCTFTrackProduc                      10*      10
    2 CaloSim              -w OscarProducer:g4                      89*      89
    3 SiStripQuality       -w SiStripRawToClus                       1*       1

* Some occurrences of this message were suppressed in all logs, due to limits.

 type    category    Examples: run/evt        run/evt          run/evt
 ---- -------------------- ---------------- ---------------- ----------------
    1 CSCTFPtMethods::Pt2StnChiSq() 1/163            1/163            1/4977
    2 CaloSim              1/90             1/116            1/5777
    3 SiStripQuality       Run: 1                            

Severity    # Occurrences   Total Occurrences
--------    -------------   -----------------
Warning               100                 100
 Message from PHCORK(MODCOR):: initialization
 MODOP=1 -- no corrections on event: DEFAULT
 Message from PHOTOS: IPHQRK(MODCOR):: (re)initialization
 MODOP=1 -- blocks emission from light quarks: DEFAULT
 Message from PHOTOS: IPHEKL(MODCOR):: (re)initialization
 MODOP=2 -- blocks emission in pi0 to gamma e+e-: DEFAULT
 MODOP=2 -- blocks emission in Kl  to gamma e+e-: DEFAULT



 ***************************************************************************
 *                         *****TAUOLA LIBRARY: VERSION 2.7 ******         *
 *                         ***********August   1995***************         *
 *                         **AUTHORS: S.JADACH, Z.WAS*************         *
 *                         **R. DECKER, M. JEZABEK, J.H.KUEHN*****         *
 *                         **AVAILABLE FROM: WASM AT CERNVM ******         *
 *                         ***** PUBLISHED IN COMP. PHYS. COMM.***         *
 *                          Physics initialization by CLEO collab          *
 *                          see Alain Weinstein www home page:             *
 *                         http://www.cithep.caltech.edu/~ajw/             *
 *                         /korb_doc.html#files                            *
 *                         *******CERN TH-6793 NOVEMBER  1992*****         *
 *                         **5 or more pi dec.: precision limited          *
 *                         *******CERN-TH-6793 NOVEMBER  1992*****         *
 *                         **5 or more pi dec.: precision limited          *
 *                         ******DEXAY ROUTINE: INITIALIZATION****         *
 *                   0     JAK1   = DECAY MODE FERMION1 (TAU+)             *
 *                   0     JAK2   = DECAY MODE FERMION2 (TAU-)             *
 ***************************************************************************


 ***************************************************************************




 ***************************************************************************
 *                         *****TAUOLA LIBRARY: VERSION 2.7 ******         *
 *                         ***********August   1995***************         *
 *                         **AUTHORS: S.JADACH, Z.WAS*************         *
 *                         **R. DECKER, M. JEZABEK, J.H.KUEHN*****         *
 *                         **AVAILABLE FROM: WASM AT CERNVM ******         *
 *                         ***** PUBLISHED IN COMP. PHYS. COMM.***         *
 *                          Physics initialization by CLEO collab          *
 *                          see Alain Weinstein www home page:             *
 *                         http://www.cithep.caltech.edu/~ajw/             *
 *                         /korb_doc.html#files                            *
 *                         *******CERN TH-6793 NOVEMBER  1992*****         *
 *                         **5 or more pi dec.: precision limited          *
 *                         ****DEKAY ROUTINE: INITIALIZATION******         *
 *                   0     JAK1   = DECAY MODE TAU+                        *
 *                   0     JAK2   = DECAY MODE TAU-                        *
 ***************************************************************************


 ***************************************************************************




 ***************************************************************************
 *                         *****TAUOLA UNIVERSAL INTERFACE: ******         *
 *                         *****VERSION 1.21, September 2005******         *
 *                         **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 ********         *
 *                                                                         *
 *                         ********** INITIALIZATION  ************         *
 *             1.00000     tau polarization switch must be 1 or 0          *
 *             1.57080     Higs scalar/pseudo mix CERN-TH/2003-166         *
 *         0               PI0 decay switch must be 1 or 0                 *
 *         0               ETA decay switch must be 1 or 0                 *
 *         0               K0S decay switch must be 1 or 0                 *
 ***************************************************************************




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

=================
 I                                  I                            I            I
 I   0 All included subprocesses    I         6000         42990 I  3.212D-11 I
 I 102 g + g -> h0                  I         4677         11690 I  2.498D-11 I
 I 123 f + f' -> f + f' + h0        I          360          8834 I  1.932D-12 I
 I 124 f + f' -> f" + f"' + h0      I          963         22466 I  5.203D-12 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 =  0.00000 *********

Using the command

cat *.stdout | grep "All included subprocesses" | awk '{print $10}' | cut -c 1-5 | awk '{ SUM += $1} END { print SUM}'
in the Res folder I got that the cross section for 150 GeV Higgs is 1.58784e-8. From DBS we see that the events we get are only 1,333,600 out of the 3 million so that mean the filter let 44.4533333% pass. This corresponds to a cross section of 7.05847808e-9. Now the units of this is mb according to the output. Petra said that this might not be right and it is perhaps nb.

9 June 2010
Output from 150 events with the filters on
     MSTU(12)       changed from              0 to          12345
WARNING: MCParticlePairFilter : size of some vectors not matching with 2!!
I love rockn Roll
WARNING: MCParticlePairFilter : size of some vectors not matching with 2!!
I love rockn Roll
----------------------------------------------
Initializing Tauola
Tauola: Polarization enabled

**********************************************************************

  PRETAUOLA_05 version for tau-lepton decays  
  VERTEX information is included
  tau-lepton decay modes used in TAUOLA package
    TAU+ force to decay into:  all decay modes              
    TAU- force to decay into:  all decay modes              

**********************************************************************

 Message from PHCORK(MODCOR):: initialization
 MODOP=1 -- no corrections on event: DEFAULT
 Message from PHOTOS: IPHQRK(MODCOR):: (re)initialization
 MODOP=1 -- blocks emission from light quarks: DEFAULT
 Message from PHOTOS: IPHEKL(MODCOR):: (re)initialization
 MODOP=2 -- blocks emission in pi0 to gamma e+e-: DEFAULT
 MODOP=2 -- blocks emission in Kl  to gamma e+e-: DEFAULT
1
 ********************************************************************************
 *                                                                              *
 *                          ==========================                          *
 *                            PHOTOS, Version:  2.15                            *
 *                            Released at:  11/10/ 5                            *
 *                          ==========================                          *
 *                                                                              *
 *                  PHOTOS QED Corrections in Particle Decays                   *
 *                                                                              *
 *         Monte Carlo Program - by E. Barberio, B. van Eijk and Z. Was         *
 *         Version 2.09  - by P. Golonka and Z.W.                               *
 *                                                                              *
 ********************************************************************************
 *                                                                              *
 *                  Internal input parameters:                                  *
 *                                                                              *
 *                  INTERF= T  ISEC= T  ITRE= F  IEXP= F  IFTOP= T   IFW= T     *
 *                  ALPHA_QED= 0.00730   XPHCUT=.100E-01                        *
 *                                                                              *
 *                  option with interference is active                          *
 *                  option with double photons is active                        *
 *                  emision in t tbar production is active                      *
 *                  correction wt in decay of W is active                       *
 *                                                                              *
 *          WARNING (1): /HEPEVT/ is not anymore the standard common block      *
 *                                                                              *
 *          PHOTOS expects /HEPEVT/ to have REAL*8 variables. To change to      *
 *          REAL*4 modify its declaration in subr. PHOTOS_GET PHOTOS_SET:       *
 *               REAL*8  d_h_phep,  d_h_vhep                                    *
 *          WARNING (2): check dims. of /hepevt/ /phoqed/ /ph_hepevt/.          *
 *          HERE:                     d_h_nmxhep=  4000 and  NMXHEP=10000       *
 *                                                                              *
 ********************************************************************************



 ***************************************************************************
 *                         *****TAUOLA LIBRARY: VERSION 2.7 ******         *
 *                         ***********August   1995***************         *
 *                         **AUTHORS: S.JADACH, Z.WAS*************         *
 *                         **R. DECKER, M. JEZABEK, J.H.KUEHN*****         *
 *                         **AVAILABLE FROM: WASM AT CERNVM ******         *
 *                         ***** PUBLISHED IN COMP. PHYS. COMM.***         *
 *                          Physics initialization by CLEO collab          *
 *                          see Alain Weinstein www home page:             *
 *                         http://www.cithep.caltech.edu/~ajw/             *
 *                         /korb_doc.html#files                            *
 *                         *******CERN TH-6793 NOVEMBER  1992*****         *
 *                         **5 or more pi dec.: precision limited          *
 *                         *******CERN-TH-6793 NOVEMBER  1992*****         *
 *                         **5 or more pi dec.: precision limited          *
 *                         ******DEXAY ROUTINE: INITIALIZATION****         *
 *                   0     JAK1   = DECAY MODE FERMION1 (TAU+)             *
 *                   0     JAK2   = DECAY MODE FERMION2 (TAU-)             *
 ***************************************************************************

  KORB: rho/rhop -> pi-pi0 mixing:
  KORB: rho   =  0.76990002      0.15120000    
  KORB: rhop  =   1.3700000      0.50999999     -0.11000000    
  KORB: K*/K*prime -> Kpi mixing:
  KORB: kstp  =   1.7000000      0.23500000     -3.79999988E-02
  KORB: a1/a1prime -> 3pi, KKpi mixing:
  KORB: a1    =   1.2750000      0.69999999    
  KORB: a1prim=   1.4610000      0.25000000       0.0000000    
  KORB: K1A/K1B -> Kpipi mixing:
  KORB: K1A   =   1.2700000      9.00000036E-02  1.00000000    
  KORB: K1B   =   1.4020000      0.17399999      0.80000001    
  KORB: rho/rhop/rhopp -> 4pi mixing:
  KORB: rho   =  0.76990002      0.15120000    
  KORB: rhopp =   1.4650000      0.31000000     -0.11000000    
  KORB: rhoppp=   1.7000000      0.23500000     -0.11000000    
  KORB: 3PI-PI0 PARAMS:   0.0000000      0.12420000      0.16040000      0.27110001      0.44430000    
  KORB: PI-3PI0 PARAMS:   0.0000000      1.00000000    
  In F3pi: add (rho-pi S-wave) + (rhop-pi S-wave) +
               (rho-pi D-wave) + (rhop-pi D-wave) +
    (f2 pi D-wave) + (sigmapi S-wave) + (f0pi S-wave)

 ***************************************************************************




 ***************************************************************************
 *                         *****TAUOLA LIBRARY: VERSION 2.7 ******         *
 *                         ***********August   1995***************         *
 *                         **AUTHORS: S.JADACH, Z.WAS*************         *
 *                         **R. DECKER, M. JEZABEK, J.H.KUEHN*****         *
 *                         **AVAILABLE FROM: WASM AT CERNVM ******         *
 *                         ***** PUBLISHED IN COMP. PHYS. COMM.***         *
 *                          Physics initialization by CLEO collab          *
 *                          see Alain Weinstein www home page:             *
 *                         http://www.cithep.caltech.edu/~ajw/             *
 *                         /korb_doc.html#files                            *
 *                         *******CERN TH-6793 NOVEMBER  1992*****         *
 *                         **5 or more pi dec.: precision limited          *
 *                         ****DEKAY ROUTINE: INITIALIZATION******         *
 *                   0     JAK1   = DECAY MODE TAU+                        *
 *                   0     JAK2   = DECAY MODE TAU-                        *
 ***************************************************************************


 ***************************************************************************




 ***************************************************************************
 *                         *****TAUOLA UNIVERSAL INTERFACE: ******         *
 *                         *****VERSION 1.21, September 2005******         *
 *                         **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 ********         *
 *                                                                         *
 *                         ********** INITIALIZATION  ************         *
 *             1.00000     tau polarization switch must be 1 or 0          *
 *             1.57080     Higs scalar/pseudo mix CERN-TH/2003-166         *
 *         0               PI0 decay switch must be 1 or 0                 *
 *         0               ETA decay switch must be 1 or 0                 *
 *         0               K0S decay switch must be 1 or 0                 *
 ***************************************************************************

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 /afs/cern.ch/cms/sw/slc5_ia32_gcc434/external/lhapdf/5.6.0-cms2/share/lhapdf/PDF
 with          1 members
 ====  initialized. ===========================
 Strong coupling at Mz for PDF is:  0.12978

 ==============================================================================
 I                                                                            I
 I              PYTHIA will be initialized for a p on p collider              I
 I                  at   7000.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   96   Semihard QCD 2 -> 2           I    9.6420D+03   I
           I  102   g + g -> h0                   I    1.0780D-10   I
           I  123   f + f' -> f + f' + h0         I    9.3532D-11   I
           I  124   f + f' -> f" + f"' + h0       I    2.3343D-10   I
           I                                      I                 I
           ==========================================================

 ****** PYMULT: initialization of multiple interactions for MSTP(82) = 4 ******
        pT0 = 2.25 GeV gives sigma(parton-parton) = 5.71D+02 mb: accepted

 ****** PYMIGN: initialization of multiple interactions for MSTP(82) = 4 ******
        pT0 = 2.25 GeV gives sigma(parton-parton) = 2.43D+02 mb: accepted

 ********************** PYINIT: initialization completed **********************
Begin processing the 1st record. Run 1, Event 1, LumiSection 1 at 09-Jun-2010 15:02:59 CEST

  PRETAUOLA: additional BR after          1 tau+- decays =   1.000000

Begin processing the 2nd record. Run 1, Event 2, LumiSection 1 at 09-Jun-2010 15:03:51 CEST
Begin processing the 3rd record. Run 1, Event 3, LumiSection 1 at 09-Jun-2010 15:03:51 CEST
Begin processing the 4th record. Run 1, Event 4, LumiSection 1 at 09-Jun-2010 15:03:51 CEST
Begin processing the 5th record. Run 1, Event 5, LumiSection 1 at 09-Jun-2010 15:03:51 CEST
Begin processing the 6th record. Run 1, Event 6, LumiSection 1 at 09-Jun-2010 15:03:51 CEST
Begin processing the 7th record. Run 1, Event 7, LumiSection 1 at 09-Jun-2010 15:03:51 CEST
Begin processing the 8th record. Run 1, Event 8, LumiSection 1 at 09-Jun-2010 15:03:51 CEST
Begin processing the 9th record. Run 1, Event 9, LumiSection 1 at 09-Jun-2010 15:04:57 CEST
Begin processing the 10th record. Run 1, Event 10, LumiSection 1 at 09-Jun-2010 15:04:57 CEST
Begin processing the 11th record. Run 1, Event 11, LumiSection 1 at 09-Jun-2010 15:04:57 CEST
Begin processing the 12th record. Run 1, Event 12, LumiSection 1 at 09-Jun-2010 15:04:57 CEST
Begin processing the 13th record. Run 1, Event 13, LumiSection 1 at 09-Jun-2010 15:06:25 CEST
Begin processing the 14th record. Run 1, Event 14, LumiSection 1 at 09-Jun-2010 15:06:25 CEST
Begin processing the 15th record. Run 1, Event 15, LumiSection 1 at 09-Jun-2010 15:07:27 CEST
Begin processing the 16th record. Run 1, Event 16, LumiSection 1 at 09-Jun-2010 15:07:27 CEST
Begin processing the 17th record. Run 1, Event 17, LumiSection 1 at 09-Jun-2010 15:07:27 CEST
Begin processing the 18th record. Run 1, Event 18, LumiSection 1 at 09-Jun-2010 15:08:53 CEST
Begin processing the 19th record. Run 1, Event 19, LumiSection 1 at 09-Jun-2010 15:09:37 CEST
Begin processing the 20th record. Run 1, Event 20, LumiSection 1 at 09-Jun-2010 15:10:15 CEST
Begin processing the 21st record. Run 1, Event 21, LumiSection 1 at 09-Jun-2010 15:11:19 CEST
Begin processing the 22nd record. Run 1, Event 22, LumiSection 1 at 09-Jun-2010 15:11:19 CEST
Begin processing the 23rd record. Run 1, Event 23, LumiSection 1 at 09-Jun-2010 15:11:19 CEST
Begin processing the 24th record. Run 1, Event 24, LumiSection 1 at 09-Jun-2010 15:11:19 CEST
Begin processing the 25th record. Run 1, Event 25, LumiSection 1 at 09-Jun-2010 15:11:59 CEST
Begin processing the 26th record. Run 1, Event 26, LumiSection 1 at 09-Jun-2010 15:12:34 CEST
Begin processing the 27th record. Run 1, Event 27, LumiSection 1 at 09-Jun-2010 15:13:28 CEST
Begin processing the 28th record. Run 1, Event 28, LumiSection 1 at 09-Jun-2010 15:13:28 CEST
Begin processing the 29th record. Run 1, Event 29, LumiSection 1 at 09-Jun-2010 15:14:28 CEST
Begin processing the 30th record. Run 1, Event 30, LumiSection 1 at 09-Jun-2010 15:14:28 CEST
Begin processing the 31st record. Run 1, Event 31, LumiSection 1 at 09-Jun-2010 15:14:28 CEST
Begin processing the 32nd record. Run 1, Event 32, LumiSection 1 at 09-Jun-2010 15:14:28 CEST
Begin processing the 33rd record. Run 1, Event 33, LumiSection 1 at 09-Jun-2010 15:15:06 CEST
Begin processing the 34th record. Run 1, Event 34, LumiSection 1 at 09-Jun-2010 15:15:38 CEST
Begin processing the 35th record. Run 1, Event 35, LumiSection 1 at 09-Jun-2010 15:15:38 CEST
Begin processing the 36th record. Run 1, Event 36, LumiSection 1 at 09-Jun-2010 15:16:18 CEST
Begin processing the 37th record. Run 1, Event 37, LumiSection 1 at 09-Jun-2010 15:16:18 CEST
Begin processing the 38th record. Run 1, Event 38, LumiSection 1 at 09-Jun-2010 15:16:18 CEST
Begin processing the 39th record. Run 1, Event 39, LumiSection 1 at 09-Jun-2010 15:16:18 CEST
Begin processing the 40th record. Run 1, Event 40, LumiSection 1 at 09-Jun-2010 15:16:18 CEST
Begin processing the 41st record. Run 1, Event 41, LumiSection 1 at 09-Jun-2010 15:16:18 CEST
Begin processing the 42nd record. Run 1, Event 42, LumiSection 1 at 09-Jun-2010 15:17:10 CEST
Begin processing the 43rd record. Run 1, Event 43, LumiSection 1 at 09-Jun-2010 15:17:10 CEST
Begin processing the 44th record. Run 1, Event 44, LumiSection 1 at 09-Jun-2010 15:17:11 CEST
Begin processing the 45th record. Run 1, Event 45, LumiSection 1 at 09-Jun-2010 15:17:46 CEST
Begin processing the 46th record. Run 1, Event 46, LumiSection 1 at 09-Jun-2010 15:17:46 CEST
Begin processing the 47th record. Run 1, Event 47, LumiSection 1 at 09-Jun-2010 15:17:46 CEST
Begin processing the 48th record. Run 1, Event 48, LumiSection 1 at 09-Jun-2010 15:18:06 CEST
Begin processing the 49th record. Run 1, Event 49, LumiSection 1 at 09-Jun-2010 15:19:35 CEST
Begin processing the 50th record. Run 1, Event 50, LumiSection 1 at 09-Jun-2010 15:19:35 CEST
Begin processing the 51st record. Run 1, Event 51, LumiSection 1 at 09-Jun-2010 15:21:11 CEST
Begin processing the 52nd record. Run 1, Event 52, LumiSection 1 at 09-Jun-2010 15:21:11 CEST
Begin processing the 53rd record. Run 1, Event 53, LumiSection 1 at 09-Jun-2010 15:21:49 CEST
Begin processing the 54th record. Run 1, Event 54, LumiSection 1 at 09-Jun-2010 15:21:49 CEST
Begin processing the 55th record. Run 1, Event 55, LumiSection 1 at 09-Jun-2010 15:22:49 CEST
Begin processing the 56th record. Run 1, Event 56, LumiSection 1 at 09-Jun-2010 15:23:30 CEST
Begin processing the 57th record. Run 1, Event 57, LumiSection 1 at 09-Jun-2010 15:23:30 CEST
Begin processing the 58th record. Run 1, Event 58, LumiSection 1 at 09-Jun-2010 15:24:51 CEST
Begin processing the 59th record. Run 1, Event 59, LumiSection 1 at 09-Jun-2010 15:26:07 CEST
Begin processing the 60th record. Run 1, Event 60, LumiSection 1 at 09-Jun-2010 15:26:52 CEST
Begin processing the 61st record. Run 1, Event 61, LumiSection 1 at 09-Jun-2010 15:26:52 CEST
Begin processing the 62nd record. Run 1, Event 62, LumiSection 1 at 09-Jun-2010 15:26:52 CEST
Begin processing the 63rd record. Run 1, Event 63, LumiSection 1 at 09-Jun-2010 15:26:53 CEST
Begin processing the 64th record. Run 1, Event 64, LumiSection 1 at 09-Jun-2010 15:26:53 CEST
Begin processing the 65th record. Run 1, Event 65, LumiSection 1 at 09-Jun-2010 15:26:53 CEST
Begin processing the 66th record. Run 1, Event 66, LumiSection 1 at 09-Jun-2010 15:26:53 CEST
Begin processing the 67th record. Run 1, Event 67, LumiSection 1 at 09-Jun-2010 15:26:53 CEST
Begin processing the 68th record. Run 1, Event 68, LumiSection 1 at 09-Jun-2010 15:26:53 CEST
Begin processing the 69th record. Run 1, Event 69, LumiSection 1 at 09-Jun-2010 15:27:52 CEST
Begin processing the 70th record. Run 1, Event 70, LumiSection 1 at 09-Jun-2010 15:29:05 CEST
Begin processing the 71st record. Run 1, Event 71, LumiSection 1 at 09-Jun-2010 15:29:05 CEST
Begin processing the 72nd record. Run 1, Event 72, LumiSection 1 at 09-Jun-2010 15:29:05 CEST
Begin processing the 73rd record. Run 1, Event 73, LumiSection 1 at 09-Jun-2010 15:30:25 CEST
Begin processing the 74th record. Run 1, Event 74, LumiSection 1 at 09-Jun-2010 15:30:25 CEST
Begin processing the 75th record. Run 1, Event 75, LumiSection 1 at 09-Jun-2010 15:31:47 CEST
Begin processing the 76th record. Run 1, Event 76, LumiSection 1 at 09-Jun-2010 15:32:47 CEST
Begin processing the 77th record. Run 1, Event 77, LumiSection 1 at 09-Jun-2010 15:32:47 CEST
Begin processing the 78th record. Run 1, Event 78, LumiSection 1 at 09-Jun-2010 15:32:47 CEST
Begin processing the 79th record. Run 1, Event 79, LumiSection 1 at 09-Jun-2010 15:32:47 CEST
Begin processing the 80th record. Run 1, Event 80, LumiSection 1 at 09-Jun-2010 15:33:34 CEST
Begin processing the 81st record. Run 1, Event 81, LumiSection 1 at 09-Jun-2010 15:33:34 CEST
Begin processing the 82nd record. Run 1, Event 82, LumiSection 1 at 09-Jun-2010 15:33:34 CEST
Begin processing the 83rd record. Run 1, Event 83, LumiSection 1 at 09-Jun-2010 15:33:34 CEST
Begin processing the 84th record. Run 1, Event 84, LumiSection 1 at 09-Jun-2010 15:33:34 CEST
Begin processing the 85th record. Run 1, Event 85, LumiSection 1 at 09-Jun-2010 15:33:34 CEST
Begin processing the 86th record. Run 1, Event 86, LumiSection 1 at 09-Jun-2010 15:34:40 CEST
Begin processing the 87th record. Run 1, Event 87, LumiSection 1 at 09-Jun-2010 15:34:40 CEST
Begin processing the 88th record. Run 1, Event 88, LumiSection 1 at 09-Jun-2010 15:35:47 CEST
Begin processing the 89th record. Run 1, Event 89, LumiSection 1 at 09-Jun-2010 15:36:36 CEST
Begin processing the 90th record. Run 1, Event 90, LumiSection 1 at 09-Jun-2010 15:36:36 CEST
Begin processing the 91st record. Run 1, Event 91, LumiSection 1 at 09-Jun-2010 15:37:02 CEST
Begin processing the 92nd record. Run 1, Event 92, LumiSection 1 at 09-Jun-2010 15:37:02 CEST
Begin processing the 93rd record. Run 1, Event 93, LumiSection 1 at 09-Jun-2010 15:37:02 CEST
Begin processing the 94th record. Run 1, Event 94, LumiSection 1 at 09-Jun-2010 15:38:25 CEST
Begin processing the 95th record. Run 1, Event 95, LumiSection 1 at 09-Jun-2010 15:39:53 CEST
Begin processing the 96th record. Run 1, Event 96, LumiSection 1 at 09-Jun-2010 15:39:54 CEST
Begin processing the 97th record. Run 1, Event 97, LumiSection 1 at 09-Jun-2010 15:39:54 CEST
Begin processing the 98th record. Run 1, Event 98, LumiSection 1 at 09-Jun-2010 15:39:54 CEST
Begin processing the 99th record. Run 1, Event 99, LumiSection 1 at 09-Jun-2010 15:39:54 CEST
Begin processing the 100th record. Run 1, Event 100, LumiSection 1 at 09-Jun-2010 15:39:54 CEST
Begin processing the 101st record. Run 1, Event 101, LumiSection 1 at 09-Jun-2010 15:39:54 CEST
Begin processing the 102nd record. Run 1, Event 102, LumiSection 1 at 09-Jun-2010 15:41:52 CEST
Begin processing the 103rd record. Run 1, Event 103, LumiSection 1 at 09-Jun-2010 15:41:52 CEST
Begin processing the 104th record. Run 1, Event 104, LumiSection 1 at 09-Jun-2010 15:43:35 CEST
Begin processing the 105th record. Run 1, Event 105, LumiSection 1 at 09-Jun-2010 15:44:22 CEST
Begin processing the 106th record. Run 1, Event 106, LumiSection 1 at 09-Jun-2010 15:44:22 CEST
Begin processing the 107th record. Run 1, Event 107, LumiSection 1 at 09-Jun-2010 15:45:01 CEST
Begin processing the 108th record. Run 1, Event 108, LumiSection 1 at 09-Jun-2010 15:45:01 CEST
Begin processing the 109th record. Run 1, Event 109, LumiSection 1 at 09-Jun-2010 15:45:01 CEST
Begin processing the 110th record. Run 1, Event 110, LumiSection 1 at 09-Jun-2010 15:45:01 CEST
Begin processing the 111th record. Run 1, Event 111, LumiSection 1 at 09-Jun-2010 15:45:01 CEST
Begin processing the 112th record. Run 1, Event 112, LumiSection 1 at 09-Jun-2010 15:45:01 CEST
Begin processing the 113th record. Run 1, Event 113, LumiSection 1 at 09-Jun-2010 15:45:29 CEST
Begin processing the 114th record. Run 1, Event 114, LumiSection 1 at 09-Jun-2010 15:45:29 CEST
Begin processing the 115th record. Run 1, Event 115, LumiSection 1 at 09-Jun-2010 15:45:29 CEST
Begin processing the 116th record. Run 1, Event 116, LumiSection 1 at 09-Jun-2010 15:45:30 CEST
Begin processing the 117th record. Run 1, Event 117, LumiSection 1 at 09-Jun-2010 15:45:30 CEST
Begin processing the 118th record. Run 1, Event 118, LumiSection 1 at 09-Jun-2010 15:45:30 CEST
Begin processing the 119th record. Run 1, Event 119, LumiSection 1 at 09-Jun-2010 15:46:20 CEST
Begin processing the 120th record. Run 1, Event 120, LumiSection 1 at 09-Jun-2010 15:47:04 CEST
Begin processing the 121st record. Run 1, Event 121, LumiSection 1 at 09-Jun-2010 15:47:51 CEST
Begin processing the 122nd record. Run 1, Event 122, LumiSection 1 at 09-Jun-2010 15:49:13 CEST
Begin processing the 123rd record. Run 1, Event 123, LumiSection 1 at 09-Jun-2010 15:49:46 CEST
Begin processing the 124th record. Run 1, Event 124, LumiSection 1 at 09-Jun-2010 15:49:46 CEST
Begin processing the 125th record. Run 1, Event 125, LumiSection 1 at 09-Jun-2010 15:50:45 CEST
Begin processing the 126th record. Run 1, Event 126, LumiSection 1 at 09-Jun-2010 15:50:45 CEST
Begin processing the 127th record. Run 1, Event 127, LumiSection 1 at 09-Jun-2010 15:50:45 CEST
Begin processing the 128th record. Run 1, Event 128, LumiSection 1 at 09-Jun-2010 15:51:17 CEST
Begin processing the 129th record. Run 1, Event 129, LumiSection 1 at 09-Jun-2010 15:52:12 CEST
Begin processing the 130th record. Run 1, Event 130, LumiSection 1 at 09-Jun-2010 15:52:12 CEST
Begin processing the 131st record. Run 1, Event 131, LumiSection 1 at 09-Jun-2010 15:52:12 CEST
Begin processing the 132nd record. Run 1, Event 132, LumiSection 1 at 09-Jun-2010 15:52:12 CEST
Begin processing the 133rd record. Run 1, Event 133, LumiSection 1 at 09-Jun-2010 15:52:12 CEST
Begin processing the 134th record. Run 1, Event 134, LumiSection 1 at 09-Jun-2010 15:53:01 CEST
Begin processing the 135th record. Run 1, Event 135, LumiSection 1 at 09-Jun-2010 15:54:03 CEST
Begin processing the 136th record. Run 1, Event 136, LumiSection 1 at 09-Jun-2010 15:55:00 CEST
Begin processing the 137th record. Run 1, Event 137, LumiSection 1 at 09-Jun-2010 15:55:00 CEST
Begin processing the 138th record. Run 1, Event 138, LumiSection 1 at 09-Jun-2010 15:55:00 CEST
Begin processing the 139th record. Run 1, Event 139, LumiSection 1 at 09-Jun-2010 15:55:25 CEST
Begin processing the 140th record. Run 1, Event 140, LumiSection 1 at 09-Jun-2010 15:56:10 CEST
Begin processing the 141st record. Run 1, Event 141, LumiSection 1 at 09-Jun-2010 15:57:04 CEST
Begin processing the 142nd record. Run 1, Event 142, LumiSection 1 at 09-Jun-2010 15:57:50 CEST
Begin processing the 143rd record. Run 1, Event 143, LumiSection 1 at 09-Jun-2010 15:57:50 CEST
Begin processing the 144th record. Run 1, Event 144, LumiSection 1 at 09-Jun-2010 15:57:50 CEST
Begin processing the 145th record. Run 1, Event 145, LumiSection 1 at 09-Jun-2010 15:57:50 CEST
Begin processing the 146th record. Run 1, Event 146, LumiSection 1 at 09-Jun-2010 15:58:36 CEST
Begin processing the 147th record. Run 1, Event 147, LumiSection 1 at 09-Jun-2010 15:59:14 CEST
Begin processing the 148th record. Run 1, Event 148, LumiSection 1 at 09-Jun-2010 15:59:15 CEST
Begin processing the 149th record. Run 1, Event 149, LumiSection 1 at 09-Jun-2010 15:59:15 CEST
Begin processing the 150th record. Run 1, Event 150, LumiSection 1 at 09-Jun-2010 16:00:00 CEST
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          150          1242 I  5.606D-11 I
 I 102 g + g -> h0                  I          109           294 I  4.029D-11 I
 I 123 f + f' -> f + f' + h0        I            9           274 I  4.331D-12 I
 I 124 f + f' -> f" + f"' + h0      I           32           674 I  1.144D-11 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 =  0.00000 *********


**********************************************************************
    PRETAUOLA_05 - summary: 
    TAU+ force to decay into:  all decay modes              
    TAU- force to decay into:  all decay modes              
    Nev =       150 N(taus) =        28 ADDITIONAL BR =   1.000000




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


=============================================

MessageLogger Summary

 type     category        sev    module        subroutine        count    total
 ---- -------------------- -- ---------------- ----------------  -----    -----
    1 CaloSim              -w OscarProducer:g4                       2*       2
    2 SiStripQuality       -w SiStripRawToClus                       1*       1

* Some occurrences of this message were suppressed in all logs, due to limits.

 type    category    Examples: run/evt        run/evt          run/evt
 ---- -------------------- ---------------- ---------------- ----------------
    1 CaloSim              1/127            1/127            
    2 SiStripQuality       Run: 1                            

Severity    # Occurrences   Total Occurrences
--------    -------------   -----------------
Warning                 3                   3
3319.372u 28.633s 1:01:11.87 91.1%      0+0k 0+0io 27052pf+0w

I got it all working so I am going to run it all today!!! Horray! I am not going to do the whole chain up to reconstruction just part of the way. I am posting my current version, just so it doesn't get lost.

# Auto generated configuration file
# using: 
# Revision: 1.177 
# Source: /cvs_server/repositories/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v 
# with command line options: Configuration/Generator/python/H200ZZ4L_Tauola_7TeV_cfi.py -s GEN,SIM,DIGI,L1,DIGI2RAW,HLT:GRun,RAW2DIGI,L1Reco --conditions auto:startup --datatier GEN-SIM-RAW --eventcontent RAWDEBUG -n 1 --no_exec
import FWCore.ParameterSet.Config as cms

process = cms.Process('HLT')

# import of standard configurations
process.load('Configuration.StandardSequences.Services_cff')
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.StandardSequences.MixingNoPileUp_cff')
process.load('Configuration.StandardSequences.GeometryDB_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('Configuration.StandardSequences.VtxSmearedEarly10TeVCollision_cff')
process.load('Configuration.StandardSequences.SimIdeal_cff')
process.load('Configuration.StandardSequences.Digi_cff')
process.load('Configuration.StandardSequences.SimL1Emulator_cff')
process.load('Configuration.StandardSequences.DigiToRaw_cff')
process.load('HLTrigger.Configuration.HLT_GRun_cff')
process.load('Configuration.StandardSequences.RawToDigi_cff')
process.load('Configuration.StandardSequences.L1Reco_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.load('Configuration.EventContent.EventContent_cff')


process.configurationMetadata = cms.untracked.PSet(
    version = cms.untracked.string('$Revision: 1.14 $'),
    annotation = cms.untracked.string('Configuration/Generator/python/H200ZZ4L_Tauola_7TeV_cfi.py nevts:10'),
    name = cms.untracked.string('PyReleaseValidation')
)
process.maxEvents = cms.untracked.PSet(
    input = cms.untracked.int32(30)
)

process.options = cms.untracked.PSet(
SkipEvent = cms.untracked.vstring('ProductNotFound')
)
# Input source
process.source = cms.Source("EmptySource")

# Output definition
process.output = cms.OutputModule("PoolOutputModule",
    splitLevel = cms.untracked.int32(0),
    outputCommands = process.RAWDEBUGEventContent.outputCommands,
    fileName = cms.untracked.string('30_skim.root'),
#H200ZZ4L_Tauola_7TeV_cfi_py_GEN_SIM_DIGI_L1_DIGI2RAW_HLT_RAW2DIGI_L1Reco.root                                  
    dataset = cms.untracked.PSet(
        dataTier = cms.untracked.string('GEN-SIM-RAW'),
        #dataTier = cms.untracked.string(''),
        filterName = cms.untracked.string('')
    ),
    SelectEvents = cms.untracked.PSet(
        SelectEvents = cms.vstring('generation_step')
    )
)

# Additional output definition

# Other statements
process.GlobalTag.globaltag = 'START37_V4::All'
process.generator = cms.EDFilter("Pythia6GeneratorFilter",
    ExternalDecays = cms.PSet(
        Tauola = cms.untracked.PSet(
            UseTauolaPolarization = cms.bool(True),
            InputCards = cms.PSet(
                mdtau = cms.int32(0),
                pjak2 = cms.int32(0),
                pjak1 = cms.int32(0)
            )
        ),
        parameterSets = cms.vstring('Tauola')
    ),
    pythiaPylistVerbosity = cms.untracked.int32(0),
    filterEfficiency = cms.untracked.double(1.0),
    pythiaHepMCVerbosity = cms.untracked.bool(False),
    comEnergy = cms.double(7000.0),
    maxEventsToPrint = cms.untracked.int32(0),
    PythiaParameters = cms.PSet(
        pythiaUESettings = cms.vstring('MSTJ(11)=3     ! Choice of the fragmentation function', 
            'MSTJ(22)=2     ! Decay those unstable particles', 
            'PARJ(71)=10 .  ! for which ctau  10 mm', 
            'MSTP(2)=1      ! which order running alphaS', 
            'MSTP(33)=0     ! no K factors in hard cross sections', 
            'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)', 
            'MSTP(52)=2     ! work with LHAPDF', 
            'MSTP(81)=1     ! multiple parton interactions 1 is Pythia default', 
            'MSTP(82)=4     ! Defines the multi-parton model', 
            'MSTU(21)=1     ! Check on possible errors during program execution', 
            'PARP(82)=1.8387   ! pt cutoff for multiparton interactions', 
            'PARP(89)=1960. ! sqrts for which PARP82 is set', 
            'PARP(83)=0.5   ! Multiple interactions: matter distrbn parameter', 
            'PARP(84)=0.4   ! Multiple interactions: matter distribution parameter', 
            'PARP(90)=0.16  ! Multiple interactions: rescaling power', 
            'PARP(67)=2.5    ! amount of initial-state radiation', 
            'PARP(85)=1.0  ! gluon prod. mechanism in MI', 
            'PARP(86)=1.0  ! gluon prod. mechanism in MI', 
            'PARP(62)=1.25   ! ', 
            'PARP(64)=0.2    ! ', 
            'MSTP(91)=1      !', 
            'PARP(91)=2.1   ! kt distribution', 
            'PARP(93)=15.0  ! '),
        processParameters = cms.vstring('PMAS(25,1)=200.0        !mass of Higgs', 
            'MSEL=0                  !(D=1) to select between full user control (0, then use MSUB) and some preprogrammed alternative: QCD hight pT processes (1, then ISUB=11, 12, 13, 28, 53, 68), QCD low pT processes (2, then ISUB=11, 12, 13, 28, 53, 68, 91, 92, 94, 95)', 
            'MSUB(102)=1             !ggH', 
            'MSUB(123)=1             !ZZ fusion to H', 
            'MSUB(124)=1             !WW fusion to H', 
            'CKIN(45)=5.             !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch', 
            'CKIN(46)=150.           !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
            'CKIN(47)=5.             !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
            'CKIN(48)=150.           !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
            'MDME(174,1)=0           !Z decay into d dbar', 
            'MDME(175,1)=0           !Z decay into u ubar', 
            'MDME(176,1)=0           !Z decay into s sbar', 
            'MDME(177,1)=0           !Z decay into c cbar', 
            'MDME(178,1)=0           !Z decay into b bbar', 
            'MDME(179,1)=0           !Z decay into t tbar', 
            'MDME(182,1)=1           !Z decay into e- e+', 
            'MDME(183,1)=1           !Z decay into nu_e nu_ebar', 
            'MDME(184,1)=1           !Z decay into mu- mu+', 
            'MDME(185,1)=1           !Z decay into nu_mu nu_mubar', 
            'MDME(186,1)=1           !Z decay into tau- tau+', 
            'MDME(187,1)=1           !Z decay into nu_tau nu_taubar', 
            'MDME(210,1)=0           !Higgs decay into dd', 
            'MDME(211,1)=0           !Higgs decay into uu', 
            'MDME(212,1)=0           !Higgs decay into ss', 
            'MDME(213,1)=0           !Higgs decay into cc', 
            'MDME(214,1)=0           !Higgs decay into bb', 
            'MDME(215,1)=0           !Higgs decay into tt', 
            'MDME(216,1)=0           !Higgs decay into', 
            'MDME(217,1)=0           !Higgs decay into Higgs decay', 
            'MDME(218,1)=0           !Higgs decay into e nu e', 
            'MDME(219,1)=0           !Higgs decay into mu nu mu', 
            'MDME(220,1)=0           !Higgs decay into tau nu tau', 
            'MDME(221,1)=0           !Higgs decay into Higgs decay', 
            'MDME(222,1)=0           !Higgs decay into g g', 
            'MDME(223,1)=0           !Higgs decay into gam gam', 
            'MDME(224,1)=0           !Higgs decay into gam Z', 
            'MDME(225,1)=1           !Higgs decay into Z Z', 
            'MDME(226,1)=0           !Higgs decay into W W'),
        parameterSets = cms.vstring('pythiaUESettings', 
            'processParameters')
    )
)

process.leptonfilter = cms.EDFilter("MCParticlePairFilter",
    Status = cms.untracked.vint32(3, 3),
    ParticleID1 = cms.untracked.vint32(11, -11, 13, -13, 15, -15),
    ParticleID2 = cms.untracked.vint32(11, -11, 13, -13, 15, -15)
)
process.neutrinofilter = cms.EDFilter("MCParticlePairFilter",
    Status = cms.untracked.vint32(3, 3),
    ParticleID1 = cms.untracked.vint32(12, -12, 14, -14, 16, -16),
    ParticleID2 = cms.untracked.vint32(12, -12, 14, -14, 16, -16)
)

process.ProductionFilterSequence = cms.Sequence(process.generator*process.leptonfilter*process.neutrinofilter)


# Path and EndPath definitions
process.generation_step = cms.Path(process.pgen)
process.simulation_step = cms.Path(process.psim)
process.digitisation_step = cms.Path(process.pdigi)
process.L1simulation_step = cms.Path(process.SimL1Emulator)
process.digi2raw_step = cms.Path(process.DigiToRaw)
#process.raw2digi_step = cms.Path(process.RawToDigi)
#process.L1Reco_step = cms.Path(process.L1Reco)
#process.RECO_step = cms.Path(process.RECO)
process.endjob_step = cms.Path(process.endOfProcess)
process.out_step = cms.EndPath(process.output)

# Schedule definition
process.schedule = cms.Schedule(process.generation_step,process.simulation_step,process.digitisation_step,process.L1simulation_step,process.digi2raw_step)
process.schedule.extend(process.HLTSchedule)
#process.schedule.extend([process.raw2digi_step,process.L1Reco_step,process.endjob_step,process.out_step])
process.schedule.extend([process.endjob_step,process.out_step])
#process.schedule.extend(process.ParticleListDrawer)
# special treatment in case of production filter sequence  
for path in process.paths:
    getattr(process,path)._seq = process.ProductionFilterSequence*getattr(process,path)._seq
#    getattr(process,path)._seq = process.generator*getattr(process,path)._seq


    

7 June 2010
So I was getting errors of this type
%MSG-s CMSException:  AfterFile 07-Jun-2010 19:39:02 CEST PostEndRun
cms::Exception caught in cmsRun
---- EventProcessorFailure BEGIN
EventProcessingStopped
---- ScheduleExecutionFailure BEGIN
ProcessingStopped
---- ProductNotFound BEGIN
getByLabel: Found zero products matching all criteria
Looking for sequence of type: reco::Candidate
Looking for module label: genParticleCandidates
Looking for productInstanceName:
cms::Exception going through module ParticleDecayDrawer/printDecay run: 1 lumi: 1 event: 1
If you wish to continue processing events after a ProductNotFound exception,
add "SkipEvent = cms.untracked.vstring('ProductNotFound')" to the "options" PSet in the configuration.
---- ProductNotFound END
Exception going through path p
---- ScheduleExecutionFailure END
an exception occurred during current event processing
cms::Exception caught in EventProcessor and rethrown
---- EventProcessorFailure END


%MSG

And I finally found the fix. You need to change the cms.InputTag from genParticlesCandidates to just genParticles. Hooray!

3 June 2010
I keep getting the error
The plugin 'HZZ4LeptonsMCGenFilter' is found in multiple files 
 'pluginGeneratorInterfacePythia6Filters.so'
 'pluginHZZ4LeptonsFilters.so'
in directory '/afs/cern.ch/user/m/mkress/scratch0/2L2NU/CMSSW_3_7_0/lib/slc5_ia32_gcc434'.
The code must be changed so the plugin only appears in one plugin file. You will need to remove the macro which registers the plugin so it only appears in one of these files.
  If none of these files register such a plugin, then the problem originates in a library to which all these files link.
The plugin registration must be removed from that library since plugins are not allowed in regular libraries.

I had accidentally my plugin twice so to fix this I had to remove the double definition in the Buildfile and then run

 scramv1 b clean 
and then recompile. I also ran the command
scramv1 b --reset 
which didn't work on it's own, but perhaps needed to be used in tandem. I am not sure, but at lease using both of them worked. (This took me two days to figure out)

2 June 2010
I needed to add to my version of HZZ4LeptonsMCGenFilter.cc at the end.
DEFINE_FWK_MODULE(HZZ4LeptonsMCGenFilter);
This finally allowed the thing to work! Hooray! Now I need to work on adding the part that lists the particles, good luck smile

I also had to add the line

SkipEvent = cms.untracked.vstring('ProductNotFound')
inside
process.options = cms.untracked.PSet(
)
in order to keep going after an event has been filtered.

1 June 2010
I am working on the Filter and found that I need to change the file HZZ4LeptonsMCGenFilter.cc to not filter neutrinos. That needed me to add those lines which I had to look up the numbers. I have attached a pdf to the main page for HtoZZto2L2NU that has those.

I am now trying to put the files from the HiggstoZZ package HZZ4LeptonsMCGenFilter.cc, HZZ4LeptonsMCGenFilter.h and the corresponding two python files in the Pythia6Generator part. I had to change in the header of HZZ4LeptonsMCGenFilter.h , FWCore/ParameterSet/interface/InputTag.h to FWCore/Utilities/interface/InputTag.h. Not a big deal (error tells us what to change it to) but took me awhile to figure it out. I think I am going to try and make the entire HiggstoZZ code compatible with this newest version CMSSW_3_7_0 seeing as it will be validated for ICHEP it is probably a good release to use. I am also thinking of using the header from the old H200ZZ4l code for the base of my H200ZZ2L2NU code instead of the random code Nicola had me generate. I will check out if it is close or the same other than the changes made.

31 May 2010
So all my jobs got errors. I am using 3.7.0 but apparently that may not be on at the fnal server I am submitting with condor. I may need to use an older version. I will be testing this out again today.

I went and talked with Nicola and he helped me make a new python file to use with cmsRun here are the steps Look up on DBS chosing the Software release that you want and the type you want for example MonteCarlo. Click on the configuration part and check out this part

# with command line options: ZMM_7TeV.cfi -s GEN:ProductionFilterSequence,SIM,DIGI,L1,DIGI2RAW,HLT:GRun,RAW2DIGI,L1Reco -n 10 --geometry DB --conditions auto:startup --no_exec
Then run cmsDriver.py with those command line options and you create the python file.

Now you need to take the part from the old python file that tells what products to have. Here is the one for 2L2NU

process.generator = cms.EDFilter("Pythia6GeneratorFilter",
    ExternalDecays = cms.PSet(
        Tauola = cms.untracked.PSet(
            UseTauolaPolarization = cms.bool(True),
            InputCards = cms.PSet(
                mdtau = cms.int32(0),
                pjak2 = cms.int32(0),
                pjak1 = cms.int32(0)
            )
        ),
        parameterSets = cms.vstring('Tauola')
    ),
    pythiaPylistVerbosity = cms.untracked.int32(0),
    filterEfficiency = cms.untracked.double(1.0),
    pythiaHepMCVerbosity = cms.untracked.bool(False),
    comEnergy = cms.double(7000.0),
    maxEventsToPrint = cms.untracked.int32(0),
    PythiaParameters = cms.PSet(
        pythiaUESettings = cms.vstring('MSTJ(11)=3     ! Choice of the fragmentation function', 
            'MSTJ(22)=2     ! Decay those unstable particles', 
            'PARJ(71)=10 .  ! for which ctau  10 mm', 
            'MSTP(2)=1      ! which order running alphaS', 
            'MSTP(33)=0     ! no K factors in hard cross sections', 
            'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)', 
            'MSTP(52)=2     ! work with LHAPDF', 
            'MSTP(81)=1     ! multiple parton interactions 1 is Pythia default', 
            'MSTP(82)=4     ! Defines the multi-parton model', 
            'MSTU(21)=1     ! Check on possible errors during program execution', 
            'PARP(82)=1.8387   ! pt cutoff for multiparton interactions', 
            'PARP(89)=1960. ! sqrts for which PARP82 is set', 
            'PARP(83)=0.5   ! Multiple interactions: matter distrbn parameter', 
            'PARP(84)=0.4   ! Multiple interactions: matter distribution parameter', 
            'PARP(90)=0.16  ! Multiple interactions: rescaling power', 
            'PARP(67)=2.5    ! amount of initial-state radiation', 
            'PARP(85)=1.0  ! gluon prod. mechanism in MI', 
            'PARP(86)=1.0  ! gluon prod. mechanism in MI', 
            'PARP(62)=1.25   ! ', 
            'PARP(64)=0.2    ! ', 
            'MSTP(91)=1      !', 
            'PARP(91)=2.1   ! kt distribution', 
            'PARP(93)=15.0  ! '),
        processParameters = cms.vstring('PMAS(25,1)=115.0        !mass of Higgs', 
            'MSEL=0                  !(D=1) to select between full user control (0, then use MSUB) and some preprogrammed alternative: QCD hight pT processes (1, then ISUB=11, 12, 13, 28, 53, 68), QCD low pT processes (2, then ISUB=11, 12, 13, 28, 53, 68, 91, 92, 94, 95)', 
            'MSUB(102)=1             !ggH', 
            'MSUB(123)=1             !ZZ fusion to H', 
            'MSUB(124)=1             !WW fusion to H', 
            'CKIN(45)=5.             !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch', 
            'CKIN(46)=150.           !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
            'CKIN(47)=5.             !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
            'CKIN(48)=150.           !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch', 
            'MDME(174,1)=0           !Z decay into d dbar', 
            'MDME(175,1)=0           !Z decay into u ubar', 
            'MDME(176,1)=0           !Z decay into s sbar', 
            'MDME(177,1)=0           !Z decay into c cbar', 
            'MDME(178,1)=0           !Z decay into b bbar', 
            'MDME(179,1)=0           !Z decay into t tbar', 
            'MDME(182,1)=1           !Z decay into e- e+', 
            'MDME(183,1)=1           !Z decay into nu_e nu_ebar', 
            'MDME(184,1)=1           !Z decay into mu- mu+', 
            'MDME(185,1)=1           !Z decay into nu_mu nu_mubar', 
            'MDME(186,1)=1           !Z decay into tau- tau+', 
            'MDME(187,1)=1           !Z decay into nu_tau nu_taubar', 
            'MDME(210,1)=0           !Higgs decay into dd', 
            'MDME(211,1)=0           !Higgs decay into uu', 
            'MDME(212,1)=0           !Higgs decay into ss', 
            'MDME(213,1)=0           !Higgs decay into cc', 
            'MDME(214,1)=0           !Higgs decay into bb', 
            'MDME(215,1)=0           !Higgs decay into tt', 
            'MDME(216,1)=0           !Higgs decay into', 
            'MDME(217,1)=0           !Higgs decay into Higgs decay', 
            'MDME(218,1)=0           !Higgs decay into e nu e', 
            'MDME(219,1)=0           !Higgs decay into mu nu mu', 
            'MDME(220,1)=0           !Higgs decay into tau nu tau', 
            'MDME(221,1)=0           !Higgs decay into Higgs decay', 
            'MDME(222,1)=0           !Higgs decay into g g', 
            'MDME(223,1)=0           !Higgs decay into gam gam', 
            'MDME(224,1)=0           !Higgs decay into gam Z', 
            'MDME(225,1)=1           !Higgs decay into Z Z', 
            'MDME(226,1)=0           !Higgs decay into W W'),
        parameterSets = cms.vstring('pythiaUESettings', 
            'processParameters')
    )
)
Then you need to take out any filters or extra garbage that is in there. This example had one.

Now I need to make a filter basing if off the Higgs package file HZZ4LeptonsMCGenFilter.cc file. Fun....

29 May 2010
So I accentually deleted my Crab files for the job I have been running. Grrrrrrr.. anyways I am attempting to submit them through fermilab and condor. We will see how that goes. I had some trouble saving to castor. I found using the following in the crab.cfg file
storage_element = srm.cern.ch
storage_path = /srm/managerv1?SFN=/castor/cern.ch/user/u/username/subdir
I get
crab: You are asking to stage out without using CMS Storage Name convention. In this case you 
    must specify both user_remote_dir and storage_path in the crab.cfg section [USER].
     For further information please visit : 
   https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideCrabHowTo#Stageout_and_publication

When I use what I used on lxplus which worked fine.

storage_element=srm-cms.cern.ch
storage_path=/srm/managerv2?SFN=/castor/cern.ch
user_remote_dir=/user/m/mkress/higgsdir
I get
crab: Problems trying remote dir check: 
   Failure while checking remote dir: 'Wrong option passed to the command'

to set up fnal

source /uscmst1/prod/sw/cms/cshrc uaf
cd ....
cmsenv
source /uscmst1/prod/grid/CRAB_2_7_2_p1/crab.csh
source /uscmst1/prod/grid/gLite_SL5_CRAB_27x.csh


---++++++24 May 2010

I had troubles writing my files to castor and has to end up giving all my folders permission 777.  I don't really like this (for obvious security reasons) but that has at least made me able to occasionally not have a write error.
Edit | Attach | Watch | Print version | History: r14 < r13 < r12 < r11 < r10 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r14 - 2010-08-17 - MatthewKress
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main 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