Changeset 5 in huonglan

Timestamp:

Aug 5, 2011, 6:20:27 PM (13 years ago)

Author:

huonglan

Message:

first full version of 2011 analysis

File:

• ANA2011/ana2011.C (modified) (4 diffs)

ANA2011/ana2011.C

@@ -740 +740 @@
 
       if (!overlap_elejet){
-        Particle q(v4j,2,1,0);
+        int flavor = -1;
+        if (!RealData)
+          flavor = abs(jet_flavor_truth_label->at(i));
+        Particle q(v4j,2,1,flavor);
         goodJets.push_back(q);
         IndGoodJets.push_back(i);
@@ -746 +749 @@
         // a small btag study
         int lab = 0;
-        int flavor = abs(jet_flavor_truth_label->at(i));
         if (flavor == 15)     lab = 1;
         else if (flavor == 4) lab = 2;
         else if (flavor == 5) lab = 3;
-        m_sv0[lab]->Fill(sv0);
-        m_ip3dsv1[lab]->Fill(jet_flavor_weight_Comb->at(i));
+        m_sv0[lab]->Fill(max(-19.99,min(sv0,39.99)));
+        float ip3d = jet_flavor_weight_Comb->at(i);
+        m_ip3dsv1[lab]->Fill(max(-19.99,min(ip3d,39.99)));
         //
         if (sv0 > 5.85){
@@ -1397 +1400 @@
   //##############         STAR CHI2       #################
   //########################################################
-  bool doChi2Ana = true;
+  bool doChi2Ana = false;
   if (!doChi2Ana) return kTRUE;
 
@@ -2753 +2756 @@
   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.;
+  b_ptjq1W=jq1W.Pt()/1000.;
+  b_ptjq2W=jq2W.Pt()/1000.;
+  b_ptjben=jben.Pt()/1000.;
   b_etajbqq=jbqq.Eta();
   b_etajq1W=jq1W.Eta();
   b_etajq2W=jq2W.Eta();
   b_etajben=jben.Eta();
+  b_phijbqq=jbqq.Phi();
+  b_phijq1W=jq1W.Phi();
+  b_phijq2W=jq2W.Phi();
+  b_phijben=jben.Phi();
+
+  b_ptLep=vlep.Pt()/1000.;
   b_etaLep=vlep.Eta();
   b_phiLep=vlep.Phi();
+  b_ptNeu=ptNeu/1000.;
   b_etaNeu=etaNeu;
   b_etaMin=etas.at(0);  
   b_etaMax=etas.at(5);
-  b_ptjbqq=jbqq.Pt()/1000.;
-  b_ptjq1W=jq1W.Pt()/1000.;
-  b_ptjq2W=jq2W.Pt()/1000.;
-  b_ptjben=jben.Pt()/1000.;
-  b_ptLep=vlep.Pt()/1000.;
-  b_ptNeu=ptNeu/1000.;
   b_dRqqW=dRqq_afterchi2;
   b_dRmin_eq=dRmineq;