Package gaiasimu.gaia

Examples of gaiasimu.gaia.Gaia

242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
      gaiasimu.universe.source.AstroSource as = sourceList.get(i);
      Random random = new Random();
      random.setSeed(i);
     
      // gogsprocessSource(as, random);
      Gaia gaia = GogSimulation.getGaia();
      double sourceMagnitude = as.getPhotometry().getGaiaMag(GaiaPhotometrySystem.photometryG,
          GogConstants.CATALOGUE_TIME, gc.isPhotometryFromSpectrum());
     
      // AstroSource as= (AstroSource)aSourcesList.get(i).getComponents()[j];
      gaia.cu1.tools.numeric.algebra.GVector3d positionICRS = as.getAstrometry().getICRSPosition(GogConstants.CATALOGUE_TIME);
      // compute all possible transits for this source, using ScanningLaw.inverseScan
      double halfAperture = GaiaParam.Satellite.FOV_AC * GaiaParam.Nature.DEGREE_RADIAN / 2;
      // This method returns the list of all the transits
      double t1 = System.currentTimeMillis();
      List<gaiasimu.gaia.spacecraft.Transit> fovTransits = gaia.attitude.inverseScan(positionICRS, halfAperture, gaia
          .getTelescope(FoV.FOV1));
      List<gaiasimu.gaia.spacecraft.Transit> fovTransits2 = gaia.attitude.inverseScan(positionICRS, halfAperture, gaia
          .getTelescope(FoV.FOV2));
      // System.out.println("HA:" + halfAperture + "HA2:" +halfAperture2 + "ZA:" + zAngle + "ZA2:" + zAngle2);
      fovTransits.addAll(fovTransits2);
      // keep here only the number of transits that have been requested
      if ((gc.getTransitNumber() > 0) && (gc.getTransitNumber() < fovTransits.size())) {
View Full Code Here
25
26
27
28
29
30
31
32
33
34
35
    gaiaInstances = new ArrayList<Gaia>();
   
    // Instantiate Gaia as often as specified by parameter
    for(int i = 0; i < NGaias; i++) {
      try {
        gaiaInstances.add(new Gaia());
      } catch (SimuException e) {
        logger.error("Failed to instantiate Gaia!");
        e.printStackTrace();
      }
    }
View Full Code Here
100
101
102
103
104
105
106
107
108
109
110
111
112
    GaiaSimuTime time = tr.getGSTime();
    GVector3d comrs = star.getAstrometry().getCoMRSPosition(time, relativity);
   
    // Calculate local scan coordinates
    double scanAngle;
    Gaia gaia = sim.getGaia();
    synchronized(gaia) {
      scanAngle = gaia.getAttitude().getScanAngle(comrs, time);
    }
   
    StellarAstrometry astro = (refAstro != null)? refAstro: astrometry;
    GVector2d lsc = Conversion.fromCoMRStoLSC(comrs, astro.getAlpha(), astro.getDelta(), scanAngle);
   
View Full Code Here
186
187
188
189
190
191
192
193
194
195
196
197
198
199
      double halfAperture = GaiaParam.Satellite.FOV_AC * GaiaParam.Nature.DEGREE_RADIAN/2;
     
      // Get transits for both fields of view
      ArrayList<Transit> fov1Transits, fov2Transits;
     
      Gaia gaia = sim.getGaia();
      synchronized(gaia) {
        fov1Transits = gaia.attitude.inverseScan(vICRS, halfAperture, gaia.getTelescope(FoV.FOV1) );
        fov2Transits = gaia.attitude.inverseScan(vICRS, halfAperture, gaia.getTelescope(FoV.FOV2) );
      }
     
      // Merge in one array
      transits.addAll(fov1Transits);
      transits.addAll(fov2Transits);
View Full Code Here
95
96
97
98
99
100
101
102
103
104
105
106
107
        fieldAngles.add(randomOffset);
      }
     
      // Calculate scanning angle
      double scanAngle;
      Gaia gaia = sim.getGaia();
      synchronized(gaia) {
        scanAngle = gaia.getAttitude().getScanAngle(tr.getDirectionCoMRS(), tr.getGSTime());
      }
     
      // Add new observed transit to array list
      obsData.addTransit(new TransitData(fieldAngles, errors, false, tr.getGSTime(), scanAngle, useAc));
     
