#! /bin/bash

date

#this is for the 2jet
signal="tb"
lepton="muon"
Lepton="muon"
NJets="2"
NTag="1"
for type in training ; do #small_training training testing yield ; do
for reco in p20 p17
do
for lepton in electron muon
do
for NJets in 2 3
do
for btag in 1 2
do
for sample in data # dilepton #lepjets zbb zcc zjj tb tq wbb wcc wjj diboson data qcd
do

#############################################################	
if [ $reco == "p17" ] ; then
    if [ $lepton == "electron" ] ; then
	if [ $NJets == "2" ] ; then
	    #p17 btag 1 2 jet electron
	    if [ $btag == "1" ] ; then
		wbb="159.96017" #"87.8644+26.1624"
		wcc="81.52595" #"43.9907+13.2086"
		wlp="182.35866" #"130.531+9.7195"
		qcd="34.4726" #"408"
		diboson="15.307" #"24.8836"
		ttbar="35.3638" #"23.4643+11.0211"
		tb="7.7957"
		tqb="16.1625"
		scalefactor="1.0"
	    #p17 btag 2 2 jet electron
	    else
		wbb="39.209604" #"11.7348+3.18438"
		wcc="5.8471635" #".915272+0.293509"
		wlp="2.9307568" #"0.27976+0.0111631"
		qcd="2.31181" #"8"
		diboson="1.87994"
		ttbar="16.92314" #"8.04585+1.82679"
		tb="5.28735"
		tqb="0.743361"
		
		scalefactor="1.0"
		
	    fi  
	else
	    #p17 btag 1 3 jet electron
	    if [ $btag == "1" ] ; then
		wbb="42.22139" #"87.8644+26.1624"
		wcc="25.06143" #"43.9907+13.2086"
		wlp="38.435908" #"130.531+9.7195"
		qcd="11.4333" #"408"
		diboson="4.81285" #"24.8836"
		ttbar="56.2416" #3.4643+11.0211"
		tb="2.56538"
		tqb="5.88742"
		scalefactor="1.0"
		acc_tqb=".6"
		acc_tqg=".4"
	    #p17 btag 2 3 jet electron
	    else
		wbb="11.806454" #"11.7348+3.18438"
		wcc="2.6030023" #".915272+0.293509"
		wlp="1.2370009" #"0.27976+0.0111631"
		qcd="1.12969" #"8"
		diboson="0.748627"
		ttbar="34.68119" #"8.04585+1.82679"
		tb="1.89689"
		tqb="1.43792"
		scalefactor="1.0"
		acc_tqb="1"
		acc_tqg="0"
	    fi
	fi
    else
    	if [ $NJets == "2" ] ; then
	    #p17 btag 1 2 jet muon
	    if [ $btag == "1" ] ; then
		wbb="175.2725" #"87.8644+26.1624"
		wcc="92.78919" #"43.9907+13.2086"
		wlp="198.98112" #"130.531+9.7195"
		qcd="54.8405" #"408"
		diboson="20.407" #"24.8836"
		ttbar="29.7739" #"23.4643+11.0211"
		tb="8.50845"
		tqb="18.6396"

		scalefactor="1.0"
		
	    #p17 btag 2 2 jet muon
	    else
		wbb="41.73643" #"11.7348+3.18438"
		wcc="6.563912" #".915272+0.293509"
		wlp="3.139028" #"0.27976+0.0111631"
		qcd="3.4136" #"8"
		diboson="2.72281"
		ttbar="15.90847" #"8.04585+1.82679"
		tb="5.95006"
		tqb="0.831591"
		
		scalefactor="1.0"
		
	    fi
	else
	    #p17 btag 1 3 jet muon
	    if [ $btag == "1" ] ; then
		wbb="68.84762" #"87.8644+26.1624"
		wcc="41.74674" #"43.9907+13.2086"
		wlp="60.31536" #"130.531+9.7195"
		qcd="18.6177" #"408"
		diboson="6.95242" #"24.8836"
		ttbar="58.397" #"23.4643+11.0211"
		tb="3.14891"
		tqb="7.44449"
		acc_tqb=".6"
		acc_tqg=".4"
		scalefactor="1.0"
		
	    #p17 btag 2 3 jet muon
	    else
		wbb="18.95941" #"11.7348+3.18438"
		wcc="4.168158" #".915272+0.293509"
		wlp="1.954356" #"0.27976+0.0111631"
		qcd="2.0661" #"8"
		diboson="1.17351"
		ttbar="38.47556" #"8.04585+1.82679"
		tb="2.37833"
		tqb="1.87256"
		scalefactor="1.0"
		acc_tqb="1"
		acc_tqg="0"
	    fi
	fi  

    fi
