Skip to content

Commit

Permalink
Enabled folding of backgrounds from random trigger. Will work as soon…
Browse files Browse the repository at this point in the history 
… as the files appear.
  • Loading branch information
@T-Britton
T-Britton committed Mar 31, 2017
1 parent 9e3ca69 commit 7503538
Show file tree
Hide file tree
Showing 4 changed files with 337 additions and 230 deletions.
221 changes: 137 additions & 84 deletions MCwrapper/MakeMC.csh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -58,21 +58,23 @@ setenv NUMTHREADS 1
endif

# PRINT INPUTS
echo "CONTEXT = $JANA_CALIB_CONTEXT"
echo "ENVIRONMENT = $ENVIRONMENT"
echo "CONFIG_FILE = $CONFIG_FILE"
echo "OUTDIR = $OUTDIR"
echo "RUN_NUMBER = $RUN_NUMBER"
echo "FILE_NUMBER = $FILE_NUMBER"
echo "NUM TO GEN = $EVT_TO_GEN"
echo "generator = $GENERATOR"
echo "generator = $GENERATOR"
echo "generation = $GENR $CLEANGENR"
echo "Geant = $GEANT $CLEANGEANT"
echo "GCONTROL = $CUSTOM_GCONTROL"
echo "BKG_FOLD = $BKGFOLDSTR"
echo "MCsmear = $SMEAR $CLEANSMEAR"
echo "Recon = $RECON $CLEANRECON"
echo "Geant = $GEANT $CLEANGEANT"
echo "GCONTROL = $CUSTOM_GCONTROL"
echo "BKG_FOLD = $BKGFOLDSTR"
echo "MCsmear = $SMEAR $CLEANSMEAR"
echo "Recon = $RECON $CLEANRECON"

echo "detected c-shell"

#printenv
#necessary to run swif, uses local directory if swif=0 is used
if ("$MCSWIF" == "1") then
Expand All @@ -91,6 +93,25 @@ else
cp $CUSTOM_GCONTROL/Gcontrol.in ./
endif


@ length_count = `echo $RUN_NUMBER | wc -c` - 1

set formatted_runNumber=""
while ( $length_count < 6 )
set formatted_runNumber = "0""$formatted_runNumber"
@ length_count = $length_count + 1
end
set formatted_runNumber=$formatted_runNumber$RUN_NUMBER

@ flength_count = `echo $FILE_NUMBER | wc -c` - 1

set formatted_fileNumber=""
while ( $flength_count < 3 )
set formatted_fileNumber = "0""$formatted_fileNumber"
@ flength_count = $flength_count + 1
end
set formatted_fileNumber=$formatted_fileNumber$FILE_NUMBER

if ("$GENR" != "0") then
if ("$GENERATOR" != "genr8" && "$GENERATOR" != "bggen" && "$GENERATOR" != "genEtaRegge" && "$GENERATOR" != "gen_2pi_amp" && "$GENERATOR" != "gen_pi0" && "$GENERATOR" != "gen_2pi_primakoff" ) then
echo "NO VALID GENERATOR GIVEN"
Expand All @@ -107,13 +128,13 @@ if ("$GENR" != "0") then

if ("$GENERATOR" == "genr8") then
echo "configuring genr8"
cp $CONFIG_FILE ./genr8\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
cp $CONFIG_FILE ./genr8\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
else if ("$GENERATOR" == "bggen") then
echo "configuring bggen"
cp $MCWRAPPER_CENTRAL/Generators/bggen/particle.dat ./
cp $MCWRAPPER_CENTRAL/Generators/bggen/pythia.dat ./
cp $MCWRAPPER_CENTRAL/Generators/bggen/pythia-geant.map ./
cp $CONFIG_FILE ./bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
cp $CONFIG_FILE ./bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
if ( `echo $eBEAM_ENERGY | grep -o "\." | wc -l` == 0) then
set eBEAM_ENERGY = $eBEAM_ENERGY\.
endif
Expand All @@ -129,129 +150,161 @@ if ("$GENR" != "0") then

else if ("$GENERATOR" == "genEtaRegge") then
echo "configuring genEtaRegge"
cp $CONFIG_FILE ./genEtaRegge\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
cp $CONFIG_FILE ./genEtaRegge\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
else if ("$GENERATOR" == "gen_2pi_amp") then
echo "configuring gen_2pi_amp"
cp $CONFIG_FILE ./gen_2pi_amp\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
cp $CONFIG_FILE ./gen_2pi_amp\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
else if ("$GENERATOR" == "gen_2pi_primakoff") then
echo "configuring gen_2pi_primakoff"
cp $CONFIG_FILE ./gen_2pi_primakoff\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
cp $CONFIG_FILE ./gen_2pi_primakoff\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
else if ("$GENERATOR" == "gen_pi0") then
echo "configuring gen_pi0"
cp $CONFIG_FILE ./gen_pi0\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
cp $CONFIG_FILE ./gen_pi0\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
endif
set config_file_name=`basename "$CONFIG_FILE"`
echo $config_file_name

