Information about the flat ntuple

Latest Tag: flatnt_59719c

59718 Tag: flatnt_59718a

Please email Collin Wolfe or Stephen Levy if you have any problems.

59717 Tag: flatnt_59717d

Please email Collin Wolfe or Stephen Levy if you have any problems.

59716 Tag: flatnt_59716

Please email Collin Wolfe if you have any problems.

Where to find the flat ntupler (directories are specific to Chicago systems):

  
  PLEASE NOTE:
    March 09, 2004
    This section has not been updated for awhile.  The ntuple script
    uses the current stable version of the UCNtuple.

The current "stable" version is in /cdf/data32a/cwolfe/59716_ntuple/. A sample .tcl file for running the ntupler is /cdf/data32a/cwolfe/59716_ntuple/UCNtuple/run.tcl. To run the ntupler, change directory to /cdf/data32a/cwolfe/59716_ntuple/, and run 

export FLOUTPUT=outputfile_for_ntuple
ucntuple-stable.exe multi.tcl -i filename_to_ntuple
(assuming you're using bash, of course; if you use tcsh say setenv FLOUTPUT outputfile_for_ntuple). The flat ntupler is linked statically, so in fact you can run this executable as-is from anywhere, copy it to anywhere and run it, etc.

Where to find flat ntuples (on Chicago systems):

  
  PLEASE NOTE:
    March 09, 2004
    This section has not been updated for awhile.
    A current listing of the most used RunII ucntuples can be found here: DataSets

For stream A, made from the super-express stream: /cdf/data26b/flatntuples

For stream A, made from the aexp03 dataset (production 4.5.2): /cdf/data26c/flatntuples_aexp03

For stream B high-pt electrons (bhel03 dataset, production 4.5.2): /cdf/data26c/flatntuples_bhel03

For stream B high-pt muons (bhmu03 dataset, production 4.5.2): /cdf/data26c/flatntuples_bhmu03

For the stream E tau+met sample (eexo03 dataset, production 4.5.2): /cdf/data26c/flatntuples_eexo03

Notes for Z->b bbar operation:

To obtain the new blocks, it is necessary to set the talk-to option "addZbbBlocks" to true (by default none of the new blocks are made). In addition, "simulatedData" must be set to true to get the HEPG, OBSP, OBSV, and Zbb blocks, and options in the "zbbMenu" control creation of the Level 3 blocks.

An example analysis script with the flat ntuple:

The following three files are an analysis script to find central-central Z->ee with the flat ntuple:

Disclaimer:  These files will not necessarily provide results consistent with the accepted values. There are no known macros providing results consistent with the accpeted values. Please find some and post them here in place of these, non-fully-debugged macros.

To use:

  1. Download these files and place them in the same directory.
  2. From that directory, type ln -s /cdf/data32a/cwolfe/59716_ntuple/include/UCNtuple . (this is to get the include files to work properly)
  3. Enter ROOT: root
  4. Create a TChain: (for example, to get runs 1456xx) 
    TChain tc("ucntuple")
tc.Add("/cdf/data26b/flatntuples/flat_4.5.0_1456*.root")

    OR, instead, you can type 

    .x load.C 

    to load a file that is already set up to load the data files
  5. Load the loop code: .L zee_tight.C
  6. Start the loop: zsearch(&tc)
  7. Some of the plots generated are zee and runnumber. You can plot them by typing zee.Draw() and runnumber.Draw().

Additional sample macros (ALSO NOT FULLY DEBUGGED, PLEASE REPLACE WITH MACROS THAT ARE!): 




How to Loop Through Data in ROOT:


    

    You can find information about looping through data in ROOT here.
Or you can view the file by typing 
     xdvi ~psgibbs/RootLoop.dvi
     at your terminal.






Talk-tos:


AC++> mod talk UCNtupleModule

UCNtupleModule> help
Command(s) available in the "UCNtupleModule" module:
      exit		Leave the current menu, module, or process
      help		Bring up help text for the current context.
      show		Display the value of any parameters or statistics.
      parmSetName	general parameter of type string, set or list:
         Name to give generated parameter set at exit.
         Usage: parmSetName set 
      readParmSet	Read and fill the value of the params of this module.
      echo		Send text argument to stdout (useful in scripts).
      action		Command used to control module actions
      verbose		general parameter of type bool, set or list:
         Turn on verbose screen output. Usage: verbose set t/f
      production	general parameter of type bool, set or list:
         Suppress all screen output. Usage: production set t/f
      root		Type root to get to ROOT prompt 
      call		define name of the interpreted function to be executed for this module
      splitMode		general parameter of type int, set or list
      bufferSize	general parameter of type int, set or list
      compressionLevel	general parameter of type int, set or list
      histfile		
      maxSize		
      processName	general parameter of type string, set or list:
	Choose which process name to use to select objects.
      jetCollName	general parameter of type string, set or list:
	Choose which jet collection to use in making the first jet block.
      jetColl2Name	general parameter of type string, set or list:
	Choose which jet collection to use in making the second jet block.
      jetColl3Name	general parameter of type string, set or list:
	Choose which jet collection to use in making the third jet block.
      tagCollName	general parameter of type string, set or list:
	Choose which SecVtxTrk collection to use in making the first tagging block.
      tagColl2Name	general parameter of type string, set or list:
	Choose which SecVtxTrk collection to use in making the second tagging block.
      tagColl3Name	general parameter of type string, set or list:
	Choose which SecVtxTrk collection to use in making the third tagging block.
      tagTrkName	general parameter of type string, set or list
      tagTrk2Name	general parameter of type string, set or list
      tagTrk3Name	general parameter of type string, set or list
      tagCollJetBlock	general parameter of type string, set or list:
	Choose which Jet block to match jets in the first tagging block.
	Should match secvtx and jet block talk-tos!
      tagColl2JetBlock	general parameter of type string, set or list:
	Choose which Jet block to match jets in the 2nd tagging block.
	Should match secvtx and jet block talk-tos!
      tagColl3JetBlock	general parameter of type string, set or list:
	Choose which Jet block to match jets in the 3rd tagging block.
	Should match secvtx and jet block talk-tos!
      vertexCollType	general parameter of type string, set or list:
	Choose which type of vertex collection to use
\tValid options are VertexColl and ZVertexColl
      vertexCollName	general parameter of type string, set or list:
	Choose which vertex collection to use.
      eventVertexStrategy	general parameter of type int, set or list:
	Choose what strategy is used for finding the event vertex
	(vertex.vzev).
		0 = highest sum-pt vertex, else zero;
		1 = vertex closest to z0 of highest pt good lepton.
      revertexMet	general parameter of type bool, set or list:
	Recompute missing Et with event vertex.
      useRefitTracks	general parameter of type bool, set or list:
	Use tracks that have been refit with TrackRefitModule.
      tracksRefitNow	general parameter of type bool, set or list:
	Set to true if TrackRefitModule is being run in this executable, or
	to false if refit tracks are already present in the data being read in.
      simulatedData	general parameter of type bool, set or list:
	Set this to true to include HEPG, OBSP, and OBSV (and Zbb) blocks.
      addZbbBlocks	general parameter of type bool, set or list:
	Set this to true to include Z->b bbar analysis specific blocks.
      addTauBlock	general parameter of type bool, set or list:
	Set this to true to include Tau analysis specific blocks.
      topSummaryName	general parameter of type string, set or list:
	Name of the TopSummary Class (Set to Top).
      topSummaryDesc	general parameter of type string, set or list:
	Description of the TopSummary object (Set to newTopCollectionPrefix value in TopEventModule)
      zbbMenu		Invokes the zbbMenu submenu


UCNtupleModule> zbbMenu

UCNtupleModule/zbbMenu> help
Commands belonging to the "Z->b bbar options" Menu:
      exit
      help
      show
      level3Info        general parameter of type bool, set or list:
        Set this to true to include Level3 blocks.
      l3jet04CollName   general parameter of type string, set or list:
        Choose name of the jet collection to use in making level3 jets
        (JetCluModule-cone0.4 to use default) and process name
        specificator (L3 is a default) separated by semicolon, like
        :.
      l3jet07CollName   general parameter of type string, set or list:
        Choose name of the jet collection to use in making level3 jets
        (JetCluModule-cone0.7 to use default) and process name
        specificator (L3 is a default) separated by semicolon, like
        :.
      l3trkCollName	general parameter of type string, set or list:
	name of the track view to match with the svt tracks (defTracks to
	 use default) and process name specificator (curent process is a 
	default) separated by semicolomn, like :
      beamStrategy      general parameter of type int, set or list:
        Please document me! :-)
      cot2svt		Invokes the cot2svt submenu

