inputTrees for >=1lepton analysis

V12 production for paper

D4PD for Hard1Lep and ZMET

  • Sample name corresponds to lists in submit/FileList_p1328.
Status Sample Resource Name Version DoSyst Comment (naming convention etc)
smile Egamma Local Yuichi, Naoko 12-01-03 False  
smile Muons Local Yuichi, Naoko 12-01-03 False  
smile Muons (QCD, mt2l) Local Naoko 12-01-10 False BUGFIX. Need to use ntupleKits-00-12-01-10 for mt2lMuo QCD tree production
smile PowhegPythiaTTbar Local Yuichi, Naoko 12-01-03 True pdfUncertainty should be computed using the sample produced with ntupleKits-00-12-01-09.
smile SingleTop Local Yuichi, Naoko 12-01-03 True  
smile ttbarV Local Naoko 12-01-03 True *Use MakeScriptsFromList.py in >=12-01-05
smile SherpaWMassiveBC Grid Manuel 12-01-03 True "group.phys-susy.schreyer.*_*.Strong1Lep.*.v12_1_trial01/"
smile SherpaZMassiveBC Grid Manuel 12-01-03 True "group.phys-susy.schreyer.*_*.Strong1Lep.*.v12_1_trial01/"
smile SherpaWMassiveBC Local Naoko 12-01-03 True  
smile SherpaZMassiveBC Local Naoko 12-01-03 True  
smile SherpaDibosons Local Yuichi 12-01-03 True  
smile Dibosons Local Naoko 12-01-03 True PowhegPythiaDibosons
smile SherpaDibosonsMassiveBC Grid Manuel 12-01-03 True

!FileList is available in trunk

"group.phys-susy.schreyer.*_*.Strong1Lep.*.v12_1_trial01/"

smile SherpaDY Grid Manuel 12-01-03 True "group.phys-susy.schreyer.*_*.Strong1Lep.*.v12_1_trial01/"
smile PowhegPythiaTTbarSyst Local Yuichi, Naoko 12-01-04 False  
smile PowhegJimmyTTbar Local Yuichi, Naoko 12-01-04 False  
  AlpgenWgammaLeptonFilter Grid Pedro 12-01-08 True  
smile AlpgenWgamma Local Naoko 12-01-08 True (no lepton filter)
smile McAtNloTTbar Local Yuichi, Naoko 12-01-04 False  
smile AlpgenTTbar Local Yuichi, Naoko 12-01-03 False  
smile SingleTopSyst Local Yuichi, Naoko 12-01-04 False  
smile ttbarVSyst Local Yuichi, Naoko 12-01-03 False  
smile HiggsSugra Local Yuichi 12-01-03 True  
smile SMGGoneStepCC Local Yuichi 12-01-03 True  
smile SMSSoneStepCC Local Yuichi 12-01-03 True  
smile SMGG2WWZZ Local Yuichi, Naoko 12-01-03 True  
smile SMSS2WWZZ Local Yuichi, Naoko 12-01-03 True  
smile SMGG2CNsl Local Yuichi 12-01-03 True  
smile SMSS2CNsl Local Yuichi 12-01-03 True  
smile SbottomTChar Local Yuichi 12-01-03 True  
smile SbottomTChar Grid Pedro 12-01-03 True  
smile Gtt Grid Pedro 12-01-03 True  
smile GttCharm Grid Perdo 12-01-03 True  
smile NUHMG Grid Pedro 12-01-03 True  
smile nGM Grid Pedro 12-01-03 12-01-03 True  
  bRPV Grid Pedro 12-01-06 True extend SignalFamily to EWkino
  GGM Grid Pedro 12-01-06 True extend SignalFamily to EWkino
smile bRPV Local Naoko 12-01-06 True extend SignalFamily to EWkino
smile GGM Local Elena 12-01-06 True extend SignalFamily to EWkino

tree for Hard1Lep

  • Sample name is one of arguments given to treeMaker.
  • FullSyst =True : 59 systematic treenames. See here
  • FullSyst = False : NoSys only
  • Warning, important please use GridManager.py which will be provided by Pedro and Adrian privately later [10 Dec].
