Ignore:
Timestamp:
Nov 21, 2011, 3:35:31 PM (13 years ago)
Author:
campagne
Message:

improve ON et OFF sur OFF filtre (jec)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • BAORadio/AmasNancay/trunk/mergeAnaFiles.cc

    r606 r607  
    304304                              ,"onoffRun01420side","onoffRun11420side"
    305305                              ,"onoffCycle01420side","onoffCycle11420side"
    306                               ,"onoffRaw0f14001420","onoffRaw1f14001420"
    307                               ,"offRaw0f14001420","offRaw1f14001420"
    308                               ,"onRaw0f14001420","onRaw1f14001420"
     306                              ,"onoffRaw0f14101415","onoffRaw1f14101415"
     307                              ,"offRaw0f14101415","offRaw1f14101415"
     308                              ,"onRaw0f14101415","onRaw1f14101415"
    309309  };
    310310  NTuple onoffevolution(NINFO,onffTupleName);
     
    327327  sa_size_t ch1400 = freqToChan(1400);
    328328  sa_size_t ch1405 = freqToChan(1400);
     329  sa_size_t ch1410 = freqToChan(1410);
    329330  sa_size_t ch1415 = freqToChan(1415);
    330331  sa_size_t ch1420 = freqToChan(1420);
     
    758759      aSpecOffovOff.Div(aSpecOffFitltered);
    759760
    760       TVector<r_4> meanInRange_1400a1420Freq_OFF_noCalib(NUMBER_OF_CHANNELS);
    761       meanInRange(aSpecOffovOff,ch1400,ch1420,meanInRange_1400a1420Freq_OFF_noCalib);
    762 
    763       xnt[21] = meanInRange_1400a1420Freq_OFF_noCalib(0);
    764       xnt[22] = meanInRange_1400a1420Freq_OFF_noCalib(1);
    765 
    766       //      TMatrix<r_4> aSpecOnFitltered(NUMBER_OF_CHANNELS,NUMBER_OF_FREQ);
    767       //      medianFiltering(aSpecOn,halfWidth,aSpecOnFitltered);
     761      TVector<r_4> meanInRange_1410a1415Freq_OFF_noCalib(NUMBER_OF_CHANNELS);
     762      meanInRange(aSpecOffovOff,ch1410,ch1415,meanInRange_1410a1415Freq_OFF_noCalib);
     763
     764      xnt[21] = meanInRange_1410a1415Freq_OFF_noCalib(0);
     765      xnt[22] = meanInRange_1410a1415Freq_OFF_noCalib(1);
     766
    768767      TMatrix<r_4> aSpecOnovOff(aSpecOn,false);
    769768      aSpecOnovOff.Div(aSpecOffFitltered);
    770769
    771       TVector<r_4> meanInRange_1400a1420Freq_ON_noCalib(NUMBER_OF_CHANNELS);
    772       meanInRange(aSpecOnovOff,ch1400,ch1420,meanInRange_1400a1420Freq_ON_noCalib);
    773 
    774       xnt[23] = meanInRange_1400a1420Freq_ON_noCalib(0);
    775       xnt[24] = meanInRange_1400a1420Freq_ON_noCalib(1);
     770      TVector<r_4> meanInRange_1410a1415Freq_ON_noCalib(NUMBER_OF_CHANNELS);
     771      meanInRange(aSpecOnovOff,ch1410,ch1415,meanInRange_1410a1415Freq_ON_noCalib);
     772
     773      xnt[23] = meanInRange_1410a1415Freq_ON_noCalib(0);
     774      xnt[24] = meanInRange_1410a1415Freq_ON_noCalib(1);
    776775     
    777776      //store infos to Ntuple
     
    852851//Used like:
    853852//
    854 void meanRawDiffOnOffCycles() throw(string) {
     853void
     854meanRawDiffOnOffCycles() throw(string) {
    855855  list<string> listOfFiles;
    856856  string directoryName;
Note: See TracChangeset for help on using the changeset viewer.