if ("$GENERATOR" == "genr8") then
echo "RUNNING GENR8"
set RUNNUM = $RUN_NUMBER+$FILE_NUMBER
sed -i 's/TEMPCOHERENT/'$COHERENT_PEAK'/' genr8\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
set RUNNUM = $formatted_runNumber+$formatted_fileNumber
sed -i 's/TEMPCOHERENT/'$COHERENT_PEAK'/' genr8\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
# RUN genr8 and convert
genr8 -r$RUN_NUMBER -M$EVT_TO_GEN -A$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.ascii < genr8\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf #$config_file_name
genr8_2_hddm $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.ascii
genr8 -r$formatted_runNumber -M$EVT_TO_GEN -A$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.ascii < genr8\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf #$config_file_name
genr8_2_hddm $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.ascii
else if ("$GENERATOR" == "bggen") then
set colsize=`rcnd $RUN_NUMBER collimator_diameter | awk '{print $1}' | sed -r 's/.{2}$//' | sed -e 's/\.//g'`
if ("$colsize" == "B" || "$colsize" == "R" ) then
if ("$colsize" == "B" || "$colsize" == "R" || "$JANA_CALIB_CONTEXT" != "variation=mc") then
set colsize = "34"
endif
set RANDOM=$$
echo $RANDOM
sed -i 's/TEMPTRIG/'$EVT_TO_GEN'/' bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPRUNNO/'$RUN_NUMBER'/' bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPCOLD/'0.00$colsize'/' bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPRAND/'$RANDOM'/' bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPELECE/'$eBEAM_ENERGY'/' bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPCOHERENT/'$COHERENT_PEAK'/' bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPMINGENE/'$GEN_MIN_ENERGY'/' bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPMAXGENE/'$GEN_MAX_ENERGY'/' bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPTRIG/'$EVT_TO_GEN'/' bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPRUNNO/'$RUN_NUMBER'/' bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPCOLD/'0.00$colsize'/' bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPRAND/'$RANDOM'/' bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPELECE/'$eBEAM_ENERGY'/' bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPCOHERENT/'$COHERENT_PEAK'/' bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPMINGENE/'$GEN_MIN_ENERGY'/' bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPMAXGENE/'$GEN_MAX_ENERGY'/' bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf

ln -s bggen\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf fort.15
ln -s bggen\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf fort.15
bggen
mv bggen.hddm $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm
mv bggen.hddm $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm
else if ("$GENERATOR" == "genEtaRegge") then
echo "RUNNING GENETAREGGE"
set colsize=`rcnd $RUN_NUMBER collimator_diameter | awk '{print $1}' | sed -r 's/.{2}$//' | sed -e 's/\.//g'`
if ("$colsize" == "B" || "$colsize" == "R" ) then
if ("$colsize" == "B" || "$colsize" == "R" || "$JANA_CALIB_CONTEXT" != "variation=mc") then
set colsize = "34"
endif
sed -i 's/TEMPCOLD/'0.00$colsize'/' genEtaRegge\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPELECE/'$eBEAM_ENERGY'/' genEtaRegge\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPCOHERENT/'$COHERENT_PEAK'/' genEtaRegge\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPMINGENE/'$GEN_MIN_ENERGY'/' genEtaRegge\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPMAXGENE/'$GEN_MAX_ENERGY'/' genEtaRegge\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
genEtaRegge -N$EVT_TO_GEN -O$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm -I'genEtaRegge'\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf
sed -i 's/TEMPCOLD/'0.00$colsize'/' genEtaRegge\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPELECE/'$eBEAM_ENERGY'/' genEtaRegge\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPCOHERENT/'$COHERENT_PEAK'/' genEtaRegge\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPMINGENE/'$GEN_MIN_ENERGY'/' genEtaRegge\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
sed -i 's/TEMPMAXGENE/'$GEN_MAX_ENERGY'/' genEtaRegge\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
genEtaRegge -N$EVT_TO_GEN -O$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm -I'genEtaRegge'\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf
else if ("$GENERATOR" == "gen_2pi_amp") then
echo "RUNNING GEN_2PI_AMP"
set optionals_line = `head -n 1 $config_file_name | sed -r 's/.//'`
echo $optionals_line
echo gen_2pi_amp -c gen_2pi_amp\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf -o $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm -hd $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY -b $GEN_MAX_ENERGY $optionals_line
gen_2pi_amp -c gen_2pi_amp\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf -hd $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm -o $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY - b $GEN_MAX_ENERGY $optionals_line
echo gen_2pi_amp -c gen_2pi_amp\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf -o $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm -hd $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY -b $GEN_MAX_ENERGY $optionals_line
gen_2pi_amp -c gen_2pi_amp\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf -hd $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm -o $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY - b $GEN_MAX_ENERGY $optionals_line
else if ("$GENERATOR" == "gen_2pi_primakoff") then
echo "RUNNING GEN_2PI_PRIMAKOFF"
set optionals_line = `head -n 1 $config_file_name | sed -r 's/.//'`
echo $optionals_line
echo gen_2pi_primakoff -c gen_2pi_primakoff\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf -o $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm -hd $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY -b $GEN_MAX_ENERGY $optionals_line
gen_2pi_primakoff -c gen_2pi_primakoff\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf -hd $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm -o $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY - b $GEN_MAX_ENERGY $optionals_line
echo gen_2pi_primakoff -c gen_2pi_primakoff\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf -o $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm -hd $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY -b $GEN_MAX_ENERGY $optionals_line
gen_2pi_primakoff -c gen_2pi_primakoff\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf -hd $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm -o $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY - b $GEN_MAX_ENERGY $optionals_line
else if ("$GENERATOR" == "gen_pi0") then
echo "RUNNING GEN_PI0"
set optionals_line = `head -n 1 $config_file_name | sed -r 's/.//'`
echo $optionals_line
gen_pi0 -c gen_pi0\_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.conf -hd $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm -o $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY -b $GEN_MAX_ENERGY -p $COHERENT_PEAK -s $FILE_NUMBER $optionals_line -m $eBEAM_ENERGY
gen_pi0 -c gen_pi0\_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.conf -hd $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm -o $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.root -n $EVT_TO_GEN -r $RUN_NUMBER -a $GEN_MIN_ENERGY -b $GEN_MAX_ENERGY -p $COHERENT_PEAK -s $formatted_fileNumber $optionals_line -m $eBEAM_ENERGY
endif

