moon_position_table_new3 Namespace Reference

Functions
def	SetupAnalysisRegion (AnalysisRegions)
def	AnalRegion (a, theta)
def	FillAnal (a, region, time, resolution)
def	tl (date)
def	computeDayBoundaryMoon (AnalysisRegions, day, theta, resolution)

Variables
	f = TFile('semfe_tmp.root','RECREATE')
	ntp_tr = TTree('ntp','My test tree')
	s = MyStruct()
string	version = "1.0"
string	usage = "usage: %prog [options]"
	parser = OptionParser(usage=usage,version=version)
	dest
	default
	help
	options
	args
	resolution = float(options.res)
	utc = pytz.utc
	central = timezone('US/Central')
	moon = ephem.Moon()
	sun = ephem.Sun()
	nova = ephem.Observer()
	lat
	lon
int	minute = 60
int	second = 1
int	hour = 60
float	oneMinute = 1.0
float	oneSecond = 1.0
	start_day = options.startdate
	numdays = int(options.numdays)
dictionary	AnalysisRegionsSun = {}
dictionary	AnalysisRegionsMoon = {}
int	last_region = 999
int	theta = 10
	day = ephem.Date(start_day)
	Loop over the days we care about. More...

Function Documentation

def moon_position_table_new3.AnalRegion	(		a,
			theta
	)

Definition at line 104 of file moon_position_table_new3.py.

Referenced by computeDayBoundaryMoon().

def AnalRegion(a, theta) :
    dayBoundary   = theta* math.pi / 180.0 
    nightBoundary = -1 *theta * math.pi / 180.0
    type = 0

        
    if a < nightBoundary :
        type=1
    if a > dayBoundary :
        type=2
    if (a < dayBoundary) and (a > nightBoundary) :
        type=3

#    print '%f %f, %f %d' %(dayBoundary, nightBoundary, a, type)

    return type

def moon_position_table_new3.computeDayBoundaryMoon	(		AnalysisRegions,
			day,
			theta,
			resolution
	)

Definition at line 149 of file moon_position_table_new3.py.

References AnalRegion(), FillAnal(), PandAna.Demos.demo1.range, and submit_syst.str.

def computeDayBoundaryMoon(AnalysisRegions, day, theta, resolution):
    global oneMinute
    AnalysisRegions['night_time']=0
    AnalysisRegions['day_time']=0
    AnalysisRegions['twilight_time']=0
#    for i in xrange(1,int(86400/resolution)):
#        nova.date = day + (oneSecond*resolution * i)
    for i in xrange(1,1440):
        nova.date = day + oneMinute * i
        #moon.compute(nova)
        #region = AnalRegion(moon.alt, theta)
        sun.compute(nova)
        region = AnalRegion(sun.alt, theta)
        FillAnal(AnalysisRegions, region, nova.date, resolution)

    # Print the transition times in UTC
    print '%20s, %20s, %20s, %20s' % ((AnalysisRegions['n2t_time_boundary']),(AnalysisRegions['t2d_time_boundary']),
                                      (AnalysisRegions['d2t_time_boundary']),(AnalysisRegions['t2n_time_boundary']))

    tmp_1 = str(AnalysisRegions['n2t_time_boundary']) + "UTC"
    proc_1 = subprocess.Popen(["/grid/fermiapp/products/nova/externals/novadaq/v04_05_00/slf6.x86_64.e9.debug/bin/NovaTimeConvert",
                               tmp_1], stdout=subprocess.PIPE)
    semfe_1 = proc_1.stdout.read()
    proc_1.wait()

    tmp_2 = str(AnalysisRegions['t2d_time_boundary']) + "UTC"
    proc_2 = subprocess.Popen(["/grid/fermiapp/products/nova/externals/novadaq/v04_05_00/slf6.x86_64.e9.debug/bin/NovaTimeConvert",
                               tmp_2], stdout=subprocess.PIPE)
    semfe_2 = proc_2.stdout.read()
    proc_2.wait()

    tmp_3 = str(AnalysisRegions['d2t_time_boundary']) + "UTC"
    proc_3 = subprocess.Popen(["/grid/fermiapp/products/nova/externals/novadaq/v04_05_00/slf6.x86_64.e9.debug/bin/NovaTimeConvert",
                               tmp_3], stdout=subprocess.PIPE)
    semfe_3 = proc_3.stdout.read()
    proc_3.wait()

    tmp_4 = str(AnalysisRegions['t2n_time_boundary']) + "UTC"
    proc_4 = subprocess.Popen(["/grid/fermiapp/products/nova/externals/novadaq/v04_05_00/slf6.x86_64.e9.debug/bin/NovaTimeConvert",
                               tmp_4], stdout=subprocess.PIPE)
    semfe_4 = proc_4.stdout.read()
    proc_4.wait()

    print semfe_1.split()[13], semfe_2.split()[13], semfe_3.split()[13], semfe_4.split()[13]

    s.Night_Twilight = long(semfe_1.split()[13])
    s.Twilight_Day = long(semfe_2.split()[13])
    s.Day_Twilight = long(semfe_3.split()[13])
    s.Twilight_Night = long(semfe_4.split()[13])
    ntp_tr.Fill()    
    
