Skip to content

Commit

Permalink
Merge branch 'master' of https://github.com/cms-edbr/ExoDiBosonResona…
Browse files Browse the repository at this point in the history 
…nces
  • Loading branch information
@santanas
santanas committed Oct 14, 2013
2 parents 93d770c + ad8ed1b commit 6165f8f
Show file tree
Hide file tree
Showing 3 changed files with 77 additions and 32 deletions.
68 changes: 42 additions & 26 deletions EDBRCommon/test/analyzerEDBR_zz2l2j_tpl.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -40,32 +40,41 @@
)


if "WW"=="<ANALYSIS>" :
process.ANEDBR.VType = cms.string("W")
else :
process.ANEDBR.VType = cms.string("Z")



### if false, use the default collections
### in ExoDiBosonResonances.EDBRCommon.analyzerEDBR_cfi
### (i.e. all the cands passing pre-selection cuts)
processFullSel=True

if processFullSel :
# process.ANEDBR.EDBREEJJColl=cms.InputTag("BestSidebandSelectorEle:doubleJet")
# process.ANEDBR.EDBRMMJJColl=cms.InputTag("BestSidebandSelectorMu:doubleJet")
# process.ANEDBR.EDBREEJColl=cms.InputTag("BestSidebandSelectorEle:singleJet")
# process.ANEDBR.EDBRMMJColl=cms.InputTag("BestSidebandSelectorMu:singleJet")
# process.ANEDBR.EDBREEJJColl=cms.InputTag("BestCandSelectorEle:doubleJet")
# process.ANEDBR.EDBRMMJJColl=cms.InputTag("BestCandSelectorMu:doubleJet")
# process.ANEDBR.EDBREEJColl=cms.InputTag("BestCandSelectorEle:singleJet")
# process.ANEDBR.EDBRMMJColl=cms.InputTag("BestCandSelectorMu:singleJet")

# process.ANEDBR.EDBREEJJColl=cms.InputTag("BestFullRangeSelectorEle:doubleJet")
# process.ANEDBR.EDBRMMJJColl=cms.InputTag("BestFullRangeSelectorMu:doubleJet")
# process.ANEDBR.EDBREEJColl=cms.InputTag("BestFullRangeSelectorEle:singleJet")
# process.ANEDBR.EDBRMMJColl=cms.InputTag("BestFullRangeSelectorMu:singleJet")

processFullSel='<PROCESS>'

if processFullSel == "fullsb" :
process.ANEDBR.EDBREEJJColl=cms.InputTag("BestSidebandSelectorEle:doubleJet")
process.ANEDBR.EDBRMMJJColl=cms.InputTag("BestSidebandSelectorMu:doubleJet")
process.ANEDBR.EDBREEJColl=cms.InputTag("BestSidebandSelectorEle:singleJet")
process.ANEDBR.EDBRMMJColl=cms.InputTag("BestSidebandSelectorMu:singleJet")
elif processFullSel == "fullsig" :
process.ANEDBR.EDBREEJJColl=cms.InputTag("BestCandSelectorEle:doubleJet")
process.ANEDBR.EDBRMMJJColl=cms.InputTag("BestCandSelectorMu:doubleJet")
process.ANEDBR.EDBREEJColl=cms.InputTag("BestCandSelectorEle:singleJet")
process.ANEDBR.EDBRMMJColl=cms.InputTag("BestCandSelectorMu:singleJet")
elif processFullSel == "fullrange" :
process.ANEDBR.EDBREEJJColl=cms.InputTag("BestFullRangeSelectorEle:doubleJet")
process.ANEDBR.EDBRMMJJColl=cms.InputTag("BestFullRangeSelectorMu:doubleJet")
process.ANEDBR.EDBREEJColl=cms.InputTag("BestFullRangeSelectorEle:singleJet")
process.ANEDBR.EDBRMMJColl=cms.InputTag("BestFullRangeSelectorMu:singleJet")
elif processFullSel == "ttbar" :
process.ANEDBR.EDBREEJJColl=cms.InputTag("BestTTBarSelectorEle:doubleJet")
process.ANEDBR.EDBRMMJJColl=cms.InputTag("BestTTBarSelectorMu:doubleJet")
process.ANEDBR.EDBREEJColl=cms.InputTag("BestTTBarSelectorEle:singleJet")
process.ANEDBR.EDBRMMJColl=cms.InputTag("BestTTBarSelectorMu:singleJet")
else :
print 'Processing preselected collections (default)'