#GEANT/smearing
#modify TEMPIN/TEMPOUT/TEMPTRIG/TOSMEAR in control.in

if ("$GEANT" != "0") then
echo "RUNNING GEANT"$GEANTVER
set colsize=`rcnd $RUN_NUMBER collimator_diameter | awk '{print $1}' | sed -r 's/.{2}$//' | sed -e 's/\.//g'`
if ("$colsize" == "B" || "$colsize" == "R" ) then
set colsize = "34"
endif
if ("$GEANT" != "0") then
echo "RUNNING GEANT"$GEANTVER
set colsize=`rcnd $RUN_NUMBER collimator_diameter | awk '{print $1}' | sed -r 's/.{2}$//' | sed -e 's/\.//g'`
if ("$colsize" == "B" || "$colsize" == "R" || "$JANA_CALIB_CONTEXT" != "variation=mc") then
set colsize = "34"
endif

if ( `echo $eBEAM_ENERGY | grep -o "\." | wc -l` == 0) then
set eBEAM_ENERGY = $eBEAM_ENERGY\.
endif
if ( `echo $COHERENT_PEAK | grep -o "\." | wc -l` == 0) then
set COHERENT_PEAK = $COHERENT_PEAK\.
endif
if ( `echo $eBEAM_ENERGY | grep -o "\." | wc -l` == 0) then
set eBEAM_ENERGY = $eBEAM_ENERGY\.
endif
if ( `echo $COHERENT_PEAK | grep -o "\." | wc -l` == 0) then
set COHERENT_PEAK = $COHERENT_PEAK\.
endif

set inputfile=$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber
cp Gcontrol.in $PWD/control'_'$formatted_runNumber'_'$formatted_fileNumber.in
chmod 777 $PWD/control'_'$formatted_runNumber'_'$formatted_fileNumber.in
sed -i 's/TEMPELECE/'$eBEAM_ENERGY'/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in
sed -i 's/TEMPCOHERENT/'$COHERENT_PEAK'/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in
sed -i 's/TEMPIN/'$inputfile.hddm'/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in
sed -i 's/TEMPRUNG/'$RUN_NUMBER'/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in
sed -i 's/TEMPOUT/'$inputfile'_geant.hddm/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in
sed -i 's/TEMPTRIG/'$EVT_TO_GEN'/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in
sed -i 's/TEMPCOLD/'0.00$colsize'/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in

if ( "$BKGFOLDSTR" != "BeamPhotons" ) then
echo "removing Beam Photon background from geant simulation"
sed -i 's/BGRATE/cBGRATE/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in
sed -i 's/BGGATE/cBGGATE/' control'_'$formatted_runNumber'_'$formatted_fileNumber.in
endif

