Changeset fe3908f in git

Timestamp:

Jun 27, 2016, 4:40:50 PM (8 years ago)

Author:

Alexandre Mertens <alexandre.mertens@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

a3a556c

Parents:

734b267

Message:

working with calos

Files:

• cards/delphes_card_CMS.tcl (modified) (4 diffs)
• examples/Validation.C (modified) (16 diffs)

cards/delphes_card_CMS.tcl

@@ -35 +35 @@
   
   FastJetFinder
+  FastCaloJetFinder
 
   JetEnergyScale
@@ -210 +211 @@
   set HCalEnergyMin 1.0
 
-  set ECalEnergySignificanceMin 1.0
-  set HCalEnergySignificanceMin 1.0
+  set ECalEnergySignificanceMin 2.0
+  set HCalEnergySignificanceMin 2.0
 
   set SmearTowerCenter true
@@ -497 +498 @@
 }
 
+################
+# caloJet finder
+################
+
+module FastJetFinder FastCaloJetFinder {
+  set InputArray Calorimeter/towers
+
+  set OutputArray jets
+
+  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
+  set JetAlgorithm 6
+  set ParameterR 0.5
+
+  set JetPTMin 20.0
+}
+
 ##################
 # Jet Energy Scale
@@ -610 +627 @@
 
   add Branch UniqueObjectFinder/jets Jet Jet
+  add Branch FastCaloJetFinder/jets CaloJet Jet
   add Branch UniqueObjectFinder/electrons Electron Electron
   add Branch UniqueObjectFinder/photons Photon Photon

examples/Validation.C

@@ -57 +57 @@
 double ptrangemin = 10;
 double ptrangemax = 10000;
-static const int Nbins = 10;
+static const int Nbins = 20;
 
 int objStyle = 1;
@@ -101 +101 @@
     obj = object;
 
-    cenResolHist = new TH1D(obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_cen", obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_cen", 500, -0.4, 0.4);
+    cenResolHist = new TH1D(obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_cen", obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_cen", 500, -1, 1);
     fwdResolHist = new TH1D(obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_fwd", obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_fwd", 500, 0.4, 0.4);
 
@@ -317 +317 @@
   Long64_t allEntries = treeReader->GetEntries();
 
-cout << "** Computing resolution of " << branchJet->GetName() << " induced by " << branchGenJet->GetName() << endl;
+  cout << "** Computing resolution of " << branchJet->GetName() << " induced by " << branchGenJet->GetName() << endl;
 
   Jet *jet, *genjet;
@@ -338 +338 @@
 
     // Loop over all reconstructed jets in event
-    for(i = 0; i < branchJet->GetEntriesFast(); ++i)
+    for(i = 0; i < TMath::Min(2,branchJet->GetEntriesFast()); ++i) //branchJet->GetEntriesFast(); ++i)
     {
 
@@ -347 +347 @@
 
      // Loop over all hard partons in event
-     for(j = 0; j < branchGenJet->GetEntriesFast(); ++j)
+     for(j = 0; j < TMath::Min(2,branchGenJet->GetEntriesFast()); ++j)
      {
         genjet = (Jet*) branchGenJet->At(j);
@@ -365 +365 @@
       }
 
-      if(deltaR < 0.1)
+      if(deltaR < 0.3)
       {
         pt  = genJetMomentum.Pt();
@@ -372 +372 @@
         for (bin = 0; bin < Nbins; bin++)
         {
-            if(pt > histos->at(bin).ptmin && pt < histos->at(bin).ptmax && eta > 0.0 && eta < 0.5) 
+            if(pt > histos->at(bin).ptmin && pt < histos->at(bin).ptmax && eta > 0.0 && eta < 2.5) 
             {
                 histos->at(bin).cenResolHist->Fill((bestGenJetMomentum.E()-jetMomentum.E())/bestGenJetMomentum.E());
@@ -389 +389 @@
     delete f1;
     return make_pair (sig, sigErr);
+    //return make_pair (hist->GetRMS(), hist->GetRMSError());
 }
 
@@ -395 +396 @@
 {
     Int_t bin;
-    TGraphErrors gr = TGraphErrors(Nbins);
+    Int_t count = 0;
+    TGraphErrors gr = TGraphErrors(Nbins/2);
     Double_t sig = 0;
     Double_t sigErr = 0;
@@ -406 +408 @@
             std::cout << " entries : " << histos->at(bin).cenResolHist->GetEntries() << std::endl;
             std::pair<Double_t, Double_t> sigvalues = GausFit(histos->at(bin).cenResolHist);
-            gr.SetPoint(bin,(histos->at(bin).ptmin+histos->at(bin).ptmax)/2.0, sigvalues.first);
-            gr.SetPointError(bin,0, sigvalues.second);
+            gr.SetPoint(count,(histos->at(bin).ptmin+histos->at(bin).ptmax)/2.0, sigvalues.first);
+            gr.SetPointError(count,0, sigvalues.second);
+            count++;
         }
         /*
@@ -516 +519 @@
   gPad->SetBottomMargin(0.2);
   gPad->SetLeftMargin(0.2);
+  gPad->Modified();
+  gPad->Update();
   
 }
@@ -653 +658 @@
   C_el2->SaveAs(elRes+".eps");
 
+  gDirectory->cd(0);
 
 /*
@@ -802 +808 @@
   TGraphErrors gr_pfjets = EresGraph(&plots_pfjets, true);
   gr_pfjets.SetName("pfJet");
+
+
+  // PFJets Energy Resolution
+  std::vector<resolPlot> plots_calojets;
+  HistogramsCollection(&plots_calojets, TMath::Log10(ptrangemin), TMath::Log10(ptrangemax), "CaloJet");
+  GetJetsEres( &plots_calojets, branchCaloJet, branchGenJet, treeReader);
+  TGraphErrors gr_calojets = EresGraph(&plots_calojets, true);
+  gr_calojets.SetName("caloJet");
+
+
+  TString jetRes = "jetERes";
+  TCanvas *C_jet = new TCanvas(jetRes,jetRes, 1000, 500);
+  TMultiGraph *mg_jet = new TMultiGraph(jetRes,jetRes);
+  TLegend *leg_jet = new TLegend(0.52,0.7,0.9,0.9);
+
+  addGraph(mg_jet, &gr_calojets, leg_jet, 3);
+  addGraph(mg_jet, &gr_pfjets, leg_jet, 1);
+
+  mg_jet->Draw("ACX");
+  leg_jet->Draw();
+
+  DrawAxis(mg_jet, leg_jet, 0.25);
+
+  C_jet->SaveAs(jetRes+".eps");
 
   
@@ -820 +850 @@
   {  
       plots_pfjets.at(bin).cenResolHist->Write();
+      plots_calojets.at(bin).cenResolHist->Write();
       plots_el.at(bin).cenResolHist->Write();
       plots_eltrack.at(bin).cenResolHist->Write();
@@ -848 +879 @@
   C_el1->Write();
   C_el2->Write();
+  C_jet->Write();
 /*
   C_mu->Write();
@@ -853 +885 @@
 */
   gr_pfjets.Write();
-  //gr_calojets.Write();
+  gr_calojets.Write();
 
   fout->Write();