##### set Ngen and xsect values ofr MC samples; xsect in pb !!!
if "TTBAR"=="<SAMPLE>":
Expand Down Expand Up @@ -697,15 +706,22 @@
process.ANEDBR.VTaggingScaleFactorLP=cms.double(1.1)

process.filterFinalSelPath = cms.EDFilter("HLTHighLevel",
TriggerResultsTag = cms.InputTag("TriggerResults","","CMG"),
HLTPaths = cms.vstring("cmgEDBRWWEle","cmgEDBRWWMu"),
#HLTPaths = cms.vstring("cmgEDBRZZEle","cmgEDBRZZMu"),
eventSetupPathsKey = cms.string(''),
andOr = cms.bool(True), # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
throw = cms.bool(True) # throw exception on unknown path names
)

if processFullSel :
TriggerResultsTag = cms.InputTag("TriggerResults","","CMG"),
HLTPaths = cms.vstring("FullPathDUMMYELE","FullPathDUMMYMU"),
#HLTPaths = cms.vstring("cmgEDBRZZEle","cmgEDBRZZMu"),
eventSetupPathsKey = cms.string(''),
andOr = cms.bool(True), # how to deal with multiple triggers: True (OR) accept if ANY is true, False (AND) accept if ALL are true
throw = cms.bool(True) # throw exception on unknown path names
)

if "WW"=="<ANALYSIS>" :
process.filterFinalSelPath.HLTPaths = cms.vstring("cmgEDBRWWEle","cmgEDBRWWMu")
else :
process.filterFinalSelPath.HLTPaths = cms.vstring("cmgEDBRZZEle","cmgEDBRZZMu")


#AB: should we add also "ttbar" ?
if processFullSel=="fullsb" or processFullSel=="fullsig" or processFullSel=="fullrange" :
process.p=cms.Path(process.filterFinalSelPath+process.ANEDBR)
else :
process.p=cms.Path(process.ANEDBR)
32 changes: 29 additions & 3 deletions EDBRCommon/test/run_AnalyzerEDBR.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -15,23 +15,49 @@ else
exit 1
fi

