StRoot: StSpinPool/StWalgoB2009/St2009W

00001 // $Id: St2009W_histo.cxx,v 1.20 2011/09/14 14:23:21 stevens4 Exp $
00002 //
00003 //*-- Author : Jan Balewski, MIT
00004 
00005 #include <TH2.h>
00006 #include <TList.h>
00007 #include <TLine.h>
00008 #include <TMath.h>
00009 
00010 #include "St2009WMaker.h"
00011 
00012 //________________________________________________
00013 //________________________________________________
00014 void
00015 St2009WMaker::initHistos(){
00016   const float PI=TMath::Pi();
00017 
00018   //...... data histograms
00019   memset(hA,0,sizeof(hA));
00020   TList *Lx;  TLine *ln;TH1 *h; 
00021   char txt[1000], txt0[100];
00022   int nCase=18;
00023   hA[0]=h=new TH1F("muStatEve","W-algo: event count",nCase,0,nCase);
00024   h->GetXaxis()->SetTitleOffset(0.4);  h->GetXaxis()->SetLabelSize(0.06);  h->GetXaxis()->SetTitleSize(0.05); h->SetMinimum(0.8);
00025   h->SetLineColor(kBlue);h->SetLineWidth(2);
00026   h->SetMarkerSize(2);//<-- large text
00027 
00028   char key[][200]={"inp","BHT3Id","L2wId","L2wBits","L2wET","L2wRnd","tpcOn","primVert","vertZ","Pt10",
00029                    "B-in","B200","TrB","Tr2Cl","eta1","goldW"};
00030   for(int i=0;i<16;i++) h->Fill(key[i],0.); // preset the order of keys
00031   
00032   hA[1]=h=new TH1F("muInTrg","muW input triggers, WARN: scrambled if manyruns are combined by hadd.C; trigID (random order)",nCase,0,nCase);
00033   h->GetXaxis()->SetLabelSize(0.06);
00034 
00035   hA[2]=h=new TH1F("mubX48","L2W-ET events vs. bXing; bXing= raw bx48",128,-0.5,127.5);
00036   h->SetFillColor(kGreen);
00037 
00038   hA[3]=h=new TH1F("mubX7","L2W-ET events vs. bXing; bXing= raw bx7",128,-0.5,127.5);
00039   h->SetFillColor(kBlue);
00040 
00041   hA[4]=new TH1F("mubX48v","L2W-ET & primVertex  vs. bXing; bXing= raw bx48",128,-0.5,127.5);
00042   hA[5]=h=new TH1F("mubX7v","L2W-ET & primVertex; bXing= raw bx7",128,-0.5,127.5);
00043   h->SetFillColor(kBlue);
00044 
00045   //... DMS data
00046   hA[6]=h=new TH1F("muDsm1","L2W-ET events DMS spectrum; DSM value",64,-0.5,63.5); 
00047   h->SetMinimum(0.8);
00048   hA[7]=h=new TH1F("muDsm2","L2W-Rnd events DMS spectrum; DSM value",64,-0.5,63.5); 
00049   h->SetMinimum(0.8);
00050   sprintf(txt,"L2W-ET events w/ DMS>%d vs.BTOW TP ID bXing; Hanks' TP ID",par_DsmThres);
00051   hA[8]=new TH1F("muDsm3",txt,300,-0.5,299.5);
00052   sprintf(txt,"L2W-ET events w/ DMS>%d & primVertexvs.BTOW TP ID bXing; Hanks' TP ID",par_DsmThres);
00053   hA[9]=h=new TH1F("muDsm4",txt,300,-0.5,299.5);
00054   h->SetFillColor(kBlue); h->SetLineColor(kBlue);
00055 
00056 
00057   //.... vertex histograms .....
00058   hA[10]=h=new TH1F("muVRf","PPV Vertex rank, funny X-axis; X=Log10(rank)+offset", 150, -9,25);
00059   Lx=h->GetListOfFunctions();
00060   ln=new TLine(0,0,0,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00061 
00062 
00063   hA[11]=h=new TH1F("muZv","Z of any vertex w/ rank>0;Z-vertex (cm)",100,-200,200);
00064   Lx=h->GetListOfFunctions();
00065   ln=new TLine(par_vertexZ,0,par_vertexZ,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00066   ln=new TLine(-par_vertexZ,0,-par_vertexZ,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00067 
00068   hA[12]=new TH1F("muNV","#  vertices per event, rank>0 & Z in range; # of vertices",10,0,10);
00069 
00070   //..... Tracks....
00071   hA[20]=h=new TH1F("muStatTrk","W-algo:  track  count; cases",nCase,0,nCase);
00072   h->GetXaxis()->SetTitleOffset(0.4);  h->GetXaxis()->SetLabelSize(0.06);  h->GetXaxis()->SetTitleSize(0.05); h->SetMinimum(0.8);
00073   h->SetLineColor(kGreen); h->SetLineWidth(2);
00074   char keyT[][200]={"101","pt1","nHit","Hfrac","Rin","Rout","ptOK","@B","CL","#Delta R","fr24",
00075                     "noNear","noAway","goldW"};
00076              
00077   for(int i=0;i<14;i++) h->Fill(keyT[i],0.); // preset the order of keys
00078 
00079 
00080   hA[21]=h=new TH1F("muTrNfit","primary track  in-selection & vertexZ; nFitPoints",50,0,50);
00081   Lx=h->GetListOfFunctions();
00082   ln=new TLine(par_nFitPts,0,par_nFitPts,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00083 
00084   hA[22]=h=new TH1F("muTrFitFrac","primary track in-selection & vertexZ; nFit/nPoss ",50,0,1.1);
00085   Lx=h->GetListOfFunctions();
00086   ln=new TLine(par_nHitFrac,0,par_nHitFrac,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00087 
00088   hA[23]=h=new TH1F("muTrRxyIn","primary track first hit  in-selection & vertexZ; Rxy (cm)",60,50,170.);
00089   Lx=h->GetListOfFunctions();
00090   ln=new TLine(par_trackRin,0,par_trackRin,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00091 
00092   hA[24]=h=new TH1F("muTrRxyOut","primary track last hit  in-selection & vertexZ; Rxy (cm)",60,100,220.);
00093   Lx=h->GetListOfFunctions();
00094   ln=new TLine(par_trackRout,0,par_trackRout,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00095   
00096   hA[25]=h=new TH1F("muTrPt1","global track PT ; track PT (GeV/c)",160,0,80);
00097   Lx=h->GetListOfFunctions();
00098   ln=new TLine(par_trackPt,0,par_trackPt,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln); h->SetFillColor(kYellow);
00099 
00100 
00101   hA[26]=h=new TH2F("muTr2D1","lastHit on track; detector eta ; detector phi (rad)",100,-1.1,1.1,200,-3.2,3.2);
00102 
00103   hA[27]=h=new TH1F("muTrPt1N","global NEGATIVE track PT; track PT (GeV/c)",160,0,80);
00104   h->SetFillColor(7);
00105 
00106   hA[28]=h=new TH2F("muTrdEdX"," dEdX vs. momentum; track P (GeV); dE/dx (keV)",20,0,10,100,0,10);
00107 
00108 
00109   hA[29]=h=new TH1F("muTrPt1Pr","primary  track PT; track PT (GeV/c)",160,0,80);
00110   hA[30]=h=new TH1F("muTrPt1NPr","primary  NEGATIVE track PT; track PT (GeV/c)",160,0,80);
00111 
00112 
00113   //..... BTOW .....
00114 
00115   hA[31]=h=new TH1F("muBmaxAdc","BTOW maxADC in event, in-selection; max tower ADC",200,0,5500);
00116   hA[32]=h=new TH1F("muBtotAdc","BTOW sum of ADC>thres , in-selection;ADC sum/event", 120,0,12000.);
00117 
00118   hA[33]=h=new TH1F("muBclET","matched BTOW 2x2 cluster ET  ;cluster  ET (GeV)",100,0,100);
00119   Lx=h->GetListOfFunctions();
00120   ln=new TLine(par_clustET,0,par_clustET,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00121 
00122   hA[34]=h=new TH2F("muBclAdcPt"," matched  TPC PT vs. 2x2 cluster ADC sum ; cluster (ADC sum);TPC  PT (GeV)",50,0,5000,75,0,150);
00123 
00124   hA[35]=new TH1F("muTrch2","track glob chi2/dof X-Y",100,0,5);
00125   hA[36]=new TH2F("muTrch2b","track glob chi2/dof; chi2/dof  X-Y; last hit eta",30,0,5.,30,-1,1);
00126 
00127   
00128   //.... 4x4 cluster
00129   hA[37]=h=new TH1F("muBclET24","matched BTOW 4x4 cluster ET ;cluster 4x4  ET (GeV)",100,0,100);
00130    hA[38]=h=new TH2F("muBclE242D","Excess energy in 4x4 cluster vs. 2x2 E;2x2 cluster E (GeV); E(4x4)-E(2x2)  E (GeV)",50,0,80,50,0,60);
00131 
00132    hA[39]=h=new TH1F("muBclET24R"," ratio (2x2/4x4) cluster ET ; fraction: cluster ET 2x2/ 4x4 ET",100,0,1.2);
00133    Lx=h->GetListOfFunctions();
00134   ln=new TLine(par_clustFrac24,0,par_clustFrac24,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00135 
00136    //..... jet energy , fraction
00137    
00138    hA[40]=h=new TH1F("muBjetET"," near 'EM jet' ET ; 'EM jet'  ET (GeV)",100,0,100);
00139 
00140    hA[41]=h=new TH2F("muBclEjetE2D","Excess nearCone ET  vs. 2x2 E;2x2 cluster ET (GeV); ET(cone-2x2) (GeV)",50,0,80,50,0,60);
00141 
00142    hA[42]=h=new TH1F("muBjetETR"," ratio (2x2/nearCone) ET ; cluster ET/ near cone ET",100,0,1.2);
00143   Lx=h->GetListOfFunctions();
00144   ln=new TLine(par_nearTotEtFrac,0,par_nearTotEtFrac,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00145   
00146 
00147   // .... track-EMC distance cuts
00148   hA[43]=h=new TH2F("muBdist1","3D Distance(track-cluster) vs. 2x2 E;2x2 cluster E (GeV); | distance | (cm)",50,0,80,50,0,25);
00149   Lx=h->GetListOfFunctions();
00150   ln=new TLine(0,par_delR3D,1.e6,par_delR3D);  ln->SetLineColor(kRed);  Lx->Add(ln);
00151   hA[44]=h=new TH2F("muBdist2","#Delta Z   (track-cluster) vs.Z-clust; Z-cluster (cm); #Delta Z (cm)",100,-300,300,40,-20,20);
00152   hA[45]=h=new TH2F("muBdist3","R#Delta #phi   (track-cluster) vs. 2x2 E;2x2 cluster E (GeV); R#Delta #phi (cm)",50,0,80,80,-20,20);
00153   hA[46]=h=new TH1F("muBdist4","3D Distance(track-cluster) vs. 2x2 E;| 3D distance |   (cm)",100,0,50);
00154   Lx=h->GetListOfFunctions();
00155   ln=new TLine(par_delR3D,0,par_delR3D,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00156 
00157   //.. continue same side jet veto
00158   hA[47]=h=new TH1F("muTjetET"," quenched near 'TPC jet' PT ; 'TPC jet'  PT (GeV)",100,0,100);
00159   hA[48]=h=new TH2F("muTjetBjet2D","quenched  near 'TPC jet' vs.  near 'EMC jet' ; EMC jet ET (GeV); TPC PT  (GeV/c)",50,0,80,50,0,60);
00160   hA[49]=h=new TH1F("muTBjetET"," near 'TPC+EMC jet' ET ; jet ET (GeV), no double counting",100,0,100);
00161 
00162   //.. away-side jet  veto
00163   hA[50]=h=new TH1F("muTwayET"," TPC away-cone PT sum;   PT (GeV)",100,0,100);
00164   hA[51]=h=new TH1F("muBwayET"," BTOW away-cone ET sum;   ET (GeV)",100,0,100);
00165   hA[52]=h=new TH2F("muTotwayET2D"," away TPC+EMC ET sum  vs. 2x2 cluster ET;  2x2 ET (GeV); away ET (GeV)",50,0,100,150,0,100);
00166 
00167 
00168   hA[53]=h=new TH2F("muAwayET2Db"," away EMC ET sum  vs. 2x2 cluster ET;  2x2 ET (GeV); away EMC ET (GeV)",50,0,100,80,0,80);
00169 
00170 
00171   hA[54]=h=new TH1F("muAwayTotEt"," away-cone TPC+EMC ET sum ; away ET (GeV)",200,0,100);
00172 
00173   hA[55]=h=new TH1F("muEwayET"," ETOW away-cone ET sum;   ET (GeV)",100,0,100); // away side energy  
00174   
00175   hA[56]=h=new TH1F("muSmallNearTpcPT",Form("TPC PT in #Delta R =%.1f from lepton candidate; PT (GeV)",par_smallNearDeltaR),100,0,100);
00176   hA[57]=h=new TH2F("muTr2D1pt5","lastHit on track (pt > 5); detector eta ; detector phi (rad)",100,-1.1,1.1,240,-PI,PI);
00177 
00178 
00179   hA[58]=new TH1F("muTrch2West","track glob chi2/dof  West TPC ; chi2/dof",100,0,5);
00180   hA[59]=new TH1F("muTrch2East","track glob chi2/dof  East TPC ; chi2/dof",100,0,5);
00181 
00182   //... final golden plots ....
00183 
00184   hA[60]=h=new TH2F("muBclETPt"," TPC PT vs.  isolated cluster 2x2 ET, matched;2x2 cluster ET (GeV) ; TPC  PT (GeV)",50,0,100,75,0,150);
00185   
00186   hA[61]=new TH1F("mubX7bht3","L2W-BHT3-rnd  events vs. bXing; bXing= raw bx7",128,-0.5,127.5); // filled on input
00187   hA[62]=h=new TH2F("muEne_Deta","final W: cluster energy vs. detector eta; barrel eta bin; 2x2 Energy (GeV)",40,0,40,50,0,100);
00188 
00189   // free 63-69
00190 
00191   //..... BSMD ......reserve 2x10= [70:e.... 80:p.... 89]
00192   const char cPlane[ mxBSmd]={'E','P'};
00193   for(int iep=0;iep<mxBSmd;iep++){
00194     sprintf(txt0,"muS%cadc1",cPlane[iep]);
00195     sprintf(txt,"BSMD-%c whole plane ADC; ADC",cPlane[iep]);
00196     hA[70+10*iep]=h=new TH1F(txt0,txt,275,0,1100);
00197   }// end of E-,P-planes
00198 
00199   //.... final Ws
00200   hA[90]=h=new TH1F("muWET","  Final selection; 2x2 cluster ET (GeV)", 100,0,100);
00201   Lx=h->GetListOfFunctions();
00202   ln=new TLine(par_highET,0,par_highET,1.e6);  ln->SetLineColor(kRed);  Lx->Add(ln);
00203 
00204   sprintf(txt," Final W selection, ET>%.0f GeV 'goldenW'; detector eta ; detector phi (rad)",par_highET);
00205   hA[91]=new TH2F("muW2D1",txt,10,-1.0,1.0,24,-PI,PI);
00206 
00207   hA[92]=new TH2F("muWdedx","Track dEdx, final W selection; 2x2 ET (GeV); dEdx (keV)",100,0,100,100,0,10);
00208   hA[93]=new TH2F("muWglDca","Track glob vertex abs(DCA), final W ; 2x2 ET (GeV); |DCA| (cm)",100,0,100,100,0,5);
00209   hA[94]=new TH2F("muWglDcaSP","Track prim POSITIVE glob signed DCA, final W; 2x2 ET (GeV); sDCA (cm)",100,0,100,100,-5,5);
00210   hA[95]=new TH2F("muWglDcaSN","Track prim NEGATIVE glob signed DCA, final W ; 2x2 ET (GeV); sDCA (cm)",100,0,100,100,-5,5);
00211 
00212   sprintf(txt,"Vertex ID, final W selection, 2x2 ET>%.0f GeV; vertex ID",par_highET);
00213   hA[96]=h=new TH1F("muWcar1",txt,10,-0.5,9.5);
00214   h->SetLineColor(kBlue);h->SetLineWidth(3);
00215 
00216   sprintf(txt,"Vertex 'funny' rank, final W selection, 2x2 ET>%.0f GeV;  X=Log10(rank)+offset",par_highET);
00217   hA[97]=new TH1F("muWcar2",txt, 150, -9,25);
00218   
00219   sprintf(txt,"Vertex Z , final W selection, 2x2 ET>%.0f GeV; Z(cm)",par_highET);
00220   hA[98]=new TH1F("muWcar3",txt, 100, -200,200);
00221   hA[99]=h=new TH1F("muWeta","final Ws ; lepton eta",100, -1.5,1.5);
00222 
00223   // free 100-109
00224 
00225   //..... series of electron ET plots after succesive cuts
00226   char tt2[][200]={"max 2x2","track matched","2x2 / 4x4","no near ET","no away ET"};
00227   for(int i=0;i<5;i++){
00228     sprintf(txt,"electron candidate, cut=%s; 2x2 ET (GeV)",tt2[i]);
00229     sprintf(txt0,"muETlive%d",i);
00230     hA[110+i]=h=new TH1F(txt0,txt, 100,0,100);
00231   }
00232   
00233   //free 115-131
00234 
00235   hA[117]=h=new TH2F("mujetQAeta_phi","Input Jet phi vs eta ;  eta ; phi ",50,-3,3,63,-PI,PI);
00236   hA[118]=h=new TH1F("mujetQApt","Input Jet pt; pt;",100,0,100);
00237 
00238   hA[132]=h=new TH2F("muptBalance_clust","ptBalance vs cluster ET; 2x2 Cluster ET; ptBalance",100,0,100,100,0,100);
00239   hA[133]=h=new TH2F("muptBalance_awayTot","ptBalance vs awayside PT; awayside PT; ptBalance",100,0,100,100,0,100);
00240 
00241   hA[134]=h=new TH2F("musPtBalance_clust","sPtBalance vs cluster ET; 2x2 Cluster ET (GeV); signed Pt balance (GeV)",100,0,100,100,-100,100);
00242   Lx=h->GetListOfFunctions();
00243   ln=new TLine(0,par_ptBalance,100,par_ptBalance);  ln->SetLineColor(kRed);  Lx->Add(ln);
00244 
00245   hA[135]=h=new TH2F("musPtBalance_awayTot","sPtBalance vs awayside PT; awayside PT; sPtBalance",100,0,100,100,-100,100);
00246   Lx=h->GetListOfFunctions();
00247   ln=new TLine(0,par_ptBalance,100,par_ptBalance);  ln->SetLineColor(kRed);  Lx->Add(ln);
00248 
00249 
00250   hA[136]=h=new TH1F("muclustPtBal",Form("PT Balance > %.1f ; 2x2 Cluster ET",par_ptBalance),100,0,100);
00251   hA[137]=h=new TH1F("muclustPtBal_bckgrd",Form("PT Balance < %.1f ; 2x2 Cluster ET",par_ptBalance),100,0,100);
00252   hA[138]=h=new TH1F("muclustAwayPt","AwaySide PT < 8 ; 2x2 Cluster ET",100,0,100);
00253   hA[139]=h=new TH1F("muclustAwayPt_bckgrd","AwaySide PT > 8 ; 2x2 Cluster ET",100,0,100);
00254   hA[140]=h=new TH1F("muclustPtBalnoE",Form("sPT Balance > %.1f (EEMC not included); 2x2 Cluster ET",par_ptBalance),100,0,100);
00255   hA[141]=h=new TH1F("muclustAwayPtnoE","AwaySide PT < 8 (EEMC not included); 2x2 Cluster ET",100,0,100);
00256 
00257   // Histograms added for background subtraction and systematic
00258   char str[200];
00259   for (int i=0; i<=20; i++) {
00260     sprintf(str,"neg_failAwaySide_Awayside_pt_bin_%d",i);
00261     hA[142+i] = new TH2F(str,str,100,0,100,21,0,21);
00262   }
00263 
00264   for (int i=0; i<=20; i++) {
00265     sprintf(str,"pos_failAwaySide_Awayside_pt_bin_%d",i);
00266     hA[163+i] = new TH2F(str,str,100,0,100,21,0,21);
00267   }
00268 
00269   hA[184+2] = new TH1F("pos_muclustpTbal_wE","pos_muclustpTbal_wE",100,0,100);
00270   hA[184+1] = new TH1F("neg_muclustpTbal_wE","neg_muclustpTbal_wE",100,0,100);
00271   hA[184+4] = new TH1F("pos_muclustpTbal_noE","pos_muclustpTbal_noE",100,0,100);
00272   hA[184+3] = new TH1F("neg_muclustpTbal_noE","neg_muclustpTbal_noE",100,0,100);
00273   hA[184+6] = new TH1F("pos_muclustpTbal_back","pos_muclustpTbal_back",100,0,100);
00274   hA[184+5] = new TH1F("neg_muclustpTbal_back","neg_muclustpTbal_back",100,0,100);
00275 
00276   hA[190] = new TH1F("muZvReweight","Reweight Z vertex to match data",100,-200,200);
00277 
00278   //algo efficiency ET dependence
00279   hA[191]=h=new TH2F("muBclET24R_ET","ratio (2x2/4x4) cluster ET vs 2x2 cluster ET ; 2x2 cluster ET (GeV); fraction: cluster ET 2x2/ 4x4 ET",100,0,100,100,0,1.2);
00280   Lx=h->GetListOfFunctions();
00281   ln=new TLine(0,par_clustFrac24,1.e6,par_clustFrac24);  ln->SetLineColor(kRed);  Lx->Add(ln);
00282   hA[192]=h=new TH2F("muBclEjetE2D_ET","ratio (2x2/nearCone) ET vs. 2x2 cluster ET; 2x2 cluster ET (GeV); ET(cone-2x2) (GeV)",100,0,100,100,0,1.2);
00283   Lx=h->GetListOfFunctions();
00284   ln=new TLine(0,par_nearTotEtFrac,1.e6,par_nearTotEtFrac);  ln->SetLineColor(kRed);  Lx->Add(ln);
00285 
00286   //eta dependent signal and background
00287   hA[200+2] = new TH2F("pos_muclustpTbal_wE_etaBin","pos_muclustpTbal_wE_etaBin",100,-1,1,100,0,100);
00288   hA[200+1] = new TH2F("neg_muclustpTbal_wE_etaBin","neg_muclustpTbal_wE_etaBin",100,-1,1,100,0,100);
00289   hA[200+4] = new TH2F("pos_muclustpTbal_noE_etaBin","pos_muclustpTbal_noE_etaBin",100,-1,1,100,0,100);
00290   hA[200+3] = new TH2F("neg_muclustpTbal_noE_etaBin","neg_muclustpTbal_noE_etaBin",100,-1,1,100,0,100);
00291   hA[200+6] = new TH2F("pos_muclustpTbal_back_etaBin","pos_muclustpTbal_back_etaBin",100,-1,1,100,0,100);
00292   hA[200+5] = new TH2F("neg_muclustpTbal_back_etaBin","neg_muclustpTbal_back_etaBin",100,-1,1,100,0,100);
00293 
00294   //eta dependent histos for background subtraction systematic
00295   for (int i=0; i<=20; i++) {
00296     sprintf(str,"neg_failsPtBal_sPtBal_bin_%d",i);
00297     hA[210+i] = new TH2F(str,str,100,-1,1,100,0,100);
00298   }
00299   for (int i=0; i<=20; i++) {
00300     sprintf(str,"pos_failsPtBal_sPtBal_bin_%d",i);
00301     hA[231+i] = new TH2F(str,str,100,-1,1,100,0,100);
00302   }
00303 
00304   //some golden W plots for comparison to embedding
00305   hA[260]=h=new TH2F("muWbX7","L2W-ET events vs. bXing; bXing= raw bx7",100,0,200000,128,-0.5,127.5);
00306   hA[261]=h=new TH2F("muWNV","#  vertices per event, rank>0 & Z in range; # of vertices",100,0,200000,10,0,10);
00307 
00308   hA[262]=h=new TH2F("muWTrNfit","prim tr nFitP; nFitPoints",100,0,200000,50,0,50);
00309   hA[263]=h=new TH2F("muWTrFitFrac","prim tr nFitFrac; nFit/nPoss ",100,0,200000,50,0,1.1);
00310   
00311   hA[264]=h=new TH2F("muWTrRxyOut","prim tr last hit filter; Rxy (cm)",100,0,200000,60,100,220.);
00312   hA[265]=h=new TH2F("muWTrRxyIn","prim tr 1st hit  filter; Rxy (cm)",100,0,200000,60,50,170.);
00313   hA[266]=h=new TH2F("muWTrPt1","primary track PT ; track PT (GeV/c)",100,0,200000,160,0,80);
00314   hA[267]=h=new TH2F("muWTrch2","track glob chi2/dof X-Y",100,0,200000,100,0,5);
00315   hA[268]=h=new TH2F("muWTrInvPt","primary track 1/PT",100,0,200000,100,0,0.1);
00316   hA[269]=h=new TH2F("muWglDcaGold","Track glob vertex abs(DCA), final W; |DCA| (cm)",100,0,200000,100,0,5);
00317   hA[270]=h=new TH2F("muWglsDcaGold","Track glob signed DCA, final W; sDCA (cm)",100,0,200000,100,-5,5);
00318   hA[271]=h=new TH2F("muWglDcaZGold","Track glob Z DCA, final W; DCAZ (cm)",100,0,200000,100,-5,5);
00319   hA[272]=h=new TH2F("muWglPt1TrVert","Track global Pt, final W from one track vertices",100,0,200000,160,0,80);
00320   hA[273]=h=new TH2F("muWglInvPt1TrVert","Track global 1/Pt, final W from one track vertices",100,0,200000,100,0,0.1);
00321 
00322   //2 eta bins for x-section ratio
00323   hA[280+2] = new TH2F("pos_muclustpTbal_wE_etaBin2","pos_muclustpTbal_wE_etaBin2",100,0,1,100,0,100);
00324   hA[280+1] = new TH2F("neg_muclustpTbal_wE_etaBin2","neg_muclustpTbal_wE_etaBin2",100,0,1,100,0,100);
00325   hA[280+4] = new TH2F("pos_muclustpTbal_noE_etaBin2","pos_muclustpTbal_noE_etaBin2",100,0,1,100,0,100);
00326   hA[280+3] = new TH2F("neg_muclustpTbal_noE_etaBin2","neg_muclustpTbal_noE_etaBin2",100,0,1,100,0,100);
00327   hA[280+6] = new TH2F("pos_muclustpTbal_back_etaBin2","pos_muclustpTbal_back_etaBin2",100,0,1,100,0,100);
00328   hA[280+5] = new TH2F("neg_muclustpTbal_back_etaBin2","neg_muclustpTbal_back_etaBin2",100,0,1,100,0,100);
00329   //eta dependent histos for background subtraction systematic
00330   for (int i=0; i<=20; i++) {
00331     sprintf(str,"neg_failsPtBal_sPtBal_bin_%d_etaBin2",i);
00332     hA[290+i] = new TH2F(str,str,100,0,1,100,0,100);
00333   }
00334   for (int i=0; i<=20; i++) {
00335     sprintf(str,"pos_failsPtBal_sPtBal_bin_%d_etaBin2",i);
00336     hA[311+i] = new TH2F(str,str,100,0,1,100,0,100);
00337   }
00338 
00339   // add histos to the list (if provided)
00340   for(int i=0;i<mxHA;i++) {
00341     if(  hA[i]==0) continue;
00342     HList->Add( hA[i]);
00343   }
00344   
00345   //  HList->ls();
00346   LOG_INFO<<Form("%s::initHistos done",GetName())<<endm;
00347 
00348 }
00349 
00350 // $Log: St2009W_histo.cxx,v $
00351 // Revision 1.20  2011/09/14 14:23:21  stevens4
00352 // update used for cross section PRD paper
00353 //
00354 // Revision 1.19  2010/11/09 23:00:50  balewski
00355 // added chi2/dof for East & West TPC separately
00356 //
00357 // Revision 1.18  2010/05/01 01:31:44  balewski
00358 // added W->JJ code & JES calibration
00359 //
00360 // Revision 1.17  2010/04/27 16:53:45  stevens4
00361 // add code to remove events tagged as Zs from W candidates
00362 //
00363 // Revision 1.16  2010/03/23 01:31:40  seelej
00364 // Fix to the filling of the histograms for the background systematic.
00365 //
00366 // Revision 1.15  2010/03/20 18:38:34  balewski
00367 // *** empty log message ***
00368 //
00369 // Revision 1.14  2010/03/18 16:52:17  balewski
00370 // corrected sPtBalance for no-endcap
00371 //
00372 // Revision 1.13  2010/03/14 22:50:31  balewski
00373 // *** empty log message ***
00374 //
00375 // Revision 1.12  2010/02/22 15:49:34  seelej
00376 // Joe : Changes to code for inclusion of background subtraction and systematic studies
00377 //
00378 // Revision 1.11  2010/01/29 01:56:01  stevens4
00379 // disable lepton track reco in TPC sector 20
00380 //
00381 // Revision 1.10  2010/01/28 20:10:05  balewski
00382 // added eta dependent spin sorting
00383 //
00384 // Revision 1.9  2010/01/27 22:12:24  balewski
00385 // spin code matched to x-section code
00386 //
00387 // Revision 1.8  2010/01/18 03:26:15  balewski
00388 // expanded TPC track filtering, not finished
00389 //
00390 // Revision 1.7  2010/01/10 01:45:10  stevens4
00391 // fix plots w/o EEMC in veto
00392 //
00393 // Revision 1.6  2010/01/09 02:29:19  stevens4
00394 // fix histo names
00395 //
00396 // Revision 1.5  2010/01/09 00:07:16  stevens4
00397 // add jet finder
00398 //
00399 // Revision 1.4  2010/01/06 19:16:48  stevens4
00400 // track cuts now on primary component, cleanup
00401 //
00402 // Revision 1.3  2010/01/06 14:11:13  balewski
00403 // one Z-plot added
00404 //
00405 // Revision 1.2  2009/12/08 04:48:35  balewski
00406 // *** empty log message ***
00407 //
00408 // Revision 1.1  2009/11/23 23:00:18  balewski
00409 // code moved spin-pool
00410 //
StSpinPool/StWalgoB2009/St2009W_histo.cxx