Status Sample Resource Name Version FullSyst Comment (naming convention etc)
smile Data Local Yuichi 12-01-04 12-01-05 False  
smile QCD Local Yuichi 12-01-04 12-01-05 False  
smile PowhegPythiaTTbar Local Naoko 12-01-04 12-01-05 True pdfUncertainty should be computed using the sample produced with ntupleKits-00-12-01-09.
smile SingleTop Local Yuichi 12-01-04 12-01-05 True  
smile ttbarV Local Naoko 12-01-04 12-01-05 True  
smile SherpaWMassiveBC Grid Adrian 12-01-04 12-01-05 True  
smile SherpaZMassiveBC Grid Adrian 12-01-04 12-01-05 True  
smile SherpaWMassiveBC Local Naoko 12-01-04 12-01-05 True  
  SherpaZMassiveBC Local Naoko 12-01-04 12-01-05 True Missing JVF Trees?
smile SherpaDibosons Local Yuichi 12-01-04 12-01-05 True  
  SherpaDibosonsMassiveBC Grid Adrian 12-01-04 12-01-05 True FileList is available in trunk. Missing JVF Trees?
smile SherpaDY Grid Adrian 12-01-04 12-01-05 True  
  WgammaLeptonFilter Grid Adrian 12-01-08 True  
smile Wgamma Local Naoko 12-01-08 True  
  Test Local Yuichi 12-01-04 12-01-05 False PowhegPythiaTTbarSyst(105861)
smile PowhegJimmyTTbar Local Yuichi 12-01-04 12-01-05 False  
smile McAtNloTTbar Local Yuichi 12-01-04 12-01-05 False  
smile AlpgenTTbar Local Yuichi 12-01-04 12-01-05 False  
smile SingleTopSyst Local Yuichi 12-01-04 12-01-05 False  
smile ttbarVSyst Local Yuichi 12-01-04 12-01-05 False  
smile HiggsSU Local Yuichi 12-01-04 12-01-05 True  
smile SMGG1step Local Yuichi 12-01-04 12-01-05 True  
smile SMSS1step Local Yuichi 12-01-04 12-01-05 True  
smile SMGGWWZZ Local Yuichi, Naoko 12-01-04 12-01-05 True  
smile SMSS2WWZZ Local Yuichi, Naoko 12-01-04 12-01-05 True  
smile SMGG2CNsl Local Yuichi 12-01-04 12-01-05 True  
smile SMSS2CNsl Local Yuichi 12-01-04 12-01-05 True  
smile SMSbottomTChar Local Yuichi 12-01-04 12-01-05 True  
smile SMGtt Grid Valentina 12-01-04 12-01-05 True  
smile SMGttCharm Grid Valnetina 12-01-04 12-01-05 True  
  bRPV Grid Marija 12-01-04 12-01-05 True use 12-01-06 if you care about EWkino processes
smile NUHMG Grid Valentina 12-01-04 12-01-05 True  
smile nGM Grid Marija 12-01-04 12-01-05 True  

tree for ZMET

  • Sample name is one of arguments given to treeMaker.
  • FullSyst =True : 59 systematic treenames. See here
  • FullSyst = False: NoSys only
  • Warning, important please use GridManager.py which will be provided by Pedro and Adrian privately later [10 Dec].
Status Sample Resource Name Version FullSyst Comment (naming convention etc)
smile Data Local Naoko 12-01-04 12-01-07 False  
smile QCD Local Naoko 12-01-04 12-01-05 False  
smile PowhegPythiaTTbar Local Naoko 12-01-04 12-01-05 True pdfUncertainty should be computed using the sample produced with ntupleKits-00-12-01-09.
smile SingleTop Local Naoko 12-01-04 12-01-05 True  
smile ttbarV Local Naoko 12-01-04 12-01-05 True  
  SherpaWMassiveBC Grid Adrian 12-01-04 12-01-05 True  
  SherpaZMassiveBC Grid Adrian 12-01-04 12-01-05 True  
  SherpaWMassiveBC Local Naoko 12-01-04 12-01-05 True  
  SherpaZMassiveBC Local Naoko 12-01-04 12-01-05 True  
smile PowhegDibosons Local Naoko 12-01-04 12-01-05 True PowhegPythiaDibosons
  SherpaDY Grid Adrian 12-01-04 12-01-05 True  
  WgammaLeptonFilter Grid Manuel 12-01-08 True  
  Wgamma Local Naoko 12-01-08 True  
