ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/plot_diff.C
(Generate patch)

Comparing UserCode/LJMet/MultivariateAnalysis/root/plot_diff.C (file contents):
Revision 1.3 by jindal, Mon Feb 9 22:31:47 2009 UTC vs.
Revision 1.4 by kukartse, Fri Feb 13 00:53:03 2009 UTC

# Line 4 | Line 4
4    const int nvar =14 ;
5    int color;
6    TString var[nvar]={"Mt","ht","M123inv","centrality","sphericity","htPlusLepton","getHt3","ht2p","ktMinPrime","DphiJMET","aplanarity", "dPhiLMet","W_MT","metHtPlusLepton"};
7 <  TString lepton_ID = "eidRobustTight";
8 <  TString  tmp = "eidRobustTight.root";
7 >  TString lepton_ID = "TM2DCompatibilityLoose_caloIso";
8 >  TString  tmp = "TM2DCompatibilityLoose_caloIso.root";
9    TFile * f1 = new TFile(tmp);
10  
11    TH1F *hnomsig[nvar];
# Line 85 | Line 85 | for(j=0;j<nvar;j++){
85          leg->SetLineColor(0);
86          leg->Draw();
87          tmp= var[j]+"_"+lepton_ID+".png";
88 <        c1[j]->Print(tmp);
88 >        // c1[j]->Print(tmp);
89          int nbins = hnomsig[j]->GetNbinsX();
90 +        double chi2 = 0;
91 +        double chi = 0;
92 +        int ndof=0;
93          for(int i = 1; i < nbins+1; i++){
94                  SigNom = hnomsig[j]->GetBinContent(i);
95                  SigRev = hrevsig[j]->GetBinContent(i);
96 <                Diff = SigNom - SigRev;
96 >                Diff = SigRev - SigNom;
97                  if( N1[i][j] >0 &&  N2[i][j] > 0)
98 <                        Err  =(1./ sqrt(N1[i][j]))*SigNom + (1./ sqrt(N2[i][j]))*SigRev;
98 >                  {    
99 >                    Err  =(1./ sqrt(N1[i][j]))*SigNom + (1./ sqrt(N2[i][j]))*SigRev;
100 >                  }
101                  else
102 <                        Err = 0;
102 >                  {
103 >                    //  std::cout << "diff = " << Diff <<std::endl;
104 >                    Err = 0;}
105 >
106 >                if(SigNom > 0)
107 >                  {
108 >                    chi = (Diff*Diff)/SigNom;
109 >                    chi2 = chi2 + chi;
110 >                    ndof++;
111 >                  }
112 >        
113                  hdiff[j]->SetBinContent(i,Diff);
114 <                hdiff[j]->SetBinError(i,Err);}
114 >                hdiff[j]->SetBinError(i,Err);
115 >        }
116 >        std::cout << "Variable = " << var[j] << " Chi2 = " << chi2<< " ndof = " << ndof-1 <<std::endl;
117                  
118          c2[j]->cd();
119          hdiff[j]->SetTitle(var[j]+" difference");
# Line 105 | Line 122 | for(j=0;j<nvar;j++){
122          hdiff[j]->Draw();
123          tmp= var[j]+"_"+lepton_ID+"_diff.png";
124          
125 <        c2[j]->Print(tmp);
125 >        //c2[j]->Print(tmp);
126  
127  
128  

Diff Legend

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