[ViewVC] Diff of: cvsroot/UserCode/cbrown/AnalysisFramework/Plotting/Modules/DustyModule.C

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/DustyModule.C (file contents):
Revision 1.4 by buchmann, Wed Jul 20 08:51:33 2011 UTC vs.
Revision 1.5 by buchmann, Wed Aug 31 08:08:33 2011 UTC

-+
-+
-+
-+
-+
-+
-+
-+
-+
-+
-+
-+
-+
+void alpha_scan_susy_space(string mcjzb, string datajzb) {
-+
+  TCanvas *c3 = new TCanvas("c3","c3");
-+
+  vector<float> binning;
-+
+  binning=allsamples.get_optimal_binsize(mcjzb,cutmass&&cutOSSF&&cutnJets,20,50,800);
-+
+  float arrbinning[binning.size()];
-+
+  for(int i=0;i<binning.size();i++) arrbinning[i]=binning[i];
-+
+  TH1F *puredata   = allsamples.Draw("puredata",  datajzb,binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,data,luminosity);
-+
+  puredata->SetMarkerSize(DataMarkerSize);
-+
+  TH1F *allbgs   = allsamples.Draw("allbgs",  "-"+datajzb,binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,data,luminosity);
-+
+  TH1F *allbgsb   = allsamples.Draw("allbgsb",  "-"+datajzb,binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,data,luminosity);
-+
+  TH1F *allbgsc   = allsamples.Draw("allbgsc",  datajzb,binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,data,luminosity);
-+
+  allbgs->Add(allbgsb,-1);
-+
+  allbgs->Add(allbgsc);
-+
+  int ndata=puredata->Integral();
-+
+  ofstream myfile;
-+
+  myfile.open ("susyscan_log.txt");
-+
+  TFile *susyscanfile = new TFile("/scratch/fronga/SMS/T5z_SqSqToQZQZ_38xFall10.root");
-+
+  TTree *suevents = (TTree*)susyscanfile->Get("events");
-+
+  TH2F *exclusionmap = new TH2F("exclusionmap","",20,0,500,20,0,1000);
-+
+  TH2F *exclusionmap1s = new TH2F("exclusionmap1s","",20,0,500,20,0,1000);
-+
+  TH2F *exclusionmap2s = new TH2F("exclusionmap2s","",20,0,500,20,0,1000);
-+
+  TH2F *exclusionmap3s = new TH2F("exclusionmap3s","",20,0,500,20,0,1000);
-+
-+
+  susy_scan_axis_labeling(exclusionmap);
-+
+  susy_scan_axis_labeling(exclusionmap1s);
-+
+  susy_scan_axis_labeling(exclusionmap2s);
-+
+  susy_scan_axis_labeling(exclusionmap3s);
-+
-+
+  Int_t MyPalette[100];
-+
+  Double_t r[]    = {0., 0.0, 1.0, 1.0, 1.0};
-+
+  Double_t g[]    = {0., 0.0, 0.0, 1.0, 1.0};
-+
+  Double_t b[]    = {0., 1.0, 0.0, 0.0, 1.0};
-+
+  Double_t stop[] = {0., .25, .50, .75, 1.0};
-+
+  Int_t FI = TColor::CreateGradientColorTable(5, stop, r, g, b, 100);
-+
+  for (int i=0;i<100;i++) MyPalette[i] = FI+i;
-+
-+
+  gStyle->SetPalette(100, MyPalette);
-+
-+
+  for(int m23=50;m23<500;m23+=25) {
-+
+    for (int m0=(2*(m23-50)+150);m0<=1000;m0+=50)
-+
+    {
-+
+      c3->cd();
-+
+      stringstream drawcondition;
-+
+      drawcondition << "pfJetGoodNum>=3&&(TMath::Abs(masses[0]-"<<m0<<")<10&&TMath::Abs(masses[2]-masses[3]-"<<m23<<")<10)&&mll>5&&id1==id2";
-+
+      TH1F *puresignal = new TH1F("puresignal","puresignal",binning.size()-1,arrbinning);
-+
+      TH1F *puresignall= new TH1F("puresignall","puresignal",binning.size()-1,arrbinning);
-+
+      stringstream drawvar,drawvar2;
-+
+      drawvar<<mcjzb<<">>puresignal";
-+
+      drawvar2<<"-"<<mcjzb<<">>puresignall";
-+
+      suevents->Draw(drawvar.str().c_str(),drawcondition.str().c_str());
-+
+      suevents->Draw(drawvar2.str().c_str(),drawcondition.str().c_str());
-+
+      if(puresignal->Integral()<60) {
-+
+        delete puresignal;
-+
+        continue;
-+
+      }
-+
+      puresignal->Add(puresignall,-1);//we need to correct for the signal contamination - we effectively only see (JZB>0)-(JZB<0) !!
-+
+      puresignal->Scale(ndata/(20*puresignal->Integral()));//normalizing it to 5% of the data
-+
+      stringstream saveas;
-+
+      saveas<<"Model_Scan/m0_"<<m0<<"__m23_"<<m23;
-+
+      dout << "PLEASE KEEP IN MIND THAT SIGNAL CONTAMINATION IS NOT REALLY TAKEN CARE OF YET DUE TO LOW STATISTICS! SHOULD BE SOMETHING LIKE THIS : "<< endl;
-+
+//        TH1F *signalpredlo   = allsamples.Draw("signalpredlo",  "-"+mcjzb,  binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,mc,  luminosity,allsamples.FindSample("LM4"));
-+
+//        TH1F *signalpredro   = allsamples.Draw("signalpredro",      mcjzb,  binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,mc,  luminosity,allsamples.FindSample("LM4"));
-+
+//        TH1F *puredata       = allsamples.Draw("puredata",          datajzb,binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,data,luminosity);
-+
+//        signalpred->Add(signalpredlo,-1);
-+
+//        signalpred->Add(signalpredro);
-+
+//        puresignal->Add(signalpred,-1);//subtracting signal contamination
-+
+//---------------------
-+
+//      dout << "(m0,m23)=("<<m0<<","<<m23<<") contains " << puresignal->Integral() << endl;
-+
+//    TH1F *puresignal = allsamples.Draw("puresignal",mcjzb,  binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,mc,  luminosity,allsamples.FindSample("LM4"));
-+
+      vector<float> results=establish_upper_limits(puredata,allbgs,puresignal,saveas.str(),myfile);
-+
+      if(results.size()==0) {
-+
+        delete puresignal;
-+
+        continue;
-+
+      }
-+
+      exclusionmap->Fill(m23,m0,results[0]);
-+
+      exclusionmap1s->Fill(m23,m0,results[1]);
-+
+      exclusionmap2s->Fill(m23,m0,results[2]);
-+
+      exclusionmap3s->Fill(m23,m0,results[3]);
-+
+      delete puresignal;
-+
+      dout << "(m0,m23)=("<<m0<<","<<m23<<") : 3 sigma at " << results[3] << endl;
-+
+    }
-+
+  }//end of model scan for loop
-+
-+
+  dout << "Exclusion Map contains" << exclusionmap->Integral() << " (integral) and entries: " << exclusionmap->GetEntries() << endl;
-+
+  c3->cd();
-+
+  exclusionmap->Draw("CONTZ");
-+
+  CompleteSave(c3,"Model_Scan/CONT/Model_Scan_Mean_values");
-+
+  exclusionmap->Draw("COLZ");
-+
+  CompleteSave(c3,"Model_Scan/COL/Model_Scan_Mean_values");
-+
-+
+  exclusionmap1s->Draw("CONTZ");
-+
+  CompleteSave(c3,"Model_Scan/CONT/Model_Scan_1sigma_values");
-+
+  exclusionmap1s->Draw("COLZ");
-+
+  CompleteSave(c3,"Model_Scan/COL/Model_Scan_1sigma_values");
-+
-+
+  exclusionmap2s->Draw("CONTZ");
-+
+  CompleteSave(c3,"Model_Scan/CONT/Model_Scan_2sigma_values");
-+
+  exclusionmap2s->Draw("COLZ");
-+
+  CompleteSave(c3,"Model_Scan/COL/Model_Scan_2sigma_values");
-+
-+
+  exclusionmap3s->Draw("CONTZ");
-+
+  CompleteSave(c3,"Model_Scan/CONT/Model_Scan_3sigma_values");
-+
+  exclusionmap3s->Draw("COLZ");
-+
+  CompleteSave(c3,"Model_Scan/COL/Model_Scan_3sigma_values");
-+
-+
+  TFile *exclusion_limits = new TFile("exclusion_limits.root","RECREATE");
-+
+  exclusionmap->Write();
-+
+  exclusionmap1s->Write();
-+
+  exclusionmap2s->Write();
-+
+  exclusionmap3s->Write();
-+
+  exclusion_limits->Close();
-+
+  susyscanfile->Close();
-+
-+
+  myfile.close();
-+
+}

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/DustyModule.C (file contents): Revision 1.4 by buchmann, Wed Jul 20 08:51:33 2011 UTC vs. Revision 1.5 by buchmann, Wed Aug 31 08:08:33 2011 UTC

Diff Legend

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/DustyModule.C (file contents):
Revision 1.4 by buchmann, Wed Jul 20 08:51:33 2011 UTC vs.
Revision 1.5 by buchmann, Wed Aug 31 08:08:33 2011 UTC