33 #include "TGraphErrors.h" 133 size_t cpf(
const netcluster &,
bool=
false,
int=0);
141 inline void clean(
int cID=0);
147 inline size_t size() {
return pList.size(); }
149 inline size_t capacity() {
return pList.capacity(); }
151 inline size_t csize() {
return cList.size(); }
155 for(
size_t i=0;
i<sCuts.size();
i++) {
156 if(
k<2 && sCuts[
i]==
k) n++;
157 if(
k>1 && sCuts[
i] <1) n++;
165 for(
size_t i=0;
i<sCuts.size();
i++) {
166 if(
k>1 && sCuts[
i] <1) n+=cList[
i].size();
167 if(
k<2 && sCuts[
i]==
k) n+=cList[
i].size();
175 inline netpixel* getPixel(
size_t n,
size_t i);
185 size_t setcore(
bool core,
int id=0);
191 if(n>0 && n<=sCuts.size()) sCuts[n-1] = 1;
197 for(
size_t i=0; i<sCuts.size(); i++)
198 if(sCuts[i]!=1) sCuts[
i] =
n;
204 virtual size_t cleanhalo(
bool=
false);
209 size_t addhalo(
int=0);
222 virtual size_t delink();
236 virtual size_t supercluster(
char atype,
double S,
bool core);
237 virtual size_t supercluster(
char atype,
double S,
double gap,
bool core, TH1F* = NULL);
242 virtual size_t defragment(
double T,
double F, TH2F* = NULL);
262 virtual size_t cluster(
int kt,
int kf);
295 size_t write(
const char* file,
int app=0);
302 size_t write(FILE* fp,
int app=0);
319 size_t read(
const char *);
328 size_t read(FILE* file,
int maxPix);
343 std::vector<int> read(TFile* froot,
TString tdir,
TString tname,
int nmax=0,
int cycle=0,
int rate=0,
int cID=0);
357 size_t loadTDampSSE(
network& net,
char c,
size_t BATCH=10000,
size_t LOUD=0);
364 double mchirp(
int ID,
double=2.5,
double=1.e20,
double=0);
367 double mchirpTEST(
int ID);
370 void chirpDraw(
int id);
414 if(!n)
return &pList[
i];
416 printf(
"getPixel(): size=%d, ID=%d\n",(
int)cList.size(),(
int)n);
419 if(cList[n-1].
size()<=i) {
420 printf(
"getPixel(): size=%d, index=%d\n",(
int)cList[n-1].
size(),(
int)i);
423 return &pList[cList[n-1][
i]];
428 while(!pList.empty()) pList.pop_back();
429 pList.clear(); std::vector<netpixel>().swap(pList);
430 while(!cList.empty()) cList.pop_back();
431 cList.clear(); std::vector<vector_int>().swap(cList);
432 while(!cData.empty()) cData.pop_back();
433 cData.clear(); std::vector<clusterdata>().swap(cData);
434 while(!sCuts.empty()) sCuts.pop_back();
435 sCuts.clear(); std::vector<int>().swap(sCuts);
436 while(!cRate.empty()) cRate.pop_back();
437 cRate.clear(); std::vector<vector_int>().swap(cRate);
439 cTime.clear(); std::vector<float>().swap(
cTime);
441 cFreq.clear(); std::vector<float>().swap(
cFreq);
442 while(!sArea.empty()) sArea.pop_back();
443 sArea.clear(); std::vector<vector_float>().swap(sArea);
444 while(!p_Map.empty()) p_Map.pop_back();
445 p_Map.clear(); std::vector<vector_float>().swap(p_Map);
446 while(!p_Ind.empty()) p_Ind.pop_back();
447 p_Ind.clear(); std::vector<vector_int>().swap(p_Ind);
452 for(
int i=0;
i<(
int)this->cList.size(); ++
i) {
454 if((cID!=0)&&((
int)pList[vint[0]].clusterID!=cID))
continue;
455 for(
int j=0;
j<(
int)vint.size();
j++) {
456 if(pList[vint[
j]].tdAmp.size()) pList[vint[
j]].clean();
461 #endif // NETCLUSTER_HH
std::vector< int > vector_int
std::vector< vector_int > cRate
int compare_PIX(const void *, const void *)
std::vector< pixel > cluster
std::vector< vector_float > sArea
std::vector< int > vector_int
void ignore(size_t n)
param: cluster ID number return void
std::vector< vector_int > cList
cout<< "SNR "<< xsnr<< endl;wavearray< double > f
std::vector< vector_float > p_Map
network ** net
NOISE_MDC_SIMULATION.
std::vector< float > vector_float
std::vector< netpixel > pList
printf("total live time: non-zero lags = %10.1f \, liveTot)
std::vector< vector_int > p_Ind
std::vector< float > cFreq
wavearray< float > mchpdf
void setcuts(int n=0)
param: sCuts flag
virtual void Browse(TBrowser *b)
netpixel * getPixel(size_t n, size_t i)
TGraphErrors chirp
chirp fit parameters (don't remove ! fix crash when exit from CINT)
std::vector< float > cTime
Meyer< double > S(1024, 2)
std::vector< clusterdata > cData
std::vector< vector_int > nTofF