**************************************************************

EXAMPLE TO GENERATE B^0 -> Jpsi K* MC:

Control cards:
/work/naia-clued0/bruce1/Jpsi_phi/MC/BH_bb_Bd-jpsi-k0star_bertram.py
/work/naia-clued0/bruce1/Jpsi_phi/MC/user.dec_b0-jpsi-k0star
/work/naia-clued0/bruce1/Jpsi_phi/MC/d0_mess_b0-jpsi-k0star.rcp

cd /work/naia-clued0/bruce1/Jpsi_phi/MC
setup D0RunII p17.09.06
setup mc_runjob
d0setwa
setenv MC_RUNJOB_CARD_PATH /work/naia-clued0/bruce1/Jpsi_phi/MC
mc_runjob -macro=BH_bb_Bd-jpsi-k0star_bertram.py

To convert to AADST's:

In naia-clued0:/work/naia-clued0/bruce1/Jpsi_phi/MC/AAProd2/p17.09.06-prod>
ls /prj_root/3001/ckm_write/bruce1/MC_Bd_Jpsi_K0star/d0reco/* > Bd_Jpsi_K0star.list
> setup D0RunII p17.09.06
> setup sam
> srt_setup SRT_QUAL=maxopt
> setup d0tools -t
> d0setwa
runreco -localfwkrcp -localbuild -localrcp -maxopt -format=anaTMB -filelist=Bd_Jpsi_K0star.list -fwkparams
 -tracking htfAA

cp d0reco/tmb_p17.09.06_NumEv-0_b+b-incl_sm.n_Bphysics_00001_555915727352988407209114405
/prj_root/3001/ckm_write/bruce1/MC_Bd_Jpsi_K0star/d0reco

To analyze, go to Jpsi_phi_20>.
/work/naia-clued0/bruce1/Jpsi_phi/MC/AAProd2/p17.09.06-prod/D0reco_x-p17.09.06-maxopt-Linux-Bd_Jpsi_K0star.list/dstOut_000
 > files

events = Events(2,2030,100)
so 2030 pythia events generated to get 100 through d0mess
*********************************************************************************
EXAMPLE TO PRODUCE AADST'S FROM MC (OR DATA) IN SAM:

PRODUCE AADST'S FOR ID 83892
B+ ->Jpsi K pipi

cd /work/naia-clued0/bruce1/Jpsi_phi/p20.09.02-prod
0.  setup D0RunII p20.09.02
    kinit -f
    setup d0cvs
    setup sam
    setup d0tools

sam translate constraints --dim="GLOBAL.REQUESTID 83892 and data_tier thumbnail"
sam translate constraints --dim="GLOBAL.REQUESTID 83892 and data_tier thumbnail and file_name\
 tmb_p20.09.02_NumEv-7500_b+b-incl_sm.n_mcc99_mcp20_osg-gw-2.t2.ucsd.edu_83892-12007323407445-JIM_MERGE_0"

sam create dataset definition --defname=MCrequestID83892_12007323407445 \
--group=dzero --dim="GLOBAL.REQUESTID 83892 \
and data_tier thumbnail and file_name \
tmb_p20.09.02_NumEv-7500_b+b-incl_sm.n_mcc99_mcp20_osg-gw-2.t2.ucsd.edu_83892-12007323407445-JIM_MERGE_0"
DatasetDefinition saved with definitionId = 1192735

    runreco -localfwkrcp -localbuild -localrcp -maxopt \
    -format=ana -userpostscript=post.csh \
    -initscript=initchunkanalyze.sh -sourcescript=chunkanalyzesource.sh \
    -defname=MCrequestID83892_12007323407445 -group=dzero -nofpe \
    -cabsrv2 -cabtar=tarball_20.09.02_2.tar -jobs=1 \
    -scratch=/prj_root/3001/ckm_write/bruce1/scratch \
    -caboutpath=/prj_root/3001/ckm_write/bruce1/B0s_MC

*************************************************

Pointers to location of MC generation files for dimu-asymmetry
are in D0 note 4996, page 9:

run-mc_bck100604
bruce.macro_bck100604
bruce.rcp_bck100604
runTMB_bck100604

in

/work/naia-clued0/bruce1/dimu_asym_2004/MC

***************************************************

copied from
/D0/dist/packages/d0_mess/devel/src/RCPReader.cpp :

     // register values with cut maker
     // momentum/energy cuts
     _cutMaker.registerDoubleValue("Pt", _ptValue) ;
     _cutMaker.registerDoubleValue("E", _eValue) ;
     // phi cuts
     _cutMaker.registerDoubleValue("Phi", _phiValue) ;
     // eta cuts
     _cutMaker.registerDoubleValue("Eta", _etaValue) ;
     _cutMaker.registerDoubleValue("AbsEta", _absetaValue) ;
     // Id cuts
     _cutMaker.registerIntValue("PdgId", _pdgValue) ;
     _cutMaker.registerIntValue("AbsPdgId", _abspdgValue) ;
     _cutMaker.registerIntValue("ParentId", _parentIdValue) ;
     _cutMaker.registerIntValue("FlavorId", _flavorIdValue) ;
     _cutMaker.registerIntValue("AbsFlavorId", _absflavorIdValue) ;
     _cutMaker.registerIntValue("GrandParentId", _gparentIdValue) ;
     _cutMaker.registerIntValue("ParentMesonFlavor",_parMesFlavValue) ;
     _cutMaker.registerIntValue("ParentBaryonFlavor",_parBarFlavValue) ;
     _cutMaker.registerIntValue("ParentHadronFlavor",_parHadFlavValue) ;
     _cutMaker.registerIntValue("ParentBosonId",_parBosIdValue) ;
     _cutMaker.registerIntValue("AbsParentId",_absParentIdValue) ;
     _cutMaker.registerIntValue("AbsGrandParentId", _absgparentIdValue) ;

***********************************************