CreateUniqueLags3ifos.C

Go to the documentation of this file.

/*
# Copyright (C) 2019 Gabriele Vedovato
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program.  If not, see <https://www.gnu.org/licenses/>.
*/
 
 
void CreateUniqueLags3ifos(TString odir, int nLags=599) {
 
  if(nLags%2==0) {
    cout << "CreateUniqueLags3ifos - Error : nLags must be odd !!!" << endl;
    exit(1);
  } 
 
  char ofname[1024];
  sprintf(ofname,"%s/Unique_n%d.lags",odir.Data(),nLags);
 
  ofstream out;
  out.open(ofname,ios::out);
  if(!out.good()) {cout << "CreateUniqueLags3ifos - Error : Opening File : " << ofname << endl;gSystem->Exit(1);}
 
  out << 0 << "\t" << 0 << "\t" << 0 << "\t" << 0 << endl;
  for(int i=1;i<nLags;i++) {
    int lag0 = 0;
    int lag1 = i%2 ? i : nLags-i+1;
    int lag2 = i%2 ? nLags-i : i-1;
    out << i << "\t" << lag0 << "\t" << lag1 << "\t" << lag2 << endl;
  }
 
  out.close();
 
  cout << endl << "Created file : " << ofname << endl << endl;
 
  exit(0);
}