#    theTime = datetime.fromtimestamp(time.time())

    unixtimes = []
    del unixtimes[:]
    for key in ['n2t_time_boundary','t2d_time_boundary','d2t_time_boundary','t2n_time_boundary']:
        t = ephem.Date(AnalysisRegions[key]).datetime()
        # unixtime = time.mktime(theTime.timetuple())
        unixtime = 0
        try:
            unixtime = time.mktime(t.timetuple())
        except:
            print 'Value out of range'
        unixtimes.append(unixtime)
    # print '%20s, %20s, %20s, %20s' % (unixtimes[0], unixtimes[1], unixtimes[2], unixtimes[3] ) 
    sys.stdout.flush()

def moon_position_table_new3.FillAnal	(		a,
			region,
			time,
			resolution
	)

Definition at line 121 of file moon_position_table_new3.py.

Referenced by computeDayBoundaryMoon().

def FillAnal(a, region, time, resolution) :
    global last_region
    #
    # Added a resolution step to the counters
    stepSize = oneSecond*resolution
    if region == 1 :
        a['night_time']+= stepSize
        a['total_night_time']+= stepSize
    if region == 2 :
         a['day_time']+= stepSize
         a['total_day_time']+= stepSize
    if region == 3 :
         a['twilight_time']+= stepSize            
         a['total_twilight_time']+= stepSize

    if (region == 1) and (last_region == 3) :
        a['t2n_time_boundary'] = time
    if (region == 3) and (last_region == 1) :
        a['n2t_time_boundary'] = time
    if region == 2 and last_region == 3 :
        a['t2d_time_boundary'] = time
    if region == 3 and last_region == 2 :
        a['d2t_time_boundary'] = time
    last_region = region
    

def moon_position_table_new3.SetupAnalysisRegion

(

AnalysisRegions

)

Definition at line 80 of file moon_position_table_new3.py.

def SetupAnalysisRegion( AnalysisRegions ) :
    AnalysisRegions['day_time']              = 0
    AnalysisRegions['night_time']            = 0
    AnalysisRegions['twilight_time']         = 0
    AnalysisRegions['total_day_time']              = 0
    AnalysisRegions['total_night_time']            = 0
    AnalysisRegions['total_twilight_time']         = 0
    AnalysisRegions['n2t_time_boundary']     = 0
    AnalysisRegions['t2d_time_boundary']     = 0
    AnalysisRegions['d2t_time_boundary']     = 0
    AnalysisRegions['t2n_time_boundary']     = 0

def moon_position_table_new3.tl

(

date

)

Definition at line 146 of file moon_position_table_new3.py.

Referenced by genie::SKKinematicsGenerator.CalculateKin_AtharSingleKaon(), genie::DFRKinematicsGenerator.ComputeMaxXSec(), genie::SKKinematicsGenerator.ComputeMaxXSec(), DrawDet(), om::GUIMain.GUIMain(), genie::COHXSec.Integrate(), genie::DFRXSec.Integrate(), genie::KPhaseSpace.IsAllowed(), art::SourceHelper.makePtr(), ana.Nue2017FourBinDivisions(), ana.Nue2018ThreeBinDivisions(), Nue2018ThreeBinDivisions_Reduced(), ana.Nue2018ThreeBinDivisions_Reduced(), genie::DFRKinematicsGenerator.ProcessEventRecord(), and genie::KPhaseSpace.TLim().