set inputfile=$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER
cp Gcontrol.in $PWD/control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
chmod 777 $PWD/control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
sed -i 's/TEMPELECE/'$eBEAM_ENERGY'/' control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
sed -i 's/TEMPCOHERENT/'$COHERENT_PEAK'/' control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
sed -i 's/TEMPIN/'$inputfile.hddm'/' control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
sed -i 's/TEMPRUNG/'$RUN_NUMBER'/' control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
sed -i 's/TEMPOUT/'$inputfile'_geant.hddm/' control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
sed -i 's/TEMPTRIG/'$EVT_TO_GEN'/' control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
#sed -i 's/TOSMEAR/'$SMEAR'/' control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
sed -i 's/TEMPCOLD/'0.00$colsize'/' control'_'$RUN_NUMBER'_'$FILE_NUMBER.in
mv $PWD/control'_'$RUN_NUMBER'_'$FILE_NUMBER.in $PWD/control.in
mv $PWD/control'_'$formatted_runNumber'_'$formatted_fileNumber.in $PWD/control.in

if ("$GEANTVER" == "3") then
hdgeant
else if ("$GEANTVER" == "4") then
#make run.mac then call it below
rm -f run.mac
echo "/run/beamOn $EVT_TO_GEN" > run.mac
echo "exit" >>! run.mac
hdgeant4 -t$NUMTHREADS run.mac
rm run.mac
else
echo "INVALID GEANT VERSION"
exit
endif
if ("$GEANTVER" == "3") then
hdgeant
else if ("$GEANTVER" == "4") then
#make run.mac then call it below
rm -f run.mac
echo "/run/beamOn $EVT_TO_GEN" > run.mac
echo "exit" >>! run.mac
hdgeant4 -t$NUMTHREADS run.mac
rm run.mac
else
echo "INVALID GEANT VERSION"
exit
endif

if ("$SMEAR" != "0") then
echo "RUNNING MCSMEAR"

if ("$BKGFOLDSTR" == "0") then
mcsmear -o$inputfile'_geant_smeared.hddm' $inputfile'_geant.hddm'
if ("$BKGFOLDSTR" == "BeamPhotons" || "$BKGFOLDSTR" == "None") then
echo "running MCsmear without folding in random background"
mcsmear -o$inputfile'_geant_smeared.hddm' $inputfile'_geant.hddm'
else
echo 'mcsmear -o'$inputfile'_geant_smeared.hddm'' '$inputfile'_geant.hddm'' '$BKGFOLDSTR
mcsmear -o$inputfile'_geant_smeared.hddm' $inputfile'_geant.hddm' $BKGFOLDSTR
if("$BKGFOLDSTR" == "DEFAULT") then
#find file and run:1
echo "Finding the right file to fold in"
set runperiod="RunPeriod-2017-01"

if ( $RUN_NUMBER > 40000) then
#set runperiod="RunPeriod-2017-10"
endif

if ( $RUN_NUMBER < 30000 ) then
echo "Warning: random triggers did not exist by this point"
endif
set bkglocstring="/cache/halld/""$runperiod""/sim/random_triggers/""$formatted_runNumber"".hddm"
#set bkglocstring="/w/halld-scifs1a/home/tbritton/converted.hddm"

if (! -f $bkglocstring) then
echo "Could not find mix-in file "$bkglocstring
exit
else
echo "mcsmear -o$inputfile"\_"geant"\_"smeared.hddm $inputfile"\_"geant.hddm $bkglocstring"\:"1"
mcsmear -o$inputfile\_geant\_smeared.hddm $inputfile\_geant.hddm $bkglocstring\:1
endif

else
#trust the user and use their string
echo 'mcsmear -o'$inputfile'_geant_smeared.hddm'' '$inputfile'_geant.hddm'' '$BKGFOLDSTR
mcsmear -o$inputfile'_geant_smeared.hddm' $inputfile'_geant.hddm' $BKGFOLDSTR
endif
endif
#run reconstruction
if ("$CLEANGENR" == "1") then
Expand All @@ -264,14 +317,14 @@ if ("$GENR" != "0") then
unlink fort.15
endif

rm $GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm
rm $GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm
endif

if ("$RECON" != "0") then
echo "RUNNING RECONSTRUCTION"
hd_root $inputfile'_geant_smeared.hddm' -PPLUGINS=danarest,monitoring_hists -PNTHREADS=$NUMTHREADS
mv hd_root.root hd_root_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.root
mv dana_rest.hddm dana_rest_$GEN_NAME\_$RUN_NUMBER\_$FILE_NUMBER.hddm
mv hd_root.root hd_root_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.root
mv dana_rest.hddm dana_rest_$GEN_NAME\_$formatted_runNumber\_$formatted_fileNumber.hddm

if ("$CLEANGEANT" == "1") then
rm *_geant.hddm
Expand Down
Loading

0 comments on commit 7503538

Please sign in to comment.