lcmath - Subtracts or adds two binned lightcurves.
lcmath infile bgfile outfile multi multb
This program takes as input two binned lightcurves and calculates a result lightcurve either subtracting (second from the first) or adding the two input files. The two input files should cover the same time interval. In particular the start time and the stop time of the first lightcurve should be included within the start and stop time of the second lightcurve. Error are prompt if the start and/or the stop are outside of the time interval of the second file. The time resolution of the two files can be different, but the integration time of the first input file should be always less than or equal to that in the seconf input file. The task does not work for input files containing "EVENT" table. Constant scaling and additive factors can be separately input for both lightcurves. By default, the task will apply any vignetting (or collimator) and deadtime correction information contained in the input file (header keywords VIGNET and DEADC) to the result lightcurve. The user may cancel this correction setting the parameter "docor". The columns names recognized by the task are TIME, COUNT, RATE and ERROR. The units (TUNITn keyword) are expected to be 'count' for column name 'COUNT' and 'count/s' for column name RATE. The input file structure is a FITS bintable with data stored in the column COUNT or RATE either as a single element or a vector of element in each row. For the last case the task works correctly only if the column n contains a vector of counts in different channel (the value of the column keyword 1CTYPn should be 'CHANNEL'). The paramaters "emin" and "emax" are applicable in this case. The structure of the output file will have always data stored as a single element in each row.
If scaling and additive factors are input for both file1 and file2 the derived counts are calculated as:
count1=M1*x1 + C1 count2=M2*x2 + C2
where M1 and M2 are the input scaling values for the first and second lightcurve respectively; C1 and C2 are the additive input factors in counts (rescaling for the fractional exposure) and x1 and x2 are the value always in counts read from the files (if the file contains rate x1 and x2 are the value after rescaling for the integration and exposure time). The final results in the output file will be:
count3= count1 -/+ count2
divided for the integration time. The error for count3 is calculated as:
err3= sqrt( err1**2 + err2**2)
where err3, err2, and err1 are the error values for the count3, count2 and count1 respectively. To calculate the individual errors (i.e. err1 for the quantity count1 and err2 for the quantity count2) the user can choose one of the 6 different methods listed below.
If the file contains an error column, sx1 is the value in the file rescaled for the integration time and exposure, otherwise sx1 is set as the sqrt(x1) or sqrt(x2). If x1 or x2 are 0 the error (sx1 or sx2) is set to 1. If the scaling factor is 1 and the additive factor is 0, the methods 1,2,3,4 are equivalent.
Report problems to firstname.lastname@example.org, email@example.com. Provide a detailed description of the problem (possible a log file).