WebSVN – f9daq – Diff – /belle2/masterclass/th1fit.php


<?php
 
function send_message($id, $message, $progress) {
    $d = array('message' => $message , 'progress' => $progress); //prepare json
    $x = array('id' => $id, 'data' => $d );
    echo json_encode($x) . PHP_EOL;
/*
   if (PHP_SAPI !== 'cli') {
     echo PHP_EOL;
     ob_flush();
   }
   flush();
*/
}
 
 
function executeCmd($commandLine) {
 
  $pipe = popen("$commandLine" , 'r');
 
  if (!$pipe) {
    print "pipe failed.";
    return "";
  }
  $output = '';
  $out    = '';
  while(!feof($pipe)) {
    $contents = fread($pipe, 12);
    $unpacked = unpack("l*",$contents);
    $len0      = sizeof($unpacked);
    if ($len0 > 2) {
      $id       =  $unpacked[1];
      $len      =  $unpacked[2];
      $progress =  $unpacked[3];
      if ($len>0){
        $out= fread($pipe, $len);
        send_message($id,$out,$progress);
      }
      $output .=$out;
    }
  }
  pclose($pipe);
  return $output;
}
 
 if (isset($_POST["min"])){
   $min =  $_POST["min"];
 } else {
   $min =  0;
 }
 
 if (isset($_POST["max"])){
   $max =  $_POST["max"];
 } else {
   $max =  10;
 }
 
 if (isset($_POST["fitfun"])){
   $fitfun =  $_POST["fitfun"];
 } else {
   $fitfun =  "gaus";
 }
 
 if (isset($_POST["prm"])){
   $prm =  $_POST["prm"];
 } else {
   $prm =  "9,11,,99,3,,7";
 }
 
 if (isset($_POST["name"])){
   $name =  $_POST["name"];
 } else {
   $name = "h101";
 }
 
 if (isset($_POST["data"])){
   $data =  $_POST["data"];
 } else {
   $data = '{ "_typename" : "TH1F", "fUniqueID" : 0, "fBits" : 50331656, "fName" : "h100", "fTitle" : "Number of particles in the event", "fLineColor" : 602, "fLineStyle" : 1, "fLineWidth" : 1, "fFillColor" : 0, "fFillStyle" : 1001, "fMarkerColor" : 1, "fMarkerStyle" : 1, "fMarkerSize" : 1, "fNcells" : 52, "fXaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "xaxis", "fTitle" : "N particles", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 1, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 50, "fXmin" : -0.5, "fXmax" : 49.5, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fYaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "yaxis", "fTitle" : "N events", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 0, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 1, "fXmin" : 0, "fXmax" : 1, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fZaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "zaxis", "fTitle" : "", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 1, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 1, "fXmin" : 0, "fXmax" : 1, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fBarOffset" : 0, "fBarWidth" : 1000, "fEntries" : 10000, "fTsumw" : 9983, "fTsumw2" : 9983, "fTsumwx" : 175205, "fTsumwx2" : 3620667, "fMaximum" : -1111, "fMinimum" : -1111, "fNormFactor" : 0, "fContour" : [], "fSumw2" : [], "fOption" : "", "fFunctions" : { "_typename" : "TList", "name" : "TList", "arr" : [], "opt" : [] }, "fBufferSize" : 0, "fBuffer" : [], "fBinStatErrOpt" : 0, "fArray" : [0, 0, 0, 0, 86, 122, 117, 135, 176, 263, 306, 367, 412, 508, 596, 643, 618, 566, 524, 552, 503, 439, 418, 374, 333, 311, 229, 222, 188, 172, 123, 120, 92, 82, 81, 53, 47, 32, 36, 28, 27, 17, 14, 9, 10, 6, 8, 6, 4, 5, 3, 17] }';
 }
 
$json  = json_decode($data, true);
$class = $json['_typename'];
$hname = $json['fName'];
$htitle= $json['fTitle'];
$xaxis = $json['fXaxis'];
$xtitle= $xaxis['fTitle'];
$xbins = $xaxis['fNbins'];
$xmin  = $xaxis['fXmin'];
$xmax  = $xaxis['fXmax'];
$yaxis = $json['fYaxis'];
$ytitle= $yaxis['fTitle'];
 
 
$hdata='';
foreach( $json['fArray'] as $key => $value){
        if($key != "0"){
            $hdata .= ", ";
        }
        $hdata .= $value;
}
 
