Python HistsGroomed.book1Fの例

コード例 #1

ファイル: smpjs_tuple_analyzer_fwlite.py プロジェクト: pvmulder/usercode

# Mjj and <Mjet> versus GENERATED/RECO Mjj and <Mjet>. Only for MC. 
hists.book3F('histAK7MjjResponseVsPtAvg', 'AK7 m_{jj} Response;m_{jj} (GeV);p_{T}^{AVG} (GeV)',
              nx=70, x1=0., x2=7000., ny=80, y1=0.8, y2=1.2, nz=10, z1=0., z2=1000.)
hists.book3F('histAK7MjetResponseVsPtAvg', 'AK7 <m_{jet}> Versus p_{T}^{AVG} ;<m_{jet}> (GeV);p_{T}^{AVG} (GeV)',
              nx=30, x1=0., x2=300., ny=20, y1=0.0, y2=2.0, nz=10, z1=0., z2=1000.)

############## Responses ##############


# Pt and eta response. Only for MC
hists.book2F('histAK7PtResponse',  ';p_{T}^{GEN} (GeV);p_{T}^{RECO} / p_{T}^{GEN}', nx=50, x1=0., x2=500., ny=20, y1=0.5, y2=1.5)
hists.book2F('histAK7EtaResponse', ';#eta^{GEN} (GeV);p_{T}^{RECO} / p_{T}^{GEN}', nx=50, x1=-5.0, x2=5.0, ny=20, y1=0.5, y2=1.5)

hists.book1F('recoEffNum', 'reco efficiency numerator;p_{T}^{GEN} (GeV)', nx=30, x1=0., x2=300.)
hists.book1F('recoEffDen', 'reco efficiency denominator;p_{T}^{GEN} (GeV)', nx=30, x1=0., x2=300.)
hists.book1F('genEffNum', 'gen efficiency numerator;p_{T}^{RECO} (GeV)', nx=30, x1=0., x2=300.)
hists.book1F('genEffDen', 'gen efficiency denominator;p_{T}^{RECO} (GeV)', nx=30, x1=0., x2=300.)


############## Basic distributions ##############

hists.book1F('histAK7DeltaPhi',  ';#Delta #phi', nx=25,x1=0,x2=ROOT.TMath.Pi())
hists.book1F('histAK7DeltaY',    ';#Delta y', nx=25,x1=0,x2=ROOT.TMath.Pi()*2.0)

hists.book2F('histAK7PtAvgVsNvtx',  ';N_{VTX};p_{T}^{RECO} (GeV)',   nx=25,x1=0,x2=25, ny=280, y1=0, y2=7000)
hists.book2F('histAK7MjetVsNvtx',   ';N_{VTX};m_{jet}^{RECO} (GeV)', nx=25,x1=0,x2=25, ny=60, y1=0, y2=300)
hists.book2F('histAK7PtAvgVsMjetGroomOverReco',   ';Reconstructed m_{jet}^{GROOM}/m_{jet}^{UNGROOM};p_{T}^{AVG}', nx=51,x1=0.0,x2=1.02, ny=len(ptBins)-1, ybins=ptBins)
hists.book2F('histAK7PtAvgVsMjetGroomOverRecoTrue', ';True m_{jet}^{GROOM}/m_{jet}^{UNGROOM};p_{T}^{AVG}', nx=51,x1=0.0,x2=1.02, ny=len(ptBins)-1, ybins=ptBins)

コード例 #2

ファイル: smpjs_tuple_analyzer_fwlite_data.py プロジェクト: pvmulder/usercode

## hists.book2F('histAK7MjetVsNvtx',
##              'AK7 <m_{jet}> Versus p_{T}^{AVG} ;<m_{jet}> (GeV);p_{T}^{AVG} (GeV)',
##              nx=30, x1=0., x2=300., ny=len(ptBins)-1, ybins=ptBins)

## # Mjj and <Mjet> versus ptAvg for the different triggers. Only data. #
## for trig in trigHelper.trigsToKeep :
##     hists.book2F('histAK7PtAvgVsNvtx_' + trig,
##                  'AK7 m_{jj} Versus p_{T}^{AVG} ' + trig +';m_{jj} (GeV);p_{T}^{AVG} (GeV)',
##                  nx=70, x1=0., x2=7000., ny=len(ptBins)-1, ybins=ptBins)
##     hists.book2F('histAK7MjetVsNvtx_' + trig,
##                  'AK7 <m_{jet}> Versus p_{T}^{AVG} ;<m_{jet}> (GeV)' + trig + ';p_{T}^{AVG} (GeV)',
##                  nx=30, x1=0., x2=300., ny=len(ptBins)-1, ybins=ptBins)


############## Basic distributions ##############
hists.book1F('histAK7DeltaPhi',  ';#Delta #phi', nx=25,x1=0,x2=ROOT.TMath.Pi())
hists.book1F('histAK7DeltaY',    ';#Delta y', nx=25,x1=0,x2=ROOT.TMath.Pi()*2.0)

hists.book2F('histAK7PtAvgVsNvtx',  ';N_{VTX};p_{T}^{RECO} (GeV)',   nx=25,x1=0,x2=25, ny=280, y1=0, y2=7000)
hists.book2F('histAK7MjetVsNvtx',   ';N_{VTX};m_{jet}^{RECO} (GeV)', nx=25,x1=0,x2=25, ny=30, y1=0, y2=300)
hists.book2F('histAK7PtAvgVsMjetGroomOverReco',   ';m_{jet}^{GROOM}/m_{jet}^{RECO};p_{T}^{AVG}', nx=51,x1=0.0,x2=1.02, ny=len(ptBins)-1, ybins=ptBins)

# Mjj and <Mjet> versus ptAvg for the different triggers. Only data. #
for trig in trigHelper.trigsToKeep :
    hists.book1F('histAK7DeltaPhi_' + trig,  trig + ';#Delta #phi', nx=25,x1=0,x2=ROOT.TMath.Pi())
    hists.book1F('histAK7DeltaY_' + trig,    trig + ';#Delta y', nx=25,x1=0,x2=ROOT.TMath.Pi())
    hists.book2F('histAK7PtAvgVsNvtx_' + trig,
                 trig +';N_{VTX};p_{T}^{RECO} (GeV)',
                 nx=25,x1=0,x2=25, ny=280, y1=0, y2=7000)
    hists.book2F('histAK7MjetVsNvtx_' + trig,
                 trig +';N_{VTX};p_{T}^{RECO} (GeV)',

コード例 #3

ファイル: smpjs_tuple_analyzer_fwlite_ptJet.py プロジェクト: pvmulder/usercode

hists.book3F('histAK7MjetGenVsRecoVsPtJetHighPt', 'AK7 <m_{jet}> Response Matrix ;<m_{jet}>^{GEN} (GeV);<m_{jet}>^{RECO} (GeV);p_{T}^{AVG} (GeV)',
              nx=60, x1=0., x2=300., ny=60, y1=0., y2=300., nz=14, z1=200., z2=1600.)


# Mjet versus GENERATED/RECO Mjet. Only for MC. 
hists.book3F('histAK7MjetResponseVsPtJet', 'AK7 <m_{jet}> Versus p_{T}^{AVG} ;<m_{jet}> (GeV);p_{T}^{AVG} (GeV)',
              nx=60, x1=0., x2=300., ny=20, y1=0.0, y2=2.0, nz=10, z1=0., z2=1000.)

############## Responses ##############


# Pt and eta response. Only for MC
hists.book2F('histAK7PtResponse',  ';p_{T}^{GEN} (GeV);p_{T}^{RECO} / p_{T}^{GEN}', nx=50, x1=0., x2=500., ny=20, y1=0.5, y2=1.5)
hists.book2F('histAK7EtaResponse', ';#eta^{GEN} (GeV);p_{T}^{RECO} / p_{T}^{GEN}', nx=50, x1=-5.0, x2=5.0, ny=20, y1=0.5, y2=1.5)

hists.book1F('recoEffNum', 'reco efficiency numerator;p_{T}^{GEN} (GeV)', nx=60, x1=0., x2=300.)
hists.book1F('recoEffDen', 'reco efficiency denominator;p_{T}^{GEN} (GeV)', nx=60, x1=0., x2=300.)
hists.book1F('genEffNum', 'gen efficiency numerator;p_{T}^{RECO} (GeV)', nx=60, x1=0., x2=300.)
hists.book1F('genEffDen', 'gen efficiency denominator;p_{T}^{RECO} (GeV)', nx=60, x1=0., x2=300.)


hists.book1F('ptAsymmetry', '(p_{T}^{1} - p_{T}^{2}) / (p_{T}^{1} + p_{T}^{2})', nx=100, x1=-1, x2=1)

############## Basic distributions ##############

hists.book2F('histAK7PtJetVsNvtx',  ';N_{VTX};p_{T}^{RECO} (GeV)',   nx=25,x1=0,x2=25, ny=280, y1=0, y2=7000)
hists.book2F('histAK7MjetVsNvtx',   ';N_{VTX};m_{jet}^{RECO} (GeV)', nx=25,x1=0,x2=25, ny=60, y1=0, y2=300)
hists.book2F('histAK7PtJetVsMjetGroomOverReco',   ';m_{jet}^{GROOM}/m_{jet}^{RECO};p_{T}^{AVG}', nx=51,x1=0.0,x2=1.02, ny=len(ptBins)-1, ybins=ptBins)

vetoPtCut = 30.0

コード例 #4

ファイル: smpjs_tuple_analyzer_fwlite.py プロジェクト: pvmulder/usercode

# Mjj and <Mjet> versus GENERATED/RECO Mjj and <Mjet>. Only for MC. 
hists.book3F('histAK7MjjResponseVsPtAvg', 'AK7 m_{jj} Response;m_{jj} (GeV);p_{T}^{AVG} (GeV)',
              nx=70, x1=0., x2=7000., ny=80, y1=0.8, y2=1.2, nz=10, z1=0., z2=1000.)
hists.book3F('histAK7MjetResponseVsPtAvg', 'AK7 <m_{jet}> Versus p_{T}^{AVG} ;<m_{jet}> (GeV);p_{T}^{AVG} (GeV)',
              nx=60, x1=0., x2=300., ny=20, y1=0.0, y2=2.0, nz=10, z1=0., z2=1000.)

############## Responses ##############


# Pt and eta response. Only for MC
hists.book2F('histAK7PtResponse',  ';p_{T}^{GEN} (GeV);p_{T}^{RECO} / p_{T}^{GEN}', nx=50, x1=0., x2=500., ny=20, y1=0.5, y2=1.5)
hists.book2F('histAK7EtaResponse', ';#eta^{GEN} (GeV);p_{T}^{RECO} / p_{T}^{GEN}', nx=50, x1=-5.0, x2=5.0, ny=20, y1=0.5, y2=1.5)

hists.book1F('recoEffNum', 'reco efficiency numerator;p_{T}^{GEN} (GeV)', nx=60, x1=0., x2=300.)
hists.book1F('recoEffDen', 'reco efficiency denominator;p_{T}^{GEN} (GeV)', nx=60, x1=0., x2=300.)
hists.book1F('genEffNum', 'gen efficiency numerator;p_{T}^{RECO} (GeV)', nx=60, x1=0., x2=300.)
hists.book1F('genEffDen', 'gen efficiency denominator;p_{T}^{RECO} (GeV)', nx=60, x1=0., x2=300.)


mjjPtCut = 50.0
mjjEtaCut = 2.0

## for ibin in xrange(len(ptBins)-1) :
##     resp = ROOT.RooUnfoldResponse(60, 0., 300., 60, 0., 300.)
##     resp.SetName( 'response' + str(ibin))
##     responses.append(resp)

responses = []