smile Test Local Naoko 12-01-04 12-01-05 False PowhegPythiaTTbarSyst(105861)
smile PowhegJimmyTTbar Local Naoko 12-01-04 12-01-05 False  
smile McAtNloTTbar Local Naoko 12-01-04 12-01-05 False  
smile AlpgenTTbar Local Naoko 12-01-04 12-01-05 False  
smile SingleTopSyst Local Naoko 12-01-04 12-01-05 False --> will produce SinlgeTopDiLeptSyst
smile ttbarVSyst Local Naoko 12-01-04 12-01-05 False  
smile SMGGWWZZ Local Naoko 12-01-04 12-01-05 True  
smile SMSS2WWZZ Local Naoko 12-01-04 12-01-05 True  
smile bRPV Grid Naoko 12-01-04 12-01-06 True Extend SignalFamily to EWKino processes /eos/atlas/user/n/nkanaya/trees/Moriond/p1328/mt2l/v12_1_6/
smile GGM Grid Elena 12-01-04 12-01-06 True Extend SignalFamily to EWKino processes. /eosatlas//eos/atlas/user/e/eromeroa/trees/v12_1_6/

D4PD for Soft1Lep, Soft2Lep

  • Sample name corresponds to lists in submit/FileList_p1328.
Status Sample Resource Name Version DoSyst Comment (Grid location etc)
smile JetTauEtmiss Local Takuya 12-01-03 False  
smile PowhegPythiaTTbar Local Takuya 12-01-03 True pdfUncertainty should be computed using the sample produced with ntupleKits-00-12-01-09.
smile SingleTop Local Takuya 12-01-03 True  
smile ttbarV Local Takuya 12-01-03 True *Use MakeScriptsFromList.py in >=12-01-05
smile SherpaWMassiveBC Grid Pedro 12-01-03 True  
smile SherpaZMassiveBC Grid Pedro 12-01-03 True  
smile SherpaDibosons Local Takuya 12-01-03 True  
smile Dibosons Local Takuya 12-01-03 True PowhegPythiaDibosons
smile SherpaDibosonMassiveBC Grid Pedro 12-01-03 True two samples (183591,183592) are still missing. FileList is available in trunk
smile SherpaDY Grid Pedro 12-01-03 True  
smile AlpgenDY Local Takuya 12-01-03 True  
smile AlpgenWgamma Local Takuya 12-01-08 True  
smile PowhegPythiaTTbarSyst Local Takuya 12-01-04 False  
smile PowhegJimmyTTbar Local Takuya 12-01-04 False  
smile McAtNloTTbar Local Takuya 12-01-04 False  
smile AlpgenTTbar Local Takuya 12-01-03 False  
smile SingleTopSyst Local Takuya 12-01-04 False  
smile ttbarVSyst Local Takuya 12-01-03 False  
smile SMGGoneStepCC Local Takuya 12-01-03 True  
smile SMSSoneStepCC Local Takuya 12-01-03 True  
smile mUED Local Takuya 12-01-03 True JETMET filter
smile mUED2Lfilter Local Takuya 12-01-03 True  
smile StopBCharDeg Local Takuya 12-01-03 True  
smile StopBCharDegMET250 Local Takuya 12-01-03 True  
smile StopBCharDegFixed Local Takuya 12-01-03 True  
smile SMGG2CNsl Local Takuya 12-01-03 True Low Priority
smile SMSS2CNsl Local Takuya 12-01-03 True Low Priority

tree for Soft1Lep, Soft2Lep

  • Sample name is one of arguments given to treeMaker.
  • FullSyst =True : means 59 systematic treenames. See here
  • FullSyst = False : NoSys only.
  • Warning, important please use GridManager.py which will be provided by Pedro and Adrian privately later [10 Dec].
Status Sample Resource Name Version FullSyst Comment (naming convention etc)
smile Data Local Takuya 12-01-04 12-01-05 (12-01-07 for Soft2Lep) False  
smile QCD Local Takuya 12-01-04 12-01-05 False  
smile PowhegPythiaTTbar Local Takuya 12-01-04 12-01-05 True pdfUncertainty should be computed using the sample produced with ntupleKits-00-12-01-09.
smile SingleTop Local Takuya 12-01-0412-01-05 True  
smile ttbarV Local Takuya 12-01-04 12-01-05 True  
smile SherpaWMassiveBC Grid Da 12-01-04 12-01-05 True  
smile SherpaZMassiveBC Grid Da 12-01-04 12-01-05 True group.phys-susy.urrejola.*.Strong1Lep.SoftLep.v12_1_3_trial2
smile PowhegDibosons Local Takuya 12-01-04 12-01-05 True  
smile SherpaDibosons Local Takuya 12-01-04 12-01-05 True  
smile SherpaDibosonsMassiveBC Grid Takuya 12-01-04 12-01-05 True

two samples (183591,183592) are still missing. FileList is available in trunk.