$cmd0  = ".L th1fit.cc+". PHP_EOL;
$cmd0 .= "double data[]={". $hdata . "};" .PHP_EOL;
$cmd0 .= "th1fit(data, \"$hname\",\"$htitle;$xtitle;$ytitle\", $xbins , $xmin, $xmax, $min ,$max ,\"$fitfun\",\"$prm\");". PHP_EOL;
 
if (file_exists('/opt/root/bin/thisroot.sh')){
  $profile ='/opt/root/bin/thisroot.sh';
}
$cmd  ="bash -c 'source $profile ; root -b -l <<EOL" . PHP_EOL;
$cmd .= $cmd0. PHP_EOL;
$cmd .= "EOL'". PHP_EOL;
executeCmd($cmd);
//echo $cmd0;
//echo "name=$name min=$min max=$max func=$fitfun data=$data";
?>
 
 

Rev 268	Rev 269
1	<?php	1	<?php
2		2
3	function send_message($id, $message, $progress) {	3	function send_message($id, $message, $progress) {
4	$d = array('message' => $message , 'progress' => $progress); //prepare json	4	$d = array('message' => $message , 'progress' => $progress); //prepare json
5	$x = array('id' => $id, 'data' => $d );	5	$x = array('id' => $id, 'data' => $d );
6	echo json_encode($x) . PHP_EOL;	6	echo json_encode($x) . PHP_EOL;
7	/*	7	/*
8	if (PHP_SAPI !== 'cli') {	8	if (PHP_SAPI !== 'cli') {
9	echo PHP_EOL;	9	echo PHP_EOL;
10	ob_flush();	10	ob_flush();
11	}	11	}
12	flush();	12	flush();
13	*/	13	*/
14	}	14	}
15		15
16		16
17	function executeCmd($commandLine) {	17	function executeCmd($commandLine) {
18		18
19	$pipe = popen("$commandLine" , 'r');	19	$pipe = popen("$commandLine" , 'r');
20		20
21	if (!$pipe) {	21	if (!$pipe) {
22	print "pipe failed.";	22	print "pipe failed.";
23	return "";	23	return "";
24	}	24	}
25	$output = '';	25	$output = '';
26	$out = '';	26	$out = '';
27	while(!feof($pipe)) {	27	while(!feof($pipe)) {
28	$contents = fread($pipe, 12);	28	$contents = fread($pipe, 12);
29	$unpacked = unpack("l*",$contents);	29	$unpacked = unpack("l*",$contents);
30	$len0 = sizeof($unpacked);	30	$len0 = sizeof($unpacked);
31	if ($len0 > 2) {	31	if ($len0 > 2) {
32	$id = $unpacked[1];	32	$id = $unpacked[1];
33	$len = $unpacked[2];	33	$len = $unpacked[2];
34	$progress = $unpacked[3];	34	$progress = $unpacked[3];
35	if ($len>0){	35	if ($len>0){
36	$out= fread($pipe, $len);	36	$out= fread($pipe, $len);
37	send_message($id,$out,$progress);	37	send_message($id,$out,$progress);
38	}	38	}
39	$output .=$out;	39	$output .=$out;
40	}	40	}
41	}	41	}
42	pclose($pipe);	42	pclose($pipe);
43	return $output;	43	return $output;
44	}	44	}
45		45
46	if (isset($_POST["min"])){	46	if (isset($_POST["min"])){
47	$min = $_POST["min"];	47	$min = $_POST["min"];
48	} else {	48	} else {
49	$min = 0;	49	$min = 0;
50	}	50	}
51		51
52	if (isset($_POST["max"])){	52	if (isset($_POST["max"])){
53	$max = $_POST["max"];	53	$max = $_POST["max"];
54	} else {	54	} else {
55	$max = 10;	55	$max = 10;
56	}	56	}
57		57
58	if (isset($_POST["fitfun"])){	58	if (isset($_POST["fitfun"])){
59	$fitfun = $_POST["fitfun"];	59	$fitfun = $_POST["fitfun"];
60	} else {	60	} else {
61	$fitfun = "gaus";	61	$fitfun = "gaus";
62	}	62	}
63		63
64	if (isset($_POST["prm"])){	64	if (isset($_POST["prm"])){
65	$prm = $_POST["prm"];	65	$prm = $_POST["prm"];
66	} else {	66	} else {
67	$prm = "~~0;1;2~~";	67	$prm = "9,11,,99,3,,7";
68	}	68	}
69		69
70	if (isset($_POST["name"])){	70	if (isset($_POST["name"])){
71	$name = $_POST["name"];	71	$name = $_POST["name"];
72	} else {	72	} else {
73	$name = "h101";	73	$name = "h101";
74	}	74	}
75		75
76	if (isset($_POST["data"])){	76	if (isset($_POST["data"])){
77	$data = $_POST["data"];	77	$data = $_POST["data"];
78	} else {	78	} else {
79	$data = '{ "_typename" : "TH1F", "fUniqueID" : 0, "fBits" : 50331656, "fName" : "h100", "fTitle" : "Number of particles in the event", "fLineColor" : 602, "fLineStyle" : 1, "fLineWidth" : 1, "fFillColor" : 0, "fFillStyle" : 1001, "fMarkerColor" : 1, "fMarkerStyle" : 1, "fMarkerSize" : 1, "fNcells" : 52, "fXaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "xaxis", "fTitle" : "N particles", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 1, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 50, "fXmin" : -0.5, "fXmax" : 49.5, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fYaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "yaxis", "fTitle" : "N events", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 0, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 1, "fXmin" : 0, "fXmax" : 1, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fZaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "zaxis", "fTitle" : "", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 1, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 1, "fXmin" : 0, "fXmax" : 1, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fBarOffset" : 0, "fBarWidth" : 1000, "fEntries" : 10000, "fTsumw" : 9983, "fTsumw2" : 9983, "fTsumwx" : 175205, "fTsumwx2" : 3620667, "fMaximum" : -1111, "fMinimum" : -1111, "fNormFactor" : 0, "fContour" : [], "fSumw2" : [], "fOption" : "", "fFunctions" : { "_typename" : "TList", "name" : "TList", "arr" : [], "opt" : [] }, "fBufferSize" : 0, "fBuffer" : [], "fBinStatErrOpt" : 0, "fArray" : [0, 0, 0, 0, 86, 122, 117, 135, 176, 263, 306, 367, 412, 508, 596, 643, 618, 566, 524, 552, 503, 439, 418, 374, 333, 311, 229, 222, 188, 172, 123, 120, 92, 82, 81, 53, 47, 32, 36, 28, 27, 17, 14, 9, 10, 6, 8, 6, 4, 5, 3, 17] }';	79	$data = '{ "_typename" : "TH1F", "fUniqueID" : 0, "fBits" : 50331656, "fName" : "h100", "fTitle" : "Number of particles in the event", "fLineColor" : 602, "fLineStyle" : 1, "fLineWidth" : 1, "fFillColor" : 0, "fFillStyle" : 1001, "fMarkerColor" : 1, "fMarkerStyle" : 1, "fMarkerSize" : 1, "fNcells" : 52, "fXaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "xaxis", "fTitle" : "N particles", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 1, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 50, "fXmin" : -0.5, "fXmax" : 49.5, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fYaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "yaxis", "fTitle" : "N events", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 0, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 1, "fXmin" : 0, "fXmax" : 1, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fZaxis" : { "_typename" : "TAxis", "fUniqueID" : 0, "fBits" : 50331648, "fName" : "zaxis", "fTitle" : "", "fNdivisions" : 510, "fAxisColor" : 1, "fLabelColor" : 1, "fLabelFont" : 42, "fLabelOffset" : 0.005, "fLabelSize" : 0.035, "fTickLength" : 0.03, "fTitleOffset" : 1, "fTitleSize" : 0.035, "fTitleColor" : 1, "fTitleFont" : 42, "fNbins" : 1, "fXmin" : 0, "fXmax" : 1, "fXbins" : [], "fFirst" : 0, "fLast" : 0, "fBits2" : 0, "fTimeDisplay" : false, "fTimeFormat" : "", "fLabels" : null, "fModLabs" : null }, "fBarOffset" : 0, "fBarWidth" : 1000, "fEntries" : 10000, "fTsumw" : 9983, "fTsumw2" : 9983, "fTsumwx" : 175205, "fTsumwx2" : 3620667, "fMaximum" : -1111, "fMinimum" : -1111, "fNormFactor" : 0, "fContour" : [], "fSumw2" : [], "fOption" : "", "fFunctions" : { "_typename" : "TList", "name" : "TList", "arr" : [], "opt" : [] }, "fBufferSize" : 0, "fBuffer" : [], "fBinStatErrOpt" : 0, "fArray" : [0, 0, 0, 0, 86, 122, 117, 135, 176, 263, 306, 367, 412, 508, 596, 643, 618, 566, 524, 552, 503, 439, 418, 374, 333, 311, 229, 222, 188, 172, 123, 120, 92, 82, 81, 53, 47, 32, 36, 28, 27, 17, 14, 9, 10, 6, 8, 6, 4, 5, 3, 17] }';
80	}	80	}
81		81
82	$json = json_decode($data, true);	82	$json = json_decode($data, true);
83	$class = $json['_typename'];	83	$class = $json['_typename'];
84	$hname = $json['fName'];	84	$hname = $json['fName'];
85	$htitle= $json['fTitle'];	85	$htitle= $json['fTitle'];
86	$xaxis = $json['fXaxis'];	86	$xaxis = $json['fXaxis'];
87	$xtitle= $xaxis['fTitle'];	87	$xtitle= $xaxis['fTitle'];
88	$xbins = $xaxis['fNbins'];	88	$xbins = $xaxis['fNbins'];
89	$xmin = $xaxis['fXmin'];	89	$xmin = $xaxis['fXmin'];
90	$xmax = $xaxis['fXmax'];	90	$xmax = $xaxis['fXmax'];
91	$yaxis = $json['fYaxis'];	91	$yaxis = $json['fYaxis'];
92	$ytitle= $yaxis['fTitle'];	92	$ytitle= $yaxis['fTitle'];
93		93
94		94
95	$hdata='';	95	$hdata='';
96	foreach( $json['fArray'] as $key => $value){	96	foreach( $json['fArray'] as $key => $value){
97	if($key != "0"){	97	if($key != "0"){
98	$hdata .= ", ";	98	$hdata .= ", ";
99	}	99	}
100	$hdata .= $value;	100	$hdata .= $value;
101	}	101	}
102		102
103	$cmd0 = ".L th1fit.cc+". PHP_EOL;	103	$cmd0 = ".L th1fit.cc+". PHP_EOL;
104	$cmd0 .= "double data[]={". $hdata . "};" .PHP_EOL;	104	$cmd0 .= "double data[]={". $hdata . "};" .PHP_EOL;
105	$cmd0 .= "th1fit(data, \"$hname\",\"$htitle;$xtitle;$ytitle\", $xbins , $xmin, $xmax, $min ,$max ,\"$fitfun\");". PHP_EOL;	105	$cmd0 .= "th1fit(data, \"$hname\",\"$htitle;$xtitle;$ytitle\", $xbins , $xmin, $xmax, $min ,$max ,\"$fitfun\",\"$prm\");". PHP_EOL;
106		106
107	if (file_exists('/opt/root/bin/thisroot.sh')){	107	if (file_exists('/opt/root/bin/thisroot.sh')){
108	$profile ='/opt/root/bin/thisroot.sh';	108	$profile ='/opt/root/bin/thisroot.sh';
109	}	109	}
110	$cmd ="bash -c 'source $profile ; root -b -l <<EOL" . PHP_EOL;	110	$cmd ="bash -c 'source $profile ; root -b -l <<EOL" . PHP_EOL;
111	$cmd .= $cmd0. PHP_EOL;	111	$cmd .= $cmd0. PHP_EOL;
112	$cmd .= "EOL'". PHP_EOL;	112	$cmd .= "EOL'". PHP_EOL;
113	executeCmd($cmd);	113	executeCmd($cmd);
114	//echo $cmd0;	114	//echo $cmd0;
115	//echo "name=$name min=$min max=$max func=$fitfun data=$data";	115	//echo "name=$name min=$min max=$max func=$fitfun data=$data";
116	?>	116	?>
117		117
118		118

Subversion Repositories f9daq

(root)/belle2/masterclass/th1fit.php @ 329 – Rev 268 → 269