#------------------------------------------------------------------------------
# etaped.tcl
# for data
# Upsilon -> mu+mu- pi+ pi-
#------------------------------------------------------------------------------
# always source the error logger early in your main tcl script
sourceFoundFile ErrLogger/ErrLog.tcl
sourceFoundFile FrameScripts/FwkCfgVar.tcl
sourceFoundFile FrameScripts/talkto.tcl
# Disable the use of envvars
set ProdTclOnly true

# set the error logging level to 'warning'.  If you encounter a configuration
# error you can get more information using 'trace'
ErrLoggingLevel warning

## allowed values of BetaMiniReadPersistence are (currently) "Kan", "Bdb"
##
FwkCfgVar BetaMiniReadPersistence Kan

## allowed (non-expert) values of levelOfDetail are "micro", "cache", "extend" 
## or "refit"
##
FwkCfgVar levelOfDetail "cache"

## allowed values of ConfigPatch are "Run1", "Run2" or "MC".  This MUST be set 
## consistent ## with your input data type or you will get INCONSISTENT OR 
## INCORRECT RESULTS
##
#FwkCfgVar ConfigPatch "Run1"
FwkCfgVar ConfigPatch "Run2"
#FwkCfgVar ConfigPatch "MC"

##
## Set the number of events to run. If this isn't set, all events in the
## input collections will be processed.
##
#FwkCfgVar NEvent

##  input collections
echo UserTcl
if [ info exists env(UserTcl) ] {  
    echo source $env(UserTcl)      
      source ./$env(UserTcl)       
    } 

## choose the flavor of ntuple to write (hbook or root) and the file name
##
set BetaMiniTuple "hbook"
#set BetaMiniTuple "root"
set histFileName "$env(UserHbook)"

sourceFoundFile BetaMiniUser/btaMiniPhysics.tcl
sourceFoundFile CompositionSequences/CompPsiInitSequence.tcl
seq append BetaMiniPhysicsSequence -a PidSequence CompPsiInitSequence

talkto PidLHElectronsGTL {
   unrefinedListName set eMicroLoose
}

talkto eBremRecoPID {
   inputList0 set eMicroLoose
}


mod clone TagFilterByValue TagPi0
seq append BetaMiniReadSequence -a KanEventUpdateTag TagPi0
talkto TagPi0 {
  TagValueName  set nGoodTrkLoose
  TagValueType  set Integer
  TagValueMinimum set 2
  TagValueMaximum set 15
}

sequence create MyAnalysis

mod clone SmpMakerDefiner MyPsi
sequence append MyAnalysis MyPsi
catch { setProduction MyPsi }
talkto MyPsi {
  decayMode          set "eta' -> mu+ mu-"
  daughterListNames  set "muCombinedLoose"
  daughterListNames  set "muCombinedLoose"
  fittingAlgorithm   set "Add4"
  fitConstraints     set "PrimaryVertex"
  postFitSelectors   set "Mass 0.85:1.05"
}

mod clone SmpSubListerDefiner MyGamma                                
sequence append MyAnalysis MyGamma                                   
catch { setProduction MyGamma }                                      
talkto MyGamma {                                                     
   unrefinedListName  set "GoodPhotonLoose"                           
   selectors          set "LabE 0.050:10.0"                           
}                                          

mod clone SmpMakerDefiner MyEta2g                                      
sequence append MyAnalysis MyEta2g                                     
catch { setProduction MyEta2g }                                        
talkto MyEta2g {                                                       
   decayMode          set "eta -> gamma gamma"                        
   daughterListNames  set "MyGamma"                                   
   daughterListNames  set "MyGamma"                                   
   preFitSelectors    set "Mass 0.45:0.65"                          
   fittingAlgorithm   set "Add4"                                      
   fitConstraints     set "PrimaryVertex"                             
   postFitSelectors   set "Mass 0.5:0.6"                          
}

mod clone SmpRefitterDefiner MyEta2gMass                               
sequence append MyAnalysis MyEta2gMass                                 
catch { setProduction MyEta2gMass }                                    
talkto MyEta2gMass {                                                   
   unrefinedListName  set MyEta2g                                       
   fittingAlgorithm   set "Add4"                                      
   fitConstraints     set "Mass"                                      
   fitConstraints     set "Momentum"                                  
}                         
        
mod clone SmpMakerDefiner MyEtap                                     
sequence append MyAnalysis MyEtap                                    
catch { setProduction MyEtap }                                       
talkto MyEtap {                                                      
   decayMode          set "eta' -> pi+ pi- eta"                       
   daughterListNames  set "piCombinedTight"                           
   daughterListNames  set "piCombinedTight"                           
   daughterListNames  set "MyEta2gMass"                                 
   preFitSelectors    set "Mass 0.7:1.3"                          
   fittingAlgorithm   set "Add4"                                      
   fitConstraints     set "PrimaryVertex"                             
   postFitSelectors   set "Mass 0.85:1.05"                           
}                          



mod clone SmpMergerDefiner MyJpsiall
sequence append MyAnalysis MyJpsiall
catch { setProduction MyJpsiall }
talkto MyJpsiall {
   inputListNames    set "MyPsi"
   inputListNames    set "MyEtap"
   disableCloneCheck set true
}

      
path append Everything MyAnalysis

#..Use BtuTupleMaker to write out ntuples for SimpleComposition job
path append Everything BtuTupleMaker
talkto BtuTupleMaker {
  #writeEveryEvent set true
  listToDump set MyJpsiall

  fillMC set true
  mcBlockContents    set "Mass CMMomentum Momentum" 
  
  eventBlockContents set "EventID CMp4"
  eventTagsInt       set "nTracks nGoodTrkLoose"
  eventTagsBool      set "BGFMultiHadron BGFNeutralHadron BGFMuMu BGFTau BGFTwoProng BGFPhiGamma BGFAllNeutralTwoPhoton BGFIsr BGFRadTwoProng BGFHighMassHadron BGFTwoPhotonTwoTrack DigiFDchEmcPreVeto DigiFL1Open DigiFL3Open L3OutDch L3OutEmc"

  ntpBlockConfigs set "eta' etaprime   3  50"
  ntpBlockConfigs set "mu+    mu   0  50"
  ntpBlockConfigs set "eta  eta   2  50"
  ntpBlockConfigs set "pi+  pion   0  50"

  ntpBlockConfigs set "gamma gamma 0  50"
  ntpBlockToTrk   set "mu"
  ntpBlockToTrk   set "pion"

  ntpBlockContents set "etaprime:  Mass Momentum CMMomentum"
  ntpBlockContents set "mu:  Momentum CMMomentum DocaXY"
  ntpBlockContents set "eta:  Momentum CMMomentum DocaXY"
  ntpAuxListContents set "eta  : MyEta2g : r : Mass"
  ntpBlockContents set "pion:  Momentum CMMomentum DocaXY"
  fillAllCandsInList set "TRK ChargedTracks"
  ntpBlockContents set "TRK  : Momentum CMMomentum DocaXY"
  fillAllCandsInList set "gamma CalorNeutral"
  ntpBlockContents set "gamma  : Momentum CMMomentum"
  trkExtraContents set BitMap:piSelectorsMap,eSelectorsMap,KSelectorsMap,muSelectorsMap
}

echo mod talk evtcounter
mod talk EvtCounter
 printFreq set 10000
exit

action enable HistoDirAction
action list

echo Starting event loop

path list

#ev begin -nev 100000
ev begin  
exit 
#ErrMsg trace "completed OK"
#exit