iceCube.uhe.analysis
public class I3ParticleAnalysisFactory extends java.lang.Object
Modifier and Type | Field and Description |
---|---|
static int[][] |
badRunID
Bad Run Numbers..
|
protected boolean |
filterOutBadRunData
Flag for filtering out the bad run data
|
protected double[][] |
histCosZenithFgQuality |
protected double[][] |
histLogEnergyCosZenith |
protected double[][] |
histLogNpeCosZenith |
protected double[][] |
histLogNpeFgQuality |
protected double[][] |
histLogNpeLogEnergy |
protected double[][] |
histLogNpeNDOM |
protected static double |
ln10 |
protected static double |
maxLogEnergy |
protected static double |
maxLogNpe |
protected double |
mcArea
Area where juliet particles are thrown in the MC simulation .
|
protected double |
mcCorsikaArea
Area of the Corsika MMC particles are thrown
|
protected double |
mcCorsikaOmega
Solid angle of the Corsika mmc particles
|
protected double |
mcOmega
Solid angle where juliet particles are thrown in
the MC simulation
|
protected static double |
minLogEnergy |
protected static double |
minLogNpe |
protected static double |
minLogNPEToAnalize
The lowest level Requirement of NPEs.
|
protected double |
minLogNumberOfEvents |
protected static int |
minNDOMsToAnalize
The lowest level Requirement of NDOMs.
|
protected AtmMuonBundleFlux |
muonBundle
AtmMuonBundleFlux object for calculating the primary cosmic ray energy
|
protected int |
numberOfFilledI3Particles |
protected double |
observationTime
Observation Time [sec]
|
protected static double |
xCenterOfIC9 |
protected static double |
yCenterOfIC9 |
Constructor and Description |
---|
I3ParticleAnalysisFactory(java.io.InputStream in)
Default Constructor
|
I3ParticleAnalysisFactory(java.io.InputStream in,
boolean filterOutBadRunData)
Constructor
|
Modifier and Type | Method and Description |
---|---|
void |
autoScale(boolean flag) |
static IComparisonResult |
calcChi2(IHistogram1D mcHistogram,
IHistogram1D dataHistogram)
Chi^2 comparison - analysis utility method.
|
static IComparisonResult |
calcChi2(IHistogram1D mcHistogram,
IHistogram1D dataHistogram,
double xFirst,
double xLast)
Chi^2 comparison - analysis utility method.
|
static IComparisonResult |
calcChi2(IHistogram1D mcHistogram,
IHistogram1D dataHistogram,
int binFirst,
int binLast)
Chi^2 comparison - analysis utility method
|
void |
changeAOmega()
Change the MC area and solid angle to the Corska configuration
|
double |
drawEnergyDistributionOnXfig()
Draw Energy distribution with slices of cos(zenith) with its bin width
|
void |
drawEventsWithAtmMuonWeights(java.lang.String modelName)
draw atmospheric muon weighted-events
|
void |
drawEventsWithGZKWeights(java.lang.String modelName)
draw GZK flux weighted-events
|
void |
drawEventsWithNoWeights()
draw unweighted-events
|
double |
drawFirstGuessQualityDistributionOnXfig()
Draw First Guess Quality distribution with slices of Npe
|
double |
drawNpeDistributionOnXfig()
Draw Npe distribution with slices of cos(zenith) with its bin width
|
double |
drawNpeDistributionWithNDOMsliceOnXfig()
Draw NPE distribution with slices of NDOM with its bin width
|
double |
drawNpeDistributionWithSliceOfEnergyOnXfig()
Draw Npe distribution with slices of Log(Energy) with its bin width
|
double |
drawZenithAngleDistributionOnXfig()
Draw cos(zenith) distribution with slices of Npe with its bin width
|
double |
drawZenithAngleDistributionWithFGsliceOnXfig()
Draw cos(zenith) distribution with slices of
First Guess quality with its bin width
|
protected static double |
getDNDLogE(double powerLaw,
double energy,
double logEnergyMinimum,
double logEnergyMaximum)
Calculate the MC Primary Spectrum dN/dLogE
except normalization constatnt
determined by the total numner of MC events.
|
protected I3Particle |
getI3ParticleLite(I3Particle iceParticle)
make I3Particle as a temporal dataclass for the analysis.
|
java.util.ListIterator |
getParticleIterator()
Return the Particle Iterator loping over I3Particles
|
protected boolean |
isBadRunData(I3Particle iceParticle)
Judge if this event has to be excluded because of
the bad run.
|
void |
makeHistogram()
Making historgram for drawing by "grafig" graphics package on Xfig.
|
protected IHistogram1D |
makeJaida1DHistogram(java.lang.String option,
IHistogramFactory jaidaHistoFactory)
Histogram name is assinged by the default - same as "option"
|
protected IHistogram1D |
makeJaida1DHistogram(java.lang.String option,
java.lang.String histName,
boolean bootstrap,
IHistogramFactory jaidaHistoFactory)
Make and Fill IHistogram1D by reading out the variables of I3Particles
|
protected IHistogram2D |
makeJaida2DHistogram(java.lang.String option,
IHistogramFactory jaidaHistoFactory)
Histogram name is assinged by the default - same as "option"
|
protected IHistogram2D |
makeJaida2DHistogram(java.lang.String option,
java.lang.String histName,
boolean bootstrap,
IHistogramFactory jaidaHistoFactory)
Make and Fill IHistogram2D by reading out the variables of I3Particles
|
void |
plotByJointLine()
Plotting histogram by jointed-line
|
void |
plotByPointsWithErrorBars()
Plotting histogram by points with error bars
|
protected void |
readI3Particles(java.io.InputStream in)
reading the I3Particle objects
|
void |
setBinSize(double deltaLogNpe,
double deltaCosZenith) |
void |
setBinSize(double deltaLogNpe,
double deltaCosZenith,
double deltaFgQuality) |
void |
setBinSize(double deltaLogEnergy,
double deltaLogNpe,
double deltaCosZenith,
double deltaFgQuality)
Set historgram bin size
|
void |
setCriteria(Criteria criteria)
Sets the criteria on making Histogram
|
void |
setMaxRangeOfNDOMsInDrawing(int nDOMs)
Set Max Number Of DOMs in Drawing
|
void |
setMCArea(double area)
Set the MC area [cm^2]
|
void |
setMCSolidAngle(double omega)
Set the MC solid angle [sr]
|
void |
setObservationTime(double time)
Set the MC solid angle [sec]
|
void |
switchToMCTruth()
Switch to parameters concerned with MCtrue.
|
void |
switchToReco()
Switch to parameters concerned with Reco results.
|
protected static final double ln10
protected static double xCenterOfIC9
protected static double yCenterOfIC9
protected double minLogNumberOfEvents
protected static double maxLogNpe
protected static double minLogNpe
protected static double maxLogEnergy
protected static double minLogEnergy
protected static int minNDOMsToAnalize
protected static double minLogNPEToAnalize
protected int numberOfFilledI3Particles
protected double mcArea
protected double mcOmega
protected double mcCorsikaArea
protected double mcCorsikaOmega
protected double observationTime
protected AtmMuonBundleFlux muonBundle
public static final int[][] badRunID
protected boolean filterOutBadRunData
protected double[][] histLogNpeCosZenith
protected double[][] histLogNpeNDOM
protected double[][] histLogNpeFgQuality
protected double[][] histCosZenithFgQuality
protected double[][] histLogNpeLogEnergy
protected double[][] histLogEnergyCosZenith
public I3ParticleAnalysisFactory(java.io.InputStream in) throws java.io.IOException
InputStream in : Stream to readout a series of I3Particles that is subject to your analysisIt reads all I3Particles from InputStream and hold them in form of List.
java.io.IOException
public I3ParticleAnalysisFactory(java.io.InputStream in, boolean filterOutBadRunData) throws java.io.IOException
InputStream in : Stream to readout a series of I3Particles that is subject to your analysis filterOutBadRunData : do not use events in the bad run. The real data analysis only.It reads all I3Particles from InputStream and hold them in form of List.
java.io.IOException
protected void readI3Particles(java.io.InputStream in) throws java.io.IOException
java.io.IOException
protected I3Particle getI3ParticleLite(I3Particle iceParticle)
public java.util.ListIterator getParticleIterator()
public void setMCArea(double area)
public void setMCSolidAngle(double omega)
public void changeAOmega()
public void setObservationTime(double time)
protected boolean isBadRunData(I3Particle iceParticle)
public void autoScale(boolean flag)
protected static double getDNDLogE(double powerLaw, double energy, double logEnergyMinimum, double logEnergyMaximum)
double powerLaw : dN/dE = E**(-powerLaw) double energy : Energy of this I3Particle [GeV] double logEnergyMinimum : log(Energy Minimum [GeV]) in the spectral range double logEnergyMaximum : log(Energy Maximum [GeV]) in the spectral range
public void setBinSize(double deltaLogEnergy, double deltaLogNpe, double deltaCosZenith, double deltaFgQuality)
double deltaLogEnergy bin size of logEnergy double deltaLogNpe bin size of logNpe double deltaCosZenith bin size of cos(zenith angle) double deltaFgQuality bin size of "First Guess Quality index" in unit of beta
public void setBinSize(double deltaLogNpe, double deltaCosZenith, double deltaFgQuality)
public void setBinSize(double deltaLogNpe, double deltaCosZenith)
public void setMaxRangeOfNDOMsInDrawing(int nDOMs)
public void switchToMCTruth()
public void switchToReco()
public void plotByJointLine()
public void plotByPointsWithErrorBars()
public void drawEventsWithAtmMuonWeights(java.lang.String modelName)
public void drawEventsWithGZKWeights(java.lang.String modelName)
public void drawEventsWithNoWeights()
public void setCriteria(Criteria criteria)
public void makeHistogram()
public double drawNpeDistributionOnXfig()
public double drawNpeDistributionWithSliceOfEnergyOnXfig()
public double drawEnergyDistributionOnXfig()
public double drawZenithAngleDistributionOnXfig()
public double drawZenithAngleDistributionWithFGsliceOnXfig()
public double drawNpeDistributionWithNDOMsliceOnXfig()
public double drawFirstGuessQualityDistributionOnXfig()
protected IHistogram1D makeJaida1DHistogram(java.lang.String option, java.lang.String histName, boolean bootstrap, IHistogramFactory jaidaHistoFactory)
String option "logEnergy" : plot logE "logRecoEnergy" : plot logRecoEnergy "logNpe" : plot logNpe "cosZenith" : plot cos(Zenith Angle) "nDOMs" : plot numberOfDOMs "firstGuess" : plot beta - first guess quality boolean bootstrap false : default. No bootstraping true : resampling the (weighted) I3Particles following the poisson distribution.
protected IHistogram1D makeJaida1DHistogram(java.lang.String option, IHistogramFactory jaidaHistoFactory)
protected IHistogram2D makeJaida2DHistogram(java.lang.String option, java.lang.String histName, boolean bootstrap, IHistogramFactory jaidaHistoFactory)
String option "logE-Npe" : plot logE-logNpe "logECR-Npe" : plot log(Primary Cosmic Ray E)-logNpe "logECR-logE" : plot log(Primary Cosmic Ray E)-logE "logRecoE-Npe" : plot logRecoE-logNpe "logE-cosZenith" : plot logE-cosZenith "logNpe-cosZenith": plot logNpe-cosZenith "logNpe-CobZ" : plot logNpe-ConterOfBrightnessZ "CobR-CobZ" : plot CenterOfBrightness(sqrt(x^2+y^2))-ConterOfBrightnessZ boolean bootstrap false : default. No bootsraping true : resampling the (weighted) I3Particles following the poisson distribution.
protected IHistogram2D makeJaida2DHistogram(java.lang.String option, IHistogramFactory jaidaHistoFactory)
public static IComparisonResult calcChi2(IHistogram1D mcHistogram, IHistogram1D dataHistogram, int binFirst, int binLast)
public static IComparisonResult calcChi2(IHistogram1D mcHistogram, IHistogram1D dataHistogram, double xFirst, double xLast)
public static IComparisonResult calcChi2(IHistogram1D mcHistogram, IHistogram1D dataHistogram)