Merge pull request #9 from mstoye/origin/master

added conmtrol region

Author: mstoye

MuonTrigControl.C

@@ -0,0 +1,75 @@
+#define MuonTrigControl_cxx
+#include "MuonTrigControl.h"
+#include <TH2.h>
+#include <TStyle.h>
+#include <TCanvas.h>
+#include  <vector>
+#include <TF1.h>
+#include <TMath.h>
+
+
+TFile* MuonTrigControl::Loop()
+{
+  TFile* myFile = new TFile(outputFileName,"recreate");
+  if (fChain == 0) return myFile;
+  Long64_t nentries = fChain->GetEntriesFast();
+
+  TH1D *h_CutFlow  =new TH1D("h_CutFlow","h_CutFlow", 25, -0.5, 24.5 ); 
+  h_CutFlow->Sumw2();
+ 
+  
+  Long64_t nbytes = 0, nb = 0;
+  for (Long64_t jentry=0; jentry<nentries;jentry++) {
+    
+    Long64_t ientry = LoadTree(jentry);
+    if (ientry < 0) break;
+    nb = fChain->GetEntry(jentry);   nbytes += nb;
+
+    float HT30 =0.;
+    int Nj30=0;
+    int Nb=0;
+   
+    for(int i=0;i<nJet;i++){//Loop Over Selected Jets, measure Nj, Nb, HT;
+      if(Jet_pt[i]>30.&& fabs(Jet_eta[i])<4.5 ) {
+	Nj30++;
+	HT30=HT30+Jet_pt[i]; 
+	if(Jet_btagCSV[i]>0.679/*0.244*/)Nb++;  //CSV[M] b-tagging
+      } 
+    }
+    if( Nb>0                )continue;  h_CutFlow->Fill(11);
+
+
+    
+     if( MLSP!=0 ){//-----(Cuts for scans to select point)---------- 
+      
+      if( (GenSusyMStop<=(Mstop-12)) || (GenSusyMStop>=(Mstop+12)) )     continue; //<---[M_stop, M_LSP: set by by the last 2 arguments of the SampleInfo objects]
+      if( (GenSusyMNeutralino<=(MLSP - 4)) || (GenSusyMNeutralino>=(MLSP  +4)) )     continue; 
+    }//-------[ CUTS for Non-Scan samples ]-------------------------
+    if( nJet==0 &&  nLepGood<2                          )continue;  h_CutFlow->Fill(1);
+    if( fabs(LepGood_pdgId[0]*LepGood_pdgId[1])!=121  )continue;  h_CutFlow->Fill(2);
+    if( LepGood_pt[0]<5 || fabs(LepGood_eta[0])>1.5   )continue;  h_CutFlow->Fill(3);
+    if( LepGood_pt[1]<3 || fabs(LepGood_eta[1])>1.5   )continue;  h_CutFlow->Fill(4);
+    if( fabs(LepGood_pdgId[1])==11 && LepGood_pt[1]<5 )continue;  h_CutFlow->Fill(5); //Threashold for electrons.
+    if( fabs(LepGood_pdgId[0])==11 && LepGood_pt[0]<7 )continue;  h_CutFlow->Fill(6); //Threashold for electrons.
+    if( LepGood_relIso[1]>.5 || LepGood_relIso[1]*LepGood_pt[1]>5 )continue;  h_CutFlow->Fill(7); //2nd lep: *Universal* iso requirements!
+    if( fabs(LepGood_dz[0])>0.02 || fabs(LepGood_dxy[0])>0.02 )    continue; 
+    if( fabs(LepGood_dz[1])>0.02 || fabs(LepGood_dxy[1])>0.02 )    continue;  h_CutFlow->Fill(8);
+    if( Jet_pt[0]<150||fabs(Jet_eta[0])>2.4           )continue;  h_CutFlow->Fill(9); //1st Jet requirements.
+    if( nJet>2 && Jet_pt[2]>60                          )continue;  h_CutFlow->Fill(10);//2nd Jet requirements.
+    for(int i=0;i<nJet;i++){//Loop Over Selected Jets, measure Nj, Nb, HT;
+      if(Jet_pt[i]>30.&& fabs(Jet_eta[i])<4.5 ) {
+	Nj30++;
+	HT30=HT30+Jet_pt[i]; 
+	if(Jet_btagCSV[i]>0.679/*0.244*/)Nb++;  //CSV[M] b-tagging
+      } 
+    }
+    if( Nb>0                )continue;  h_CutFlow->Fill(11);
+    
+  }
+  
+  myFile->Write();
+  myFile->Close();
+  
+  return (myFile);
+      // if (Cut(ientry) < 0) continue;
+}

MuonTrigControl.h

@@ -0,0 +1,13 @@
+#ifndef MuonTrigControl_h
+#define MuonTrigControl_h
+
+#include "SUSYLooperHistsSoftBase.h"
+#include "SampleInfo.h" 
+
+class MuonTrigControl : public SUSYLooperHistsSoftBase
+{
+   public:
+  MuonTrigControl(SampleInfo mySample, float Mstop_ = 0, float MLSP_ = 0) {setTree(mySample, Mstop_, MLSP_);};
+  TFile*   Loop();
+};
+#endif 

SUSYLooperHistsSoftBase.C

@@ -0,0 +1,67 @@
+#define SUSYLooperHistsSoftBase_cxx
+#include "SUSYLooperHistsSoftBase.h"
+#include <TH2.h>
+#include <TStyle.h>
+#include <TCanvas.h>
+#include  <vector>
+
+#include <TF1.h>
+#include <TMath.h>
+#include <TMatrixF.h>
+#include <TVectorF.h>
+
+TFile* SUSYLooperHistsSoftBase::Loop()
+{
+  TFile* myFile = new TFile(outputFileName,"recreate");
+   if (fChain == 0) return myFile;
+   Long64_t nentries = fChain->GetEntriesFast();
+  
+ 
+  
+   Long64_t nbytes = 0, nb = 0;
+   for (Long64_t jentry=0; jentry<nentries;jentry++) {
+
+      Long64_t ientry = LoadTree(jentry);
+      if (ientry < 0) break;
+      nb = fChain->GetEntry(jentry);   nbytes += nb;
+
+    
+   }
+
+
+  
+   return (myFile);
+      // if (Cut(ientry) < 0) continue;
+}
+
+
+double SUSYLooperHistsSoftBase::DiTau_InvMass( TLorentzVector Met, TLorentzVector L1, TLorentzVector L2, float Al ){   
+  TLorentzVector T1,T2;
+  double DiTauMass; 
+  TMatrixF A(2,2); 
+  TVectorF C(2),X(2); 
+  A(0,0)=L1.Px();
+  A(0,1)=L2.Px();
+  A(1,0)=L1.Py();
+  A(1,1)=L2.Py();
+  A=A.Invert();
+  C(0)=(Met+L1+L2).Px();
+  C(1)=(Met+L1+L2).Py();
+  X=A*C;// double X0i=X(0), X1i=X(1);
+  //---------------[ MET ReAlignement subsection ]------------------------------
+  if(X(0)<0||X(1)<0){ 
+    if     ( fabs(L1.DeltaPhi(Met))>Al && fabs(L2.DeltaPhi(Met))>Al                                                  ) {/*DO NOTHING just normaly a non-Z event!*/}
+    else if( fabs(L1.DeltaPhi(Met))<Al && fabs(L2.DeltaPhi(Met))>Al                                                  ) Met.SetPtEtaPhiM(Met.Pt(),0,L1.Phi(),0);
+    else if( fabs(L2.DeltaPhi(Met))<Al && fabs(L1.DeltaPhi(Met))>Al                                                  ) Met.SetPtEtaPhiM(Met.Pt(),0,L2.Phi(),0);
+    else if( fabs(L1.DeltaPhi(Met))<Al && fabs(L2.DeltaPhi(Met))<Al && fabs(L1.DeltaPhi(Met))<fabs(L2.DeltaPhi(Met)) ) Met.SetPtEtaPhiM(Met.Pt(),0,L1.Phi(),0);
+    else if( fabs(L1.DeltaPhi(Met))<Al && fabs(L2.DeltaPhi(Met))<Al && fabs(L1.DeltaPhi(Met))>fabs(L2.DeltaPhi(Met)) ) Met.SetPtEtaPhiM(Met.Pt(),0,L2.Phi(),0);
+  }//---------------------------------------------------------------------------
+  C(0)=(Met+L1+L2).Px(); C(1)=(Met+L1+L2).Py(); X=A*C;
+  T1.SetPxPyPzE( L1.Px()*X(0), L1.Py()*X(0), L1.Pz()*X(0), sqrt( 3.1571 +L1.P()*L1.P()*X(0)*X(0) ) );
+  T2.SetPxPyPzE( L2.Px()*X(1), L2.Py()*X(1), L2.Pz()*X(1), sqrt( 3.1571 +L2.P()*L2.P()*X(1)*X(1) ) );
+  if( X(0)>0 && X(1)>0 ) DiTauMass=(T1+T2).M();  else DiTauMass=-(T1+T2).M();  return DiTauMass;
+  //if((X(0)!=X0i||X(1)!=X1i))std::cout<<X(0)<<" "<<X(1)<<" <--"<<X0i<<" "<<X1i<<" RMETal.phi="<<(T1+T2-L1-L2).Phi()<<" RMETal.eta"<<(T1+T2-L1-L2).Eta()<<" MZ="<<DiTauMass<<endl; 
+}
+
+ 
+