Changeset 8 in huonglan

Timestamp:

Aug 22, 2011, 2:04:28 PM (13 years ago)

Author:

huonglan

Message:

22August : changing b-tagging & before mu in b-jet treat

Location:

ANA2011

Files:

• InitHistos.h (modified) (1 diff)
• MyAdd.h (modified) (2 diffs)
• SetBranchTreeAna.h (modified) (1 diff)
• ana2011.C (modified) (8 diffs)

ANA2011/InitHistos.h

@@ -26 +26 @@
   TString Match[2] = {"_match_allevents","_match_chi2"};
   for (int i = 0; i < 2; i++) {
-    name = "category"; name+= Match[i];
+    name = "category"; name += Match[i];
     m_category[i] = new TH1F(name, "", 50, 0, 50);
     fOutput->Add(m_category[i]);

ANA2011/MyAdd.h

@@ -35 +35 @@
 pair<unsigned int, vector<TLorentzVector> > matchQstat2Jets(std::vector<Particle>, std::vector<Particle>, bool);
 void TruthAna(double,vector<Particle>,vector<Particle>,bool,Particle,int,int,vector<TLorentzVector>,vector<TLorentzVector>,vector<TLorentzVector>,TLorentzVector,TLorentzVector,double,double,double);
+void FillFlavourChi2(vector<TLorentzVector>,vector<Particle>,TLorentzVector,TLorentzVector,TLorentzVector,TLorentzVector);
+
 void TagTop172(vector<Particle>,vector<TLorentzVector>,vector<TLorentzVector>,bool);
 bool emuOverlap(vector<int>);
@@ -239 +241 @@
   b_LeadingJetsInChi2,
   b_catAll,     
-  b_catChi2; 
+  b_catChi2,
+  b_flavjbqq,
+  b_flavjq1W,
+  b_flavjq2W,
+  b_flavjben,
+  b_NRealbChi2;
+
 Double_t b_dsigmaGG, b_dsigmaQQ, b_dsigmatot,
   b_chi2minmin, b_chi2second, b_chi2New,

ANA2011/SetBranchTreeAna.h

@@ -123 +123 @@
 
 tree_ana->Branch("catChi2",                  &b_catChi2,                  "catChi2/I");
+tree_ana->Branch("flavjbqq",                 &b_flavjbqq,             "flavjbqq/I");
+tree_ana->Branch("flavjq1W",                 &b_flavjq1W,             "flavjq1W/I");
+tree_ana->Branch("flavjq2W",                 &b_flavjq2W,             "flavjq2W/I");
+tree_ana->Branch("flavjben",                 &b_flavjben,             "flavjben/I");
+tree_ana->Branch("NRealbChi2",               &b_NRealbChi2,         "NRealbChi2/I");

ANA2011/ana2011.C

@@ -68 +68 @@
 #include "MatchTrigger.C"
 #include "ana2011AddFunctions.C"
+#include "FunctionsFlavourChi2.C"
 #include "emuOverlap.C"
 #include "SF.C"
@@ -755 +756 @@
         float ip3d = jet_flavor_weight_Comb->at(i);
         m_ip3dsv1[lab]->Fill(max(-19.99,min(ip3d,39.99)));
+        float combnn = jet_flavor_weight_JetFitterCOMBNN->at(i);
         //
-        if (sv0 > 5.85){
+        //if (sv0 > 5.85){
+        if (combnn > 2.0){
           hasbtagged = true;
           bjets.push_back(v4j);   
@@ -774 +777 @@
   double Ptmax = 1e-30;
   vector<TLorentzVector> vGoodJets;
+  vector<int> vGoodJetsFlav;
 
   unsigned int nGoodJets = goodJets.size();
@@ -782 +786 @@
     //for (unsigned int i = 0; i < 4; i++){
     Particle parj = goodJets.at(i);
+    vGoodJetsFlav.push_back(parj.POL());
     TLorentzVector v4 = parj.vector4();
     vGoodJets.push_back(v4);
@@ -1400 +1405 @@
   //##############         STAR CHI2       #################
   //########################################################
-  bool doChi2Ana = false;
+  bool doChi2Ana = true;
   if (!doChi2Ana) return kTRUE;
 
@@ -1769 +1774 @@
   FourJetsChi2.push_back(jq2W);
   FourJetsChi2.push_back(jben);
+
+  //Flavour Chi2
+  FillFlavourChi2(vGoodJets,goodJets,jbqq,jq1W,jq2W,jben);
 
   //Matrix Element
@@ -2756 +2764 @@
   b_pTtoplep=vtlminmin.Pt()/1000;
   b_pZtophad=vthminmin.Z()/1000;
-  b_pZttbar=vttpairchi2.Z()/1000
+  b_pZttbar=vttpairchi2.Z()/1000;
 
   b_ptjbqq=jbqq.Pt()/1000.;
@@ -3021 +3029 @@
     b_jet_flavor_weight_SV0->GetEntry(entry);
     b_jet_flavor_weight_Comb->GetEntry(entry);
+    b_jet_flavor_weight_JetFitterCOMBNN->GetEntry(entry);
     if (!RealData)
       b_jet_flavor_truth_label->GetEntry(entry);