EMCal-analysis/macro/Likelihood.sh

0001 #!/bin/tcsh -f
0002
0003 if ($#argv < 2) then
0004         echo "Missing parameters. Example usages:"
0005         echo "$0 eta0.90 /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/ "
0006         echo "$0 eta0 /phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/ "
0007
0008         exit;
0009 endif
0010
0011 echo "$0 $*"
0012
0013 set base_dir =  $argv[2]
0014 # set base_dir = '/phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/emcstudies/pidstudies/spacal2d/fieldmap/';
0015 # set base_dir = '/phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal2d/fieldmap/';
0016 #set base_dir = '/phenix/sim02/phnxreco/ePHENIX/jinhuang/sPHENIX_work/production_analysis/embedding/emcstudies/pidstudies/spacal1d/fieldmap/';
0017
0018 #set kine_config = 'eta0';
0019 #set kine_config = 'eta0.90';
0020 set kine_config = $argv[1];
0021
0022 # set pids = ('e-' 'pi-');
0023 #set pids = ('e-' 'pi-' 'kaon-' 'anti_proton');
0024 set pids = ('e-' 'pi-' 'kaon-' 'anti_proton' 'e+' 'pi+' 'kaon+' 'proton');
0025
0026 set energys = (2GeV 4GeV 8GeV);
0027 # set energys = (2GeV 4GeV );
0028 # set energys = (8GeV );
0029
0030 set q = '"';
0031
0032
0033 foreach energy ($energys)
0034         foreach pid ($pids)
0035
0036                 mkdir ${base_dir}/G4Hits_sPHENIX_${pid}_${kine_config}_${energy}-ALL.root_Ana.root.lst.old/;
0037
0038                 mv -fv ${base_dir}/G4Hits_sPHENIX_${pid}_${kine_config}_${energy}-ALL.root_Ana.root.lst* ${base_dir}/G4Hits_sPHENIX_${pid}_${kine_config}_${energy}-ALL.root_Ana.root.lst.old/;
0039
0040                 /bin/ls -1 ${base_dir}/G4Hits_sPHENIX_${pid}_${kine_config}_${energy}-[0-9][0-9][0-9][0-9].root_Ana.root   > ${base_dir}/G4Hits_sPHENIX_${pid}_${kine_config}_${energy}-ALL.root_Ana.root.lst;
0041                 wc -l ${base_dir}/G4Hits_sPHENIX_${pid}_${kine_config}_${energy}-ALL.root_Ana.root.lst;
0042
0043         end
0044 end
0045
0046 # exit
0047
0048 foreach energy ($energys)
0049         foreach pid ($pids)
0050
0051                 echo root -b -q "Fun4All_EMCalLikelihood.C(1000, ${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q})";
0052                 root -b -q "Fun4All_EMCalLikelihood.C(1000, ${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q})" & ;
0053 # echo root -b -q "Fun4All_EMCalLikelihood.C(10000000, ${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q})";
0054 #root -b -q "Fun4All_EMCalLikelihood.C(10000000, ${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q})" & ;
0055         end
0056 end
0057                 wait;
0058
0059
0060 foreach energy ($energys)
0061         foreach pid ($pids)
0062
0063                 echo root -b -q "DrawEcal_Likelihood.C(${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q},0 )";
0064                 root -b -q "DrawEcal_Likelihood.C(${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q},0 )" & ;
0065
0066         end
0067 end
0068                 wait;
0069
0070
0071 foreach energy ($energys)
0072         foreach pid ($pids)
0073
0074                 echo root -b -q "Fun4All_EMCalLikelihood.C(1000, ${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q})";
0075
0076                 root -b -q "Fun4All_EMCalLikelihood.C(1000, ${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q})" & ;
0077
0078         end
0079 end
0080                 wait;
0081
0082 foreach energy ($energys)
0083         foreach pid ($pids)
0084
0085                 echo root -b -q "DrawEcal_Likelihood.C(${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q},1 )"  ;
0086                 root -b -q "DrawEcal_Likelihood.C(${q}${base_dir}${q}, ${q}${pid}${q}, ${q}${kine_config}_${energy}${q},1 )" & ;
0087
0088         end
0089 end
0090
0091                 wait;
0092
0093 echo done;
0094 date;
0095
0096 exit;
0097
0098
0099
0100
0101