UCNtupleModule/zbbMenu> cot2svt

UCNtupleModule/cot2svt> help
Commands belonging to the "Cot2SvtMatch parameters" Menu:
      exit		
      help		
      show		
      chi2Max		general parameter of type float, set or list:
			 Maximum allowed CHI2 of SVT track
      dphiMax		general parameter of type float, set or list:
			 Maximum allowed difference in PHI between COT and SVT tracks
      dcurvMax		general parameter of type float, set or list:
			 Maximum allowed difference in curvature between COT and SVT tracks


UCNtupleModule/zbbMenu> exit

UCNtupleModule> show
Current value of item(s) in the "UCNtupleModule" module:

      Value of parmSetName for this module  is IGNORE
      Value of verbose for this module  is false
      Value of production for this module  is false
                  Value of splitMode for this module  is 99
      Value of bufferSize for this module  is 64000
      Value of compressionLevel for this module  is 1
                  Value of processName for this module  is PROD
      Value of jetCollName for this module  is JetCluModule-cone0.4
      Value of jetColl2Name for this module  is JetCluModule-cone0.7
      Value of jetColl3Name for this module  is Top_JETCLU0.4_ReclusterJetColl
      Value of tagCollName for this module  is SecVtxColl
      Value of tagColl2Name for this module  is SecVtxColl
      Value of tagColl3Name for this module  is SecVtxColl
      Value of tagTrkName for this module  is SecVtxTracks
      Value of tagTrk2Name for this module  is SecVtxTracks
      Value of tagTrk3Name for this module  is SecVtxTracks
      Value of tagCollJetBlock for this module  is jet
      Value of tagColl2JetBlock for this module  is jet2
      Value of tagColl3JetBlock for this module  is jet3
      Value of vertexCollType for this module  is VertexColl
      Value of vertexCollName for this module  is DEFAULT_VXPRIM_VERTICES
      Value of eventVertexStrategy for this module  is 0
      Value of revertexMet for this module  is true
      Value of useRefitTracks for this module  is true
      Value of tracksRefitNow for this module  is true
      Value of simulatedData for this module  is false
      Value of addZbbBlocks for this module  is false
      Value of addTauBlock for this module  is true
      Value of topSummaryName for this module  is Top
      Value of topSummaryDesc for this module  is Top
      Current value of items on the "zbbMenu" menu:
            Value of level3Info for this module  is false
            Value of l3jet04CollName for this module  is JetCluModule-cone0.4:L3
            Value of l3jet07CollName for this module  is JetCluModule-cone0.7:L3
            Value of l3trkCollName for this module  is L3
            Value of beamStrategy for this module  is 0
            Current value of items on the "cot2svt" menu:
            Value of chi2Max for this module  is 25
            Value of dphiMax for this module  is 0.015
            Value of dcurvMax for this module  is 0.00015






