validation Namespace Reference

Class to run in order to validate ChannelInfo. More...

Classes
class	NoiseClusterAnalysis
class	ValidateChannelInfo
class	ValidateDAQ2RawDigit

Functions
def	setup_nova ()
def	setup_rootlogin ()
def	sendemail ()
def	mcvalidate (release1='S11.04.30', release2='S11.06.14', process='genie', horncurrent='nhc', swapping='none', det='ndos')
def	recovalidate (release1='S11.04.30', release2='S11.06.14', det='NDOS', beam='numi', directory='recohist')

Detailed Description

Class to run in order to validate ChannelInfo.

DAQ2RawDigit validation.

Author

nowak.nosp@m.@umn.nosp@m..edu

Jarek Nowak

Date

Function Documentation

def validation.mcvalidate	(		release1 = 'S11.04.30',
			release2 = 'S11.06.14',
			process = 'genie',
			horncurrent = 'nhc',
			swapping = 'none',
			det = 'ndos'
	)

Definition at line 121 of file validation.py.

References SetTitle().

Referenced by recovalidate().

def mcvalidate(release1='S11.04.30', release2='S11.06.14', process='genie',horncurrent='nhc', swapping='none', det='ndos'):

        histos = ('fVertexX',
                  'fVertexY',
                  'fVertexZ',
                  'fDCosX',
                  'fDCosY',
                  'fDCosZ',
                  'fMuMomentum',
                  'fMuDCosX',
                  'fMuDCosY',
                  'fMuDCosZ',
                  'fEMomentum',
                  'fEDCosX',
                  'fEDCosY',
                  'fEDCosZ',
                  'fECons',
                  'fVertexHadronDecayX',
                  'fVertexHadronDecayY',
                  'fVertexHadronDecayZ')
        
        gROOT.ForceStyle()
        c1 = TCanvas("c1","c1",10,10,800,600)
        gPad.SetLeftMargin(0.15)
        c1.Clear()
        c1.Divide(1,2)
        gStyle.SetPalette(1)
        gStyle.SetOptStat(000)
        gStyle.SetFillColor(10)
        
        c1.Print("%s_%s_%s_%s_%s_%s.pdf["%(release1,release2,process,horncurrent,swapping,det)) 

        output1 = TFileMerger()
        output1.OutputFile('%s.root'%release1)
        
        path = '/nova/data/mc/%s/%s/ndos/mccheckout'%(release1,process)
        n = 0
        for f in glob.glob(os.path.join(path,'*%s_*_%s_%s*.root'%(horncurrent,swapping,det))):
                n = n + 1
                if n <= 50:
                        output1.AddFile("%s"%f)
                else:
                        pass
        output1.Merge()
        path = '/nova/data/mc/%s/%s/ndos/mccheckout'%(release2,process)
        output2 = TFileMerger()
        output2.OutputFile('%s.root'%release2)
        n = 0
        for f in glob.glob(os.path.join(path,'*%s_*_%s_%s*.root'%(horncurrent,swapping,det))):
                n = n + 1
                if n <= 20:
                        output2.AddFile("%s"%f)
                else:
                        pass
        output2.Merge()

        for histo in histos:
        
                f_tag1 = TFile('%s.root'%release1)
                histo1 = f_tag1.Get("NeutrinoAna/%s"%histo)

                f_tag2 = TFile('%s.root'%release2)
                histo2 = f_tag2.Get("NeutrinoAna/%s"%histo)
                
                scale = histo2.Integral()/histo1.Integral()
                
                c1.cd(1)
                histo1.Scale(scale)
                histo1.SetLineColor(kBlue)
                histo1.SetFillColor(kBlue)
                histo1.SetFillStyle(3004)
                histo1.GetXaxis().SetTitle("%s"%histo)
                histo1.GetYaxis().SetTitleOffset(1.2)    
                histo1.GetYaxis().SetTitle("Events")
                histo1.DrawCopy()
                histo2.SetLineColor(kRed)
                histo2.DrawCopy("same")
                kolmogorov = histo1.KolmogorovTest(histo2)
                print kolmogorov
                
                legend1 = TLegend(0.70,0.80,0.89,0.89,"")
                legend1.AddEntry(histo1,"%s"%release1,"f")
                legend1.AddEntry(histo2,"%s"%release2,"l")
                legend1.Draw()
                
                c1.cd(2)
                histo1.Sumw2()
                histo2.Sumw2()
                ratio = histo1.Clone("ratio");
                ratio.Divide(ratio,histo2,1,1,"")
                ratio.SetLineColor(kBlue)
                ratio.GetXaxis().SetTitle("%s"%histo)
                ratio.GetYaxis().SetTitle("Ratio (%s / %s)"%(release2,release1))
                ratio.DrawCopy()
                
                legend2 = TLegend(0.70,0.80,0.89,0.89,"")
                legend2.AddEntry(ratio,"Ratio","l")
                legend2.Draw()
                c1.Print("%s_%s_%s_%s_%s_%s.pdf"%(release1,release2,process,horncurrent,swapping,det))  

        c1.Print("%s_%s_%s_%s_%s_%s.pdf]"%(release1,release2,process,horncurrent,swapping,det))                 