def tl(date):
    return ephem.Date(date - 5 * ephem.hour)

Variable Documentation

dictionary moon_position_table_new3.AnalysisRegionsMoon = {}

Definition at line 94 of file moon_position_table_new3.py.

dictionary moon_position_table_new3.AnalysisRegionsSun = {}

Definition at line 93 of file moon_position_table_new3.py.

moon_position_table_new3.args

Definition at line 46 of file moon_position_table_new3.py.

moon_position_table_new3.central = timezone('US/Central')

Definition at line 52 of file moon_position_table_new3.py.

moon_position_table_new3.day = ephem.Date(start_day)

Loop over the days we care about.

Compute the positions

Definition at line 233 of file moon_position_table_new3.py.

moon_position_table_new3.default

Definition at line 39 of file moon_position_table_new3.py.

moon_position_table_new3.dest

Definition at line 39 of file moon_position_table_new3.py.

moon_position_table_new3.f = TFile('semfe_tmp.root','RECREATE')

Definition at line 14 of file moon_position_table_new3.py.

moon_position_table_new3.help

Definition at line 40 of file moon_position_table_new3.py.

int moon_position_table_new3.hour = 60

Definition at line 71 of file moon_position_table_new3.py.

int moon_position_table_new3.last_region = 999

Definition at line 100 of file moon_position_table_new3.py.

moon_position_table_new3.lat

Definition at line 65 of file moon_position_table_new3.py.

moon_position_table_new3.lon

Definition at line 66 of file moon_position_table_new3.py.

int moon_position_table_new3.minute = 60

Definition at line 69 of file moon_position_table_new3.py.

moon_position_table_new3.moon = ephem.Moon()

Definition at line 56 of file moon_position_table_new3.py.

moon_position_table_new3.nova = ephem.Observer()

Definition at line 58 of file moon_position_table_new3.py.

moon_position_table_new3.ntp_tr = TTree('ntp','My test tree')

Definition at line 15 of file moon_position_table_new3.py.

moon_position_table_new3.numdays = int(options.numdays)

Definition at line 78 of file moon_position_table_new3.py.

float moon_position_table_new3.oneMinute = 1.0

Definition at line 72 of file moon_position_table_new3.py.

float moon_position_table_new3.oneSecond = 1.0

Definition at line 73 of file moon_position_table_new3.py.

moon_position_table_new3.options

Definition at line 46 of file moon_position_table_new3.py.

moon_position_table_new3.parser = OptionParser(usage=usage,version=version)

Definition at line 37 of file moon_position_table_new3.py.

moon_position_table_new3.resolution = float(options.res)

Definition at line 48 of file moon_position_table_new3.py.

moon_position_table_new3.s = MyStruct()

Definition at line 29 of file moon_position_table_new3.py.

int moon_position_table_new3.second = 1

Definition at line 70 of file moon_position_table_new3.py.