How to recompile the flat ntupler for yourself:


Please follow these steps to recompile the flat ntupler version
59719c tagged on 13 Oct 2004:


    

  1. Create a new test release of offline version 5.3.3_nt.
This is the only release that the current version of the flat
ntupler is known for a fact to work properly with.  Setup the
necessary environment, then execute

    newrel -t 5.3.3_nt ucntuple; cd ucntuple

    
at the command line.

    2. 

  2. addpkg FlatNtuple flatnt_59719c

    3. 

  3. If you want to use shared libraries (and we recommend that you do):
    
setenv USESHLIBS 1 
    4. 

  4. Execute the following:

    sh FlatNtuple/build.sh

    

    5. 

  5. You will now have the executable!  The recommended tcl file to use
in 5.3.X is FlatNtuple/tcl/53X_multi.tcl.

    6. 






Release History:


  static const Int_t NTVERSION = 59719;

  // 59719:  Add conversionpartner to track block; Update
  //         myelectron.convword to use top group conversion code; Depricate
  //         myelectron.convrad; Re-vertex met from x-y beam position; add slope info
  //         header block; Add caltower.emTdcNhits, caltower.emTdcTime,
  //         and caltower.emTdcWidth to take advantage of new em TDCs;
  //         add 4-momentum variables to phoenix block.
  //         Add muon region (i.e., BlueBeam, Miniskirt...) info to muon block;
  //         Update muon.fiddet to make fiducial cuts.
  // 59718 : Use FlatNtuple/FlatNtuple/convert.pl to create
  //         declarestructs instead of writing it by hand.  This script
  //         also makes the FlatNtuple/src/my(object).cc files which
  //         construct the block for each object.
  // 59717d: Make D0FinderBlock tcl settable and turn it off by
  //         default.  Takes an order of mag longer cpu sec / event
  //         than any other module we run.  Run cosmic modules for 491 MC.
  // 59717c: Fix typo (called ConstructCalTowerBlock instead of
  //         ComputeCalTowerBlock).
  // 59717b: Add CdfMetModule and tcl file for remaking 4.8.4 processed
  //         data in 4.11.2 a la top group.
  // 59717 : Remove duplicate code structure in TagBlock (we were
  //       : cutting & pasting code for tagblock, tagblock2, and tagblock3 and
  //       : it looks like tagblock2 was == tagblock3, see additional comments
  //       : in code.)  Change structure of EleBlock.  First we find the tight
  //       : electrons and put the information needed by the VertexBlock to
  //       : find the best vertex into myelectron (etcorr, detector, ztrk,
  //       : ztrerr).  We go back to the electrons and set the zvtx for plug
  //       : electrons to the event zVertex.  Then we get the myelectron
  //       : information for all remaining variables.  Added tagblock2
  //       : and tagblock3 to declarestructs.  Also remove calls to
  //       : unused functions in UCNtuple so it's easier to see what is
  //       : actually happening.
  // 59716 : Fix up Vertex algo. Correct Jet info for vertex.
  //       : Add third tagging block. Add phoenix electrons.
  //       : Update to 4.11.2
  // 59715 : Add Zbb group changes & tau code from M. Reece.
  // 59714 : Compensate for change to HADTDC information, update jet 
  //       : corrections, update to 4.11.1
  // 59713 : Integrate Collin's top block, update to 4.10.4
  // 59712 : Fix bug that caused cot-only tracks not to be beam 
  //       : constrained; add variable for number of class 12 vertices;
  //       : use updated jet corrections code; add ability to remake 
  //       : CES,CPR clusters if desired; add new pi0 stuff in tau block
  //       : (from Matt R.);
  // 59711 : Fix Zbb declaration problem; bring muon and electron id words
  //       : into substantial compliance with ewk group cuts; always store
  //       : CosmicFinderModule bitword; add et out of time from 
  //       : CosmicRayInfo; add ntracksinit variable for photons;
  // 59710 : Fix jet-met problem; fix beamlines for monte carlo
  // 59709 : fix level 3 trigger bank to deal with changes in 4.8.4 L3
  // 59708 : Add ability to use vertex closest to highest-pt good lepton as
  //       : event vertex
  // 59707 : Fix jet-to-track link slowness; add new Level3 block changes;
  //       : add Beate's (major) fixes to get correct calorimetry information
  // 59706 : Fix zbb block multi-creation bug, let jet2 roam free, prevent
  //       : jet troubles from terminating jet2 creation, include Si hit
  //       : map in track block, update to 4.9.1_v02 tags for btagging,
  //       : include track-to-jet links, fix zbb block ordering
  // 59705 : Fix track isolation computation bug
  // 59704 : Add 4.9.1 track refits
  // 59703 : Add CMIO corrections to met block
  // 59702 : Update to 4.9.1
  // 59701 : Fix bug in jetTracksDefault HLO routine
  // 59700 : Z->bbbar group blocks merge; fix so empty events don't cause
  //       : crashes in met computation
  // 59601 : Fix bug in eastbhocc beam halo variable
  // 59600 : Add stubless muons; add "official" fidele information;
  //       : initial track-pointing-at-photon information; electron-to-jet
  //       : cross-reference; per-superlayer track hit count bitmap;
  //       : photon-to-jet cross-reference;
  //       : changes in photon timing code; add COT-only information for tracks;
  //       : add dof count for tracks; add small corrections to ces delta x
  //       : and delta z; add ces chisq fix for electrons; change electron
  //       : energies to two-tower sums; fix significant muon finding-
  //       : and track association-bugs; include Prereq module; use
  //       : processName talk-to for object finding;
  // 59503 : Fix bug in electron id word
  // 59502 : Fix stupid bug in muon met correction
  // 59501 : Add ability to use ZVertex for vertex block
  // 59500 : Add beam location info; add beam constrained track parameters;
  //       : add secvtx tagging (unstable); add d0 corrected for beam offset;
  //       : add corrections to met (unstable); add electron-photon
  //       : cross-reference; add electron conversion information;
  //       : add BMU muons; add individual PJW corrections; add more
  //       : photon beam halo variables; add Matt's tau variables;
  // 59415 : Add version information for production and offline
  // 59414 : Fix bugs in photon cosmic stub count (duplicate stubs) and
  //       : crashes in reverse tracking on pathologies
  // 59413 : 7/26/02 Fix bugs at end of jet block
  // 59412 : Ashutosh's reverse-tracking works!
  // 59411 : Made photon costub actually work; track extrapolated detector
  //       : eta is done from plug if track would intersect that plane before
  //       : the CES radius; added number of towers used in PEM 3x3 chisq fit
  //       : the CES radius; added number of towers used in PEM 3x3 chisq fit
  //       : to electron block; added "5 over 9" poor-man's PES chisquare to
  //       : electron block; photon CES strip association made only when matches
  //       : to centroid within 25 cm in z, and wire must be in same
  //       : half-chamber as centroid when centroid is more than 25 cm away
  //       : from break
  // 59410 : Fixed stupidity in tau block; added specific isolation
  //       : corrections for photons
  // 59409 : New tau code; electron PEM 3x3 chisq eta and phi added;
  //       : first attempt at integrating reverse-tracking cosmic filter
  // 59408 : Met is now vertexed from vertex.vzev
  // 59407 : electron detector eta, four-vectors fixed
  // 59406 : HEPG block added
  // 59405 : jet-to-vertex links made, anti-beam-halo variables added to
  //       : photons
  // 59404 : tau track links added
  // 59403 : electron likelihood added
  // 59402 : A whole bunch of (mostly plug) photon variables added
  // 59401 : Include Anyes Taffard's cosmic ray detection (muon.cosmic);
  //       : muon HADTDC timing; instantaneous luminosity;
  //       : per-track tof information
  // 59400 : Change to CES cluster for computing electron detector eta
  // 59301 : Turn on compression=2 for trigger block
  // 59300 : Overhaul trigger block, add tau and simulation code
  // 59205 : Get rid of fBits and fUniqueId in trigger branches
  // 59204 : Fix a few error logger severities
  // 59203 : Fix bounds checking for jets, tracks, muons, electrons, vertices
  // 59202 : Fill jet corrections
  // 59201 : Fix photon block crashes when too many EMObjects
  // 59200 : Bring back muon vertexing after incorporating private changes
  //        to 4.5.0 HLO muon code