smile SherpaDY Grid Takuya 12-01-04 12-01-05 True  
smile Wgamma Local Takuya 12-01-08 True  
smile Test Local Takuya 12-01-04 12-01-05 False PowhegPythiaTTbarSyst(105861)
smile PowhegJimmyTTbar Local Takuya 12-01-04 12-01-05 False  
smile McAtNloTTbar Local Takuya 12-01-04 12-01-05 False  
smile AlpgenTTbar Local Takuya 12-01-04 12-01-05 False  
smile SingleTopSyst Local Takuya 12-01-04 12-01-05 False  
smile ttbarVSyst Local Takuya 12-01-04 12-01-05 False  
smile SMGG1step Local Takuya 12-01-04 12-01-05 True  
smile SMSS1step Local Takuya 12-01-04 12-01-05 True  
smile mUEDMETfilter Local Takuya 12-01-04 12-01-05 True JETMET filter
smile mUED2Lfilter Local Takuya 12-01-04 12-01-05 True  
smile StopBCharDeg Local Takuya 12-01-04 12-01-05 True  
smile SMGG2CNsl Local Takuya 12-01-04 12-01-05 True Lower priority
smile SMSS2CNsl Local Takuya 12-01-04 12-01-05 True Lower priority

V12 theory systematic uncertainty for paper

Files

TFile *file = TFile::Open("root://eosatlas//eos/atlas/user/y/ysasaki/trees/Moriond/v4_4_6/<FileName>.root");

Suffix comment
HardEle(Muo) >=1lepton. Leading lepton is Electron(Muon)
mt2lEle(Muo) More than 2leptons. The leading lepton is Electron(Muon)
Soft1Lep ==1 soft lepton. The branch named "lep1Flavor" indicates the flavor.
Soft2Lep >=2 soft leptons. The branch named "lep1Flavor" indicates the flavor.

file contents comment
datatree_*.root Data and QCD (run=200842-212272) Unprescaled luminosity (after GRL) : 14.3ifb
sigtree_*_SM_{SS,GG}1step.root Simplfiied Model GG(SS) onestep Decay  
sigtree_*_SM_{SS,GG}2CNsl.root Simplfiied Model GG(SS) twostep Decay ()  
sigtree_*_HiggsSU.root mSugra ( maximum mixing i.e. Higgs-aware samples )  
sigtree_*_pMSSM.root pMSSM grid (LO xsec)  
sigtree_*_GGM.root GGM grid  
bkgtree_*.root Backgrounds  

bkgtree contents
SherpaW3jets 3jet filtered W sample
SherpaWMassiveB Massive B samples ( Pt0_BFilter + Pt0_BVeto + Pt140_BFilter + Pt140_BVeto ). C is massless.
AlpgenW Alpgen W sample (HF is included)
SherpaZMassiveB Massive B+C samples ( LightJets + HeavyJets )
AlpgenZ AlpgenZ sample (HF is included)
SherpaTTbar Sherpa ttbar sample
AlpgenTTbar Alpgen ttbar sample
PowhegPythiaTTbar PowhegPythia ttbar sample
SingleTop SingleTop samples
Dibosons Diboson samples
ttbarV ttbar+V samples

You should choose one sample from 3jets/MassiveB for your analysis. My recommendation is MassiveB sample as it has the best agreement with data.

SherpaZ ZnunuLightJets+ZnunuHeavyJets+Zll
SherpaZMassiveB (Recommend) ZnunuLightJets+ZnunuHeavyJets+ZllLightJets+ZllHeavyJets
SherpaWMassiveB (Recommend) WlnuMassiveBPt0_BFilter(pt<140)+WlnuMassiveBPt0_VetoB(pt<140)+WlnuMassiveBPt140_BFilter(pt>140)+WlnuMassiveB_VetoB(pt>140)
SherpaW3jets Wlnu3jets

* SherpaWincl is deleted to reduce file size.

Skimming

AnalysisType Skim
HardEle(Muo) Trigger, jet1Pt>80, met>100
mt2lEle(Muo) Trigger
Soft1Lep Trigger
Soft2Lep Trigger

Object Definition

  • No big changes from the previous definitions. Please see the diff for detailed changes.
  • Trigger : See AnalysisType.
  • Jet : |eta|<2.5.

Caveats

  • Please be careful for the trigger bias and evgen-level filters!
  • Event weight for Btag SF has changed a lot. Please see below for further information.

Contents

