Commit 03fe27b8 authored by Yizhou Cai's avatar Yizhou Cai
Browse files

update

parent 26ef6f04
CXX = g++ -std=c++11
FC = g77
ROOTINC= -I/fsa/home/sjc_caiyz/delphes_install/external/ExRootAnalysis/
ROOTLIB= -L/fsa/home/sjc_caiyz/delphes_install/external/ExRootAnalysis/ -lExRootAnalysis \
-L/fsa/home/sjc_caiyz/delphes_install/external/ExRootAnalysis/
ROOTFLAG=${shell root-config --cflags --libs}
##FJFLAG=${shell fastjet-config --cxxflags --libs --plugins}
ROOTINC2= -I/fsa/home/sjc_caiyz/delphes_install/external \
-I/fsa/home/sjc_caiyz/delphes_install/display \
-I/fsa/home/sjc_caiyz/root_install/include \
-I/fsa/home/sjc_caiyz/delphes_install
ROOTLIB2= -L/fsa/home/sjc_caiyz/delphes_install -lDelphes
CPPFLAGS=${ROOTFLAG} ${ROOTINC2} ${FJFLAG}
LDFLAGS= ${ROOTFLAG} ${ROOTLIB2} ${FJFLAG}
#FOBJS := $(patsubst %.f,%.o,$(wildcard *.f))
CCOBJS := $(patsubst %.cpp,%.o,$(wildcard *.cpp))
mainAll: ${CCOBJS}
$(CXX) *.o ${LDFLAGS} -o $@
%.o : %.cpp
$(CXX) -c $(CPPFLAGS) $& $<
clean:
rm *.o mainAll
#include <cmath>
#include <fstream>
#include <iostream>
#include <string>
#include <cstdlib>
#include <iomanip>
#include <fstream>
#include <map>
#include <complex>
#include <stdexcept>
#include <vector>
#include <stdio.h>
#include <stdlib.h>
#include <sstream>
#include <iomanip>
#include <iostream>
#include <fstream>
#include <math.h>
#include <exception>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <stdlib.h>
#include <sstream>
#include <time.h>
#include <unistd.h>
#include "TROOT.h"
#include "TClonesArray.h"
#include "TRint.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TF1.h"
#include "TCanvas.h"
#include "TChain.h"
#include "TVector3.h"
#include "TMath.h"
#include "TVector3.h"
#include "TMath.h"
#include "TFile.h"
#include "TSystem.h"
#include "TGStatusBar.h"
#include "TSystem.h"
#include "TXMLEngine.h"
#include "TTree.h"
#include "TLorentzVector.h"
#include "TNtuple.h"
#include "TStyle.h"
#include "TGraphErrors.h"
#include "TGraph.h"
#include "classes/ClassesLinkDef.h"
#include "classes/DelphesClasses.h"
#include "classes/DelphesFactory.h"
#include "classes/DelphesModule.h"
#include "classes/DelphesStream.h"
#include "ExRootAnalysis/ExRootTreeReader.h"
#include "ExRootAnalysis/ExRootAnalysisLinkDef.h"
#include "ExRootAnalysis/ExRootProgressBar.h"
#include "ExRootAnalysis/ExRootTask.h"
#include "ExRootAnalysis/ExRootTreeBranch.h"
#include "fastjet/PseudoJet.hh"
#include "fastjet/ClusterSequence.hh"
#include "fastjet/internal/BasicRandom.hh"
#include "fastjet/JetDefinition.hh"
#include "fastjet/tools/MassDropTagger.hh"
#include "fastjet/tools/Filter.hh" // needed for HEPTopTagger
#include "fastjet/tools/Pruner.hh" // needed for HEPTopTagger
#include "fastjet/tools/JHTopTagger.hh"
using namespace std;
using namespace TMath;
using namespace fastjet;
bool sortPt(const TLorentzVector a, const TLorentzVector b) {return (a.Pt() < b.Pt());}
double GetdPhi(double phi1, double phi2){
double math_pi = 3.14159265358979;
double dphi = phi1 - phi2;
if(fabs(dphi) > math_pi) dphi = (fabs(dphi)-2*math_pi);
return fabs(dphi)/math_pi;
}
double Imb(double pt1, double pt2) {return (fabs(pt1-pt2)/(pt1+pt2));}
vector<string> split(const string &str, const char pattern){
vector<string> res;
stringstream input(str);
string temp;
while(getline(input, temp, pattern))
{
res.push_back(temp);
}
return res;
}
//Mode should be selected from ('_300','_500')
vector<double> BBMET(TString inputFileName, double value){
// Create chain of root trees
TChain chain("Delphes");
string name = split(split(string(inputFileName),'/').back(),'.')[0];
if (string(inputFileName).find("root") == string::npos){
chain.Add(inputFileName+"/Events/run_01/tag_1_delphes_events.root");
chain.Add(inputFileName+"/Events/run_02/tag_1_delphes_events.root");
chain.Add(inputFileName+"/Events/run_03/tag_1_delphes_events.root");
}
else {chain.Add(inputFileName);}
// Create object of class ExRootTreeReader
ExRootTreeReader *treeReader = new ExRootTreeReader(&chain);
Long64_t numberOfEntries = treeReader->GetEntries();
TClonesArray *branchEvent = treeReader->UseBranch("Event");
TClonesArray *branchJet = treeReader->UseBranch("Jet");
TClonesArray *branchMissingET = treeReader->UseBranch("MissingET");
TClonesArray *branchPhoton = treeReader->UseBranch("Photon");
TClonesArray *branchElectron = treeReader->UseBranch("Electron");
TClonesArray *branchMuon = treeReader->UseBranch("Muon");
Jet* jet;
Photon* photon;
Electron* electron;
Muon* muon;
MissingET* met;
Track* track;
string path1=("./data/left_"), path2=(".dat");
string path=(path1.append(name)).append(path2);
string path3="./total.dat";
ofstream total;
total.open(path3, ios::app);
double cut_=0.0, cut_nlp=0.0, cut_njet=0.0, cut_mh=0.0, cut_met=0.0, cut_pth=0.0, cut_imb=0.0;
double total_weight = 0.0;
double weight_nlp=0.0, weight_njet=0.0, weight_mh=0.0, weight_met=0.0, weight_pth=0.0, weight_imb=0.0;
for(int entry = 0; entry < numberOfEntries; ++entry){
treeReader->ReadEntry(entry);
int fl_nlp=0, fl_njet=0, fl_mh=0, fl_met=0, fl_pth=0, fl_imb=0;
//Event weight
HepMCEvent *event = (HepMCEvent*) branchEvent -> At(0);
float weight = event->Weight;
total_weight += weight;
//Containers filling
vector<TLorentzVector> All_jet, All_bjet, All_photon, All_lepton;
TLorentzVector MET, TrkMiss;
for(int i=0; i<branchJet->GetEntriesFast(); i++){
jet = (Jet*) branchJet -> At(i);
TLorentzVector tjet;
tjet = jet->P4();
if (jet->TauTag) All_lepton.push_back(tjet);
if (jet->PT > 20. && fabs(jet->Eta) < 2.5){
All_jet.push_back(tjet);
if(jet->BTag) All_bjet.push_back(tjet);
}
}
sort(All_bjet.begin(), All_bjet.end(), sortPt);
for(int i=0; i<branchPhoton->GetEntriesFast(); i++){
photon = (Photon*) branchPhoton -> At(i);
if (photon->PT > 10. && fabs(photon->Eta) < 2.5){
TLorentzVector tphoton;
tphoton = photon->P4();
All_photon.push_back(tphoton);
}
}
for(int i=0; i<branchElectron->GetEntriesFast(); i++){
electron = (Electron*) branchElectron -> At(i);
if (electron->PT > 10. && fabs(electron->Eta) < 2.5){
TLorentzVector telectron;
telectron = electron->P4();
All_lepton.push_back(telectron);
}
}
for(int i=0; i<branchMuon->GetEntriesFast(); i++){
muon = (Muon*) branchMuon -> At(i);
if (muon->PT > 10. && fabs(muon->Eta) < 2.4){
TLorentzVector tmuon;
tmuon = muon->P4();
All_lepton.push_back(tmuon);
}
}
if(branchMissingET->GetEntriesFast() > 0){
met = (MissingET*) branchMissingET->At(0);
MET = met->P4();
}
float deta = 0;
float dPhi = 0;
if (All_bjet.size() == 2){
deta = fabs(All_bjet[0].Eta()-All_bjet[1].Eta());
dPhi = fabs(All_bjet[0].DeltaPhi(All_bjet[1]));
}
double dRmin = 99999.;
if (All_jet.size() > 2) {
for (unsigned int ir=0;ir<All_jet.size();ir++)
for (unsigned int jr=ir+1;jr<All_jet.size();jr++)
if(All_jet[ir].DeltaR(All_jet[jr])<dRmin){
dRmin=All_jet[ir].DeltaR(All_jet[jr]);
}
}
if (All_lepton.size() == 0 && All_photon.size() == 0) fl_nlp = 1;
if (All_bjet.size() == 2 && All_jet.size() < 4) fl_njet = 1;
if (All_bjet.size() == 2 && ((All_bjet[0] + All_bjet[1]).M() < 140.) && ((All_bjet[0] + All_bjet[1]).M() > 105.)) fl_mh = 1;
if (MET.Pt() > 40.) fl_met = 1;
if ((deta>0.5)&&(dPhi>2.2)) fl_imb = 1;
if (dRmin>2.8) fl_pth = 1;
//fl_met = 1;fl_imb = 1;fl_pth = 1;
cut_+=1.0;
if(fl_nlp==1){cut_nlp+=1.0;weight_nlp+=weight;}
if(fl_nlp==1 && fl_njet==1){cut_njet+=1.0;weight_njet+=weight;}
if(fl_nlp==1 && fl_njet==1 && fl_mh==1){cut_mh+=1.0;weight_mh+=weight;}
if(fl_nlp==1 && fl_njet==1 && fl_mh==1 && fl_met==1){cut_met+=1.0;weight_met+=weight;}
if(fl_nlp==1 && fl_njet==1 && fl_mh==1 && fl_met==1 && fl_imb==1){cut_imb+=1.0;weight_imb+=weight;}
if(fl_nlp==1 && fl_njet==1 && fl_mh==1 && fl_met==1 && fl_imb==1 && fl_pth==1){cut_pth+=1.0;weight_pth+=weight;}
}
vector<double> Ratio;
//double R_nlp=cut_nlp/cut_, R_njet=cut_njet/cut_, R_mh=cut_mh/cut_, R_met=cut_met/cut_, R_imb=cut_imb/cut_, R_pth=cut_pth/cut_;
double R_nlp=weight_nlp/total_weight, R_njet=weight_njet/total_weight, R_mh=weight_mh/total_weight, R_met=weight_met/total_weight, R_imb=weight_imb/total_weight, R_pth=weight_pth/total_weight;
Ratio.push_back(R_nlp); Ratio.push_back(R_njet); Ratio.push_back(R_mh); Ratio.push_back(R_met); Ratio.push_back(R_imb); Ratio.push_back(R_pth);
total<<name<<" EventNumbers="<<numberOfEntries<<" TotalWeight="<<total_weight<<endl;
total<<"Ratio_N(l+p) = "<<Ratio[0]<<" "<<Ratio[0]<<" "<<cut_nlp<<endl;
total<<"Ratio_N(jet) = "<<Ratio[1]<<" "<<Ratio[1]/Ratio[0]<<" "<<cut_njet<<endl;
total<<"Ratio_M(bb) = "<<Ratio[2]<<" "<<Ratio[2]/Ratio[1]<<" "<<cut_mh<<endl;
total<<"Ratio_MET = "<<Ratio[3]<<" "<<Ratio[3]/Ratio[2]<<" "<<cut_met<<endl;
total<<"Ratio_bJetSplit = "<<Ratio[4]<<" "<<Ratio[4]/Ratio[3]<<" "<<cut_imb<<endl;
total<<"Ratio_dRmin = "<<Ratio[5]<<" "<<Ratio[5]/Ratio[4]<<" "<<cut_pth<<endl;
total<<endl;
return Ratio;
}
void round(TString inputFileName, double value){
vector<double> Ratio_Cutflow_s,Ratio_Cutflow_1,Ratio_Cutflow_2,Ratio_Cutflow_3,Ratio_Cutflow_4,Ratio_Cutflow_5,Ratio_Cutflow_6,Ratio_Cutflow_7,Ratio_Cutflow_8;
Ratio_Cutflow_s=BBMET(inputFileName,value);
Ratio_Cutflow_1=BBMET("/fsa/home/sjc_caiyz/samples/ttbar/pp_ttbar.root",value);
Ratio_Cutflow_2=BBMET("/fsa/home/sjc_caiyz/samples/VH/VH_W_lv.root",value);
Ratio_Cutflow_3=BBMET("/fsa/home/sjc_caiyz/samples/VH/VH_Z_vv.root",value);
Ratio_Cutflow_4=BBMET("/fsa/home/sjc_caiyz/samples/Vjj/Wjj.root",value);
Ratio_Cutflow_5=BBMET("/fsa/home/sjc_caiyz/samples/Vjj/Zjj.root",value);
Ratio_Cutflow_6=BBMET("/fsa/home/sjc_caiyz/samples/WZ/pp_WZ_bblv.root",value);
Ratio_Cutflow_7=BBMET("/fsa/home/sjc_caiyz/samples/WZ/pp_WZ_bjvv.root",value);
Ratio_Cutflow_8=BBMET("/fsa/home/sjc_caiyz/samples/WZ/pp_ZZ_bbvv.root",value);
double total_bkg = 3000000*(Ratio_Cutflow_1.back()*89+Ratio_Cutflow_2.back()*0.3+Ratio_Cutflow_3.back()*0.2+Ratio_Cutflow_4.back()*45871+Ratio_Cutflow_5.back()*12810+Ratio_Cutflow_6.back()*1.58+Ratio_Cutflow_7.back()*3.02+Ratio_Cutflow_8.back()*4.2);
double cls = 1.96*pow(total_bkg,0.5)/(Ratio_Cutflow_s.back()*3000000);
cout<<value+20<<" "<<total_bkg<<" "<<cls<<endl;
Ratio_Cutflow_s.pop_back();Ratio_Cutflow_1.pop_back();Ratio_Cutflow_2.pop_back();Ratio_Cutflow_3.pop_back();Ratio_Cutflow_4.pop_back();Ratio_Cutflow_5.pop_back();Ratio_Cutflow_6.pop_back();Ratio_Cutflow_7.pop_back();Ratio_Cutflow_8.pop_back();
total_bkg = 3000000*(Ratio_Cutflow_1.back()*89+Ratio_Cutflow_2.back()*0.3+Ratio_Cutflow_3.back()*0.2+Ratio_Cutflow_4.back()*45871+Ratio_Cutflow_5.back()*12810+Ratio_Cutflow_6.back()*1.58+Ratio_Cutflow_7.back()*3.02+Ratio_Cutflow_8.back()*4.2);cls = 1.96*pow(total_bkg,0.5)/(Ratio_Cutflow_s.back()*3000000);
cout<<value+10<<" "<<total_bkg<<" "<<cls<<endl;
Ratio_Cutflow_s.pop_back();Ratio_Cutflow_1.pop_back();Ratio_Cutflow_2.pop_back();Ratio_Cutflow_3.pop_back();Ratio_Cutflow_4.pop_back();Ratio_Cutflow_5.pop_back();Ratio_Cutflow_6.pop_back();Ratio_Cutflow_7.pop_back();Ratio_Cutflow_8.pop_back();
total_bkg = 3000000*(Ratio_Cutflow_1.back()*89+Ratio_Cutflow_2.back()*0.3+Ratio_Cutflow_3.back()*0.2+Ratio_Cutflow_4.back()*45871+Ratio_Cutflow_5.back()*12810+Ratio_Cutflow_6.back()*1.58+Ratio_Cutflow_7.back()*3.02+Ratio_Cutflow_8.back()*4.2);cls = 1.96*pow(total_bkg,0.5)/(Ratio_Cutflow_s.back()*3000000);
cout<<value<<" "<<total_bkg<<" "<<cls<<endl;
/*for(int i=0; i<Ratio_Cutflow_s.size(); i++){
double total_bkg = 3000000*(Ratio_Cutflow_1[i]*89+Ratio_Cutflow_2[i]*0.3+Ratio_Cutflow_3[i]*0.2+Ratio_Cutflow_4[i]*45871+Ratio_Cutflow_5[i]*12810+Ratio_Cutflow_6[i]*1.58+Ratio_Cutflow_7[i]*3.02+Ratio_Cutflow_8[i]*4.2);
double cls = 1.96*pow(total_bkg,0.5)/(Ratio_Cutflow_s[i]*3000000);
cout<<i+1<<": "<<Ratio_Cutflow_s[i]<<" "<<total_bkg<<" "<<cls<<endl;
}*/
}
int main(int argc, char *argv[])
{
TString inputFileName(argv[1]);
double value;
vector<double> Ratio_Cutflow_s,Ratio_Cutflow_1,Ratio_Cutflow_2,Ratio_Cutflow_3,Ratio_Cutflow_4,Ratio_Cutflow_5,Ratio_Cutflow_6,Ratio_Cutflow_7,Ratio_Cutflow_8;
cout<<"1#"<<endl;
value=10.; round(inputFileName,value);
cout<<"2#"<<endl;
value=40.; round(inputFileName,value);
cout<<"3#"<<endl;
value=70.; round(inputFileName,value);
cout<<"4#"<<endl;
value=100.; round(inputFileName,value);
cout<<"5#"<<endl;
value=130.; round(inputFileName,value);
cout<<"6#"<<endl;
value=160.; round(inputFileName,value);
cout<<"7#"<<endl;
value=190.; round(inputFileName,value);
cout<<"8#"<<endl;
value=220.; round(inputFileName,value);
cout<<"9#"<<endl;
value=250.; round(inputFileName,value);
cout<<"10#"<<endl;
value=280.; round(inputFileName,value);
return 0;
}
This diff is collapsed.
This diff is collapsed.
#BSUB -q e74820v4
#BSUB -n 3
#BSUB -J ex1
#BSUB -o out_1
#BSUB -e err
source /fsa/home/sjc_caiyz/samples/exroot/setup.sh
OMP_NUM_THREADS="$LSB_DJOB_NUMPROC"
cd /fsa/home/sjc_caiyz/samples/exroot/total
./mainAll 1
#BSUB -q e74820v4
#BSUB -n 3
#BSUB -J ex2
#BSUB -o out_2
#BSUB -e err
source /fsa/home/sjc_caiyz/samples/exroot/setup.sh
OMP_NUM_THREADS="$LSB_DJOB_NUMPROC"
cd /fsa/home/sjc_caiyz/samples/exroot/total
./mainAll 2
#BSUB -q e74820v4
#BSUB -n 3
#BSUB -J ex3
#BSUB -o out_3
#BSUB -e err
source /fsa/home/sjc_caiyz/samples/exroot/setup.sh
OMP_NUM_THREADS="$LSB_DJOB_NUMPROC"
cd /fsa/home/sjc_caiyz/samples/exroot/total
./mainAll 3
#BSUB -q e74820v4
#BSUB -n 3
#BSUB -J ex4
#BSUB -o out_4
#BSUB -e err
source /fsa/home/sjc_caiyz/samples/exroot/setup.sh
OMP_NUM_THREADS="$LSB_DJOB_NUMPROC"
cd /fsa/home/sjc_caiyz/samples/exroot/total
./mainAll 4
#BSUB -q e74820v4
#BSUB -n 3
#BSUB -J ex5
#BSUB -o out_5
#BSUB -e err
source /fsa/home/sjc_caiyz/samples/exroot/setup.sh
OMP_NUM_THREADS="$LSB_DJOB_NUMPROC"
cd /fsa/home/sjc_caiyz/samples/exroot/total
./mainAll 5
#BSUB -q e74820v4
#BSUB -n 3
#BSUB -J ex6
#BSUB -o out_6
#BSUB -e err
source /fsa/home/sjc_caiyz/samples/exroot/setup.sh
OMP_NUM_THREADS="$LSB_DJOB_NUMPROC"
cd /fsa/home/sjc_caiyz/samples/exroot/total
./mainAll 6
1#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.1 6 1.20238e+06 0.017556
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
2#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.2 6 857120 0.021951
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
3#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.3 6 532268 0.0313002
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
4#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.4 6 259428 0.0536897
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
5#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.5 6 105659 0.0971637
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
6#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.6 6 36260.5 0.201092
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
7#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.7 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
8#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.8 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
9#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
0.9 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
10#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
1 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
11#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
1.1 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
12#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
1.2 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
13#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
1.3 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
14#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
1.4 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
15#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
1.5 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
16#
/fsa/home/sjc_caiyz/samples/signal/s_300_2
1.6 6 0 -nan
5 2.25336e+06 0.0107995
4 2.29889e+06 0.0108532
3 2.44907e+06 0.0109468
2 1.45091e+07 0.0152664
1 1.02264e+11 0.226474
------------------------------------------------------------
Sender: LSF System <lsfadmin@r012>
Subject: Job 35720998: <ex1> in cluster <hpc> Done
Job <ex1> was submitted from host <login1> by user <sjc_caiyz> in cluster <hpc> at Thu Dec 16 09:58:22 2021
Job was executed on host(s) <3*r012>, in queue <e74820v4>, as user <sjc_caiyz> in cluster <hpc> at Thu Dec 16 09:58:24 2021
</fsa/home/sjc_caiyz> was used as the home directory.
</fsa/home/sjc_caiyz/samples/exroot> was used as the working directory.
Started at Thu Dec 16 09:58:24 2021
Terminated at Fri Dec 17 05:47:32 2021
Results reported at Fri Dec 17 05:47:32 2021
Your job looked like:
------------------------------------------------------------
# LSBATCH: User input
#BSUB -q e74820v4
#BSUB -n 3
#BSUB -J ex1
#BSUB -o out_1
#BSUB -e err
source /fsa/home/sjc_caiyz/samples/exroot/setup.sh
OMP_NUM_THREADS="$LSB_DJOB_NUMPROC"
cd /fsa/home/sjc_caiyz/samples/exroot/total
./mainAll 1
------------------------------------------------------------
Successfully completed.
Resource usage summary:
CPU time : 30812.76 sec.
Max Memory : 305 MB
Average Memory : 246.85 MB
Total Requested Memory : -
Delta Memory : -
Max Swap : -
Max Processes : 4
Max Threads : 5
Run time : 71348 sec.
Turnaround time : 71350 sec.
The output (if any) is above this job summary.
PS:
Read file <err> for stderr output of this job.
1#
/fsa/home/sjc_caiyz/samples/signal/s_400_2
0.1 6 29428.6 0.00575338
5 38605.1 0.00497448
4 39897.3 0.00503079
3 2.44907e+06 0.0104598
2 1.45091e+07 0.0149312
1 1.02264e+11 0.227218
2#
/fsa/home/sjc_caiyz/samples/signal/s_400_2
0.2 6 26241.6 0.00618847
5 38605.1 0.00497448
4 39897.3 0.00503079
3 2.44907e+06 0.0104598
2 1.45091e+07 0.0149312
1 1.02264e+11 0.227218
3#
/fsa/home/sjc_caiyz/samples/signal/s_400_2
0.3 6 21101.8 0.00677369
5 38605.1 0.00497448
4 39897.3 0.00503079
3 2.44907e+06 0.0104598
2 1.45091e+07 0.0149312
1 1.02264e+11 0.227218
4#
/fsa/home/sjc_caiyz/samples/signal/s_400_2
0.4 6 16045.2 0.00807941
5 38605.1 0.00497448
4 39897.3 0.00503079
3 2.44907e+06 0.0104598
2 1.45091e+07 0.0149312
1 1.02264e+11 0.227218