Changeset 416 in Sophya for trunk/Poubelle/archTOI.old
- Timestamp:
- Sep 23, 1999, 11:08:43 AM (26 years ago)
- Location:
- trunk/Poubelle/archTOI.old
- Files:
-
- 1 added
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Poubelle/archTOI.old/archeopsfile.cc
r415 r416 886 886 } 887 887 888 long ArcheopsFile::getRawBoloCN(int ibolo, int imesure) { 889 // Si pas bloc comprime -> 0 890 if (lastBoloComp() == NULL) return 0; 891 if (numero_block(lastBoloComp()) != numero_block(lastBolo())) return 0; 892 block_type_reglage* reglage = lastReglage(); 893 if (!reglage) return 0; 894 895 int nb_coups= reglage->reglage.horloge.nb_mesures/2 - reglage->reglage.horloge.temp_mort; 896 unsigned int4* data = lastBoloComp()->data_bolo[ibolo]; 897 // Les deux premieres valeurs sont codees directement... 898 if (imesure<2) return 0; 899 int iExp = (imesure-2)/4 + 1; 900 int expo = data[iExp] & 0xf; 901 int noise = 1 << expo; 902 return (noise << 1)/nb_coups; 903 } 904 905 906 888 907 def_gains 889 908 … … 899 918 //return ((1e5*y)/(65536.*gain_ampli(reglage->reglage.bolo[ibolo]))); 900 919 } 920 921 double ArcheopsFile::getMuVBoloCN(int ibolo, int imesure) { // microvolts, sur valeur unique 922 double y = getRawBoloCN(ibolo, imesure); 923 block_type_reglage* reglage = lastReglage(); 924 block_type_param* param = lastParam(); 925 if (!reglage) return 0; 926 if (!param) return 0; 927 return bol_micro_volt(y,(double)param->param.bolo[ibolo].bolo_gain*gain_ampli(reglage->reglage.bolo[ibolo])); 928 } 929 901 930 902 931 double ArcheopsFile::getMuVBolo2(int ibolo, int imesure) { // microvolts, filtre … … 1001 1030 1002 1031 // Attention, on ne transmet pas les canaux 0 et 4.... 1032 if (ichannel == 0 || ichannel == 4) return 0; 1003 1033 int i = ichannel - 1; 1004 1034 if (i >= 4) i--; -
trunk/Poubelle/archTOI.old/archeopsfile.h
r413 r416 131 131 int getBoloCompBlockNum(); 132 132 long getRawBolo(int ibolo, int imesure); // donnee brute, avec seulement soustraction offset 133 long getRawBoloCN(int ibolo, int imesure); // bruit de compression, sur donnee brute 134 double getMuVBoloCN(int ibolo, int imesure); // bruit de compression, en microVolts 133 135 double getMuVBolo(int ibolo, int imesure); // microvolts, filtre avec filtre carre 134 136 double getMuVBolo2(int ibolo, int imesure); // microvolts, filtre +sioux -
trunk/Poubelle/archTOI.old/archtoi.cc
r408 r416 12 12 #include "archtoi.h" 13 13 #include "archparam.h" 14 #include "archvers.h" 14 15 15 16 using namespace std; -
trunk/Poubelle/archTOI.old/archtoi.h
r407 r416 67 67 }; 68 68 69 #define ARCHTOI_VERS "1.1"70 #define ARCHTOI_TAG "notag"71 72 69 #endif -
trunk/Poubelle/archTOI.old/auxinterpgps.cc
r363 r416 4 4 #include "auxinterpgps.h" 5 5 #include "fitsio.h" 6 #include "arch toi.h"6 #include "archvers.h" 7 7 8 8 AuxInterpGPS::AuxInterpGPS() { -
trunk/Poubelle/archTOI.old/ssthandler.cc
r405 r416 15 15 #include "pisteetoile.h" 16 16 17 #ifdef SST_DUMPLOG 17 18 ofstream SSTHandler::sstchass("SSTChassLogFile"); 18 19 bool SSTHandler::sstchassinit=false; 20 #endif 19 21 20 22 // diodpermut[i] = channel de la diode i … … 255 257 LastBlockStarVec.push_back((*StarHistReIter).second); 256 258 257 //#ifdef SST_DEBUG259 #ifdef SST_DUMPLOG 258 260 if(NbStarInBlock>0) { 259 261 … … 265 267 } 266 268 267 //#endif269 #endif 268 270 269 271 return NbStarInBlock; … … 274 276 for(int i=0; i<NbPhotDiodBarette; i++) PisteBar[i]=new PisteEtoile(i); 275 277 StarHistoryMap.clear(); 278 #ifdef SST_DUMPLOG 276 279 if (!sstchassinit) { 277 280 sstchass<<LastStar.printHeader()<<endl; 278 281 sstchassinit = true; 279 282 } 283 #endif 280 284 281 285 #ifdef SST_DEBUG -
trunk/Poubelle/archTOI.old/ssthandler.h
r412 r416 24 24 25 25 //#define SST_DEBUG 26 //#define SST_DUMPLOG 26 27 #define NbChanDump (10) 27 28 #define NbPhotDiodBarette (46) -
trunk/Poubelle/archTOI.old/toiiter.cc
r413 r416 215 215 case boloRaw: 216 216 case boloRawCN: 217 case boloTensCN: 217 218 trigMask |= block_bolo_mask; 218 219 break; … … 444 445 if (imes==0 && file->llastBolo()==NULL) return false; 445 446 return file->lastBolo() != NULL; 447 case boloTensCN: 446 448 case boloRaw: 447 449 return file->lastBolo() != NULL; … … 533 535 case boloRaw: 534 536 return file->getRawBolo(index, imes); 537 case boloRawCN: 538 return file->getRawBoloCN(index, imes); 539 case boloTensCN: 540 return file->getMuVBoloCN(index, imes); 535 541 case boloGainAmpli: 536 542 return file->getGainAmpli(index); … … 613 619 return file->blockNum() == file->getBoloBlockNum(); 614 620 case boloRaw: 621 case boloRawCN: 622 case boloTensCN: 615 623 return file->blockNum() == file->getBoloBlockNum(); 616 624 case boloGainAmpli: -
trunk/Poubelle/archTOI.old/toiiter.h
r410 r416 26 26 boloDACI, // courant (microAmps) 27 27 boloRawCN, // compression noise on bolo raw 28 boloTensCN, // compression noise on bolo (microVolts), no filter 28 29 boloTens2T, // total, microVolts, avec consigne DAC 29 30 boloRes, // Resistance, ohms -
trunk/Poubelle/archTOI.old/toisvr.cc
r414 r416 8 8 9 9 TOISvr::TOISvr() 10 {} 10 { 11 defaultInclude(); 12 } 11 13 12 14 // To avoid special copy constructor handling, we will not … … 109 111 else tsttoi(boloTens2T) 110 112 else tsttoi(boloRawCN) 113 else tsttoi(boloTensCN) 111 114 else tsttoi(dilDAC) 112 115 else tsttoi(dilSwitch) … … 215 218 return true; 216 219 } 220 221 void TOISvr::defaultInclude() { 222 processRequest("#REQVERSION V_230999"); 223 processRequest("#MJD0 1376.8358818"); 224 processRequest("#PERECH 0.005836818076"); 225 processRequest("#COMMENT ***WARNING***"); 226 processRequest("#COMMENT ***SOME TOI'S ARE PRELIMINARY***"); 227 processRequest("#COMMENT gyroSpeed is not calibrated"); 228 processRequest("#COMMENT azimut/alpha/delta use galaxy crossings"); 229 processRequest("#COMMENT and assume no pendulation"); 230 231 } -
trunk/Poubelle/archTOI.old/toisvr.h
r407 r416 42 42 virtual bool processOption(string keyw, string args); 43 43 44 void defaultInclude(); 45 44 46 protected: 45 47 TOIIter iter;
Note:
See TracChangeset
for help on using the changeset viewer.