AnalysisType Value Trigger
HardEle 1 EF_e24vh_medium1_EFxe35_tclcw OR EF_e60_medium1
HardMuo 2 EF_mu24_j65_a4tchad_EFxe40_tclcw
mt2lEle 3 EF_e24vh_medium1_e7_medium1 OR EF_2e12Tvh_loose1 OR EF_e12Tvh_medium1_mu8
mt2lMuo 4 EF_mu18_tight_mu8_EFFS OR EF_2mu13 OR EF_mu18_tight_e7_medium1
Soft1/2Lep 5 EF_xe80_tclcw

DatasetNumber corresponds to mc_channel_number in D3PD.

  • lepPt, lepSignal, etc...
 vector<Float_t> lepPt;
 vector<Int_t> lepCharge;
 vector<Int_t> lepFlavor; // -1 : Muon, +1 : Electron
 vector<Bool_t> lepSignal; // True for signal lepton. False for baseline lepton.
 vector<Bool_t> lepSignalExpLoose; // True for ExpLoose lepton, False for baseline lepton.
 vector<Bool_t> lepSignalExpMedium; // True for ExpMedium lepton, False for baseline lepton.

Example :

- You can draw the flavor of 2nd lepton by tree->Draw("lepFlavor[1]")
- You can access the number of leptons by tree->Draw("@lepPt.size()")
- More examples can be found in http://root.cern.ch/root/html/TTree.html#TTree:Draw@2

*isNoCrackElectron True if there's no crack electron. The electron used here is 'baseline electron'.

Number of baseline/signal leptons.

  • nTruthJetXX
Number of Truth (anti-kt) jets with pt>XXGeV.

  • mll, Zpt
Dilepton mass and Zpt calculated from leading 2 signal leptons.

  • Wpt
Wpt calculated from leading lepton and reconstructed MET.

  • Vpt
Truth pt of 'V'ector boson.

  • mt2