def validation.recovalidate	(		release1 = 'S11.04.30',
			release2 = 'S11.06.14',
			det = 'NDOS',
			beam = 'numi',
			directory = 'recohist'
	)

Definition at line 224 of file validation.py.

References BlessedPlots.main(), mcvalidate(), sendemail(), SetTitle(), setup_nova(), and setup_rootlogin().

def recovalidate(release1='S11.04.30', release2='S11.06.14', det='NDOS', beam='numi', directory='recohist'):
        
        histos = ('dcosX',
                  'dcosY',
                  'dcosZ',
                  'startX',
                  'startY',
                  'startZ',
                  'endX',
                  'endY',
                  'endZ',
                  'diffX',
                  'diffY',
                  'diffZ',
                  'dirX',
                  'dirY',
                  'fractionHitsUsed',
                  'fractionSigUsed',
                  'cellHitSig',
                  'endDistToEdge',
                  'vtxDistToEdge',
                  'hitXPosition',
                  'hitYPosition',
                  'hitZPosition',
                  'trackLengthInCell',
                  'trackLength',
                  'coscosmic',
                  'xRecoHitTNS',
                  'yRecoHitTNS')        

        gROOT.ForceStyle()
        c1 = TCanvas("c1","c1",10,10,800,600)
        gPad.SetLeftMargin(0.15)
        c1.Clear()
        c1.Divide(1,2)
        gStyle.SetPalette(1)
        gStyle.SetOptStat(000)
        gStyle.SetFillColor(10)
        
        c1.Print("%s_%s_%s_%s.pdf["%(release1,release2,det,directory))  

        output1 = TFileMerger()
        output1.OutputFile('%s.root'%release1)
        
        path = '/nova/data/novareco/%s/%s/%s'%(det,release1,directory)
        n = 0
        for f in glob.glob(os.path.join(path,'*%s*.root'%(beam))):
                n = n + 1
                if n <= 10:
                        output1.AddFile("%s"%f)
                else:
                        pass
        output1.Merge()
        path = '/nova/data/novareco/%s/%s/%s'%(det,release2,directory)
        output2 = TFileMerger()
        output2.OutputFile('%s.root'%release2)
        n = 0
        for f in glob.glob(os.path.join(path,'*%s*.root'%(beam))):
                n = n + 1
                if n <= 10:
                        output2.AddFile("%s"%f)
                else:
                        pass
        output2.Merge()

        for histo in histos:
        
                f_tag1 = TFile('%s.root'%release1)
                histo1 = f_tag1.Get("NeutrinoAna/%s"%histo)

                f_tag2 = TFile('%s.root'%release2)
                histo2 = f_tag2.Get("NeutrinoAna/%s"%histo)
                
                scale = histo2.Integral()/histo1.Integral()
                
                c1.cd(1)
                histo1.Scale(scale)
                histo1.SetLineColor(kBlue)
                histo1.SetFillColor(kBlue)
                histo1.SetFillStyle(3004)
                histo1.GetXaxis().SetTitle("%s"%histo)
                histo1.GetYaxis().SetTitleOffset(1.5)    
                histo1.GetYaxis().SetTitle("Events")
                histo1.DrawCopy()
                histo2.SetLineColor(kRed)
                histo2.DrawCopy("same")
                kolmogorov = histo1.KolmogorovTest(histo2)
                print kolmogorov
                
                legend1 = TLegend(0.70,0.80,0.89,0.89,"")
                legend1.AddEntry(histo1,"%s"%release1,"f")
                legend1.AddEntry(histo2,"%s"%release2,"l")
                legend1.Draw()
                
                c1.cd(2)
                histo1.Sumw2()
                histo2.Sumw2()
                ratio = histo1.Clone("ratio");
                ratio.Divide(ratio,histo2,1,1,"")
                ratio.SetLineColor(kBlue)
                ratio.GetXaxis().SetTitle("%s"%histo)
                ratio.GetYaxis().SetTitle("Ratio (%s / %s)"%(release2,release1))
                ratio.DrawCopy()
                
                legend2 = TLegend(0.70,0.80,0.89,0.89,"")
                legend2.AddEntry(ratio,"Ratio","l")
                legend2.Draw()
                c1.Print("%s_%s_%s_%s_%s.pdf"%(release1,release2,horncurrent,swapping,det))     

        c1.Print("%s_%s_%s_%s_%s.pdf]"%(release1,release2,horncurrent,swapping,det))   