View Full Code Here
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
      gaiasimu.universe.source.AstroSource as = sourceList.get(i);
      Random random = new Random();
      random.setSeed(i);
     
      // gogsprocessSource(as, random);
      Gaia gaia = GogSimulation.getGaia();
     
      // AstroSource as= (AstroSource)aSourcesList.get(i).getComponents()[j];
      gaia.cu1.tools.numeric.algebra.GVector3d positionICRS = as.getAstrometry().getICRSPosition(GogConstants.CATALOGUE_TIME);
     
      // compute all possible transits for this source, using ScanningLaw.inverseScan
      double halfAperture = GaiaParam.Satellite.FOV_AC * GaiaParam.Nature.DEGREE_RADIAN / 2;
     
      // This method returns the list of all the transits
      List<gaiasimu.gaia.spacecraft.Transit> fovTransits = gaia.attitude.inverseScan(positionICRS, halfAperture, gaia
          .getTelescope(FoV.FOV1));
      List<gaiasimu.gaia.spacecraft.Transit> fovTransits2 = gaia.attitude.inverseScan(positionICRS, halfAperture, gaia
          .getTelescope(FoV.FOV2));
     
      // System.out.println("HA:" + halfAperture + "HA2:" +halfAperture2 + "ZA:" + zAngle + "ZA2:" + zAngle2);
      fovTransits.addAll(fovTransits2);
     
View Full Code Here
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
   
    double halfAperture = GaiaParam.Satellite.FOV_AC*GaiaParam.Nature.DEGREE_RADIAN/2;
   
    // Create and initialize a Gaia instance with default values

    Gaia gaia = new Gaia();
   
    // This method returns the list of all the transits
   
    ArrayList<Transit> fov1Transits = gaia.attitude.inverseScan(vICRS, halfAperture, gaia.getTelescope(FoV.FOV1));
    ArrayList<Transit> fov2Transits = gaia.attitude.inverseScan(vICRS, halfAperture, gaia.getTelescope(FoV.FOV2));
    System.out.println("GJ 876 will transit "+fov1Transits.size()+" times in FOV1 during the mission.");
    System.out.println("GJ 876 will transit "+fov2Transits.size()+" times in FOV2 during the mission.");
   
    // Merge arrays
   
    ArrayList<Transit> allTransits = new ArrayList<Transit>();
    allTransits.addAll(fov1Transits);
    allTransits.addAll(fov2Transits);
   
    // Print astrometry information
   
    System.out.println("\n-- astrometric data -- ");
   
    for ( Transit tr : allTransits ) {
     
      final boolean relativity = false; // whether to include relativity
      final double lambda = 600.0; // wavelength in nm

      GaiaSimuTime time = tr.getGSTime();
      Astrometry ast = star.getAstrometry();
      double trAlpha = ast.getCoMRSAstrometricParam(time, relativity).getAlpha();
      double trDelta = ast.getCoMRSAstrometricParam(time, relativity).getDelta();
     
      GVector3d comrs = ast.getCoMRSPosition(time, relativity);
      System.out.println(comrs);
      System.out.println(time);
      System.out.println(tr.getFov());
      GVector3d fprs = gaia.fromCoMRStoFPRS(comrs, time, tr.getFov(), GaiaSimuEnums.FpField.AF, lambda);
      GVector3d fieldAngles = gaia.fromFPRStoFieldAngles(fprs, tr.getFov());
     
      System.out.printf("%.15e,\t%.15e,\t%.15e,\t%.15e,\t%.15e\n", time.getJD(),
          fieldAngles.getX(), fieldAngles.getY(), trAlpha, trDelta);
     
    }
View Full Code Here
TOP

Related Classes of gaiasimu.gaia.Gaia

  • newExamples.archive.GenerateSimulateOnePlanet
  • newExamples.archive.GenerateSimulateOnePlanet2
  • newExamples.archive.GSInvScanLaw
  • newExamples.DataGenerator
  • newExamples.GaiaWrapper
  • newExamples.Model
    • 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.