Changeset 16 in huonglan

Timestamp:

Aug 27, 2011, 12:20:55 AM (13 years ago)

Author:

huonglan

Message:

Add a muon to jet - first try

Location:

ANA2011

Files:

• AddMuonJetFunction.C (added)
• MyAdd.h (modified) (1 diff)
• ana2011.C (modified) (9 diffs)

ANA2011/MyAdd.h

@@ -10 +10 @@
 //bool GoodElectron(int,bool);
 bool GoodMuon(int);
+bool GoodMuonJet(int,bool,bool,vector<int>);//this function will be used to choose muon and add to jets
+std::vector<std::pair<TLorentzVector,TLorentzVector> > AddMuonJet(bool,bool,std::vector<int>);
 bool GoodJets(int);
 bool MatchTrigger(int,vector<TLorentzVector>,vector<TLorentzVector>);

ANA2011/ana2011.C

@@ -71 +71 @@
 #include "emuOverlap.C"
 #include "SF.C"
+#include "AddMuonJetFunction.C"
 
 bool sortdecreaseSV0(TVector2 v1, TVector2 v2)
@@ -504 +505 @@
 
   //###### SELECTION OF MUON ######
+  vector<int> indGoodMuon;
   for (unsigned int i = 0; i < mu_pt->size(); i++) {
     if(GoodMuon(i)) {
+      indGoodMuon.push_back(i);
       TLorentzVector v4m;
 
@@ -718 +721 @@
   m_loose->Fill(10,weight);
   if (ncate==12||ncate==13||ncate==21||ncate==31) m_looseSig->Fill(10,weight);
+
+  //Before selecting jets, look if there is a muon coming from jet
+  std::vector<std::pair<TLorentzVector,TLorentzVector> > MuonJetPairs = AddMuonJet(RealData,m_smearMC,indGoodMuon);
+  //cout << "Number of pair muon and jet that overlap " << MuonJetPairs.size() << endl;
+
+  vector<TLorentzVector> vMuon, vJet;
+  for (unsigned int i = 0; i < MuonJetPairs.size(); i++){
+    pair<TLorentzVector,TLorentzVector> MuonJet = MuonJetPairs.at(i);
+    vMuon.push_back(MuonJet.first);
+    vJet.push_back(MuonJet.second);
+  }
 
   //SELECTION OF JETS
@@ -729 +743 @@
     TLorentzVector v4j;
     v4j.SetPtEtaPhiM(jet_pt->at(i),jet_eta->at(i),jet_phi->at(i),jet_m->at(i));
+
+    if (vJet.size()>0){
+      for (unsigned int iJ = 0; iJ < vJet.size(); iJ++){
+        TLorentzVector vJ = vJet.at(iJ);
+        if (v4j==vJ) v4j += vMuon.at(iJ);
+      }
+    }
+
+    jet_pt->at(i) = v4j.Pt();
+    jet_E->at(i)  = v4j.E();
+
     if(GoodJets(i)) {
       double sv0 = jet_flavor_weight_SV0->at(i);
@@ -1069 +1094 @@
     m_mjj3->Fill(mjj3,weight);  
 
+  /*
   bool isMuonJet = false;
   TLorentzVector vMuonJet;
@@ -1115 +1141 @@
     }
   }
+  */ // comment for 26 August
 
   // #####################    START THE FIRST ANALYSIS    ####################### 
@@ -1479 +1506 @@
       if (!bjetusedforchi2) continue;
 
-      //cout << "enter here???" << endl;
-      
-      //cout << "vGoodJets size = " << vGoodJets.size() << endl;
-      //cout << "bjetsInd size = " << bjetsInd.size() << endl;
-      //cout << "index of bjets in vGoodJets : ";
-      /*for (unsigned int i = 0; i < bjetsInd.size(); i++) {
-        cout << bjetsInd.at(i) << " " ;
-        }*/
-      //cout << endl;
-
       if (vGoodJets.size()==4 && bjetsInd.size()==3){
         threebjets4GoodJets = true;
@@ -1511 +1528 @@
     }//end of (askbtagged)
     
-    for (unsigned int j=0; j < index.size(); j++){
+    /*for (unsigned int j=0; j < index.size(); j++){
       if (PRINT) cout << "index : " << index.at(j) << endl;
       TLorentzVector vp4 = vGoodJets.at(index.at(j));
@@ -1519 +1536 @@
     }
     vecpars2.push_back(goodLepton.at(0).vector4());
+    */ //comment for 26 August
+    
+    for (unsigned int j=0; j < index.size(); j++){
+      if (PRINT) cout << "index : " << index.at(j) << endl;
+      TLorentzVector vp4 = vGoodJets.at(index.at(j));
+      vecpars2.push_back(vp4);
+    }
+    vecpars2.push_back(goodLepton.at(0).vector4());   
 
     //##############################################################