Changeset 18 in huonglan

Timestamp:

Aug 29, 2011, 12:59:40 AM (13 years ago)

Author:

huonglan

Message:

adding muon to b-jet

Location:

ANA2011

Files:

• CutFlow.C (modified) (1 diff)
• ana2011.C (modified) (13 diffs)

ANA2011/CutFlow.C

@@ -1 +1 @@
 {
-  TFile *f = TFile::Open("res.root");
+  TFile *f = TFile::Open("./CUTFLOW/res180710Muons_satisfyCutFlow.root");
 
   TString s[6] = {"","E","M","Sig","ESig","MSig"};

ANA2011/ana2011.C

@@ -723 +723 @@
 
   //Before selecting jets, look if there is a muon coming from jet
+  cout << "1" << endl;
   std::vector<std::pair<TLorentzVector,TLorentzVector> > MuonJetPairs = AddMuonJet(RealData,m_smearMC,indGoodMuon);
-  //cout << "Number of pair muon and jet that overlap " << MuonJetPairs.size() << endl;
+  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);
-  }
+  if (MuonJetPairs.size() > 0){
+    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);
+    }
+  }
+
+  cout << "MuonJetPairs.size() " << MuonJetPairs.size() << endl;
 
   //SELECTION OF JETS
@@ -747 +752 @@
       for (unsigned int iJ = 0; iJ < vJet.size(); iJ++){
         TLorentzVector vJ = vJet.at(iJ);
-        if (v4j==vJ) v4j += vMuon.at(iJ);
+        if (v4j==vJ){
+          cout << "jet_pt before adding " << jet_pt->at(i) << endl;
+          cout << "vJ pt " << vJ.Pt() << endl;
+          cout << "vM pt " << vMuon.at(iJ).Pt() << endl;
+          v4j += vMuon.at(iJ);
+        }
       }
+      jet_pt->at(i) = v4j.Pt();
+      jet_E->at(i)  = v4j.E();
+
+      //cout << "jet_pt after adding " << jet_pt->at(i) << endl;
     }
 
-    jet_pt->at(i) = v4j.Pt();
-    jet_E->at(i)  = v4j.E();
+    cout << "before GoodJets" << endl;
 
     if(GoodJets(i)) {
@@ -798 +811 @@
   m_weiSV0Max[0]->Fill(sv0Max,weight);
 
+  cout << "2" << endl;
+
   bool askbtagged = false;
 
@@ -831 +846 @@
   }
   HT += topMET_et;
+
+  cout << "3" << endl;
 
   //pt leading jet : C15 : add for u4
@@ -1146 +1163 @@
   // #### The Naive analysis : no fit, test all combinations,         ###########
   // #### take the one minimizing the difference between the two tops ###########  
+
+  cout << "begin first analysis" << endl;
 
   //###### RECONSTRUCT the LEPTONIC W ######
@@ -1435 +1454 @@
   if (!doChi2Ana) return kTRUE;
 
+  cout << "begin chi2" << endl;
+
   //Here I have to add this condition cause my chi2 method is happy 
   //just in the case I have one Lepton (e or mu)
@@ -1476 +1497 @@
     vector<TLorentzVector> vecpars2;
     vector<int> index = indjets.at(id);
-    //cout << "index size : " << index.size() << endl;
+    cout << "index size : " << index.size() << endl;
 
     vector<TLorentzVector> vJetsUseForChi2;
@@ -1484 +1505 @@
     }
     
+    cout << "vJetsUseForChi2 " << vJetsUseForChi2.size() << endl;
+
     //in case we ask for b-tagging
     //we have to force the bjets to be in the jets used to calculate chi2 first
@@ -1506 +1529 @@
       if (!bjetusedforchi2) continue;
 
+      cout << "before case 3 b-jet in 4 good jets" << endl;
+
       if (vGoodJets.size()==4 && bjetsInd.size()==3){
+        cout << "threebjets4GoodJets" << endl;
+
         threebjets4GoodJets = true;
         TVector2 IndSV0;
@@ -1513 +1540 @@
           vj1.SetPtEtaPhiM(jet_pt->at(inew),jet_eta->at(inew),jet_phi->at(inew),jet_m->at(inew));
           double SV0j = jet_flavor_weight_SV0->at(inew);
+          if (vJet.size()>0){
+            cout << "back to add muon prob" << endl;
+            for (unsigned int iJ = 0; iJ < vJet.size(); iJ++){
+              TLorentzVector vJ = vJet.at(iJ);
+              if (vj1==vJ){
+                cout << "find one bjet" << endl;
+                vj1 += vMuon.at(iJ);
+              }
+            }
+          }
           for (unsigned int jnew = 0; jnew < vJetsUseForChi2.size(); jnew++){
             TLorentzVector vj2 = vJetsUseForChi2.at(jnew);
@@ -1522 +1559 @@
           }
         }
+        
+        cout << "now here?" << endl;
+
         sort(IndSV0jets34.begin(),IndSV0jets34.end(),sortdecreaseSV0);
+        cout << "IndSV0jets34 size " << IndSV0jets34.size() << endl;
+
         bjetsInd34.push_back((int)IndSV0jets34.at(0).X());
+        cout << "IndSV0jets34.at(0).X() " << IndSV0jets34.at(0).X() << endl;
+
         bjetsInd34.push_back((int)IndSV0jets34.at(1).X());
+        cout << "IndSV0jets34.at(1).X() " << IndSV0jets34.at(1).X() << endl;
       }
     }//end of (askbtagged)
     
-    /*for (unsigned int j=0; j < index.size(); j++){
-      if (PRINT) cout << "index : " << index.at(j) << endl;
-      TLorentzVector vp4 = vGoodJets.at(index.at(j));
-      if (isMuonJet && vp4 == vJetOverlapNotGoodMuon) 
-        vp4 += vMuonJet;
-      vecpars2.push_back(vp4);
-    }
-    vecpars2.push_back(goodLepton.at(0).vector4());
-    */ //comment for 26 August
-    
+    cout << "come here ? " << endl;
+
     for (unsigned int j=0; j < index.size(); j++){
       if (PRINT) cout << "index : " << index.at(j) << endl;
@@ -1544 +1581 @@
     }
     vecpars2.push_back(goodLepton.at(0).vector4());   
+
+    cout << "vecpars2 size " << vecpars2.size() << endl;
 
     //##############################################################
@@ -3326 +3365 @@
 bool ana2011::GoodJets(int ij) 
 {
+  //cout << "jet_pt in GoodJet loop " << jet_pt->at(ij) << endl;
 
   // Acceptance (remove crack also):