Changeset 400 in Sophya


Ignore:
Timestamp:
Sep 2, 1999, 11:24:37 AM (26 years ago)
Author:
ansari
Message:

findstar DY qui sort des etoiles

Location:
trunk/Poubelle/archTOI.old
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/Poubelle/archTOI.old/Makefile

    r397 r400  
    11CXX=cxx
    22CFITSIODIR=../cfitsio
    3 #CPPFLAGS=-O -I$(CFITSIODIR)
    4 CPPFLAGS=-g -I$(CFITSIODIR)
     3CPPFLAGS=-O -I$(CFITSIODIR)
     4#CPPFLAGS=-g -I$(CFITSIODIR)
    55LDLIBS=-L$(CFITSIODIR) -lcfitsio
    66archtoi: archtoimain.o archtoi.o archeopsfile.o decompress.o toiiter.o toisvr.o archeops.o \

  • trunk/Poubelle/archTOI.old/archtoi.cc

    r394 r400  
    7373  else tsttoi(sstDiode)
    7474  else tsttoi(sstChannel)
     75  else tsttoi(sstStarCnt)
    7576  else tsttoi(sstStarZ)
    7677  else tsttoi(sstStarF)

  • trunk/Poubelle/archTOI.old/pisteetoile.cc

    r399 r400  
    116116//      if(!isAlive()) return;
    117117        indexDebutPiste=InDebutPiste;
    118         //      cout << "fill " << indexDebutPiste << "\n";
    119118        int* PhDiod=PhDiodArray;
    120119        int* pDataC=pData;
    … …  
    131130        for(int i=0; i<nbdata; i++, PPhDiod++,pDataC++) *PPhDiod=(*pDataC);
    132131        indexDebutPiste+=nbdata;
    133         //      cout << "push " << nbdata << " " << indexDebutPiste << "\n";
    134132        return;
    135133}
    … …  
    155153                if(testEtoile) {
    156154                // On emballe une etoile.
    157                   //              cout << "star " << setprecision(10) << TIndex << "\n";
    158155                        LastEtoile.FitForme=SuccesAFine;
    159156                        LastEtoile.X2Calc=X2Min;
    … …  
    275272//      On stocke les resultats
    276273                TIndex=indexDebutPiste+MaxIndex+tmaxfit;
    277                 //              cout << "tindex " << TIndex << " " <<indexDebutPiste<<" " <<Prepulselong << " " <<MaxIndex << " " << tmaxfit << "\n";
    278274                X2Min=0.;
    279275                SuccesAFine=false;

  • trunk/Poubelle/archTOI.old/ssthandler.cc

    r399 r400  
    184184        int jd0 = diodeHistLength-NbSampleBlock;
    185185
    186         //      cout << "find " << NoFirstSpInBlock << "\n";
    187186        bool FlagLBlockPerdu=false;
    188187        int ThisBlockNumb=numero_block(blk);
    … …  
    203202
    204203                if (FlagLBlockPerdu)  {
    205                         pPiste->fill(Diodedata+offseDataDiod,FirstSplNb(blk));
     204                        pPiste->fill(Diodedata+offseDataDiod,FirstSplNb(blk),PhDiodTabLong);
    206205                        offseDataDiod+=PhDiodTabLong;
    207206                }
    … …  
    232231                        }
    233232                // Est-on en bout de Piste?
    234                     int npouss = Pousslong;
    235                     if (offseDataDiod+npouss > NbSampleBlock) npouss = NbSampleBlock-offseDataDiod;
    236                     if (npouss<=0) break;
    237                         pPiste->push(Diodedata+offseDataDiod,npouss);
    238                         offseDataDiod+=npouss;
     233                    if (offseDataDiod+Pousslong > NbSampleBlock) break;
     234                    pPiste->push(Diodedata+offseDataDiod,Pousslong);
     235                    offseDataDiod+=Pousslong;
    239236                }
    240237        }

  • trunk/Poubelle/archTOI.old/toiiter.cc

    r394 r400  
    206206    if ((*i).triggering) {
    207207      switch ((*i).kind) {
    208         case boloTens:
    209         case boloRaw:
    210           trigMask |= block_bolo_mask;
    211           break;
    212         case gpsTime:
    213         case longitude:
    214         case latitude:
    215         case altitude:
    216           trigMask |= block_gps_mask;
    217           break;
    218         case azimut:
    219           file->needSSTProcessMask(SSTHandler::findPeriod);
    220           trigMask |= block_sst_mask;
    221           break;
    222         case sstStarZ:
    223         case sstStarF:
    224         case sstStarT:
    225           file->needSSTProcessMask(SSTHandler::findStars);
    226           trigMask |= block_sst_mask;
    227           break;
    228         case sstDiode:
    229           file->needSSTProcessMask(SSTHandler::permDiode);
    230           trigMask |= block_sst_mask;
    231           break;
    232         case sstChannel:
    233           trigMask |= block_sst_mask;
    234           break;
    235         case gyroRaw:
    236           trigMask |= block_gyro_mask;
    237           break;
    238         case alphaAxis:
    239         case deltaAxis:
    240         case alphaBolo:
    241         case deltaBolo:
    242           file->needSSTProcessMask(SSTHandler::findAxis);
    243           trigMask |= block_sst_mask;
    244           break;
     208      case boloTens:
     209      case boloRaw:
     210        trigMask |= block_bolo_mask;
     211        break;
     212      case gpsTime:
     213      case longitude:
     214      case latitude:
     215      case altitude:
     216        trigMask |= block_gps_mask;
     217        break;
     218      case azimut:
     219        file->needSSTProcessMask(SSTHandler::findPeriod);
     220        trigMask |= block_sst_mask;
     221        break;
     222      case sstStarCnt:
     223      case sstStarZ:
     224      case sstStarF:
     225      case sstStarT:
     226        file->needSSTProcessMask(SSTHandler::findStars);
     227        trigMask |= block_sst_mask;
     228        break;
     229      case sstDiode:
     230        file->needSSTProcessMask(SSTHandler::permDiode);
     231        trigMask |= block_sst_mask;
     232        break;
     233      case sstChannel:
     234        trigMask |= block_sst_mask;
     235        break;
     236      case gyroRaw:
     237        trigMask |= block_gyro_mask;
     238        break;
     239      case alphaAxis:
     240      case deltaAxis:
     241      case alphaBolo:
     242      case deltaBolo:
     243        file->needSSTProcessMask(SSTHandler::findAxis);
     244        trigMask |= block_sst_mask;
     245        break;
    245246      }
    246247    }
    … …  
    404405   int index = infos[column].index;
    405406   switch (kind) {
    406      case sampleNum:
    407      case internalTime:
    408      case mjd:
    409        return true;
    410      case boloTens:
    411        if (imes==0 && file->llastBolo()==NULL) return false;
    412        return file->lastBolo() != NULL;
    413      case boloRaw:
    414        return file->lastBolo() != NULL;
    415      case sstDiode:
    416      case sstChannel:
    417        return file->lastSST() != NULL;
    418      case sstStarZ:
    419      case sstStarF:
    420      case sstStarT:{
    421        if (file->lastSST() == NULL) return false;
    422        int n = file->getNumbStar(imes);
    423        return (n > 0 && imes < n);
    424        }
    425      case gyroRaw:
    426        return (file->lastGyro() != NULL);
    427      case gpsTime:
    428        return file->hasGPSTime();
    429      case longitude:
    430      case latitude:
    431        return file->hasGPSPos();
    432      case altitude:
    433        return file->hasGPSAlt();
    434      case tsid:
    435        return file->hasGPSPos();
    436      case azimut:
    437      case alphaAxis:
    438      case deltaAxis:
    439        return false;
    440        //return (file->lastGPS() != NULL && file->lastSST() != NULL);
    441      case alphaSst:
    442      case deltaSst:
    443      case alphaBolo:
    444      case deltaBolo:
    445        return false;
    446        
    447      case boloTemp:
    448        return false;
     407   case sampleNum:
     408   case internalTime:
     409   case mjd:
     410     return true;
     411   case boloTens:
     412     if (imes==0 && file->llastBolo()==NULL) return false;
     413     return file->lastBolo() != NULL;
     414   case boloRaw:
     415     return file->lastBolo() != NULL;
     416   case sstDiode:
     417   case sstChannel:
     418     return file->lastSST() != NULL;
     419   case sstStarCnt:
     420   case sstStarZ:
     421   case sstStarF:
     422   case sstStarT:{
     423     if (file->lastSST() == NULL) return false;
     424     int n = file->getNumbStar(imes);
     425     return (n > 0 && index < n);
     426   }
     427   case gyroRaw:
     428     return (file->lastGyro() != NULL);
     429   case gpsTime:
     430     return file->hasGPSTime();
     431   case longitude:
     432   case latitude:
     433     return file->hasGPSPos();
     434   case altitude:
     435     return file->hasGPSAlt();
     436   case tsid:
     437     return file->hasGPSPos();
     438   case azimut:
     439   case alphaAxis:
     440   case deltaAxis:
     441     return false;
     442     //return (file->lastGPS() != NULL && file->lastSST() != NULL);
     443   case alphaSst:
     444   case deltaSst:
     445   case alphaBolo:
     446   case deltaBolo:
     447     return false;
     448     
     449   case boloTemp:
     450     return false;
    449451   }
    450452   return false;
    451453}
    452  
     454
    453455double TOIIter::getValue(int column) {
    454456  if (!initDone) Init();
    … …  
    475477   int index = infos[column].index;
    476478   switch (kind) {
    477      case sampleNum:
    478        return getSampleIndex();
    479      case internalTime:
    480        return getSampleIndex() * archParam.acq.perEch;
    481      case mjd:
    482        return getMJD();
    483      case boloTens:
    484        return file->getMuVBolo(index, imes);
    485      case boloRaw:
    486        return file->getRawBolo(index, imes);
    487      case sstDiode:
    488        return file->getSSTSignal(index, imes);
    489      case sstChannel:
    490        return file->getSSTRawSignal(index, imes);
    491      case sstStarZ:
    492        return file->getSSTStarZ(index, imes);
    493      case sstStarF:
    494        return file->getSSTStarF(index, imes);
    495      case sstStarT:
    496        return file->getSSTStarT(index, imes);
    497      case gyroRaw:
    498        return file->getGyro(index, imes);
    499      case gpsTime:
    500        return file->getGPSUTC();
    501      case longitude:
    502        return file->getGPSLong();
    503      case latitude:
    504        return file->getGPSLat();
    505      case altitude:
    506        return file->getGPSAlt();
    507      case tsid:
    508        tSid.setLongitude(file->getGPSLong());
    509        return tSid.getLST(getMJD());
    510      case azimut:
    511        return file->getAzimut(imes);
    512      case alphaAxis:
    513        return file->getAlpha(imes);
    514      case deltaAxis:
    515        return file->getDelta(imes);
     479   case sampleNum:
     480     return getSampleIndex();
     481   case internalTime:
     482     return getSampleIndex() * archParam.acq.perEch;
     483   case mjd:
     484     return getMJD();
     485   case boloTens:
     486     return file->getMuVBolo(index, imes);
     487   case boloRaw:
     488     return file->getRawBolo(index, imes);
     489   case sstDiode:
     490     return file->getSSTSignal(index, imes);
     491   case sstChannel:
     492     return file->getSSTRawSignal(index, imes);
     493   case sstStarCnt:
     494     return file->getNumbStar(imes);
     495   case sstStarZ:
     496     return file->getSSTStarZ(index, imes);
     497   case sstStarF:
     498     return file->getSSTStarF(index, imes);
     499   case sstStarT:
     500     return file->getSSTStarT(index, imes);
     501   case gyroRaw:
     502     return file->getGyro(index, imes);
     503   case gpsTime:
     504     return file->getGPSUTC();
     505   case longitude:
     506     return file->getGPSLong();
     507   case latitude:
     508     return file->getGPSLat();
     509   case altitude:
     510     return file->getGPSAlt();
     511   case tsid:
     512     tSid.setLongitude(file->getGPSLong());
     513     return tSid.getLST(getMJD());
     514   case azimut:
     515     return file->getAzimut(imes);
     516   case alphaAxis:
     517     return file->getAlpha(imes);
     518   case deltaAxis:
     519     return file->getDelta(imes);
    516520   }
    517521   return -1;
    … …  
    541545     case sstChannel:
    542546     case sstDiode:
     547     case sstStarCnt:
    543548     case sstStarZ:
    544549     case sstStarF:

  • trunk/Poubelle/archTOI.old/toiiter.h

    r394 r400  
    2323  sstDiode,   // signal brut SST, indice = diode, 0-45
    2424  sstChannel, // signal brut SST, indice = canal, 0-47
     25  sstStarCnt, // nombre d'etoiles detectees durant ce sample
    2526  sstStarZ,   // numero de diode d''une eventuelle etoile,   i=etoile 0..n
    2627  sstStarF,   // flux d''une eventuelle etoile,             <0 si moins de i+1 etoiles
Note: See TracChangeset for help on using the changeset viewer.