Changeset 21 in huonglan

Timestamp:

Sep 1, 2011, 9:53:56 AM (13 years ago)

Author:

huonglan

Message:

Add efficiency test folder

Location:

SPLOT/NEWREALEASE

Files:

• EFFICIENCY (added)
• EFFICIENCY/Cat1.eps (added)
• EFFICIENCY/Cat2.eps (added)
• EFFICIENCY/Cat3.eps (added)
• EFFICIENCY/Cat4.eps (added)
• EFFICIENCY/Cat5.eps (added)
• EFFICIENCY/Cat6.eps (added)
• EFFICIENCY/CatEPS.tar (added)
• EFFICIENCY/EfficiencyDependOnCut.C (added)
• EFFICIENCY/EfficiencyDependOnCut_C.d (added)
• EFFICIENCY/EfficiencyDependOnCut_saved.C (added)
• EFFICIENCY/goC.C (added)
• MERGEFILES/MakeFakeData.C (modified) (9 diffs)

SPLOT/NEWREALEASE/MERGEFILES/MakeFakeData.C

@@ -39 +39 @@
 
   TString s3 = "v3.root";
-  TString ver = "NEW";
+  TString ver = "ToTestEfficiency";
 
   //set the name of 78 files
@@ -141 +141 @@
 
   //###################################
+  //************ File Data  ***********
+  //###################################
+
+  TFile *fData    = new TFile("res_Data"+ver+".root","RECREATE");
+  TTree *tWeiD   = new TTree("T","T");
+  float b_weiD  = -100;
+  float b_mtophad= -100;    
+  float b_HT= -100;         
+  float b_cos1= -100;       
+  float b_cos2= -100;       
+  float b_MinLep1= -100;    
+  short b_hfor= -100;       
+  short b_hasBadJet= -100;  
+  short b_isEleEvent= -100; 
+  int b_RunNumber= -100;  
+  tWeiD->Branch("wei",        &b_weiD,               "wei/F");
+  tWeiD->Branch("mtophad",    &b_mtophad,        "mtophad/F");
+  tWeiD->Branch("HT",         &b_HT,                  "HT/F");
+  tWeiD->Branch("cos1",       &b_cos1,              "cos1/F");
+  tWeiD->Branch("cos2",       &b_cos2,              "cos2/F");
+  tWeiD->Branch("MinLep1",    &b_MinLep1,        "MinLep1/F");
+  tWeiD->Branch("hfor",       &b_hfor,              "hfor/S");
+  tWeiD->Branch("hasBadJet",  &b_hasBadJet,    "hasBadJet/S");
+  tWeiD->Branch("isEleEvent", &b_isEleEvent,  "isEleEvent/S");
+  tWeiD->Branch("RunNumber",  &b_RunNumber,    "RunNumber/I");
+  gROOT->cd();
+  
+  int countD = 0;
+  for (int ifi = 0; ifi < tData->GetEntries(); ifi++){
+    countD++;
+    tData->GetEntry(ifi);
+    b_mtophad = tSig->GetLeaf("mtophad")->GetValue();
+    b_HT = tSig->GetLeaf("HT")->GetValue();
+
+    //if (b_mtophad<250 || b_HT<620) continue; 
+    if (ifi%1000==0) cout << "Reading event " << ifi << " from Signal tree" << endl;
+
+    //TString currentFile = tSig->GetTree()->GetCurrentFile()->GetName();
+
+    b_weiS = Wei[0]*tSig->GetLeaf("WeiEvent")->GetValue();
+    b_cos1 = tSig->GetLeaf("cos_tlepCMtpair")->GetValue();
+    b_cos2 = tSig->GetLeaf("cosWqq")->GetValue();
+    b_MinLep1 = tSig->GetLeaf("MinLep1")->GetValue();           
+    b_hfor = (short)tSig->GetLeaf("hfor")->GetValue();
+    b_isEleEvent = (short)tSig->GetLeaf("isEleEvent")->GetValue();
+    b_hasBadJet = (short)tSig->GetLeaf("hasBadJet")->GetValue();
+    b_RunNumber = (int)tSig->GetLeaf("RunNumber_ana")->GetValue();
+
+    tWeiS->Fill();
+  }
+
+  fSig->cd();
+  tWeiS->Write();
+  fSig->Close();
+
+
+  //###################################
   //********** File Signal  ***********
   //###################################
@@ -175 +232 @@
     b_HT = tSig->GetLeaf("HT")->GetValue();
 
-    if (b_mtophad<250 || b_HT<620) continue; 
+    //if (b_mtophad<250 || b_HT<620) continue; 
     if (ifi%1000==0) cout << "Reading event " << ifi << " from Signal tree" << endl;
 
@@ -235 +292 @@
     b_HT = ht;
 
-    if (mt<250 || ht<620) continue; 
+    //if (mt<250 || ht<620) continue; 
     if (ifi%1000==0) cout << "Reading event " << ifi << " from AllBkg tree" << endl;
 
@@ -301 +358 @@
     b_HT = ht;
 
-    if (mt<250 || ht<620) continue; 
+    //if (mt<250 || ht<620) continue; 
     if (ifi%1000==0) cout << "Reading event " << ifi << " from Top tree" << endl;
 
@@ -350 +407 @@
     b_HT = ht;
 
-    if (mt<250 || ht<620) continue; 
+    //if (mt<250 || ht<620) continue; 
     if (ifi%1000==0) cout << "Reading event " << ifi << " from AllW tree" << endl;
 
@@ -402 +459 @@
     b_HT = ht;
 
-    if (mt<250 || ht<620) continue; 
+    //if (mt<250 || ht<620) continue; 
     if (ifi%1000==0) cout << "Reading event " << ifi << " from SingT tree" << endl;
 
@@ -451 +508 @@
     b_HT = ht;
 
-    if (mt<250 || ht<620) continue; 
+    //if (mt<250 || ht<620) continue; 
     if (ifi%1000==0) cout << "Reading event " << ifi << " from AllZ tree" << endl;
 
@@ -503 +560 @@
     b_HT = ht;
 
-    if (mt<250 || ht<620) continue; 
+    //if (mt<250 || ht<620) continue; 
     if (ifi%1000==0) cout << "Reading event " << ifi << " from QCD tree" << endl;