Modify cuts definition for the HMS calrimeter calibration.

In hcal_calib/hcal_replay_cuts.def:
   add one_clust condition;
   activate good_beta condition (defined for run 52949);
   define 2 versions of the Reconstruct_master cut, one for initial
   calibration and another for iterative calibration.

examples/output_52949.def

@@ -12,6 +12,9 @@ TH1F hfptime1x ' HMS s1x fptime' H.hod.1x.fptime 80 0. 80.
 TH1F hfptime1y ' HMS s1y fptime' H.hod.1y.fptime 80 0. 80.
 TH1F hfptime2x ' HMS s2x fptime' H.hod.2x.fptime 80 0. 80.
 TH1F hfptime2y ' HMS s2y fptime' H.hod.2y.fptime 80 0. 80.
+# Beta
+TH1F hfpBeta      'HMS Hodo Beta'  H.hod.fpBeta      250 -0.5 2.0
+TH1F hfpBetaChisq 'HMS Hodo Chisq' H.hod.fpBetaChisq 550 -5. 50.
 # TDC hits per paddle
 TH1F hpostdc1 'HMS s1x+ TDC hits' H.hod.1x.postdchits 16 0.5 16.5
 TH1F hnegtdc1 'HMS s1x- TDC hits' H.hod.1x.negtdchits 16 0.5 16.5

examples/output_mkj.def

@@ -12,6 +12,9 @@ TH1F hfptime1x ' HMS s1x fptime' H.hod.1x.fptime 80 0. 80.
 TH1F hfptime1y ' HMS s1y fptime' H.hod.1y.fptime 80 0. 80.
 TH1F hfptime2x ' HMS s2x fptime' H.hod.2x.fptime 80 0. 80.
 TH1F hfptime2y ' HMS s2y fptime' H.hod.2y.fptime 80 0. 80.
+# Beta
+TH1F hfpBeta      'HMS Hodo Beta'  H.hod.fpBeta      250 -0.5 2.0
+TH1F hfpBetaChisq 'HMS Hodo Chisq' H.hod.fpBetaChisq 550 -5. 50.
 # TDC hits per paddle
 TH1F hpostdc1 'HMS s1x+ TDC hits' H.hod.1x.postdchits 16 0.5 16.5
 TH1F hnegtdc1 'HMS s1x- TDC hits' H.hod.1x.negtdchits 16 0.5 16.5

hcal_calib/hcal_calib.cpp

@@ -16,7 +16,7 @@ void hcal_calib(Int_t RunNumber) {
  theShowerCalib.ComposeVMs();      // Compute vectors amd matrices for calib.
  theShowerCalib.SolveAlphas();     // Solve for the calibration constants
  theShowerCalib.SaveAlphas();      // Save the constants
- // theShowerCalib.SaveRawData();  // Save raw data into file for debug purposes
+ theShowerCalib.SaveRawData();  // Save raw data into file for debug purposes
  theShowerCalib.FillHEcal();       // Fill histograms
 
  // Plot histograms

hcal_calib/hcal_replay_cuts.def

@@ -24,17 +24,23 @@ RawCoarseReconstruct !Pedestal_event
 Block: Reconstruct
 
 one_track H.tr.n==1
-#one_clust H.cal.nclust==1
+one_clust H.cal.nclust==1
 one_sh_track H.cal.ntracks==1
 in_delta  H.tr.tg_dp>-10.&&H.tr.tg_dp<10.
 good_cer H.cer.npesum>3.
-good_beta H.tr.beta>0.9&&H.tr.beta<1.1
-in_calx   H.cal.trx>-60.&&H.cal.trx<60.
-in_caly   H.cal.try>-30.&&H.cal.try<30.
+#good_beta H.tr.beta>0.9&&H.tr.beta<1.1
+good_beta H.hod.fpBeta>0.740&&H.hod.fpBeta<0.935   # 3sigma cut for run 52949
+in_calx   H.cal.trx>-65.4&&H.cal.trx<54.6    #top+5cm, bottom-5cm (see hcal.pos)
+in_caly   H.cal.try>-30.&&H.cal.try<30.      #left+5cm, right-5cm.(see hcal.pos)
 in_cal    in_calx&&in_caly
 
-# This version matches engine
+# This version is for calibration from scratch (no calibration constants exist,
+# first time calibration).
+Reconstruct_master one_track && one_clust && in_delta && good_cer && good_beta && in_cal
+
+# This version can be used for iterative calibration (improve existing
+# constants).
 #Reconstruct_master one_track && one_sh_track && in_delta && good_cer && good_beta && in_cal
 
-# Drop beta, which is not present yet
-Reconstruct_master one_track && one_sh_track && in_delta && good_cer && in_cal
+
+

hcal_calib/output_hcal_replay.def

@@ -4,7 +4,8 @@
 
 block H.cal.*
 variable H.cer.npesum
-variable H.tr.beta
+#variable H.tr.beta
+variable H.hod.fpBeta
 variable H.tr.p
 variable H.tr.tg_dp
 variable H.tr.ph	#tan(phi), wrt Y axis