File indexing completed on 2024-04-06 12:19:00
0001
0002 {
0003 gSystem->Load("drawEfficiency_C");
0004 gSystem->Load("drawEfficiencies_C");
0005
0006 TFile f("doubleTauL2_L25_Val_Z2Tau.root");
0007
0008 TString ReleaseVersion("CMSSW173");
0009 TString Scale("LinearScale");
0010
0011 TH1F* nMCetaTau= (TH1F*)f.Get("DQMData/TausAtGenLevel_isolatedL25PixelTau/eta_Tau_GenLevel");
0012 TH1F* nMCptTau= (TH1F*)f.Get("DQMData/TausAtGenLevel_isolatedL25PixelTau/pt_Tau_GenLevel");
0013
0014 TH1F* nMCetaTauJet = (TH1F*)f.Get("DQMData/TausAtGenLevel_isolatedL25PixelTau/nMC_Taus_vs_etaTauJet");
0015 TH1F* nMCptTauJet = (TH1F*)f.Get("DQMData/TausAtGenLevel_isolatedL25PixelTau/nMC_Taus_vs_ptTauJet");
0016 TH1F* nMCenergyTauJet = (TH1F*)f.Get("DQMData/TausAtGenLevel_isolatedL25PixelTau/nMC_Taus_vs_energyTauJet");
0017
0018 TH1F* nRecoJetetaTauJet = (TH1F*)f.Get("DQMData/ReconstructedJet_isolatedL25PixelTau/n_RecoJet_vs_etaTauJet");
0019 TH1F* nRecoJetptTauJet = (TH1F*)f.Get("DQMData/ReconstructedJet_isolatedL25PixelTau/n_RecoJet_vs_ptTauJet");
0020 TH1F* nRecoJetenergyTauJet = (TH1F*)f.Get("DQMData/ReconstructedJet_isolatedL25PixelTau/n_RecoJet_vs_energyTauJet");
0021 TH1F* nAssociatedTracks = (TH1F*)f.Get("DQMData/ReconstructedJet_isolatedL25PixelTau/Number_Associated_Tracks");
0022 TH1F* nSelectedTracks = (TH1F*)f.Get("DQMData/ReconstructedJet_isolatedL25PixelTau/Number_Selected_Tracks");
0023
0024 TH1F* nRecoJetLTetaTauJet = (TH1F*)f.Get("DQMData/ReconstructedJetWithLeadingTrack_isolatedL25PixelTau/n_RecoJet+LeadingTrack_vs_etaTauJet");
0025 TH1F* nRecoJetLTptTauJet = (TH1F*)f.Get("DQMData/ReconstructedJetWithLeadingTrack_isolatedL25PixelTau/n_RecoJet+LeadingTrack_vs_ptTauJet");
0026 TH1F* nRecoJetLTenergyTauJet = (TH1F*)f.Get("DQMData/ReconstructedJetWithLeadingTrack_isolatedL25PixelTau/n_RecoJet+LeadingTrack_vs_energyTauJet");
0027
0028 TH1F* nTaggedJetsetaTauJet = (TH1F*)f.Get("DQMData/TauTaggedJets_isolatedL25PixelTau/n_IsolatedTauTaggedJets_vs_etaTauJet");
0029 TH1F* nTaggedJetsptTauJet = (TH1F*)f.Get("DQMData/TauTaggedJets_isolatedL25PixelTau/n_IsolatedTauTaggedJets_vs_ptTauJet");
0030 TH1F* nTaggedJetsenergyTauJet = (TH1F*)f.Get("DQMData/TauTaggedJets_isolatedL25PixelTau/n_IsolatedTauTaggedJets_vs_energyTauJet");
0031
0032 TH1F* nEMTaggedJetsetaTauJet = (TH1F*)f.Get("DQMData/TauEMTaggedJets_isolatedL25PixelTau/n_IsolatedTauTaggedJets_vs_etaTauJet");
0033 TH1F* nEMTaggedJetsetaTauJet = (TH1F*)f.Get("DQMData/TauEMTaggedJets_ecalIsolation/n_EMIsolatedTauTaggedJets_vs_etaTauJet");
0034 TH1F* nEMTaggedJetsptTauJet = (TH1F*)f.Get("DQMData/TauEMTaggedJets_ecalIsolation/n_EMIsolatedTauTaggedJets_vs_ptTauJet");
0035 TH1F* nEMTaggedJetsenergyTauJet = (TH1F*)f.Get("DQMData/TauEMTaggedJets_ecalIsolation/n_EMIsolatedTauTaggedJets_vs_energyTauJet");
0036
0037
0038 TH1F* LeadingTrackPt_TaggedTau = (TH1F*)f.Get("DQMData/TauTaggedJets_isolatedL25PixelTau/LeadingTrackPt_After_Isolation");
0039 TH1F* DeltaR_LTandJet_TaggedTau = (TH1F*)f.Get("DQMData/TauTaggedJets_isolatedL25PixelTau/DeltaR_LT_and_Jet_After_Isolation");
0040 TH1F* SignalTracks_TaggedTau = (TH1F*)f.Get("DQMData/TauTaggedJets_isolatedL25PixelTau/Signal_Tks_After_Isolation");
0041 TH1F* SelectedTracks_TaggedTau = (TH1F*)f.Get("DQMData/TauTaggedJets_isolatedL25PixelTau/Selected_Tks_After_Isolation");
0042 TH1F* AssociatedTracks_TaggedTau = (TH1F*)f.Get("DQMData/TauTaggedJets_isolatedL25PixelTau/Associated_Tks_After_Isolation");
0043
0044 gStyle->SetOptStat("ie");
0045
0046 cout << endl<< " entries "<< nMCetaTauJet->Integral()<<endl;
0047
0048 c21 = new TCanvas("c21","", 10, 10, 500, 410);
0049 nMCetaTauJet->SetTitle("# MC Visible Taus (Visible Energy) "+ReleaseVersion);
0050 nMCetaTauJet->GetXaxis()->SetTitle("eta");
0051 nMCetaTauJet->DrawNormalized();
0052 c21->Print(TString("MCetaTauJets.gif"),"gif");
0053
0054
0055 c22 = new TCanvas("c22","", 10,10,500,410);
0056 nMCptTauJet->SetTitle("# MC Visible Taus (Visible Eneryg) "+ReleaseVersion);
0057 nMCptTauJet->GetXaxis()->SetTitle("Pt (GeV)");
0058 nMCptTauJet->DrawNormalized();
0059 c22->Print(TString("MCPtTauJets.gif"),"gif");
0060
0061 c23 = new TCanvas("c23","", 10,10,500,410);
0062 nMCenergyTauJet->SetTitle("# MC Visible Taus (Visible Energy) "+ReleaseVersion);
0063 nMCenergyTauJet->GetXaxis()->SetTitle("Energy (GeV)");
0064 nMCenergyTauJet->DrawNormalized();
0065 c23->Print(TString("MCEnergyTauJets.gif"),"gif");
0066
0067
0068
0069
0070
0071
0072
0073
0074
0075
0076
0077
0078
0079 c26 = new TCanvas("c26","", 10,10,500,410);
0080 nAssociatedTracks->SetTitle("Associated Tracks Matched Reconstructed Jet "+ReleaseVersion);
0081 nAssociatedTracks->GetXaxis()->SetTitle("# Associated Tracks");
0082 nAssociatedTracks->Draw();
0083 c26->Print(TString("AssociatedTracks.gif"),"gif");
0084
0085 c27 = new TCanvas("c27","", 10,10,500,410);
0086 nSelectedTracks->SetTitle("Selected Tracks Matched Reconstructed Jet "+ ReleaseVersion);
0087 nSelectedTracks->GetXaxis()->SetTitle("# Selected Tracks");
0088 nSelectedTracks->Draw();
0089 c27->Print(TString("SelectedTracks.gif"),"gif");
0090
0091
0092
0093
0094 c28 = new TCanvas("c28","", 110, 110, 600, 510);
0095 LeadingTrackPt_TaggedTau->SetTitle("Leading Track Pt Tagged Taus "+ ReleaseVersion);
0096 LeadingTrackPt_TaggedTau->GetXaxis()->SetTitle("Leading Track Pt (GeV)");
0097 LeadingTrackPt_TaggedTau->Draw();
0098 c28->Print(TString("LeadingTrackPt_TaggedTaus.gif"),"gif");
0099
0100 c29 = new TCanvas("c29","", 110, 110, 600, 510);
0101 DeltaR_LTandJet_TaggedTau->SetTitle("#Delta R(LeadingTrack, JetMomentum) Tagged Taus"+ ReleaseVersion);
0102 DeltaR_LTandJet_TaggedTau->GetXaxis()->SetTitle("#Delta R(LeadingTrack, JetMomentum)");
0103 DeltaR_LTandJet_TaggedTau->Draw();
0104 c29->Print(TString("DeltaR_LT_Jet_TaggedTaus.gif"),"gif");
0105
0106 c30 = new TCanvas("c30","",110,110,600,510);
0107 SignalTracks_TaggedTau->SetTitle("# Tracks inside Signal cone Tagged Taus > 1.0 GeV "+ReleaseVersion);
0108 SignalTracks_TaggedTau->GetXaxis()->SetTitle("# Signal Tracks");
0109 SignalTracks_TaggedTau->Draw();
0110 c30->Print(TString("SignalTracks_TaggedTaus.gif"),"gif");
0111
0112 c31 = new TCanvas("c31","",110,110,600,510);
0113 SelectedTracks_TaggedTau->SetTitle("# Selected Tracks Tagged Taus "+ReleaseVersion);
0114 SelectedTracks_TaggedTau->GetXaxis()->SetTitle("# Selected Tracks");
0115 SelectedTracks_TaggedTau->Draw();
0116 c31->Print(TString("SelectedTracks_TaggedTaus.gif"),"gif");
0117
0118 c32 = new TCanvas("c32","",110,110,600,510);
0119 AssociatedTracks_TaggedTau->SetTitle("# Associated Tracks Tagged Taus "+ReleaseVersion);
0120 AssociatedTracks_TaggedTau->GetXaxis()->SetTitle("# Associated Tracks");
0121 AssociatedTracks_TaggedTau->Draw();
0122 c32->Print(TString("AssociatedTracks_TaggedTaus.gif"),"gif");
0123
0124
0125
0126
0127
0128
0129
0130
0131
0132
0133
0134
0135
0136
0137
0138
0139
0140
0141
0142
0143
0144
0145
0146
0147
0148
0149
0150
0151
0152
0153
0154
0155
0156
0157
0158
0159
0160 TH1F* deltaRLeadingTrackJetP = (TH1F *)f.Get("DQMData/LeadingTrackPtAndDeltaRStudies_isolatedL25PixelTau/DeltaR_LeadingTrack_in_RecoJet");
0161 TH1F* leadingTrackPtJet = (TH1F *)f.Get("DQMData/LeadingTrackPtAndDeltaRStudies_isolatedL25PixelTau/Leading_track_pt_in_RecoJet");
0162
0163 c21 = new TCanvas("c21","",125,125,615,615);
0164 deltaRLeadingTrackJetP->Draw();
0165 c22 = new TCanvas("c22","",150,150,640,640);
0166 leadingTrackPtJet->Draw();
0167
0168
0169 TPaveText* text_ptCut = new TPaveText(0.66, 0.82, 0.92, 0.92, "brNDC");
0170 text_ptCut->SetBorderSize(0);
0171 text_ptCut->SetFillColor(0);
0172 text_ptCut->AddText(Form(" P_{T} > %.1f GeV", 5.0));
0173
0174 TGraphAsymmErrors* gr0 = new TGraphAsymmErrors(nTaggedJetsetaTauJet, nMCetaTauJet);
0175 c4 = new TCanvas("c4","IsolationTotalEfficiencyEta",30,30,520,520);
0176 drawEfficiency("Isolated/MC Visible Tau "+ReleaseVersion, gr0, "#eta (MC Vis. Tau)", (TH1F*) nMCetaTauJet->Clone(), c4, text_ptCut,Scale );
0177
0178 TGraphAsymmErrors* gr2 = new TGraphAsymmErrors(nRecoJetetaTauJet, nMCetaTauJet);
0179 c6 = new TCanvas("c6","CaloJetTotalEfficiencyEta",35,35,525,525);
0180 drawEfficiency("Jets(Matched)/MC Visible Taus "+ReleaseVersion, gr2, "#eta (MC Vis. Tau)",(TH1F*) nMCetaTauJet->Clone(), c6, text_ptCut, Scale);
0181
0182 TGraphAsymmErrors* gr100 = new TGraphAsymmErrors(nRecoJetLTetaTauJet, nMCetaTauJet);
0183 c100 = new TCanvas("c100", "EfficienciesStepByStepEta", 230, 230, 720, 720);
0184 drawEfficiencies("Different Steps/MC Visible Taus "+ReleaseVersion, gr2, gr100, gr0,gr100, "#eta (MC Vis. Tau)",(TH1F*) nMCetaTauJet->Clone(), c100, text_ptCut,"Calo", Scale);
0185
0186 TGraphAsymmErrors* gr3 = new TGraphAsymmErrors(nRecoJetLTetaTauJet, nRecoJetetaTauJet);
0187 c7 = new TCanvas("c7","FindingLeadTrackPartialEfficiencyEta",35,35,525,525);
0188 drawEfficiency("Jets+LeadTr/Jets(Matched) "+ReleaseVersion, gr3, "#eta (MC Vis. Tau)",(TH1F*) nRecoJetetaTauJet->Clone(), c7, text_ptCut, Scale);
0189
0190 TGraphAsymmErrors* gr4 = new TGraphAsymmErrors(nTaggedJetsetaTauJet, nRecoJetLTetaTauJet);
0191 c8 = new TCanvas("c8","IsolationPartialEfficencyEta",35,35,525,525);
0192 drawEfficiency("Isolated/Jets+LeadTr "+ReleaseVersion, gr4, "#eta (MC Vis. Tau)",(TH1F*) nRecoJetetaTauJet->Clone(), c8, text_ptCut, Scale);
0193
0194
0195
0196
0197 TPaveText* text_etaCut = new TPaveText(0.66, 0.82, 0.92, 0.92, "brNDC");
0198 text_etaCut->SetBorderSize(0);
0199 text_etaCut->SetFillColor(0);
0200 text_etaCut->AddText(Form(" %.1f < #eta < %.1f ", -2.5, 2.5));
0201
0202 TGraphAsymmErrors* gr50 = new TGraphAsymmErrors(nTaggedJetsptTauJet, nMCptTauJet);
0203 c40 = new TCanvas("c40","IsolationTotalEfficiencyPt",30,30,520,520);
0204 drawEfficiency("Isolated/MC Visible Tau "+ReleaseVersion, gr50, "P_{T} (GeV/c) (MC Vis. Tau)", (TH1F*) nMCptTauJet->Clone(), c40, text_etaCut, Scale);
0205
0206
0207 TGraphAsymmErrors* gr6 = new TGraphAsymmErrors(nRecoJetptTauJet, nMCptTauJet);
0208 c10 = new TCanvas("c10","CaloJetTotalEfficiencyPt",35,35,525,525);
0209 drawEfficiency("Jets(Matched)/MC Visible Taus "+ ReleaseVersion, gr6, "P_{T} (GeV/c) (MC Vis. Tau)",(TH1F*) nMCptTauJet->Clone(), c10, text_etaCut, Scale);
0210
0211 TGraphAsymmErrors* gr101 = new TGraphAsymmErrors(nRecoJetLTptTauJet, nMCptTauJet);
0212 c101 = new TCanvas("c101", "EfficienciesStepByStepPt", 230, 230, 720, 720);
0213 drawEfficiencies("Different Steps/MC Visible Taus "+ ReleaseVersion, gr6, gr101, gr50,gr50, "P_{T} (GeV/c) (MC Vis. Tau)",(TH1F*) nMCptTauJet->Clone(), c101, text_etaCut, "Calo", Scale);
0214 TGraphAsymmErrors* gr7 = new TGraphAsymmErrors(nRecoJetLTptTauJet, nRecoJetptTauJet);
0215 c11 = new TCanvas("c11","FindingLeadTrackPartialEfficiencyPt",35,35,525,525);
0216 drawEfficiency("Jets+LeadTr/Jets(Matched) "+ ReleaseVersion, gr7, "P_{T} (GeV/c) (MC Vis. Tau)",(TH1F*) nRecoJetptTauJet->Clone(), c11, text_etaCut, Scale);
0217
0218 TGraphAsymmErrors* gr8 = new TGraphAsymmErrors(nTaggedJetsptTauJet, nRecoJetLTptTauJet);
0219 c12 = new TCanvas("c12","IsolationPartialEfficencyPt",35,35,525,525);
0220 drawEfficiency("Isolated/Jets+LeadTr "+ ReleaseVersion, gr8, "P_{T} (GeV/c) (MC Vis. Tau)",(TH1F*) nRecoJetptTauJet->Clone(), c12, text_etaCut, Scale);
0221
0222
0223
0224
0225 TPaveText* text_bothCuts = new TPaveText(0.66, 0.82, 0.92, 0.92, "brNDC");
0226 text_bothCuts->SetBorderSize(0);
0227 text_bothCuts->SetFillColor(0);
0228 text_bothCuts->AddText(Form(" P_{T} > %.1f GeV", 5.0));
0229 text_bothCuts->AddText(Form(" %.1f < #eta < %.1f ", -2.5, 2.5));
0230
0231 TGraphAsymmErrors* gr90 = new TGraphAsymmErrors(nTaggedJetsenergyTauJet, nMCenergyTauJet);
0232 c90 = new TCanvas("c90","IsolationTotalEfficiencyEnergy",30,30,520,520);
0233 drawEfficiency("Isolated/MC Visible Tau "+ReleaseVersion, gr90, "Energy (GeV) (MC Vis. Tau)", (TH1F*) nMCenergyTauJet->Clone(), c90, text_bothCuts, Scale);
0234
0235
0236 TGraphAsymmErrors* gr10 = new TGraphAsymmErrors(nRecoJetenergyTauJet, nMCenergyTauJet);
0237 c14 = new TCanvas("c14","CaloJetTotalEfficiencyEnergy",35,35,525,525);
0238 drawEfficiency("Jets(Matched)/MC Visible Taus "+ ReleaseVersion, gr10, "Energy (GeV) (MC Vis. Tau)",(TH1F*) nMCenergyTauJet->Clone(), c14, text_bothCuts, Scale);
0239
0240 TGraphAsymmErrors* gr102 = new TGraphAsymmErrors(nRecoJetLTenergyTauJet, nMCenergyTauJet);
0241 c102 = new TCanvas("c102", "EfficienciesStepByStepEnergy", 230, 230, 720, 720);
0242 drawEfficiencies("Different Steps/MC Visible Taus "+ ReleaseVersion, gr10, gr102, gr90, gr90, "Energy (GeV) (MC Vis. Tau)",(TH1F*) nMCenergyTauJet->Clone(), c102, text_bothCuts, "Calo", Scale);
0243
0244 TGraphAsymmErrors* gr11 = new TGraphAsymmErrors(nRecoJetLTenergyTauJet, nRecoJetenergyTauJet);
0245 c15 = new TCanvas("c15","FindingLeadTrackPartialEfficiencyEnergy",35,35,525,525);
0246 drawEfficiency("Jets+LeadTr/Jets(Matched) "+ReleaseVersion, gr11, "Energy (GeV) (MC Vis. Tau)",(TH1F*) nRecoJetenergyTauJet->Clone(), c15, text_bothCuts, Scale);
0247
0248 TGraphAsymmErrors* gr12 = new TGraphAsymmErrors(nTaggedJetsenergyTauJet, nRecoJetLTenergyTauJet);
0249 c16 = new TCanvas("c16","IsolationPartialEfficencyEnergy",35,35,525,525);
0250 drawEfficiency("Isolated/Jets+LeadTr " +ReleaseVersion, gr12, "Energy (GeV) (MC Vis. Tau)",(TH1F*) nRecoJetenergyTauJet->Clone(), c16, text_bothCuts, Scale);
0251
0252
0253
0254
0255
0256
0257
0258
0259
0260
0261
0262
0263
0264
0265
0266
0267
0268
0269
0270
0271
0272
0273
0274
0275
0276
0277
0278
0279 }
0280
0281