# Author: Ketevi A. Assamagan 
# BNL, July 16, 2004 
# Modified by Theodota, 13.12.2004 
# job options for JPsi -> mm analysis example 
# Read in the AOD from POOL 
# this file is produced by running the AOD builders of the RDO raw data 
# The user is supposed to modify the input data and the flags as
# appropriate for his/her analysis

# import the analysis flags 
from AnalysisFlags import AnalysisFlags

# load common options and common DLL 
include( "AnalysisExamples/AnalysisCommon_jobOptions.py" )

# switch off ParticleJet, BJet and Photon 
DoMuon1 = True

# Print the flags before running the analysis code 
AnalysisFlags.Print()

# particle pre selection job options 
if DoMuon1:  
   include( "DoMuon_jobOptions1.py" )

# User analysis steering algorithm

AOD_ANA.members += [ "JPsi" ] 
JPsi = Algorithm( "JPsi" ) 
PreSelectedMuons = "MyPreSelectedMuons" 
JPsi.doJPsimm = TRUE 
NtupleLocID = "/FILE1/JPsi/700"
NtuplePrefix = "JPsi" 
OutputLevel = INFO
MuonPtCut = 3.0*GeV 
# LeptonEtaCut = 2.5

# The AOD input file 
InputCollections = [  
"AOD_JPsiset0.pool.root",  "AOD_JPsiset2.pool.root", "AOD_JPsiset4.pool.root",
"AOD_JPsiset5.pool.root",  "AOD_JPsiset7.pool.root",  "AOD_JPsiset8.pool.root", 
"AOD_JPsiset9.pool.root"  ]

# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL ) 
MessageSvc = Service( "MessageSvc" ) 
OutputLevel = INFO

# Number of Events to process 
theApp.EvtMax = 10000

# Root Ntuple output 
theApp.Dlls += [ "RootHistCnv" ] 
theApp.HistogramPersistency = "ROOT" 
NTupleSvc = Service( "NTupleSvc" ) 
NTupleSvc.Output = [ "FILE1 DATAFILE='JPsi.ntuple.root' OPT='NEW'" ] 
HistogramPersistencySvc = Service( "HistogramPersistencySvc" )
OutputFile = "JPsi.hist.root";

-- TheodotaLagouri - 06 Jan 2005