#------------------------------------------------------------------------------
# etaped_2pi0.tcl
# for data
# Upsilon -> mu+mu- pi0 p0
#------------------------------------------------------------------------------
# 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 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 6
}

sequence create MyAnalysis

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 Mypi02g                                      
sequence append MyAnalysis Mypi02g                                     
catch { setProduction Mypi02g }                                        
talkto Mypi02g {                                                       
   decayMode          set "pi0 -> gamma gamma"                        
   daughterListNames  set "MyGamma"                                   
   daughterListNames  set "MyGamma"                                   
   preFitSelectors    set "Mass 0.09:0.17"                          
   fittingAlgorithm   set "Add4"                                      
   fitConstraints     set "PrimaryVertex"                             
   postFitSelectors   set "Mass 0.09:0.17"                          
}
mod clone SmpRefitterDefiner Mypi02gMass                               
sequence append MyAnalysis Mypi02gMass                                 
catch { setProduction Mypi02gMass }                                    
talkto Mypi02gMass {                                                   
   unrefinedListName  set Mypi02g                                       
   fittingAlgorithm   set "Add4"                                      
   fitConstraints     set "Mass"                                      
   fitConstraints     set "Momentum"                                  
}  

mod clone SmpRefitterDefiner Mypi02gMass_1
sequence append MyAnalysis Mypi02gMass_1
catch { setProduction Mypi02gMass_1 }
talkto Mypi02gMass_1 {
   unrefinedListName  set Mypi02g
   fittingAlgorithm   set "Add4"
}


mod clone SmpMakerDefiner MyPsi
sequence append MyAnalysis MyPsi
catch { setProduction MyPsi }
talkto MyPsi {
  decayMode          set "J/psi -> mu+ mu-"
  daughterListNames  set "muCombinedVeryLoose"
  daughterListNames  set "muCombinedVeryLoose"
  fittingAlgorithm   set "Add4"
  fitConstraints     set "PrimaryVertex"
  postFitSelectors   set "Mass 6.0:12.0"
}


mod clone SmpRefitterDefiner MyPsiMass                               
sequence append MyAnalysis MyPsiMass                                 
catch { setProduction MyPsiMass }                                    
talkto MyPsiMass {                                                   
   unrefinedListName  set MyPsi                                       
   fittingAlgorithm   set "Add4"                               
}  


mod clone SmpMakerDefiner MyPhi
sequence append MyAnalysis MyPhi
catch { setProduction MyPhi }
talkto MyPhi {
  decayMode          set "Upsilon(3S) -> J/psi pi0 pi0"
  daughterListNames  set "MyPsiMass"
  daughterListNames  set "Mypi02gMass"
  daughterListNames  set "Mypi02gMass_1"
  fittingAlgorithm   set "Add4"
  fitConstraints     set "PrimaryVertex"
  postFitSelectors   set "Mass 9.5:11.0"
}

mod clone SmpRefitterDefiner MyPhiMass
sequence append MyAnalysis MyPhiMass
catch { setProduction MyPhiMass }
talkto MyPhiMass {
   unrefinedListName  set MyPhi
   fittingAlgorithm   set "Add4"
   fitConstraints     set "Mass"
   fitConstraints     set "Momentum"
}





      
path append Everything MyAnalysis

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

  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 "Upsilon(3S) phi 3 50"
  ntpBlockConfigs set "pi0 pi0 2 50"
  ntpBlockConfigs set "J/psi psi   2  50"
  ntpBlockConfigs set "mu+    mu    1  50"
  ntpBlockConfigs set "gamma gamma 0  50"

  ntpBlockToTrk   set "mu"
  
  ntpBlockContents set "phi:  Mass Momentum CMMomentum"
  ntpAuxListContents    set "phi : MyPhi : r : Mass"
  ntpBlockContents set "pi0:  Mass Momentum CMMomentum"
  ntpAuxListContents    set "pi0 : Mypi02g : r : Mass"
  ntpBlockContents set "psi:  Mass Momentum CMMomentum"
  ntpBlockContents set "mu:  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 10000
ev begin  
exit 
#ErrMsg trace "completed OK"
#exit