mt2 calculated from two leptons and missingEt using ChengHanBisect _Mt2_332_Calculator (http://www.hep.phy.cam.ac.uk/~lester/mt2/index.html )

  • meffXJetYY
Effective mass calculated with Leading XJets. Jets with Pt<YYGeV is not considered.

  • nJetXX
Number of signal jets with Pt>XXGeV.

  • nBXJetsYY
Number of BJets in the leading XJets. Jets with Pt<YYGeV is not considered.

Rough classification of ttbar sample.
DecayIndexTop Decay
0 Full hadronic
1 Dileptonic
2 Dileptonic (one of the decay is tau(had))
3 Semi-leptonic

Detailed classfication of ttbar sample.
DecayModeTop Decay
0 tt -> bb + lep + lep
1 tt -> bb + lep + tau(lep)
2 tt -> bb + tau(lep) + tau(lep)
3 tt -> bb + lep + had
4 tt -> bb + lep + tau(had)
5 tt -> bb + tau(lep) + had
6 tt -> bb + tau(lep) + tau(had)
7 tt -> bb + had + had
8 tt -> bb + had + tau(had)
9 tt -> bb + tau(had) + tau(had)

  • jetXdPhi
Delta Phi between X-th jet and MET.

  • bTagWeight, bTagWeightBUp(Down), bTagWeightCUp(Down), bTagWeightMUp(Down)
Defined as vector. For example, bTagWeight[1] corresponds to the bTag SF considering leading 3jets (JetPt<35GeV is not taken into account i.e. if Jet pt = (100,55,25), then the weight for the leading two jets are multiplied and stored.).

The conditions for the other indices are summarized here

Systematic variations for the bTag SF are stored in the branches with BUp(Down), CUp(Down), MUp(Down) suffixes (MUp(Down) is 'Misstag' systmatics).

  • qfacUp(Down)WeightTTbar, ktfacUp(Down)WeightTTbar, qfacUp(Down)WeightW, ktfacUp(Down)WeightW, iqopt2(3)WeightW
nParton systematic variation for k-factor on qfac, ktfac and iqopt. These values are taken from ScaleVariatioReweighter

Datasets

From this version, DecayIndexTop and DecayModeTop are filled correctly.
105861 PowhegPythia_AUET2BCT10_ttbar_LeptonFilter 238.06 * 1.000000 * 0.543

117800 Sherpa_CT10_TtbarLeptLept 9.161600 * 1.243 * 1.000000
117801 Sherpa_CT10_TtbarLeptTaulept 9.172400 * 1.243 * 0.352
117802 Sherpa_CT10_TtbarTauleptTaulept 2.259500 * 1.243 * 0.123904
117803 Sherpa_CT10_TtbarLeptHad 55.53000 * 1.243 * 1.000000
117804 Sherpa_CT10_TtbarLeptTauhad 9.258200 * 1.243 * 0.648
117805 Sherpa_CT10_TtbarTauleptHad 27.766000 * 1.243 * 0.352
117806 Sherpa_CT10_TtbarTauleptTauhad 2.310700 * 1.243 * 0.456192
117808 Sherpa_CT10_TtbarHadTauhad 27.788000 * 1.243 * 0.648
117809 Sherpa_CT10_TtbarTauhadTauhad 2.303100 * 1.243 * 0.419904

147770 Sherpa_CT10_Zee 1207.900000
147771 Sherpa_CT10_Zmumu 1207.800000
147772 Sherpa_CT10_Ztautau 1207.100000

146820 Sherpa_CT10_ZeeLightJets 1051.200000
146821 Sherpa_CT10_ZmumuLightJets 1052.300000
146822 Sherpa_CT10_ZtautauLightJets 1051.300000
128975 Sherpa_CT10_ZeeHeavyJets 58.780000
128976 Sherpa_CT10_ZmumuHeavyJets 58.786000
128977 Sherpa_CT10_ZtautauHeavyJets 58.794000

144992 Sherpa_CT10_Wenu_3jets 461.690000 * 1.137000
144993 Sherpa_CT10_Wmunu_3jets 461.940000 * 1.137000
144994 Sherpa_CT10_Wtaunu_3jets 461.090000 * 1.137000

147774 Sherpa_CT10_Wenu 11878.000000 * 1.025900 VpT<200GeV
147775 Sherpa_CT10_Wmunu 11879.000000 * 1.025800 VpT<200GeV
147776 Sherpa_CT10_Wtaunu 11872.000000 * 1.026400 VpT<200GeV
157534 Sherpa_CT10_WenuPt200 8.067200 * 1.200000 VpT>200GeV
157535 Sherpa_CT10_WmunuPt200 8.068800 * 1.200000 VpT>200GeV
157536 Sherpa_CT10_WtaunuPt200 8.060800 * 1.200000 VpT>200GeV

167150 Sherpa_CT10_WenuMassiveBPt0_BFilter 11030.00 * 1.104800 * 0.016646
167151 Sherpa_CT10_WmunuMassiveBPt0_BFilter 11030.00 * 1.104800 * 0.016635
167152 Sherpa_CT10_WtaunuMassiveBPt0_BFilter 11028.00 * 1.105000 * 0.016649
167153 Sherpa_CT10_WenuMassiveBPt0_VetoB 11028.00 * 1.105000 * 0.983310
167154 Sherpa_CT10_WmunuMassiveBPt0_VetoB 11030.00 * 1.104800 * 0.983280
167155 Sherpa_CT10_WtaunuMassiveBPt0_VetoB 11033.00 * 1.104500 * 0.983250
167156 Sherpa_CT10_WenuMassiveBPt140_BFilter 33.19900 * 1.105000 * 0.065963
167157 Sherpa_CT10_WmunuMassiveBPt140_BFilter 33.19400 * 1.105000 * 0.065895
167158 Sherpa_CT10_WtaunuMassiveBPt140_BFilter 33.19400 * 1.105000 * 0.065901
167159 Sherpa_CT10_WenuMassiveBPt140_VetoB 33.15000 * 1.105000 * 0.934290
167160 Sherpa_CT10_WmunuMassiveBPt140_VetoB 33.16700 * 1.105000 * 0.934070
167161 Sherpa_CT10_WtaunuMassiveBPt140_VetoB 33.15000 * 1.105000 * 0.934000

  • Dibosons
105921 McAtNloJimmy_AUET2CT10_WpWmenuenu 0.621550
105922 McAtNloJimmy_AUET2CT10_WpWmenumunu 0.619510
105923 McAtNloJimmy_AUET2CT10_WpWmenutaunu 0.621790
105924 McAtNloJimmy_AUET2CT10_WpWmmunuenu 0.620800
105925 McAtNloJimmy_AUET2CT10_WpWmmunumunu 0.620740
105926 McAtNloJimmy_AUET2CT10_WpWmmunutaunu 0.620400
105927 McAtNloJimmy_AUET2CT10_WpWmtaunuenu 0.620690
105928 McAtNloJimmy_AUET2CT10_WpWmtaunumunu 0.620760
105929 McAtNloJimmy_AUET2CT10_WpWmtaunutaunu 0.619370
105986 Herwig_AUET2CTEQ6L1_ZZ 4.691400 * 1.550000 * 0.211650
105987 Herwig_AUET2CTEQ6L1_WZ 12.00900 * 1.900000 * 0.305460

117360 AcerMCPythia_AUET2BCTEQ6L1_singletop_tchan_e 8.604 * 1.10
117361 AcerMCPythia_AUET2BCTEQ6L1_singletop_tchan_mu 8.604 * 1.10
117362 AcerMCPythia_AUET2BCTEQ6L1_singletop_tchan_tau 8.604 * 1.10
108343 McAtNloJimmy_CT10NLOME_AUET2CTEQ6L1MPI_SingleTopSChanWenu 0.564440 * 1.074
108344 McAtNloJimmy_CT10NLOME_AUET2CTEQ6L1MPI_SingleTopSChanWmunu 0.564260 * 1.074
108345 McAtNloJimmy_CT10NLOME_AUET2CTEQ6L1MPI_SingleTopSChanWtaunu 0.564040 * 1.074
108346 McAtNloJimmy_CT10NLOME_AUET2CTEQ6L1MPI_SingleTopWtChanIncl 20.658000 * 1.083

  • ttbar+V
119353 MadGraphPythia_AUET2BCTEQ6L1_ttbarW 0.104100 * 1.175
119354 MadGraphPythia_AUET2BCTEQ6L1_ttbarWj 0.093317 * 1.175
119355 MadGraphPythia_AUET2BCTEQ6L1_ttbarZ 0.067690 * 1.34
119356 MadGraphPythia_AUET2BCTEQ6L1_ttbarZj 0.087339 * 1.34

How to run

CAUTION! The following explanation is valid after v7_1_4!

Setup MyAna ( D3PD -> D4PD )

mkdir YOURFOLDER
svn co svn+ssh://svn.cern.ch/reps/atlas-ysasaki/ysasaki/ntupleKits/tags/ntupleKits-00-XX-XX-XX YOURFOLDER
cd YOURFOLDER/MyAna
./SetupSUSYTools.sh

source setup.sh
gmake clean
gmake

Submit Jobs to Grid

Make sure your dq2 environment is correctly setup.

Open submit/GridManager.py and edit Submit() function. What you need to change are: YOURNAME, IteratoinString(this is something like a version number. It would be better to be finished with "_trialX" suffix for the consistency.).

Choose the Submit function that best fits your needs (gm.Submit(...) for data processing and gm.SubmitAllSyst(...) for MC samples) and edit it to resemble your case (below is an example where the GridManager will process the SherpaWMassiveB for the HardEle channel on the grid).

def Submit():
    gm = GridManager(UserName='YOURNAME')
    gm.SetExcludeSite(['ANALY_INFN-NAPOLI','ANALY_INFN-QMUL','ANALY_IN2P3-CC-T2_RD','ANALY_INFN-T1','ANALY_QMUL'])
    gm.RecreateTarFile(True)

    ItrStr = 'v10_5_trial1'
    UseSystGroup=True
    gm.SubmitAllSyst(LeptonMode='HardEle',FileList=gm.ConverFileList('FileList_p1328/SherpaWMassiveB'),IterationString=ItrStr,IsOfficial=True)
    #gm.Submit('HardEle',gm.ConverFileList('FileList_p1328/EgammaGrid'),ItrStr,IsOfficial=True,SystType="NoSys",UseSystGroup=False)
    return
    ... (but what comes after return doesn't matter)
    

And modify the script to call the Submit() function.

if __name__=='__main__':
    Submit()

Then, execute the script. If you're asked to "reuse the tar file", then answer "yes" unless you have modified MyAna.

python GridManager.py

Auto Retry

Grid jobs fail quite often. The reason of it is said to be a lack of memory space in the site, which is difficult to avoid from user side. Auto retry is implemented in GridManager.py. First, open it and modify Loop() function. As shown below, you can set a list of regular expressions to be retried.

def Loop():
    gm = GridManager(UserName='YOURNAME')
    while True:
        gm.MakeDSList(outDSFilterList=[r'Strong1Lep.*v7_1_trial1/',r'Strong1Lep.*v7_1_trial1.v01_trial2/'])
        gm.PrintJobStatus()
        gm.Resubmit(noSubmit=False)

and modify the script to execute Loop() function.

if __name__=='__main__':
    Loop()

Then, execute the script.

After the first iteration, you get a list of running/failed/finished Datasets in ./gridInfo.

treeMaker Setup

cd treeMaker
source setup.sh -q -c
make clean
make

Submit treeMaker Jobs

First, let's edit submit/treeSubmitter.py. You need to change : UserName, IterationString, and templList. templList is the regular expressions for the input D4PDs.

def SubmitMC():
    sm = Submit(UserName='ysasaki')
    SystList = ['NoSys','JESCloseBypos','JESCloseByneg','JESFlavorCompUncertneg','JESFlavorCompUncertpos','JESFlavorResponseUncertneg','JESFlavorResponseUncertpos','BJespos','BJesneg','JESpos','JESLowpos','JESMediumpos','JESHighpos','JESneg','JESLowneg','JESMediumneg','JESHighneg','JER','EGZEEpos','EGZEEneg','EGMATpos','EGMATneg','EGPSpos','EGPSneg','EGLOWpos','EGLOWneg','EGRESpos','EGRESneg','MSCALEpos','MSCALEneg','MMSpos','MMSneg','MIDpos','MIDneg','RESOST','SCALESTpos','SCALESTneg','EffectiveNP1pos','EffectiveNP1neg','EffectiveNP2pos','EffectiveNP2neg','EffectiveNP3pos','EffectiveNP3neg','EffectiveNP4pos','EffectiveNP4neg','EffectiveNP5pos','EffectiveNP5neg','EffectiveNP6pos','EffectiveNP6neg','EffectiveNP7pos','EffectiveNP7neg','EffectiveNP8pos','EffectiveNP8neg','EffectiveNP9pos','EffectiveNP9neg','EffectiveNP10pos','EffectiveNP10neg','EffectiveNP11pos','EffectiveNP11neg','EffectiveNP12pos','EffectiveNP12neg','EffectiveNP13pos','EffectiveNP13neg']
    SamList = ['PowhegPythiaTTbar','PowhegJimmyTTbar','AlpgenJimmyZ','AlpgenJimmyW','AlpgenPythiaZ','AlpgenPythiaW','SherpaZMassiveBC','SherpaWMassiveB','SherpaDibosons','Dibosons','SingleTop','ttbarV']
    from SampleList import sampleList
    for chan in ['HardEle','HardMuo']:
        for sam in SamList:
            sm.MakeFileList(sample=sam,channel=chan,templList=[r'group.phys-susy*.*%(runnum)d.*.Strong1Lep.%(channel)s.v7_1_trial1/'],runList=sampleList[sam])
    for sam in SamList:
        sm.Submit(IterationString='v01_trial1',IsOfficial=True,filt='*%s*'%sam,systlist=SystList)

After modifying the script, you can execute it to submit the jobs.

Resubmit jobs

You can resubmit the jobs by re-executing the script.

Downloading

Let's download it to downrootfiles.

mkdir downrootfiles
cd downrootfiles
dq2-get -D -f "*.root*" "group.phys-susy.*_8TeV.*.v7_1_trial1.v01_trial1/"

Merging the files

You need to merge the root files. There are several steps. (1) Expanding the tar files, (2) Merging files (recalculating genWeight for the files.) (3) Sum up the divided runs into a sample (4) Sum up the samples into one "bkgtree" files.

cd submit
make clean; make;
python GridMerger.py

The script automatically expand the files in '../downrootfiles' and merge the files. Several hours needed when processing all the files.

Please check the file named 'complist' after executing the script. This file is a summary of the number of events in the merged files. If you find "0k events" sample in the file, it means that the sample is failed somewhere in the production.

Then, you have the sample trees in 'submit/final'. Let's merge them into one "bkgtree" file. Open 'treeMaker/FileMerger.py' and modify Add_Bkg() function. fileName = "results/..." -> fileName = "../submit/final/...".

    for c in Channel:
        fm = FileMerger()
        fm.SetOutputFile("merged/bkgtree_%s.root"%(c))
        for sam in Sample:
            for sys in Systematics:
                fileName = "../submit/final/%s_%s_%s.root"%(c,sam,sys)
                fm.AddInputFile(fileName)
        fm.Merge()

And execute the script.

Production Log


InputTreeLog (Written In Japanese)

-- YuichiSasaki - 11-Oct-2012

Topic attachments
I Attachment History Action Size Date Who Comment
Texttxt systematic.txt r1 manage 1.0 K 2013-11-27 - 14:25 NaokoKanaya systematic name given to treeMaker
Edit | Attach | Watch | Print version | History: r57 < r56 < r55 < r54 < r53 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r57 - 2014-02-17 - PedroUrrejola
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback