Skip to content
Snippets Groups Projects
Commit e53a157b authored by Stephen A. Wood's avatar Stephen A. Wood
Browse files

Update examples for recent hcana changes

  New required Cerenkov and Aerogel parameters in hcana.def
  cuts and output.def updates
  remove now unsupported mirror efficiencies from report.template
parent 5abda424
No related branches found
No related tags found
No related merge requests found
......@@ -82,3 +82,73 @@ sdc_fix_propcorr = 0
# Total number of PMTs in Gas Cherenkov detector.
hcer_tot_pmts = 2
#
# Aerogel and Cerenkov parameters added in 3/31/2017
#
h_aero_six_gev_data = 1
; Track matching reduced chi2 cut values
haero_red_chi2_min = 0.0
haero_red_chi2_max = 25.0
; Track matching beta cut values
haero_beta_min = 0.0
haero_beta_max = 1.2
; Track matching normailized energy (E/p) cut values
haero_enorm_min = 0.0
haero_enorm_max = 1.5
; ADC time window cut values used to select "good" ADC events
haero_adcTimeWindowMin = 500.
haero_adcTimeWindowMax = 2500.
; NPE Threshold for "good" hit (for measureing effic.)
haero_npe_thresh = 0.5
; Value of zpos for diffusion box cited by Hamlet M. 03/29/2017
; FIX ME! This value is for the SHMS!
haero_diff_box_zpos = 231.0
; Number of regions in which tracks can be matched too
haero_num_regions = 1
; phgcer_region: 8 values for each region
; central x,y,dx,dy values and x,y,dx,dy half widths.
; Defined in clomuns, region 1, region 2, region 3, region 4...
; FIX ME! These values are for the SHMS!
; Assumes rectangular volume with x = 103.0 cm and y = 113.0 cm
haero_region = 0.0,
0.0,
0.0,
0.0,
51.5,
56.5,
0.15,
0.1
hcer_red_chi2_min = 0.0
hcer_red_chi2_max = 25.0
hcer_enorm_min = 0.0
hcer_enorm_max = 1.5
hcer_adcTimeWindowMin = 1500.
hcer_adcTimeWindowMax = 3500.
; Threshold for good hit (for measureing effic.)
hcer_npe_thresh = 0.5
hcer_num_regions = 3 ; This value is from the 6 GeV era
; hcer_region: 8 values for each region (1 per mirror + sum)
; central x,y,dx,dy values and x,y,dx,dy half widths.
hcer_region = 30, -30, 0,
0, 0, 0,
0, 0, 0,
0, 0, 0,
20, 20, 60,
30, 30, 30,
.1, .1, .1,
.1, .1, .1
......@@ -124,8 +124,8 @@ RawCoarseReconstruct !Pedestal_event
Block: Physics
HMSScinGood H.hod.goodscinhit == 1
HMSScinShould HMSScinGood && H.cal.etotnorm > 0.7 && H.cer.npesum > 2.0
HMSScinShoulde HMSScinGood && H.cer.npesum > 2.0
HMSScinShould HMSScinGood && H.cal.etotnorm > 0.7 && H.cer.npeSum > 2.0
HMSScinShoulde HMSScinGood && H.cer.npeSum > 2.0
HMSScinShouldh H.hod.goodscinhit == 1 && g.evtyp == 1
HMSScinDid HMSScinShould && H.dc.ntrack > 0
......@@ -133,8 +133,8 @@ HMSScinDide HMSScinShoulde && H.dc.ntrack > 0
HMSScinDidh HMSScinShouldh && H.dc.ntrack > 0
BHMSScinShould HMSScinShould && H.hod.betanotrack >= 0.1 && g.evtyp==3
BHMSScinShoulde BHMSScinShould && H.cer.npesum > 0.2 && ( H.gold.beta > ( H.hod.betap - 0.5 ) ) && ( H.gold.beta < ( H.hod.betap + 0.5 ) )
BHMSScinShouldh BHMSScinShould && H.cer.npesum > 0.2
BHMSScinShoulde BHMSScinShould && H.cer.npeSum > 0.2 && ( H.gold.beta > ( H.hod.beta - 0.5 ) ) && ( H.gold.beta < ( H.hod.beta + 0.5 ) )
BHMSScinShouldh BHMSScinShould && H.cer.npeSum > 0.2
BHMSScinDid BHMSScinShould && H.dc.ntrack > 0
BHMSScinDide BHMSScinShoulde && H.dc.ntrack > 0
......
TdcCounter# See $ANALYZER/examples/output_example.def for examples
# See $ANALYZER/examples/output_example.def for examples
#
block H.dc.*
block H.hod.*
......
......@@ -385,12 +385,7 @@ Good Track Requirements:
4. {hcer_et_min:%6.2f} < calorimeter energy < {hcer_et_max:%6.2f}
Threshold set at {hcer_threshold:%6.2f}
Mirror 1 Efficiency = {H.cer.cerfiredcounter[0]/( H.cer.certrackcounter[0] + 0.0001 ):%5.3f}
Mirror 1 {H.cer.cerfiredcounter[0]:%6d} of {H.cer.certrackcounter[0]:%6d} good tracks.
Mirror 2 Efficiency = {H.cer.cerfiredcounter[1]/( H.cer.certrackcounter[1] + 0.0001 ):%5.3f}
Mirror 2 {H.cer.cerfiredcounter[1]:%6d} of {H.cer.certrackcounter[1]:%6d} good tracks.
Overall Efficiency = {H.cer.cerfiredcounter[2]/( H.cer.certrackcounter[2] + 0.0001 ):%5.3f}
Overall {H.cer.cerfiredcounter[2]:%6d} of {H.cer.certrackcounter[2]:%6d} good tracks.
Mirror efficiencies removed. Not supported in latest THcAerogel
****************************
* Calorimeter Efficiencies *
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment