Added stdDev plots to Analysis and summary plot method in Postprocessor

createDelphesScripts.sh

@@ -0,0 +1,60 @@
+#!/bin/bash
+
+# Create fast simulation resolutions plotting scripts and submit to ifarm
+METHOD="Ele" # Switch this to select reconstruction method from {"Ele","JB","DA"}
+script="$PWD/macro/analysis_resolution.C"
+postScript="$PWD/macro/postprocess_resolution.C"
+submitScript="$PWD/submit.sh"
+jobScript="$PWD/job.sh"
+out="$PWD/macro/"
+cd datarec
+for file in *-xm25.config
+do
+    energies=`echo $file | sed "s;.config;;g"`
+    config="${METHOD}_${energies}"
+    mkdir -p $out/$config
+    cp $script $out/$config/
+    newscript=${out}${config}/*.C
+    eleIn=`echo $file | grep -Eo "[0-9][0-9]*x[0-9]*" | sed "s;x.*;;g"`
+    beamIn=`echo $file | grep -Eo "[0-9][0-9]*x[0-9]*" | sed "s;.*x;;g"`
+    xAng=`echo $file | grep -Eo "[0-9][0-9]*x[0-9]*-x[0-9]*" | sed "s;.*-x;;g"`
+    xAngM=`echo $file | grep -Eo "[0-9][0-9]*x[0-9]*-xm[0-9]*" | sed "s;.*-xm;;g"`
+    echo "file=${file}"
+    echo "energies=${energies}"
+    echo "config=${config}"
+    echo "newscript=${newscript}"
+    echo "eleIn=${eleIn}"
+    echo "beamIn=${beamIn}"
+    echo "xAng=${xAng}"
+    echo "xAngM=${xAngM}"
+
+    # Modify analysis script
+    sed -i "s;datarec/dis-5x41;datarec/${energies};g" $newscript
+    sed -i "s;Ele_dis-5x41;${config};g" $newscript
+    sed -i "s;eleBeamEn=5;eleBeamEn=${eleIn};g" $newscript
+    sed -i "s;ionBeamEn=41;ionBeamEn=${beamIn};g" $newscript
+    if [ $xAng ]; then
+        sed -i "s;crossingAngle=25;crossingAngle=-${xAng};g" $newscript
+    fi
+    if [ $xAngM ]; then
+        sed -i "s;crossingAngle=25;crossingAngle=${xAngM};g" $newscript
+    fi
+    sed -i "s;Ele;${METHOD};g" $newscript
+
+    # Postprocessor
+    cp $postScript $out/$config
+    sed -i "s;dis-5x41;${config};g" $out/$config/postprocess*.C
+	sed -i "s;testheader;${eleIn}x${beamIn}GeV;g" $out/$config/postprocess*.C
+
+    # And job scripts
+    cp $submitScript $out/$config
+    cp $jobScript $out/$config
+    sed -i "s;dis-5x41;${config};g" $out/$config/*.sh
+
+    # And submit
+	sbatch $out/$config/submit.sh 
+    echo --------------------
+
+done
+cd ..
+echo DONE

createFullScripts.sh

@@ -0,0 +1,48 @@
+#!/bin/bash
+# Create full simulation resolution plotting scripts and submit to ifarm
+script="$PWD/macro/analysis_resolution_SD_Full.C"
+postScript="$PWD/macro/postprocess_resolution_SD.C"
+submitScript="$PWD/submit.sh"
+jobScript="$PWD/job.sh"
+out="$PWD/macro/canyonlands/"
+cd datarec/canyonlands/
+for file in */*.config
+do
+    config=`echo $file | sed "s;.config;;g" | grep -Eo 'canyonlands/.*' | sed 's;canyonlands/;;g' | sed 's;/files;;g'`
+    mkdir -p $out/$config
+    cp $script $out/$config/
+    newscript=${out}${config}/*.C
+    eleIn=`echo $file | grep -Eo "[0-9][0-9]*x[0-9]*" | sed "s;x.*;;g"`
+    beamIn=`echo $file | grep -Eo "[0-9][0-9]*x[0-9]*" | sed "s;.*x;;g"`
+    xAng=`echo $file | grep -Eo "[0-9][0-9]*x[0-9]*-x[0-9]*" | sed "s;.*-x;;g"`
+    xAngM=`echo $file | grep -Eo "[0-9][0-9]*x[0-9]*-xm[0-9]*" | sed "s;.*-xm;;g"`
+    echo "name=${name}"
+    echo "file=${file}"
+    echo "config=${config}"
+    echo "newscript=${newscript}"
+    echo "eleIn=${eleIn}"
+    echo "beamIn=${beamIn}"
+    echo "xAng=${xAng}"
+    echo "xAngM=${xAngM}"
+
+    sed -i "s;datarec/dis-5x41.config;${file};g" $newscript
+    sed -i "s;full_dis-5x41;full_${config};g" $newscript
+    sed -i "s;eleBeamEn=5;eleBeamEn=${eleIn};g" $newscript
+    sed -i "s;ionBeamEn=41;ionBeamEn=${beamIn};g" $newscript
+    if [ $xAng ]; then
+        sed -i "s;crossingAngle=25;crossingAngle=-${xAng};g" $newscript
+    fi
+    if [ $xAngM ]; then
+        sed -i "s;icrossingAngle=25;crossingAngle=${xAngM};g" $newscript
+    fi
+    cp $postScript $out/$config
+    sed -i "s;dis-5x41;${config};g" $out/$config/*.C
+	sed -i "s;testheader;${eleIn}x${beamIn}GeV;g" $out/$config/*.C
+    cp $submitScript $out/$config
+    cp $jobScript $out/$config
+    sed -i "s;dis-5x41;canyonlands/${config};g" $out/$config/*.sh
+    sbatch $out/$config/submit.sh
+    echo --------------------
+done
+cd ..
+echo DONE

datarec/dis-5x41-xm25.config

@@ -0,0 +1,3 @@
+1 4.002e-4 /lustre19/expphy/volatile/halla/solid/dcbyer/eic-data/new/dis-5x41-xm25.root            -1
+100 2.869e-7 /lustre19/expphy/volatile/halla/solid/dcbyer/eic-data/new/dis-5x41-q2min100-xm25.root -1
+

job.sh

@@ -0,0 +1,6 @@
+#!/bin/bash
+
+cd /work/clas12/users/$USER/eic3/largex-eic
+echo "root -q -b macro/dis-5x41/analysis_resolution.C" | ./container/shell.sh
+echo "root -q -b macro/dis-5x41/postprocess_resolution.C" | ./container/shell.sh
+echo DONE

macro/analysis_barak.C

@@ -0,0 +1,47 @@
+R__LOAD_LIBRARY(Largex)
+
+// make resolution plots
+void analysis_barak(
+    TString infiles="datarec/dis-18x275-xm25.config", /* list of input files */
+    Double_t eleBeamEn=18, /* electron beam energy [GeV] */
+    Double_t ionBeamEn=275, /* ion beam energy [GeV] */
+    Double_t crossingAngle=25, /* crossing angle [mrad] */
+    TString outfilePrefix="barak_dis-18x275-xm25" /* output filename prefix*/
+) {
+
+  // setup analysis ========================================
+  AnalysisDelphes *A = new AnalysisDelphes(
+      infiles,
+      eleBeamEn,
+      ionBeamEn,
+      crossingAngle,
+      outfilePrefix
+      );
+  A->NBINS = 1; // use this to set the number of bins along each axis, e.g., z binning (except resolution axes) for each overall bin in e.g. x and Q2
+  A->NBINSRES = 100; // use this to set the number of bins along the resolution axes for each overall bin in e.g. x and Q2
+  // A->RESHIGH = 1;
+  // A->RESLOW = -1;
+  //A->maxEvents = 30000; // use this to limit the number of events
+  A->SetReconMethod("Ele"); // set reconstruction method
+  // A->SetReconMethod("JB"); // set reconstruction method
+  //A->SetReconMethod("DA"); // set reconstruction method
+  A->AddFinalState("pipTrack"); // pi+ final state
+  A->AddFinalState("pimTrack"); // pi- final state
+  A->AddFinalState("KpTrack"); // K+ final state
+  A->AddFinalState("KmTrack"); // K- final state
+  A->AddFinalState("pTrack"); // proton final state
+
+  // define cuts ====================================
+  A->AddBinScheme("w");  A->BinScheme("w")->BuildBin("Min",3.0); // W > 3 GeV
+  A->AddBinScheme("y");  A->BinScheme("y")->BuildBin("Range",0.01,0.95); // 0.01 < y < 0.95
+  A->AddBinScheme("z");  A->BinScheme("z")->BuildBin("Range",0.2,0.9); // 0.2 < z < 0.9
+  A->AddBinScheme("xF"); A->BinScheme("xF")->BuildBin("Min",0.0); // xF > 0
+  A->AddBinScheme("ptLab");  A->BinScheme("ptLab")->BuildBin("Min",0.1); // pT_lab > 0.1 GeV (tracking limit)
+
+  // set binning scheme ====================================
+  A->AddBinScheme("q2"); A->BinScheme("q2")->BuildBins( 25, 0.1,    10000,  true );
+  A->AddBinScheme("x");  A->BinScheme("x")->BuildBins(  25, 0.00001, 1,    true );
+
+  // perform the analysis ==================================
+  A->Execute();
+};

