|
wave | Add (net_file_name) |
|
live | Add (liv_file_name) |
|
TB | checkFile (gSystem->Getenv("CWB_ROOTLOGON_FILE")) |
|
TB | checkFile (gSystem->Getenv("CWB_PARAMETERS_FILE")) |
|
TB | checkFile (gSystem->Getenv("CWB_UPARAMETERS_FILE")) |
|
TB | checkFile (gSystem->Getenv("CWB_PPARAMETERS_FILE")) |
|
TB | checkFile (gSystem->Getenv("CWB_UPPARAMETERS_FILE")) |
|
TB | checkFile (gSystem->Getenv("CWB_EPPARAMETERS_FILE")) |
|
c1 | Clear () |
|
lag | Divide (llag) |
|
wave | Draw ("lag[0]-lag[1]:lag[1]-lag[2]", cut,"goff") |
|
lag | Draw ("colz") |
|
FarPlot | Draw ("AP") |
|
| exit (0) |
|
| for (int i=0;i< num;i++) |
|
gROOT | ForceStyle () |
|
FarPlot | GetHistogram () -> SetXTitle("#lag distance [s]") |
|
cout<< "Number of lags : "
<< count<< endl;double liveTot=llag-> | GetSumOfWeights () |
|
lag | GetYaxis () -> SetTitle("lag[1]-lag[2] shift [s]") |
|
| if (nIFO!=3) |
|
| if (!test_veto.IsNull()) sprintf(cut |
|
TChain | live ("liveTime") |
|
c1 | Range (-0.75,-1.23433, 6.75, 8.17182) |
|
c1 | SaveAs (fname) |
|
c1 | SetBorderSize (1) |
|
gStyle | SetCanvasColor (kWhite) |
|
wave | SetEstimate (num) |
|
c1 | SetFillColor (0) |
|
gStyle | SetFrameBorderMode (0) |
|
c1 | SetFrameFillColor (0) |
|
c1 | SetGridx () |
|
c1 | SetGridy () |
|
gStyle | SetLineColor (kWhite) |
|
cout<< "Total BKG live time : "
<< liveTot<< endl;double FAR2[count];double
eFAR2[count];double N[count];double
eN[count];for(int i=0;i< count;i++){N[i]=D[i];}Int_t
*index=new Int_t[count];TMath::Sort(count,
N, index, false);for(int i=0;i
< count;i++){FAR2[i]=FAR[index[i]];eFAR2[i]=eFAR[index[i]];N[i]=D[index[i]];eN[i]=0.0;}TGraphErrors
*FarPlot=new TGraphErrors(count,
N, FAR2, eN, eFAR2);FarPlot-> | SetLineColor (kBlue) |
|
FarPlot | SetMarkerColor (kBlue) |
|
FarPlot | SetMarkerSize (1) |
|
FarPlot | SetMarkerStyle (20) |
|
gStyle | SetNumberContours (256) |
|
gStyle | SetOptStat (kFALSE) |
|
gStyle | SetPalette (1, 0) |
|
gStyle | SetStatBorderSize (1) |
|
lag | SetStats (kFALSE) |
|
FarPlot | SetTitle ("FAR vs lag distance") |
|
gStyle | SetTitleFillColor (kWhite) |
|
| sprintf (cut,"rho[1]>%f", T_cut) |
|
| sprintf (cut2,"!(lag[%d]==0&&slag[%d]==0)", nIFO, nIFO) |
|
| sprintf (fname,"%s/lag.png", netdir) |
|
| sprintf (fname,"%s/FAR_vs_lag.png", netdir) |
|
c1 | Update () |
|
TChain | wave ("waveburst") |
|
|
TCanvas * | c1 = new TCanvas("lag", "lag",32,55,750,502) |
|
int | count = 0 |
|
char | cut [512] |
|
char | cut2 [256] |
|
double | D [100000] |
|
double | eFAR [100000] |
|
double | FAR [100000] |
|
char | fname [1024] |
|
double | liveMax = llag->GetMaximum() |
|
TH2F * | llag = new TH2F("LSLAG","FAR distribution over lags",Nlag,xplotlagMin,xplotlagMax,Nlag,yplotlagMin,yplotlagMax) |
|
int | num = (int)wave.GetEntries(cut) |
|
s && | s |
|
cout<<"Number of Entries: "
<< num<< endl;double *lag01=wave.GetV1();double
*lag12=wave.GetV2();int num2=(int)
live.GetEntries(cut2);live.SetEstimate(num2);live.Draw("lag[0]-lag[1]:lag[1]-lag[2]:live",
cut2,"goff");double *llag01=live.GetV1();double
*llag12=live.GetV2();double
*Live=live.GetV3();char title[256];double
xplotlagMax=TMath::MaxElement(num,
lag01);double xplotlagMin=TMath::MinElement(num,
lag01);double yplotlagMax=TMath::MaxElement(num,
lag12);double yplotlagMin=TMath::MinElement(num,
lag12);cout<< xplotlagMin
<<" < lag[o]-lag[1]
< "<< xplotlagMax<< endl;cout
<< yplotlagMin<<" < lag[1]-lag[2]
< "<< yplotlagMax<< endl;int
Nlag=100;sprintf(title,"FAR
distribution over lags (post
cat3 & rho>%f)", T_cut);TH2F
*lag=new TH2F("SLAG", title,
Nlag, xplotlagMin, xplotlagMax,
Nlag, yplotlagMin, yplotlagMax);lag->
GetXaxis() ->SetTitle("lag[0]-lag[1]
shift [s]");lag->GetXaxis() -> | SetNdivisions (10, kFALSE) |
|
CWB::Toolbox | TB |
|
TString | test_veto = veto_not_vetoed |
|
s | veto_not_vetoed |
|
TAxis * | xaxis = lag->GetXaxis() |
|
TAxis * | yaxis = lag->GetYaxis() |
|