def validation.sendemail

(

)

Definition at line 87 of file validation.py.

Referenced by recovalidate().

def sendemail():
        import smtplib

        SERVER = "localhost"
        
        FROM = "tianxc@gmail.com"
        TO = ["tianxc@fnal.gov"] 
# must be a list
        
        SUBJECT = "Today's auto-validation results"

        today                    = datetime.date.today()
        
        TEXT = "Today's auto-validation results: http://nusoft.fnal.gov/nova/novasoft/dailyValidation/daily_validation_%s.pdf"%(today)
        TEXT2 = "With K-S or chi^2 p-value less than 0.05: http://nusoft.fnal.gov/nova/novasoft/dailyValidation/daily_low_pValue_%s.pdf"%(today)
        
        # Prepare actual message
        
        message = """\
From: %s
To: %s
Subject: %s

%s
%s
        """ % (FROM, ", ".join(TO), SUBJECT, TEXT, TEXT2)
        
        # Send the mail
        
        server = smtplib.SMTP(SERVER)
        server.sendmail(FROM, TO, message)
        server.quit()
        
        

def validation.setup_nova

(

)

Definition at line 8 of file validation.py.

Referenced by recovalidate().

def setup_nova():
    print os.getenv("PWD")
    if os.getenv("SRT_PUBLIC_CONTEXT") and os.getenv("SRT_PRIVATE_CONTEXT"):
        print "The base and test release have been setup correctly"
    else:
        print "The base and test release have not been setup correctly"

def validation.setup_rootlogin

(

)

Definition at line 15 of file validation.py.

Referenced by recovalidate().

def setup_rootlogin():
    # Defaults to classic style, but that's OK, we can fix it
    novaStyle = TStyle("novaStyle", "NOvA Style");
    
    # Centre title
    novaStyle.SetTitleAlign(22);
    novaStyle.SetTitleX(.5);
    novaStyle.SetTitleY(.95);
    novaStyle.SetTitleBorderSize(0);
    
    # No info box
    novaStyle.SetOptStat(0);
    
    # set the background color to white
    novaStyle.SetFillColor(10);
    novaStyle.SetFrameFillColor(10);
    novaStyle.SetCanvasColor(10);
    novaStyle.SetPadColor(10);
    novaStyle.SetTitleFillColor(0);
    novaStyle.SetStatColor(10);
    
    # Don't put a colored frame around the plots
    novaStyle.SetFrameBorderMode(0);
    novaStyle.SetCanvasBorderMode(0);
    novaStyle.SetPadBorderMode(0);
    
    # Set the default line color for a fit function to be red
    novaStyle.SetFuncColor(kRed);
    
    # Marker settings
    #    novaStyle.SetMarkerStyle(kFullCircle);
    
    # No border on legends
    novaStyle.SetLegendBorderSize(0);
    
    # Scientific notation on axes
    #    TGaxis::SetMaxDigits(3);

    # Axis titles
    novaStyle.SetTitleSize(.05, "xyz");
    novaStyle.SetTitleOffset(.8, "xyz");
    # More space for y-axis to avoid clashing with big numbers
    novaStyle.SetTitleOffset(.9, "y");
    # This applies the same settings to the overall plot title
    novaStyle.SetTitleSize(.05, "");
    novaStyle.SetTitleOffset(.8, "");
    # Axis labels (numbering)
    novaStyle.SetLabelSize(.04, "xyz");
    novaStyle.SetLabelOffset(.005, "xyz");
    
    # Thicker lines
    novaStyle.SetHistLineWidth(2);
    novaStyle.SetFrameLineWidth(2);
    novaStyle.SetFuncWidth(2);
    
    # Set the number of tick marks to show
    novaStyle.SetNdivisions(506, "xyz");
    
    # Set the tick mark style
    novaStyle.SetPadTickX(1);
    novaStyle.SetPadTickY(1);
    
    # Fonts
    novaStyle.SetStatFont(42);
    novaStyle.SetLabelFont(42, "xyz");
    novaStyle.SetTitleFont(42, "xyz");
    novaStyle.SetTitleFont(42, ""); # Apply same setting to plot titles
    novaStyle.SetTextFont(42);
    #  novaStyle.SetLegendFont(kNovaFont);
    
    gROOT.SetStyle("novaStyle");