---

How to run the code

Before to be able to run my likelihood code (run in DASH), I have to make sure that all the right variables are in the hbook files.


Need to make new zbs files containing
    - Polfit2 variables (need to run polfit2)
                /home/fdufour/work24-t2kk/polfit2
         /home/fdufour/work24-t2kk/mypolfit2

        In the .sh file we must have:

    For SK-I: (as in polfit2.sh)

    For SK-II: (as in sk2.sh)
 


    - Danka's variables (need to add epi0sep bank). This is a DASH code.
                /home/fdufour/work24-t2kk/make_zbs

 Once I have the new zbs files, I need to create new ntuples by running fillnt.F.
                    /home/fdufour/epi0/likelihood/ntuple/ (for the executable)
         /home/fdufour/work24-t2kk/t2kk/  (and work23 for scripts)
   
Only then can I run the likelihood code.

NB    for 100yrs of MC I need ??? of disk space
         for SK-I data I need ??? of disk space  



Likelihood code:

General structure:

Everything is located in: /work/fdufour/t2kk/t2k-like

Note: both parts of the code need precuts.F

Run the "training" code in order to create the tables needed.

    atm-tables.kumac
    atmsk2-tables.kumac
    t2k-tables.kumac

and the main code is ftables.F

Then I have to run the likelihood analysis itself.
        atm-like.kumac
    atmsk2-like.kumac
        t2k-like.kumac

and the main code is likelihood_7.F


Details to run properly:

To decide which variable to use in the likelihood, use the file
variable.txt where:
    1. efrac
    2. totpe
    3.xalong
    4.cosopen


The code is clearly not perfect right now and some paramters are hard-coded. So be careful to recompile the code, with the
parameters you want depeding on which normalization/MC sample you want. * Set SK version
      sk_version=1
      write(*,*) 'SK version =', sk_version

If running for SK MC, then choose between SK1 and SK2, if running T2K MC then chose sk_version=1


* Set normalization code:
      code=1





   

---

Official Plots & Tables

For plots going up to 10 GeV (only without CCQE and gif available) go there

Without CCQE hatched:

Atmospheric MC SK-I Go to atmplots	Atmospheric MC SK-II Go to atmsk2plots	T2K MC Go to T2K plots
Likelihood (total) (gif) (eps)	Likelihood (total) (gif) (eps)	Likelihood (total) (gif) (eps)
Likelihood (ebin) (gif) (eps)	Likelihood (ebin) (gif) (eps)	Likelihood (ebin) (gif) (eps)
nring (gif) (eps)	nring (gif) (eps)	nring (gif) (eps)
pid (gif) (eps)	pid (gif) (eps)	pid (gif) (eps)
cosnue (gif) (eps)	cosnue (gif) (eps)	cosnue (gif) (eps)
pi0 mass (gif) (eps)	pi0 mass (gif) (eps)	pi0 mass (gif) (eps)
pi0 like (gif) (eps)	pi0 like (gif) (eps)	pi0 like (gif) (eps)
energy fraction (gif) (eps)	energy fraction (gif) (eps)	energy fraction (gif) (eps)
totpe (gif) (eps)	totpe (gif) (eps)	totpe (gif) (eps)
xalong (gif) (eps)	xalong (gif) (eps)	xalong (gif) (eps)
cos(open) (gif) (eps)	cos(open) (gif) (eps)	cos(open) (gif) (eps)

With CCQE hatched:

Atmospheric MC SK-I Go to atmplots	Atmospheric MC SK-II Go to atmsk2plots	T2K MC Go to T2K plots
Likelihood (total) (gif) (eps)	Likelihood (total) (gif) (eps)	Likelihood (total) (gif) (eps)
Likelihood (ebin) (gif) (eps)	Likelihood (ebin) (gif) (eps)	Likelihood (ebin) (gif) (eps)
nring (gif) (eps)	nring (gif) (eps)	nring (gif) (eps)
pid (gif) (eps)	pid (gif) (eps)	pid (gif) (eps)
cosnue (gif) (eps)	cosnue (gif) (eps)	cosnue (gif) (eps)
pi0 mass (gif) (eps)	pi0 mass (gif) (eps)	pi0 mass (gif) (eps)
pi0 like (gif) (eps)	pi0 like (gif) (eps)	pi0 like (gif) (eps)
energy fraction (gif) (eps)	energy fraction (gif) (eps)	energy fraction (gif) (eps)
totpe (gif) (eps)	totpe (gif) (eps)	totpe (gif) (eps)
xalong (gif) (eps)	xalong (gif) (eps)	xalong (gif) (eps)
cos(open) (gif) (eps)	cos(open) (gif) (eps)	cos(open) (gif) (eps)

Efficiency tables:
   

Atmospheric MC SK-I	Atmospheric MC SK-II	T2K MC
5 variables		5 variables
5 + efrac		5 + efrac
5 + totpe		5 + totpe
5 + Xalong		5 + Xalong
5+ cospen		5+ cospen
5+ efrac+xalong+cosopen	5+ efrac+xalong+cosopen	5+ efrac+xalong+cosopen
All		All

S over sqrt (B) plots:     gif  or eps
Each points corresponds to:
    5 variables
   5 + efrac
    5+totpe
    5+xalong
    5+cosopen
    5+efrac+xalong+cosopen
    all


Updated plots: 0.35 < E < 0.85 gif or eps
Updated plots: 1.5 < E gif or eps

back to top