#-------------------------------------------------------------- #USER - you need to set these switches!! IMPORTANT 12/18/2002 #-------------------------------------------------------------- # beamoffsetMC will set correctly the beam position talk-to for TopEventModule, TopNtuple and SecVtxModule # Set to 2 for MC with no beam offset # Set to 1 for MC with realistic beam offset in 4.9.1 # Set to 0 for data set beamoffsetMC 0 #To run on caf, comment two lines below #Change $INPUTFILE to $env(INPUTFILE) and $PCDFDIR to $env(PCDFDIR) #make sure you export these symbols in your caf script like this! #export PCDFDIR=$PROJECT_DIR #export INPUTFILE=${ifile} set INPUTFILE GI0749.1 #set PCDFDIR /home/cdfsoft/dist/releases/4.11.1 creator set PROD module input FileInput talk FileInput input add $env(topfind_input_file) luminosityLog set lumi.log # # Report on every event processed in validation # statusFile set 1 #will remake these when we run VxPrim dropList set VertexColl:DEFAULT_VXPRIM_VERTICES show exit ##module input DHInput ##talk DHInput ## ## cache set DCACHE ### include dataset ttop2i book=cdfptop ### include file root://fcdfdata015.fnal.gov//cdf/scratch/cdfdata/top/results/Inclusive-ele_4104_REMAKE_$INPUTFILE.* ## include file ./$env(topfind_input_file) ## statusFile set 100 ## ###Limit set by beamline availability ## selectEvents set run<=163527 ## ###will remake these when we run VxPrim ## dropList set VertexColl:DEFAULT_VXPRIM_VERTICES ###Phoenix electrons wrong in inclusive data files - so drop and redo (enable Phoenix modules) ### ok in lepton, W, Z files so keep em (disable Phoenix module) ### dropList add CdfTrackColl:Phoenix_Tracking CdfTrackColl:PhoenixSI_Tracking ### dropList add CdfEmObjectColl:Phoenix_Electrons ## ## show ##exit #--------------------------- # Puff stuff #--------------------------- talk PuffModule #Needed for BTagging puffOnly set SiHitSet show exit set HISTFILE phoenix_test2.root if [file exists $HISTFILE] { exec rm $HISTFILE } mod talk HepRootManager histfile set $HISTFILE createHistoFile set false show exit #---------------------------- #USED-SETs for beamline only #---------------------------- module talk CalibrationManager if {$beamoffsetMC ==0} then { #FOR DATA - USED-SET beamlines ProcessName set PAD_PHYSICS_CDF PassName set 09 } if {$beamoffsetMC !=0} then { #For official 4.9.1 MC ProcessName set SIM_PHYSICS_CDF } show exit #--------------------------- #TopEventModule - Run II analysis cuts #--------------------------- module talk TopEventModule printLevel set 0 parmSetName set TopEvent #default is to look for most recent beamline #TopFilter, TopNtuple and SecVtxModule also need to know the prefix in their talk-to's newTopCollectionPrefix set WZ leptonIdCuts set 2 eventZPosition set 2 jets calDataProcessName set PROD #the default reclusterJetAlgorithmList set JETCLU reclusterJetConeRadiusList set 0.4 #the extras - add as many as you like! reclusterJetAlgorithmList add JETCLU reclusterJetConeRadiusList add 0.7 reclusterJetAlgorithmList add KTCLUS reclusterJetConeRadiusList add 0.5 reclusterJetAlgorithmList add MIDPOINT reclusterJetConeRadiusList add 0.4 exit #For MC with beam offset if {$beamoffsetMC ==1} then { BeamlineFromDB set false mcBeamXPos set -0.064 mcBeamYPos set 0.310 mcBeamXSlope set -0.00021 mcBeamYSlope set 0.00031 mcBeamXWidth set 0.0025 mcBeamYWidth set 0.0025 } show help exit #--------------------------- #TopFilter - needs to know prefix #This keeps events with 1 or more CEM leptons and RAW MET>20.0 #--------------------------- module talk TopFilter printLevel set 0 newTopCollectionPrefix set WZ nTightElectrons set 1 MissingEt set 20.0 #OR of region requirements nTightCEM set 1 # nTightCMUP set 1 # nTightCMX set 1 show exit #--------------------------- #SecVtxModule - Set the parameters for the tagging # algsToKill set 16 means kill Si standalone tracks #--------------------------- mod talk SecVtxModule main algsToKill set 16 #TopEventModule info jetcol set WZ_JETCLU0.4_ReclusterJetColl getZfrom set TOPMOD_JETZVERT topSummaryDesc set WZSummary #Normal info # beamHistCode set 19 jetetcut set 1.0 jetetacut set 10.0 #For MC with offset beam-spot only if {$beamoffsetMC ==1} then { BeamlineFromDB set false sigma_x set 0.0025 sigma_y set 0.0025 sigma_z set 28.0 pv_central_x set -0.064 pv_central_y set 0.310 pv_central_z set 2.5 pv_slope_dxdz set -0.00021 pv_slope_dydz set 0.00031 } if {$beamoffsetMC ==2} then { BeamlineFromDB set false } show exit exit #--------------------------- #VxPrim #--------------------------- module talk VxPrim verbose set f show exit #--------------------------- #JetProbModule #--------------------------- module talk JetProbModule verbose set f CutsMenu Jet_EtCut set 7.0 Jet_EtaRange set 2.8 ConeSize set 0.4 TrkPtMin set 0.5 MinNCotHitInAxSL set 6 MinNCotHitInStSL set 6 MinNumCotAxSL set 2 MinNumCotStSL set 2 MinNumCotAxHit set 20 MinNumCotStHit set 17 TrkD0Max set 0.15 TrkSvxChi2Max set 10000.0 MinNumRPhiHitLayer set 3 MinNumRPhiHitLayer_SVXII set 3 TrkDeltaZ0Max set 5.0 TrkZ0Max set 1.0 MinNum90DegHitLayer set 1 MinNumSmAngHitLayer set 0 TrkProcessName set "PROD" UseVxPrim set t # JetCollDescription set "JetCluModule-cone0.4" JetCollDescription set "WZ_JETCLU0.4_ReclusterJetColl" JetProcessName set "PROD" # JetProcessName set "" VertexCollDescription set "DEFAULT_VXPRIM_VERTICES" TrkMatchForShrHit set t RunSiExpect set t # # This is only temporary TrkWithNoL00Hit set t exit show exit #--------------------------- #TopNtuple - needs to know prefix and reco options #Saves TopEventModule collections and default electrons and muons #--------------------------- module talk TopNtuple #default is to look for most recent beamline printLevel set 0 makeNtuple set 1 # TopNtuplefile set ntuples/W-ele_4104_REMAKE_ntuple_$INPUTFILE.root TopNtuplefile set $env(topfind_output_root) #For MC with beam offset if {$beamoffsetMC ==1} then { mcBeamXPos set -0.064 mcBeamYPos set 0.310 mcBeamXSlope set -0.00021 mcBeamYSlope set 0.00031 mcBeamXWidth set 0.0025 mcBeamYWidth set 0.0025 } #Must be set to same values as for TopEventModule #Below are values used for leptonIdCuts set 2 = Run II analysis cuts eleUseCotParent set true eleNumberOfTowers set 2 eleBeamConstraint set true eleCemCorr set true muonBeamConstraint set true muonUseCotParent set true #Describe collections/views you want to read in #Replace prefix "WZ" with value of newTopCollectionPrefix in TopEventModule talk-to eleList set WZTightEmObjectView eleList add WZLooseEmObjectView eleList add defEmObjects eleList add PhoenixEmObjectView muoList set WZTightMuonView muoList add WZLooseMuonView muoList add defMuons muoList add defStublessMuons jetList set WZTightJetView jetList add WZLooseJetView jetList add WZ_JETCLU0.4_ReclusterJetColl jetList add WZ_JETCLU0.7_ReclusterJetColl jetList add WZ_KTCLUS0.5_ReclusterJetColl jetList add WZ_MIDPOINT0.4_ReclusterJetColl metDescription set WZMetCorrected summaryDescription set WZSummary otrkList set defTracks otrkList add Phoenix_Tracking otrkList add PhoenixSI_Tracking # # describe blocks you want to fill getElectron set 1 getMuon set 1 getJet set 1 getMet set 1 getOfflTrack set 1 getTopSummary set 1 getTopFit set 1 getSecvtxTrack set 1 getTrackJetAssoc set 1 getJetProbTrack set 1 getJetProbTrackJetAssoc set 1 getZVertices set 1 getPrimaryVertex set 1 getHepg set 1 getOBSP set 1 getCotHit set 0 getXFTHit set 0 getXFTPixel set 0 getXFTTrack set 1 getTrigger set 1 getTriggerName set 1 show exit #************************ #Only if these don't exist already talk TrackRefitModule dropL00 set t show exit #------------------------ #Cosmic Ray stuff #------------------------ module talk TofManager CalibratorType set Test exit module talk CT_TrackingModule NoTracking set t exit module talk CosmicFinderModule CalorimetryMenu JetCollName set JetCluModule-cone0.4 exit show exit #****Phoenix stuff starts here**** #======================= # Set up PhoenixTracking #======================= module talk PhoenixTracking debugLevel set 0 #Avoid making histograms createHistoFile set false makeHistograms set false makeNtuples set false minEt set 15.0 maxEta set 2.8 # beamHistCode set 19 if {$beamoffsetMC ==0} then { allowNoBeam set f BeamlineFromDB set t plugSubMenu #Chill noted these should be 1.0 now #June 17 - fixed PhoenixTracking talk-to for 4.11.1 -no manual alignment needed wPlugEscale set 1.0 ePlugEscale set 1.0 exit } if {$beamoffsetMC ==1} then { allowNoBeam set f BeamlineFromDB set false sigma_x set 0.0025 sigma_y set 0.0025 sigma_z set 28.0 pv_central_x set -0.064 pv_central_y set 0.310 pv_central_z set 2.5 pv_slope_dxdz set -0.00021 pv_slope_dydz set 0.00031 } if {$beamoffsetMC ==2} then { allowNoBeam set t BeamlineFromDB set false } show help exit #=============================================== #Set up SiPatternRecModule for PhoenixTracking #=============================================== module talk SiPatternRecModule #== no L00 for now missingLayer set 0 PerformKalSvxStandaloneTracking set false PerformKalOISvxStandaloneTracking set false method set GLOBAL inputTracks set Phoenix_Tracking outputTracks set PhoenixSI_Tracking show exit module talk PhoenixElectrons #Avoid making histograms and needing to enable HepRootManager createHistoFile set false show exit #****Phoenix stuff ends here**** #--------------------------- # Turn on/off modules #--------------------------- module disable ConfigManager module disable Prereq #Cosmic module disable TofManager module disable CT_TrackingModule module disable CosmicFinderModule module enable HepRootManager #Must not run TrackRefitModule over REMAKE data since it has already been done if {$beamoffsetMC ==0} then { module disable TrackRefitModule #phoenix module disable PhoenixTracking module disable SiPatternRecModule module disable PhoenixElectrons } #Must run TrackRefitModule over MC if {$beamoffsetMC !=0} then { module enable TrackRefitModule #phoenix module enable PhoenixTracking module enable SiPatternRecModule module enable PhoenixElectrons } #topevent module enable TopEventModule module disable TopFilter #silicon module enable SecVtxModule module enable VxPrim module enable JetProbModule #topntuple module enable TopNtuple #--------------------------- # setup the path #--------------------------- path disable AllPath path create Path-TopAna ManagerSequence \ VxPrim \ TofManager CT_TrackingModule CosmicFinderModule \ PhoenixTracking SiPatternRecModule PhoenixElectrons \ TopEventModule TopFilter \ TrackRefitModule \ SecVtxModule \ JetProbModule \ TopNtuple path enable Path-TopAna #--------------------------- #Comment this section if you only want the TopNtuple and not the output data file #--------------------------- source $env(SRT_PUBLIC_CONTEXT)/Production/defaultOutputKeepList.tcl talk FileOutput parmSetName set TopSel output create Stream-TopAna $env(topfind_output_file) output path Stream-TopAna Path-TopAna defaultKeepList Stream-TopAna output keepList Stream-TopAna -add CdfEmObjectView CdfMuonView CdfJetView CdfJetColl CdfMet \ TopSummary FastZVertexColl \ TopMassFitSummary SecVtxColl SecVtxTracks CosmicRayInfo CdfTrackColl JetProbColl Stream-TopAna compression set true fileSize set 500000 exit show exit show # Break ErrorLoggerManager's Trappist vows! # Make ErrorLoggerManager talk. INFO level messages are suppressed by default # We want to read INFO level messages, especially for short tests! #talk ErrorLoggerManager # severity set SUCCESS # show #exit # # talk ErrorLoggerManager # remoteDest # localFile set myerror.log # threshold set INCIDENTAL # exit # exit # # process all the event from the input file # # begin -nev 50 begin # # show timing etc in the end of job # show exit