Referenced by NdmcClient._adoptMetric(), NdmcClient._getMetric(), NdmcClient._updateGroupName(), AddExptCorrNumu(), fuzz::FuzzyKVertex.AddProng(), novaddt::DAQHeaderAna.analyze(), air::CosmicTrackRate.analyze(), comi::Leana.analyze(), bpf::BPFCVNAna.analyze(), upmuana::UpMuRecoAna.analyze(), upmuana::UpMuAnalysis.analyze(), genie::Registry.Append(), mcchk::DetSimAna.beginJob(), runh::BadDataFilter.beginSubRun(), CAF_makeCAFSensitivities_for_FNEX(), caf::Proxy< caf::SRHeader >.CheckEquals(), checkSwap(), cheat::CheckBackTracking.CheckTrackIdToParticle(), fuzz::FuzzyKMeanAlg.CheckUnique(), flat::Flat< caf::SRHeader >.Clear(), genie::Registry.CloneRegistryItem(), demoFitContours(), demoFitSlices(), NdmcClient.destroyMetrics(), ana::NumuCCIncPionTemplateFitter.doFullFit(), ana::nueccinc::NueCCIncTemplateFitter.doFullFit(), draw_track_in_one_view(), ana::SystMaker.DrawBkg(), DrawSurfacePoint(), comi::NearlineAna.endJob(), tbana::SpillStreamAna.EndSpill(), gov::fnal::cd::rms::provider.equal(), flat::Flat< caf::SRHeader >.Fill(), runh::BadDataFilter.filter(), util.find_median(), cmf::ContourFromLibrary.FindChiSqrForPoint(), genie::Algorithm.FindConfig(), fuzz::FuzzyKMeanAlg.FindUnclustered(), FourBinLabels(), genie::Algorithm.GetConfig(), getEffDiffErrPc(), ana::NumuCCIncPionTemplateFitter.getFitNormalizationAndErrors(), ana::nueccinc::NueCCIncTemplateFitter.getFitNormalizationAndErrors(), getFlatWeightedFittingSpectra(), ndreco::NDRecoFxs.getProngProximity(), getResolutionSpectra(), ndreco::NDRecoFxs.getTrackProximity(), cheat::BackTracker.HitsToParticle(), joint_fit_2017_contours(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), joint_fit_2017_slices(), joint_fit_2018_contours(), joint_fit_2018_slices(), jointsensitivity(), lastUpdated(), Legend(), cmf::Location.Location(), make_fc_mass_and_oct_nersc_2018(), make_fc_mass_and_oct_nersc_2019(), make_fc_mh_nersc_2018(), make_fc_mh_nersc_2019(), make_fc_oct_nersc_2018(), make_fc_oct_nersc_2019(), make_fc_slices_nersc_2018(), make_fc_slices_nersc_2018_stats(), make_fc_slices_nersc_2019(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), make_fc_surfaces_nersc_2019(), make_mass_and_oct_fc_input_2019(), MakeCAFSensitivities_for_FNEX(), dt::View.MakeChunkCombos(), MakeFakeExperiment(), evtsum::EventSummary.MakeOutput(), trk::WindowTrackingAlg.MakeViewTrack(), genie::supernova::Limits< double >.Max(), genie::Registry.Merge(), comi::Leana.MuonCandidate(), myBarChart(), caf::Proxy< caf::SRHeader >.operator=(), ana.PDG_cvnmax(), plot_ratios(), ana::PredictionScaleComp.PredictionScaleComp(), ana::SystMaker.ProcessSyst(), cvnprongtf::CVNProngTF.produce(), cvnneutronprongtf::CVNNeutronProngTF.produce(), comi::Cana.produce(), slid::LIDBuilder.produce(), cvntf::CVNProngEvaluatorTF.produce(), slid::SPIDBuilder.produce(), cheat::MCCheater.produce(), art::DataViewImpl.put(), ana::BayesianMarginal.QuantileThreshold(), filter::Filter.reconfigure(), dbutils::daqconfig::GeneralDataReader.remapValues(), calib::DriftCorrection.respondToCloseOutputFiles(), ana::PredictionCombinePeriods.SaveTo(), fuzz::FuzzyKMeanAlg.SeedA(), sensitivity2018(), sensitivity2020(), genie::supernova::Limits< double >.Set(), specprod_numuccinc(), SumSysts(), syst_table_fit(), syst_table_fit_new(), systematics_extrap_comp_from_pred_interp(), systematics_summary_from_pred_interp(), evd::HeaderDrawer.Text(), util.TimeSlice(), timing_peak(), ana::TreeMaker.TreeMaker(), genie::supernova::Limits< double >.UpdateMax(), ana::Nus17SystFromHist.WeightFor(), ana::NusSystFromHist.WeightFor(), ana::Nus18SystFromHist.WeightFor(), xsec_tot_uncert_optimization(), xsec_uncertainty_per_bin(), and NdmcClient.~NdmcClient().

moon_position_table_new3.start_day = options.startdate

Definition at line 77 of file moon_position_table_new3.py.

moon_position_table_new3.sun = ephem.Sun()

Definition at line 57 of file moon_position_table_new3.py.

int moon_position_table_new3.theta = 10

Definition at line 102 of file moon_position_table_new3.py.

string moon_position_table_new3.usage = "usage: %prog [options]"

Definition at line 36 of file moon_position_table_new3.py.

moon_position_table_new3.utc = pytz.utc

Definition at line 51 of file moon_position_table_new3.py.

string moon_position_table_new3.version = "1.0"

Definition at line 35 of file moon_position_table_new3.py.