Merge pull request #1 from mstoye/origin/master

made it simple with eos example

Author: mstoye

SUSYLooperHists.C

@@ -7,7 +7,6 @@
 
 #include <TF1.h>
 #include <TMath.h>
-#include <HTransformToCS.c>
 #include <TMatrixF.h>
 #include <TVectorF.h>
 
@@ -40,19 +39,24 @@ TFile* SUSYLooperHists::Loop()
 
    if (fChain == 0) return NULL;
    Long64_t nentries = fChain->GetEntriesFast();
-   HTransformToCS mtCS;
-   mtCS.SetCmsEnergy(8000);
-
-
+  
    TFile* myFile = new TFile(outputFileName,"recreate");
+   // checks the Mtt mass ++ simple plot example ++
    TH1D* LepTwoZmass= new TH1D("LepTwoZmass",";approximate Z(#tau#tau) mass [GeV];",100,0,2000);
+   LepTwoZmass->Sumw2();
+
+   // baseline plots +++++++++++++++++++++
+   // allow to fill cutflow to sychronize to others
+   TH1D* CutFlow = new TH1D("CutFlow",";CutFlow [unweighted];",20,-0.5,19.5);
+   CutFlow->Sumw2();
+
+ // a scan should be filled without weights before any cut to get the efficiency. The reason is that for a scan each points have different x-section, which are not accessable during the loop
    TH2D* scan = new TH2D("scan","scan",8,112.5,312.5,40,112.5,312.5);
+   // fill after some cuts the scans to get efficiency, i.e. ->Divide(scan) after the loop
    TH2D* scanA = new TH2D("scanA","scan",8,112.5,312.5,40,112.5,312.5);
    TH2D* scanB = new TH2D("scanB","scan",8,112.5,312.5,40,112.5,312.5);
    TH2D* scanC = new TH2D("scanC","scan",8,112.5,312.5,40,112.5,312.5);
-   // baseline plots +++++++++++++++++++++
-   TH1D* CutFlowAntonis = new TH1D("CutFlowAntonis",";P(l);",20,-0.5,19.5);
-   CutFlowAntonis->Sumw2();
+
 
    Long64_t nbytes = 0, nb = 0;
    for (Long64_t jentry=0; jentry<nentries;jentry++) {
@@ -61,7 +65,7 @@ TFile* SUSYLooperHists::Loop()
       if (ientry < 0) break;
       nb = fChain->GetEntry(jentry);   nbytes += nb;
 
-      if(jentry%10000==0) cout << "Event: "<<jentry <<" weight: " << weight  <<endl;
+      if(jentry%100000==0) cout << "Event: "<<jentry <<endl;
       scan->Fill(GenSusyMStop,GenSusyMNeutralino);
       // select lepton and put it into TLorentzvector 
 
@@ -109,28 +113,28 @@ TFile* SUSYLooperHists::Loop()
 
       // preselection syncronizd to Antonis ++++++++++++++++++++
       if( nLepGood == 0  ) continue;
-      CutFlowAntonis->Fill(0.);
+      CutFlow->Fill(0.);
       if(met.Pt()<200) continue;
       if(nLepGood<2)  continue;
       if(nJet==0)  continue;
       if(nTauGood!=0)  continue;
-      CutFlowAntonis->Fill(1.);
+      CutFlow->Fill(1.);
       if(fabs(LepGood_pdgId[0]* LepGood_pdgId[1])<122)   continue;
-      CutFlowAntonis->Fill(2.);
+      CutFlow->Fill(2.);
       if(lep.Pt()<5||fabs(lep.Eta())>1.5||lep.Pt()>60) continue;
-      CutFlowAntonis->Fill(3.);
+      CutFlow->Fill(3.);
       if(fabs(LepGood_pdgId[0])==11&&lep.Pt()<7) continue;
-      CutFlowAntonis->Fill(4.);
+      CutFlow->Fill(4.);
       // check that a cut was ahead in nLepton !!!!!!!!
       TLorentzVector secondLep;
       secondLep.SetPtEtaPhiM(LepGood_pt[1], LepGood_eta[1], LepGood_phi[1],0.1);
       if(secondLep.Pt()<3||fabs(secondLep.Eta())>1.5||secondLep.Pt()>60) continue;
-      CutFlowAntonis->Fill(5.);
+      CutFlow->Fill(5.);
       // lepton cuts
       if(LepGood_relIso[0]*lep.Pt()>10.) continue;
       if(LepGood_relIso[1]*secondLep.Pt()>5.) continue;
       if(LepGood_relIso[1]>.5) continue;
-      CutFlowAntonis->Fill(6.);
+      CutFlow->Fill(6.);
       // JET CUTS
       if(nJet>0){
 	if (Jet_pt[0] < 150) continue;
@@ -139,28 +143,30 @@ TFile* SUSYLooperHists::Loop()
       if(nJet>2){
 	if (Jet_pt[2] > 60) continue;
       }
-      CutFlowAntonis->Fill(7.);
+      CutFlow->Fill(7.);
       if(met.Pt()/HT30<2./3.) continue;
       if(nb40!=0) continue;    
-      CutFlowAntonis->Fill(9.);
-    // preselection syncronizd to Antonis ++++++++++++++++++++
+      CutFlow->Fill(9.);
+    // preselection syncronizd to  ++++++++++++++++++++
 
       if(MT>60&&MT<100) continue;
 
 
       float pairmass = DiTau_InvMass(met,lep,secondLep,0);     
       if(lep.Pt()>25||secondLep.Pt()>15) continue;
+ 
+      // define some over and underflow 
       if (pairmass>2000) pairmass=1999.;
       if (pairmass<0) pairmass=0.;
 
-      LepTwoZmass->Fill(pairmass,weight);
+      LepTwoZmass->Fill(pairmass,weight*puWeight);
 
 
 	if((LepGood_pdgId[0]* LepGood_pdgId[1]<-121))
 	  {
 
 	    if(pairmass>160||pairmass<20){
-	      CutFlowAntonis->Fill(10., weight*puWeight);
+	      CutFlow->Fill(10., weight*puWeight);
 	      scanB->Fill(GenSusyMStop,GenSusyMNeutralino);
 
 	      if(nb30L==0)
@@ -191,7 +197,7 @@ TFile* SUSYLooperHists::Loop()
    scanB->Divide(scan);
    scanC->Divide(scan);
 
-   CutFlowAntonis->DrawCopy();
+   CutFlow->DrawCopy();
 
    myFile->Write();
    return (myFile);

SUSYLooperHists.h

@@ -262,15 +262,14 @@ public :
 SUSYLooperHists::SUSYLooperHists(SampleInfo mySample, float Mstop_ , float MLSP_ )
 {
 
-  TFile *f = new TFile(mySample.SamplePath);
+  TFile *f = TFile::Open(mySample.FilePath);
   TTree* tree = (TTree*) f->Get("markusTreeProducer"); 
-  outputFileName = mySample.SampleName+".root";
+  outputFileName = mySample.OutputFileNameTag+".root";
   weight = mySample.weight();
-  cout << "weight again " << weight <<endl;;
   Mstop= Mstop_;
   MLSP=MLSP_;
-  maxEvents = mySample.TotalEvNum;
   Init(tree);
+  cout << "looping over: \""<< mySample.FilePath<< "\" to put plots into: \""<< outputFileName<<"\""<<endl;;
 
 }
 
@@ -443,6 +442,7 @@ Int_t SUSYLooperHists::Cut(Long64_t entry)
 {
 // This function may be called from Loop.
 // returns  1 if entry is accepted.
+  cout << " can be used to put preselection cuts for " << entry<<endl;
 // returns -1 otherwise.
    return 1;
 }

SampleInfo.h

@@ -0,0 +1,36 @@
+#ifndef SampleInfo_h
+#define SampleInfo_h
+
+#include "TString.h"
+
+class SampleInfo
+{
+
+ public:
+  SampleInfo( float Lumi_,float NEvents_,float CrossSection_,TString FilePath_, TString LegendLabel_,TString OutputFileNameTag_,int color_): Lumi(Lumi_),NEvents(NEvents_), CrossSection(CrossSection_),FilePath(FilePath_),LegendLabel(LegendLabel_),OutputFileNameTag(OutputFileNameTag_), color(color_)
+{;}
+    ;
+  float weight();
+
+
+ private:
+  float Lumi;
+  float NEvents; // number of events before any skimmming in MC
+  float CrossSection; // put -1 for data
+  TString FilePath; 
+  TString LegendLabel;
+  TString OutputFileNameTag;
+  int color;
+
+};
+
+
+float SampleInfo::weight()
+{
+  if(CrossSection==-1) return 1.; // for data
+  else
+    {
+      return (CrossSection/NEvents*Lumi);
+    }
+};
+#endif 

myplots.C

@@ -4,27 +4,41 @@
 
 void myplots()
 {
+  // all files are at eos, so start with these
+  SampleInfo eosZll(40.5,2655795,19.7,TString("root://eoscms.cern.ch//eos/cms/store/cmst3/group/susy/markus/preliminary_samples_11thApril_2014/DYJetsToLLPtZ100/markusTreeProducer/markusTreeProducer_tree.root"),TString("Zll"),TString("Z(ll)"),5);
 
+  // You can copy them over to your local directory to gain speed
+  SampleInfo Zll(40.5,2655795,19.7,TString("../data_test_antonios/DYJetsToLLPtZ100/markusTreeProducer/markusTreeProducer_tree.root"),TString("Zll"),TString("Z(ll)"),5);;
 
-  // 
+  // pathes below to local directories 
   SampleInfo wjet(282,12742383,19.7 ,TString("../data_test_antonios/WPT100JETS/markusTreeProducer/markusTreeProducer_tree.root"),TString("W_PT_100"),TString("W PT 100"),2);
   SampleInfo ttsemi(106.15,320./963.*24895259,19.7,TString("../data_test_antonios/myhaddedFileTsemi.root"),TString("TTbarSemiLep"),TString("tt(l)"),4);
-  SampleInfo ttlep(25.8,11947326*0.467,19.7,TString("../data_test_antonios/myhaddedFileTlep.root"),TString("TTbarFullyLep"),TString("tt(l)"),5);
-  SampleInfo Zll(40.5,2655795,19.7,TString("../data_test_antonios/DYJetsToLLPtZ100/markusTreeProducer/markusTreeProducer_tree.root"),TString("Zll"),TString("Z(ll)"),5);
+  SampleInfo ttlep(25.8,11947326*0.467,19.7,TString("../data_test_antonios/myhaddedFileTlep.root"),TString("TTbarFullyLep"),TString("tt(ll)"),5);
   SampleInfo Znunu50(435.71,4040980,19.7,TString("../data_test_antonios/ZJetsToNuNu100/markusTreeProducer/markusTreeProducer_tree.root"),TString("Znunu50"),TString("Z(#nu#nu)"),5);
   SampleInfo Znunu100(186.84,4416646,19.7,TString("../data_test_antonios/ZJetsToNuNu100/markusTreeProducer/markusTreeProducer_tree.root"),TString("Znunu100"),TString("Z(#nu#nu)"),5);
   SampleInfo Znunu200(45.60,5055885,19.7,TString("../data_test_antonios/ZJetsToNuNu100/markusTreeProducer/markusTreeProducer_tree.root"),TString("Znunu200"),TString("Z(#nu#nu)"),5);
   SampleInfo Znunu400(6.25,1006928,19.7,TString("../data_test_antonios/ZJetsToNuNu100/markusTreeProducer/markusTreeProducer_tree.root"),TString("Znunu400"),TString("Z(#nu#nu)"),5);
   SampleInfo ZZ4mu(6.25,1006928,19.7,TString("/afs/cern.ch/work/m/mstoye/CMGTools_SUSY/CMSSW_5_3_14/src/CMGTools/TTHAnalysis/cfg/Zz4mu/ZZTo4mu/markusTreeProducer/markusTreeProducer_tree.root"),TString("ZZ4m"),TString("Z(#nu#nu)"),5);
-  SampleInfo ZZ(6.25,1006928,19.7,TString("ZZpythia.root"),TString("ZZ4m"),TString("Z(#nu#nu)"),5);
+  //  SampleInfo ZZ(6.25,1006928,19.7,TString("ZZpythia.root"),TString("ZZ4m"),TString("Z(#nu#nu)"),5);
 
 
   vector <TFile*> Files;
   vector <SampleInfo> Names;
-  SUSYLooperHists myZll(Zll);
+
+  SUSYLooperHists myZll(eosZll);
   TFile* myZllfile = myZll.Loop();
-  Files.push_back(myZll);
+  Files.push_back(myZllfile);
+  Names.push_back(eosZll);
+
+  cout << endl<< "see the difference local vs. eos" <<endl <<endl;
+
+  SUSYLooperHists myLocalZll(Zll);
+  TFile* myLocalZllfile = myLocalZll.Loop();
+  Files.push_back(myLocalZllfile);
   Names.push_back(Zll);
+
+  cout <<" further samples "<<endl;
+
   SUSYLooperHists mytt2l(ttlep);
   TFile* myttlep = mytt2l.Loop();
   Files.push_back(myttlep);