#############################################################	
elif [ $reco == "p20" ] ; then
    if [ $lepton == "electron" ] ; then
	if [ $NJets == "2" ] ; then
	    #p20 btag 1 2 jet electron
	    if [ $btag == "1" ] ; then
		wbb="147.37387" #"87.8644+26.1624"
		wcc="81.7291" #"43.9907+13.2086"
		wlp="235.23745" #"130.531+9.7195"
		qcd="28.0414" #"408"
		diboson="17.7964" #"24.8836"
		ttbar="39.8798" #"23.4643+11.0211"
		tb="8.58333"
		tqb="16.7317"
		scalefactor="1.0"
	    #p20 btag 2 2 jet electron
	    else
		wbb="34.194697" #"11.7348+3.18438"
		wcc="6.119026" #".915272+0.293509"
		wlp="4.962038" #"0.27976+0.0111631"
		qcd="1.92212" #"8"
		diboson="1.92212"
		ttbar="17.56971" #"8.04585+1.82679"
		tb="5.15753"
		tqb="0.920617"
		
		scalefactor="1.0"
		
	    fi  
	else
	    #p20 btag 1 3 jet electron
	    if [ $btag == "1" ] ; then
		wbb="51.82264" #"87.8644+26.1624"
		wcc="32.30146" #"43.9907+13.2086"
		wlp="68.26887" #"130.531+9.7195"
		qcd="11.3622" #"408"
		diboson="5.74098" #"24.8836"
		ttbar="63.7187" #3.4643+11.0211"
		tb="2.98619"
		tqb="6.95447"
		scalefactor="1.0"
		acc_tqb=".6"
		acc_tqg=".4"
	    #p20 btag 2 3 jet electron
	    else
		wbb="14.370442" #"11.7348+3.18438"
		wcc="3.633907" #".915272+0.293509"
		wlp="2.838795" #"0.27976+0.0111631"
		qcd="1.4282" #"8"
		diboson="0.877271"
		ttbar="36.56662" #"8.04585+1.82679"
		tb="1.97043"
		tqb="1.83184"
		scalefactor="1.0"
		acc_tqb="1"
		acc_tqg="0"
	    fi
	fi
    else
    	if [ $NJets == "2" ] ; then
	    #p20 btag 1 2 jet muon
	    if [ $btag == "1" ] ; then
		wbb="171.5033" #"87.8644+26.1624"
		wcc="80.6821" #"43.9907+13.2086"
		wlp="268.5796" #"130.531+9.7195"
		qcd="78.897" #"408"
		diboson="23.2336" #"24.8836"
		ttbar="33.805" #"23.4643+11.0211"
		tb="9.12762"
		tqb="18.4177"

		scalefactor="1.0"
		
	    #p20 btag 2 2 jet muon
	    else
		wbb="38.7147" #"11.7348+3.18438"
		wcc="7.969853" #".915272+0.293509"
		wlp="5.632883" #"0.27976+0.0111631"
		qcd="5.0165" #"8"
		diboson="3.00602"
		ttbar="16.08076" #"8.04585+1.82679"
		tb="5.55393"
		tqb="0.974043"
		
		scalefactor="1.0"
		
	    fi
	else
	    #p20 btag 1 3 jet muon
	    if [ $btag == "1" ] ; then
		wbb="68.19498" #"87.8644+26.1624"
		wcc="38.2244" #"43.9907+13.2086"
		wlp="88.73736" #"130.531+9.7195"
		qcd="30.0048" #"408"
		diboson="8.46593" #"24.8836"
		ttbar="64.8643" #"23.4643+11.0211"
		tb="3.94465"
		tqb="8.52248"
		acc_tqb=".6"
		acc_tqg=".4"
		scalefactor="1.0"
		
	    #p20 btag 2 3 jet muon
	    else
		wbb="18.7391" #"11.7348+3.18438"
		wcc="5.475645" #".915272+0.293509"
		wlp="3.665287" #"0.27976+0.0111631"
		qcd="3.0324" #"8"
		diboson="1.355"
		ttbar="39.1317" #"8.04585+1.82679"
		tb="2.59354"
		tqb="2.17447"
		scalefactor="1.0"
		acc_tqb="1"
		acc_tqg="0"
	    fi
	fi  

    fi
fi

    scalefactor=1.0

    if [ $lepton == "electron" ] ; then
	Lep="CC"
	lep="cc"
	LEP="CC"
	lepno="1"
    else
	Lep="mu"
	lep="mu"
	LEP="MU"
	lepno="2"
	
    fi
    
    if [ $NJets == "3" ] ; then
	TopoNJet="EqThreeJet"
    else
	TopoNJet="EqTwoJet"
    fi

   
    if [ $btag == "1" ] ; then
	tagname="EqOneTag"
    else
	tagname="EqTwoTag"
    fi	
    
    if [ $sample == "QCD" ] ; then
	signal="qcd"
    elif [ $sample == "DATA" ] ; then
	signal="data"
    elif [ $sample == "wjj" ] ; then
        signal="wlp"
    elif [ $sample == "zjj" ] ; then
        signal="zlp"
    elif [ $sample == "tq" ] ; then
        signal="tqb"
    else
	signal=$sample
    fi    

    METopo="/rooms/ds3/meena/SingleTOP/TopoVars/sep14/default_disc/output_root_catall_sep11/${type}_Btagging/${sample}_${reco}_${lepton}_NJets${NJets}_Btagging_${btag}tag.root"

#    mkdir -p  /rooms/ds3/meena/SingleTOP/TopoVars/sep16/${type}

    outputfilename="/rooms/ds3/meena/SingleTOP/TopoVars/sep16/${type}/${reco}_${Lep}_${signal}_${tagname}_${TopoNJet}_zero.root"
#    outputfilename="${type}/${reco}_${Lep}_${signal}_${tagname}_${TopoNJet}_zero.root"

    echo root -l -b -q ''run_opt_disc.C'("'${METopo}'","'${outputfilename}'",'${NJets}','${wbb}','${wcc}','${wlp}','${qcd}','${diboson}','${ttbar}','${tb}','${tqb}','${scalefactor}')'
    
    root -l -b -q ''run_opt_disc.C'("'${METopo}'","'${outputfilename}'",'${NJets}','${wbb}','${wcc}','${wlp}','${qcd}','${diboson}','${ttbar}','${tb}','${tqb}','${scalefactor}')'

done
done
done
done
done
done

echo "finished with all sets!"
date