CHANNEL="ZZ"
if [ $# -gt 1 ]
then
CHANNEL=$2

if [[ ${CHANNEL} != "WW" && ${CHANNEL} != "ZZ" ]]
then
echo "Wrong analysis channel provided to run_AnalyzerEDBR.sh !!! You gave as input: ${CHANNEL}. Valid options are: WW ZZ"
exit 2
fi
fi

echo
echo "My working area is $MYWORKAREA"
echo

type="presel" #### could be: presel (default), fullsb, fullsig, fullrange, ttbar
if [ $# -gt 2 ]
then
type=$3

if [[ ${type} != "presel" && ${type} != "fullsb" && ${type} != "fullsig" && ${type} != "fullrange" && ${type} != "ttbar" ]]
then
echo "Wrong analysis channel provided to run_AnalyzerEDBR.sh !!! You gave as input: ${type}. Valid options are: presel fullsb fullsig fullrange ttbar"
exit 3
fi
fi

echo "Inputs to run_AnalyzerEDBR.sh: $TAG $CHANNEL $type"

type="fullsb" #### could be fullsb or presel or fullsig, denpends on analyzerEDBR_zz2l2j_tpl.py

CFGAREA=${MYWORKAREA}/ExoDiBosonResonances/EDBRCommon/test/tmp_cfg_${type}/
mkdir -p $CFGAREA


OUTDIR=/afs/cern.ch/user/b/bonato/work/PhysAnalysis/EXOVV_2012/analyzer_trees/productionv2d/${type}
OUTDIR=/afs/cern.ch/user/b/bonato/work/PhysAnalysis/EXOVV_2012/analyzer_trees/productionv2g/${type}
mkdir -p $OUTDIR/logs

# ok, now do the job:
sed -e "s/<SAMPLE>/${TAG}/g" < ${MYWORKAREA}/ExoDiBosonResonances/EDBRCommon/test/analyzerEDBR_zz2l2j_tpl.py > ${CFGAREA}/analyzerEDBR_zz2l2j_${type}_${TAG}_cfg.py
sed -e "s/<SAMPLE>/${TAG}/g" < ${MYWORKAREA}/ExoDiBosonResonances/EDBRCommon/test/analyzerEDBR_zz2l2j_tpl.py > ${CFGAREA}/analyzerEDBR_zz2l2j_${type}_${TAG}_TMP.py
sed -e "s/<ANALYSIS>/${CHANNEL}/g" < ${CFGAREA}/analyzerEDBR_zz2l2j_${type}_${TAG}_TMP.py > ${CFGAREA}/analyzerEDBR_zz2l2j_${type}_${TAG}_TMP2.py
sed -e "s/<PROCESS>/${type}/g" < ${CFGAREA}/analyzerEDBR_zz2l2j_${type}_${TAG}_TMP2.py > ${CFGAREA}/analyzerEDBR_zz2l2j_${type}_${TAG}_cfg.py
rm -f ${CFGAREA}/analyzerEDBR_zz2l2j_${type}_${TAG}_TMP*.py
cd $OUTDIR
cmsRun ${CFGAREA}/analyzerEDBR_zz2l2j_${type}_${TAG}_cfg.py &> "logs/log_AnalyzerPresel_${type}_${TAG}.out"

Expand Down
9 changes: 6 additions & 3 deletions EDBRCommon/test/submit_AnalyzerEDBR.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,10 @@

### BulkG_ZZ_lljj_c0p2_M800 BulkG_ZZ_lljj_c0p2_M1100 BulkG_ZZ_lljj_c0p2_M1400 BulkG_ZZ_lljj_c0p2_M1800

SampleName=( )
TYPE="fullsig" ### presel fullsb fullsig fullrange ttbar
CHANNEL="ZZ" ###"WW" "ZZ"

SampleName=( RSG_ZZ_lljj_c0p2_M1500 RSG_ZZ_lljj_c0p2_M1000 RSG_ZZ_lljj_c0p05_M1000 ) ###DoubleMu_Run2012A_22Jan2013 DoubleMuParked_Run2012B_22Jan2013 DoubleMuParked_Run2012C_22Jan2013 DoubleMuParked_Run2012D_22Jan2013 Photon_Run2012A_22Jan2013 DoublePhotonHighPt_Run2012B_22Jan2013 DoublePhotonHighPt_Run2012C_22Jan2013 DoublePhotonHighPt_Run2012D_22Jan2013 DYJetsPt70To100 DYJetsPt100 TTBARpowheg DYJetsPt70To100 DYJetsPt100 WW WZ ZZ BulkG_ZZ_lljj_c0p2_M600 BulkG_ZZ_lljj_c0p2_M700 BulkG_ZZ_lljj_c0p2_M800 BulkG_ZZ_lljj_c0p2_M900 BulkG_ZZ_lljj_c0p2_M1000 BulkG_ZZ_lljj_c0p2_M1100 BulkG_ZZ_lljj_c0p2_M1200 BulkG_ZZ_lljj_c0p2_M1300 BulkG_ZZ_lljj_c0p2_M1400 BulkG_ZZ_lljj_c0p2_M1500 BulkG_ZZ_lljj_c0p2_M1600 BulkG_ZZ_lljj_c0p2_M1700 BulkG_ZZ_lljj_c0p2_M1800 BulkG_ZZ_lljj_c0p2_M1900 BulkG_ZZ_lljj_c0p2_M2000 BulkG_ZZ_lljj_c0p2_M2100 BulkG_ZZ_lljj_c0p2_M2200 BulkG_ZZ_lljj_c0p2_M2300 BulkG_ZZ_lljj_c0p2_M2400 BulkG_ZZ_lljj_c0p2_M2500 )


#### DATA:
Expand All @@ -15,14 +18,14 @@ SampleName=( )
# BulkG_ZZ_lljj_c0p2_M600 BulkG_ZZ_lljj_c0p2_M700 BulkG_ZZ_lljj_c0p2_M800 BulkG_ZZ_lljj_c0p2_M900 BulkG_ZZ_lljj_c0p2_M1000 BulkG_ZZ_lljj_c0p2_M1100 BulkG_ZZ_lljj_c0p2_M1200 BulkG_ZZ_lljj_c0p2_M1300 BulkG_ZZ_lljj_c0p2_M1400 BulkG_ZZ_lljj_c0p2_M1500 BulkG_ZZ_lljj_c0p2_M1600 BulkG_ZZ_lljj_c0p2_M1700 BulkG_ZZ_lljj_c0p2_M1800 BulkG_ZZ_lljj_c0p2_M1900 BulkG_ZZ_lljj_c0p2_M2000 BulkG_ZZ_lljj_c0p2_M2100 BulkG_ZZ_lljj_c0p2_M2200 BulkG_ZZ_lljj_c0p2_M2300 BulkG_ZZ_lljj_c0p2_M2400 BulkG_ZZ_lljj_c0p2_M2500

#### alternative signal samples:
# BulkG_ZZ_lljj_c1p0_M600 BulkG_ZZ_lljj_c1p0_M1500 BulkG_ZZ_lljj_M700_G40 BulkG_ZZ_lljj_M700_G80 BulkG_ZZ_lljj_M700_G120 BulkG_ZZ_lljj_M2000_G40 BulkG_ZZ_lljj_M2000_G80 BulkG_ZZ_lljj_M2000_G120 BulkG_ZZ_lljj_M1000_G40 BulkG_ZZ_lljj_M1000_G80 BulkG_ZZ_lljj_M1000_G120 BulkG_ZZ_lljj_M1500_G40 BulkG_ZZ_lljj_M1500_G80 BulkG_ZZ_lljj_M1500_G120 BulkG_ZZ_lljj_c0p2_M1700-JHU-herwig BulkG_ZZ_lljj_c0p2_M1400-JHU-herwig BulkG_ZZ_lljj_c0p2_M800-JHU-herwig BulkG_ZZ_lljj_c0p2_M1100-JHU-herwig
# BulkG_ZZ_lljj_c1p0_M600 BulkG_ZZ_lljj_c1p0_M1500 BulkG_ZZ_lljj_M700_G40 BulkG_ZZ_lljj_M700_G80 BulkG_ZZ_lljj_M700_G120 BulkG_ZZ_lljj_M2000_G40 BulkG_ZZ_lljj_M2000_G80 BulkG_ZZ_lljj_M2000_G120 BulkG_ZZ_lljj_M1000_G40 BulkG_ZZ_lljj_M1000_G80 BulkG_ZZ_lljj_M1000_G120 BulkG_ZZ_lljj_M1500_G40 BulkG_ZZ_lljj_M1500_G80 BulkG_ZZ_lljj_M1500_G120 BulkG_ZZ_lljj_c0p2_M1700-JHU-herwig BulkG_ZZ_lljj_c0p2_M1400-JHU-herwig BulkG_ZZ_lljj_c0p2_M800-JHU-herwig BulkG_ZZ_lljj_c0p2_M1100-JHU-herwig RSG_ZZ_lljj_c0p2_M1500 RSG_ZZ_lljj_c0p2_M1000 RSG_ZZ_lljj_c0p05_M1000



for sample in "${SampleName[@]}"
do
echo "Submitting $sample"
bsub -q 8nh -J "treeEDBR_${sample}" run_AnalyzerEDBR.sh $sample
bsub -q 8nh -J "treeEDBR_${sample}" run_AnalyzerEDBR.sh $sample $CHANNEL $TYPE
echo
done

Expand Down

0 comments on commit 6165f8f

Please sign in to comment.