Package gaiasimu.universe.source.stellar

Examples of gaiasimu.universe.source.stellar.ComponentPhotometry

276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
    // Initialize ComponentAstrometry pointing to the primary and then correct with setCompanion()
    ComponentAstrometry compAstrom = new ComponentAstrometry(planets.size(), orbParam);
    ExoPlanetPhysicalParameters compPhys =
      new ExoPlanetPhysicalParameters(tEq, plRadius, vSini, massM2, MBol2, star.getStarPhysicalParams());
   
    ComponentPhotometry compPhot =
      new ComponentPhotometry(compAstrom, magMV2, distance, vmi2, star.getPhotometry().getAbsV(), star.getPhotometry().getExtinctionCurve());
   
    ExoPlanet planet = new ExoPlanet(star.getRootId(), compPhot, compAstrom, compPhys);
    compPhot.setTargetAstroSource(planet);
    compAstrom.setTargetAstroSource(planet);
    compPhys.setTargetAstroSource(planet);
   
    return this.addPlanet(planet);
   
View Full Code Here
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
    // Initialize ComponentAstrometry pointing to the primary and then correct with setCompanion()
    ComponentAstrometry compAstrom = new ComponentAstrometry(index, orbParam);
    ExoPlanetPhysicalParameters compPhys =
      new ExoPlanetPhysicalParameters(tEq, plRadius, vSini, massM2, MBol2, star.getStarPhysicalParams());
   
    ComponentPhotometry compPhot =
      new ComponentPhotometry(compAstrom, magMV2, distance, vmi2, star.getPhotometry().getAbsV(), star.getPhotometry().getExtinctionCurve());
   
    ExoPlanet planet = new ExoPlanet(star.getRootId(), compPhot, compAstrom, compPhys);
    compPhot.setTargetAstroSource(planet);
    compAstrom.setTargetAstroSource(planet);
    compPhys.setTargetAstroSource(planet);
   
    return planet;
   
View Full Code Here
102
103
104
105
106
107
108
109
110
111
112
    }
    StarPhysicalParameters starPhysicalParameters = new StarPhysicalParameters(mass, spectralType, magBol, pop, feH, alphaE);
   
    // Generate empty astrometry
    ComponentAstrometry primaryAstrometry = new ComponentAstrometry();
    ComponentPhotometry primaryPhotometry = new ComponentPhotometry(primaryAstrometry, magMv, 1.0e3 / astrometry.getVarpi(), vMinusI, absV);
   
    Star star = new Star(idRoot, primaryPhotometry, primaryAstrometry, starPhysicalParameters);
   
    return star;
   
View Full Code Here
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
    // init primary star

    long primaryID = star.getRootId();
    ComponentAstrometry primaryAstrometry = new ComponentAstrometry();

    ComponentPhotometry primaryPhotometry = new ComponentPhotometry(primaryAstrometry, magMv,
        distance, vMinusI, absV, star.getPhotometry().getExtinctionCurve());
    StarPhysicalParameters origPhys = star.getStarPhysicalParams();

    // TODO: some physical parameters of the primary may not be coherent anymore...
    StarPhysicalParameters primaryPhys = new StarPhysicalParameters(
        origPhys.getTeff(), origPhys.getFeH(), origPhys.getlogg(), origPhys.getAlphaElements(),
        origPhys.getMass(), origPhys.getSpectralType(), 10.0, origPhys.getMeanRadius(),
        origPhys.getStellarPop(), origPhys.getvsini());

    star = new Star(primaryID, primaryPhotometry, primaryAstrometry, primaryPhys);
    components[0] = star;

    //
    // init Planet
    //
    int compnb = 1;

    double tEqu2 = origPhys.getTeff() * Math.sqrt(origPhys.getMeanRadius() / (2 * orbParam.semiMajorAxis));
    tEqu2 *= Math.pow((1 - bondAlbedo), 1./4);
   
    // Calculating radius -- not necessary for astrometry, could be replaced by simpler model
    final double a[] = {1.2158123, -0.15080468, 0.69730554, -0.017303049, 0.00026388263, -1.2822719e-6,
      3.3135757e-9, -4.8366013e-12, 4.0017118e-15, -1.7477892e-18, 3.1347937e-22};
    double plRadius = a[0] + a[1] * Math.pow(mass, a[2]);
    for (int i = 3; i < a.length; i++) {
      plRadius += a[i] * Math.pow(mass, a[i - 2]);
    }
   
    plRadius *= GaiaParam.Nature.JUPITER_POLARRADIUS/GaiaParam.Nature.SUN_EQUATORIALRADIUS;
   
    // TODO put a correct value for absolute visual and bolometric magnitude, vsini, V-I
    double magMV2 = 32.58, MBol2 = magMV2-4; //Mag MV from Baraffe models.
    double vSini2 = 0, vmi2 = 9.9999;

    // Initialize ComponentAstrometry pointing to the primary and then correct with setCompanion()
    ComponentAstrometry compAstrom = new ComponentAstrometry(compnb, orbParam);
    ExoPlanetPhysicalParameters compPhys = new ExoPlanetPhysicalParameters(
        tEqu2, plRadius, vSini2, massM2
        / GaiaParam.Nature.SUN_MASS* GaiaParam.Nature.JUPITER_MASS, MBol2, primaryPhys) ;
   
    ComponentPhotometry compPhot = new ComponentPhotometry(compAstrom, magMV2,
        distance, vmi2, absV, star.getPhotometry().getExtinctionCurve());

    StellarSource comp = new ExoPlanet(primaryID, compPhot, compAstrom, compPhys);
    compPhot.setTargetAstroSource(comp);
    compAstrom.setTargetAstroSource(comp);
    compPhys.setTargetAstroSource(comp);

    components[compnb] = comp;

View Full Code Here
TOP

Related Classes of gaiasimu.universe.source.stellar.ComponentPhotometry

  • newExamples.archive.GSInvScanLaw
  • newExamples.Model
  • newExamples.PlanetTemplate
  • newExamples.StarTemplate
    • Copyright © 2018 www.massapicom. All rights reserved.
    All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.