macro/analysis_resolution.C

@@ -0,0 +1,45 @@
+R__LOAD_LIBRARY(Largex)
+
+// make resolution plots
+void analysis_resolution(
+    TString infiles="datarec/dis-5x41.config", /* list of input files */
+    Double_t eleBeamEn=5, /* electron beam energy [GeV] */
+    Double_t ionBeamEn=41, /* ion beam energy [GeV] */
+    Double_t crossingAngle=25, /* crossing angle [mrad] */
+    TString outfilePrefix="Ele_dis-5x41" /* output filename prefix*/
+) {
+
+  // setup analysis ========================================
+  AnalysisDelphes *A = new AnalysisDelphes(
+      infiles,
+      eleBeamEn,
+      ionBeamEn,
+      crossingAngle,
+      outfilePrefix
+      );
+  A->NBINS = 1; // use this to set the number of bins along each axis, e.g., z binning (except resolution axes) for each overall bin in e.g. x and Q2
+  A->NBINSRES = 100; // use this to set the number of bins along the resolution axes for each overall bin in e.g. x and Q2
+  // A->RESHIGH = 1;
+  // A->RESLOW = -1;
+  //A->maxEvents = 30000; // use this to limit the number of events
+  A->SetReconMethod("Ele"); // set reconstruction method
+  A->AddFinalState("pipTrack"); // pion final state
+  // A->AddFinalState("pimTrack"); // pion final state
+  // A->AddFinalState("KpTrack"); // pion final state
+  // A->AddFinalState("KmTrack"); // pion final state
+  // A->AddFinalState("pTrack"); // pion final state
+
+  // define cuts ====================================
+  A->AddBinScheme("w");  A->BinScheme("w")->BuildBin("Min",3.0); // W > 3 GeV
+  A->AddBinScheme("y");  A->BinScheme("y")->BuildBin("Range",0.01,0.95); // 0.01 < y < 0.95
+  A->AddBinScheme("z");  A->BinScheme("z")->BuildBin("Range",0.2,0.9); // 0.2 < z < 0.9
+  A->AddBinScheme("xF"); A->BinScheme("xF")->BuildBin("Min",0.0); // xF > 0
+  A->AddBinScheme("ptLab");  A->BinScheme("ptLab")->BuildBin("Min",0.1); // pT_lab > 0.1 GeV (tracking limit)
+
+  // set binning scheme ====================================
+  A->AddBinScheme("q2"); A->BinScheme("q2")->BuildBins( 4, 1,    1000,  true );
+  A->AddBinScheme("x");  A->BinScheme("x")->BuildBins(  6, 0.01, 1,    true );
+
+  // perform the analysis ==================================
+  A->Execute();
+};

macro/analysis_resolution_SD_Full.C

@@ -0,0 +1,41 @@
+R__LOAD_LIBRARY(Largex)
+
+// make resolution plots
+void analysis_resolution_SD_Full(
+    TString infiles="datarec/full_dis-5x41.config", /* list of input files */
+    Double_t eleBeamEn=5, /* electron beam energy [GeV] */
+    Double_t ionBeamEn=41, /* ion beam energy [GeV] */
+    Double_t crossingAngle=25, /* crossing angle [mrad] */
+    TString outfilePrefix="full_dis-5x41" /* output filename prefix*/
+) {
+
+  // setup analysis ========================================
+  AnalysisDD4hep *A = new AnalysisDD4hep(
+      infiles,
+      eleBeamEn,
+      ionBeamEn,
+      crossingAngle,
+      outfilePrefix
+      );
+  A->NBINS = 5; // use this to set the number of bins along each axis, e.g., z binning (except resolution axes) for each overall bin in e.g. x and Q2
+  A->NBINSRES = 100; // use this to set the number of bins along the resolution axes for each overall bin in e.g. x and Q2
+  // A->RESHIGH = 1;
+  // A->RESLOW = -1;
+  //A->maxEvents = 30000; // use this to limit the number of events
+  A->SetReconMethod("Ele"); // set reconstruction method
+  A->AddFinalState("pipTrack"); // pion final state
+
+  // define cuts ====================================
+  A->AddBinScheme("w");  A->BinScheme("w")->BuildBin("Min",3.0); // W > 3 GeV
+  A->AddBinScheme("y");  A->BinScheme("y")->BuildBin("Range",0.01,0.95); // 0.01 < y < 0.95
+  A->AddBinScheme("z");  A->BinScheme("z")->BuildBin("Range",0.2,0.9); // 0.2 < z < 0.9
+  A->AddBinScheme("xF"); A->BinScheme("xF")->BuildBin("Min",0.0); // xF > 0
+  A->AddBinScheme("ptLab");  A->BinScheme("ptLab")->BuildBin("Min",0.1); // pT_lab > 0.1 GeV (tracking limit)
+
+  // set binning scheme ====================================
+  A->AddBinScheme("q2"); A->BinScheme("q2")->BuildBins( 4, 1,    100,  true );
+  A->AddBinScheme("x");  A->BinScheme("x")->BuildBins(  6, 0.01, 1,    true );
+
+  // perform the analysis ==================================
+  A->Execute();
+};