00001
00002
00003
00004
00005 #include <TH2.h>
00006 #include <TList.h>
00007 #include <TLine.h>
00008 #include <TMath.h>
00009
00010 #include "St2011WMaker.h"
00011
00012
00013
00014 void
00015 St2011WMaker::initHistos(){
00016 const float PI=TMath::Pi();
00017
00018
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","Barrel 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);
00027
00028 char key[][200]={"inp","L2bwId","L2bwBits","L2bwET","L2bwRnd","tpcOn","primVert","vertZ","Pt10",
00029 "B-in","B200","TrB","Tr2Cl","eta1","noNear","noAway","goldW"};
00030 for(int i=0;i<17;i++) h->Fill(key[i],0.);
00031
00032 hA[1]=h=new TH1F("muInTrg","mu Barrel W 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","L2WB-ET events vs. bXing; bXing= raw bx48",128,-0.5,127.5);
00036 h->SetFillColor(kGreen);
00037
00038 hA[3]=h=new TH1F("mubX7","L2WB-ET events vs. bXing; bXing= raw bx7",128,-0.5,127.5);
00039 h->SetFillColor(kBlue);
00040
00041 hA[4]=new TH1F("mubX48v","L2WB-ET & primVertex vs. bXing; bXing= raw bx48",128,-0.5,127.5);
00042 hA[5]=h=new TH1F("mubX7v","L2WB-ET & primVertex; bXing= raw bx7",128,-0.5,127.5);
00043 h->SetFillColor(kBlue);
00044
00045
00046 hA[6]=h=new TH1F("muDsm1","L2WB-ET events DMS spectrum; DSM value",64,-0.5,63.5);
00047 h->SetMinimum(0.8);
00048 hA[7]=h=new TH1F("muDsm2","L2WB-Rnd events DMS spectrum; DSM value",64,-0.5,63.5);
00049 h->SetMinimum(0.8);
00050 sprintf(txt,"L2WB-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,"L2WB-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
00058 hA[10]=h=new TH1F("muVRf","L2WB: 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","L2WB: 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","L2WB: # vertices per event, rank>0 & Z in range; # of vertices",10,0,10);
00069
00070
00071 hA[20]=h=new TH1F("muStatTrk","Barrel 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]={"flag","pt1","nHit","Hfrac","Rin","Rout","ptOK","@B","CL","fr24",
00075 "#Delta R","eta1","noNear","noAway","goldW"};
00076
00077 for(int i=0;i<15;i++) h->Fill(keyT[i],0.);
00078
00079
00080 hA[21]=h=new TH1F("muTrNfit","Barrel: 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","Barrel: 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","Barrel: 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","Barrel: 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","Barrel: 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","Barrel: 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","Barrel: global NEGATIVE track PT; track PT (GeV/c)",160,0,80);
00104 h->SetFillColor(7);
00105
00106 hA[28]=h=new TH2F("muTrdEdX","Barrel: dEdX vs. momentum; track P (GeV); dE/dx (keV)",20,0,10,100,0,10);
00107
00108
00109 hA[29]=h=new TH1F("muTrPt1Pr","Barrel: primary track PT; track PT (GeV/c)",160,0,80);
00110 hA[30]=h=new TH1F("muTrPt1NPr","Barrel: primary NEGATIVE track PT; track PT (GeV/c)",160,0,80);
00111
00112
00113
00114
00115 hA[31]=h=new TH1F("muBmaxAdc","Barrel: BTOW maxADC in event, in-selection; max tower ADC",200,0,5500);
00116 hA[32]=h=new TH1F("muBtotAdc","Barrel: 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","Barrel: 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","Barrel: track glob chi2/dof X-Y",100,0,5);
00125 hA[36]=new TH2F("muTrch2b","Barrel: track glob chi2/dof; chi2/dof X-Y; last hit eta",30,0,5.,30,-1,1);
00126
00127
00128
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","Barrel: 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","Barrel: 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
00137
00138 hA[40]=h=new TH1F("muBjetET","Barrel: near 'EM jet' ET ; 'EM jet' ET (GeV)",100,0,100);
00139
00140 hA[41]=h=new TH2F("muBclEjetE2D","Barrel: 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","Barrel: 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
00148 hA[43]=h=new TH2F("muBdist1","Barrel: 3D Distance(track-cluster) vs. 2x2 E;2x2 cluster E (GeV); | distance | (cm)",50,0,80,50,0,25);
00149 hA[44]=h=new TH2F("muBdist2","Barrel: #Delta Z (track-cluster) vs.Z-clust; Z-cluster (cm); #Delta Z (cm)",100,-300,300,40,-20,20);
00150 hA[45]=h=new TH2F("muBdist3","Barrel: R#Delta #phi (track-cluster) vs. 2x2 E;2x2 cluster E (GeV); R#Delta #phi (cm)",50,0,80,80,-20,20);
00151 hA[46]=h=new TH1F("muBdist4","Barrel: 3D Distance(track-cluster) vs. 2x2 E;| 3D distance | (cm)",100,0,50);
00152 Lx=h->GetListOfFunctions();
00153 ln=new TLine(par_delR3D,0,par_delR3D,1.e6); ln->SetLineColor(kRed); Lx->Add(ln);
00154
00155
00156 hA[47]=h=new TH1F("muTjetET","Barrel: quenched near 'TPC jet' PT ; 'TPC jet' PT (GeV)",100,0,100);
00157 hA[48]=h=new TH2F("muTjetBjet2D","Barrel: quenched near 'TPC jet' vs. near 'EMC jet' ; EMC jet ET (GeV); TPC PT (GeV/c)",50,0,80,50,0,60);
00158 hA[49]=h=new TH1F("muTBjetET","Barrel: near 'TPC+EMC jet' ET ; jet ET (GeV), no double counting",100,0,100);
00159
00160
00161 hA[50]=h=new TH1F("muTwayET","Barrel: TPC away-cone PT sum; PT (GeV)",100,0,100);
00162 hA[51]=h=new TH1F("muBwayET","Barrel: BTOW away-cone ET sum; ET (GeV)",100,0,100);
00163 hA[52]=h=new TH2F("muTotwayET2D","Barrel: away TPC+EMC ET sum vs. 2x2 cluster ET; 2x2 ET (GeV); away ET (GeV)",50,0,100,150,0,100);
00164
00165
00166 hA[53]=h=new TH2F("muAwayET2Db","Barrel: away EMC ET sum vs. 2x2 cluster ET; 2x2 ET (GeV); away EMC ET (GeV)",50,0,100,80,0,80);
00167
00168
00169 hA[54]=h=new TH1F("muAwayTotEt","Barrel: away-cone TPC+EMC ET sum ; away ET (GeV)",200,0,100);
00170
00171 hA[55]=h=new TH1F("muEwayET","Barrel: ETOW away-cone ET sum; ET (GeV)",100,0,100);
00172
00173 hA[57]=h=new TH2F("muTr2D1pt5","Barrel: lastHit on track (pt > 5); detector eta ; detector phi (rad)",100,-1.1,1.1,240,-PI,PI);
00174
00175
00176 hA[58]=new TH1F("muTrch2West","Barrel: track glob chi2/dof West TPC ; chi2/dof",100,0,5);
00177 hA[59]=new TH1F("muTrch2East","Barrel: track glob chi2/dof East TPC ; chi2/dof",100,0,5);
00178
00179
00180
00181 hA[60]=h=new TH2F("muBclETPt","Barrel: Awayside TPC PT vs. isolated cluster 2x2 ET, matched;2x2 cluster ET (GeV) ; Awayside TPC PT (GeV)",50,0,100,75,0,150);
00182
00183 hA[61]=new TH1F("mubX7bht","L2BW-BHT-rnd events vs. bXing; bXing= raw bx7",128,-0.5,127.5);
00184 hA[62]=h=new TH2F("muEne_Deta","Barrel W: cluster energy vs. detector eta, final selection; barrel eta bin; 2x2 Energy (GeV)",40,0,40,50,0,100);
00185
00186
00187
00188
00189 const char cPlane[ mxBSmd]={'E','P'};
00190 for(int iep=0;iep<mxBSmd;iep++){
00191 sprintf(txt0,"muS%cadc1",cPlane[iep]);
00192 sprintf(txt,"BSMD-%c whole plane ADC; ADC",cPlane[iep]);
00193 hA[70+10*iep]=h=new TH1F(txt0,txt,275,0,1100);
00194 }
00195
00196
00197 hA[90]=h=new TH1F("muWET","Barrel W: Final Selection; 2x2 cluster ET (GeV)", 100,0,100);
00198 Lx=h->GetListOfFunctions();
00199 ln=new TLine(par_highET,0,par_highET,1.e6); ln->SetLineColor(kRed); Lx->Add(ln);
00200
00201 sprintf(txt,"Barrel W: Final selection, ET>%.0f GeV 'goldenW'; detector eta ; detector phi (rad)",par_highET);
00202 hA[91]=new TH2F("muW2D1",txt,10,-1.0,1.0,24,-PI,PI);
00203
00204 hA[92]=new TH2F("muWdedx","Barrel W: Track dEdx, final selection; 2x2 ET (GeV); dEdx (keV)",100,0,100,100,0,10);
00205 hA[93]=new TH2F("muWglDca","Barrel W: Track glob vertex abs(DCA), final selection ; 2x2 ET (GeV); |DCA| (cm)",100,0,100,100,0,5);
00206 hA[94]=new TH2F("muWglDcaSP","Barrel W: Track prim POSITIVE glob signed DCA, final selection; 2x2 ET (GeV); sDCA (cm)",100,0,100,100,-5,5);
00207 hA[95]=new TH2F("muWglDcaSN","Barrel W: Track prim NEGATIVE glob signed DCA, final selection ; 2x2 ET (GeV); sDCA (cm)",100,0,100,100,-5,5);
00208
00209 sprintf(txt,"Barrel W: Vertex ID, final selection, 2x2 ET>%.0f GeV; vertex ID",par_highET);
00210 hA[96]=h=new TH1F("muWcar1",txt,10,-0.5,9.5);
00211 h->SetLineColor(kBlue);h->SetLineWidth(3);
00212
00213 sprintf(txt,"Barrel W: Vertex 'funny' rank, final W selection, 2x2 ET>%.0f GeV; X=Log10(rank)+offset",par_highET);
00214 hA[97]=new TH1F("muWcar2",txt, 150, -9,25);
00215
00216 sprintf(txt,"Barrel W: Vertex Z, final selection, 2x2 ET>%.0f GeV; Z(cm)",par_highET);
00217 hA[98]=new TH1F("muWcar3",txt, 100, -200,200);
00218 hA[99]=h=new TH1F("muWeta","Barrel W: lepton eta final selection; lepton eta",100, -1.5,1.5);
00219
00220
00221
00222
00223 char tt2[][200]={"max 2x2","track matched","no near ET","no away ET"};
00224 for(int i=0;i<4;i++){
00225 sprintf(txt,"Barrel electron candidate, cut=%s; 2x2 ET (GeV)",tt2[i]);
00226 sprintf(txt0,"muETlive%d",i);
00227 hA[110+i]=h=new TH1F(txt0,txt, 100,0,100);
00228 }
00229
00230
00231
00232 hA[117]=h=new TH2F("mujetQAeta_phi","Input Jet phi vs eta ; eta ; phi ",50,-3,3,63,-PI,PI);
00233 hA[118]=h=new TH1F("mujetQApt","Input Jet pt; pt;",100,0,100);
00234
00235 hA[132]=h=new TH2F("muptBalance_clust","ptBalance vs cluster ET; 2x2 Cluster ET; ptBalance",100,0,100,100,0,100);
00236 hA[133]=h=new TH2F("muptBalance_awayTot","ptBalance vs awayside PT; awayside PT; ptBalance",100,0,100,100,0,100);
00237
00238 hA[134]=h=new TH2F("musPtBalance_clust","Barrel: sPtBalance vs cluster ET; 2x2 Cluster ET (GeV); signed Pt balance (GeV)",100,0,100,100,-100,100);
00239 Lx=h->GetListOfFunctions();
00240 ln=new TLine(0,par_ptBalance,100,par_ptBalance); ln->SetLineColor(kRed); Lx->Add(ln);
00241
00242 hA[135]=h=new TH2F("musPtBalance_awayTot","Barrel: sPtBalance vs awayside PT; awayside PT; sPtBalance",100,0,100,100,-100,100);
00243 Lx=h->GetListOfFunctions();
00244 ln=new TLine(0,par_ptBalance,100,par_ptBalance); ln->SetLineColor(kRed); Lx->Add(ln);
00245
00246
00247 hA[136]=h=new TH1F("muclustPtBal",Form("Barrel: PT Balance > %.1f ; 2x2 Cluster ET",par_ptBalance),100,0,100);
00248 hA[137]=h=new TH1F("muclustPtBal_bckgrd",Form("Barrel: PT Balance < %.1f ; 2x2 Cluster ET",par_ptBalance),100,0,100);
00249 hA[140]=h=new TH1F("muclustPtBalnoE",Form("Barrel: sPT Balance > %.1f (EEMC not included); 2x2 Cluster ET",par_ptBalance),100,0,100);
00250
00251
00252 char str[200];
00253 for (int i=0; i<=20; i++) {
00254 sprintf(str,"neg_failAwaySide_Awayside_pt_bin_%d",i);
00255 hA[142+i] = new TH2F(str,str,100,0,100,21,0,21);
00256 }
00257
00258 for (int i=0; i<=20; i++) {
00259 sprintf(str,"pos_failAwaySide_Awayside_pt_bin_%d",i);
00260 hA[163+i] = new TH2F(str,str,100,0,100,21,0,21);
00261 }
00262
00263 hA[184+2] = new TH1F("pos_muclustpTbal_wE","Barrel: pos_muclustpTbal_wE",100,0,100);
00264 hA[184+1] = new TH1F("neg_muclustpTbal_wE","Barrel: neg_muclustpTbal_wE",100,0,100);
00265 hA[184+4] = new TH1F("pos_muclustpTbal_noE","Barrel: pos_muclustpTbal_noE",100,0,100);
00266 hA[184+3] = new TH1F("neg_muclustpTbal_noE","Barrel: neg_muclustpTbal_noE",100,0,100);
00267 hA[184+6] = new TH1F("pos_muclustpTbal_back","Barrel: pos_muclustpTbal_back",100,0,100);
00268 hA[184+5] = new TH1F("neg_muclustpTbal_back","Barrel: neg_muclustpTbal_back",100,0,100);
00269
00270
00271 hA[190]=h=new TH2F("muclustEt_etaWp","Barrel W+: 2x2 Cluster ET vs. lepton eta, final selection; lepton eta in LAB; lepton 2x2 Cluster ET (GeV)",32,-2.,2.,60,0.,60.);
00272 hA[191]=h=new TH2F("muclustEt_etaWm","Barrel W-: 2x2 Cluster ET vs. lepton eta, final selection; lepton eta in LAB; lepton 2x2 Cluster ET (GeV)",32,-2.,2.,60,0.,60.);
00273
00274
00275 hA[195]=h=new TH1F("muEtaLT0_maxTowADC","max tower ADC for eta < 0",4096,0,4096);
00276 hA[196]=h=new TH1F("muEtaGT0_maxTowADC","max tower ADC for eta > 0",4096,0,4096);
00277 hA[197]=h=new TH2F("muTr2D1pt10","lastHit on track (pt > 10); detector eta ; detector phi (rad)",100,-1.1,1.1,240,-PI,PI);
00278 hA[198]=h=new TH2F("muTrPt_eta","Track PT vs lastHit eta; detector eta ; prim track pt (GeV)",100,-1.1,1.1,100,0,100);
00279 hA[199]=h=new TH2F("muBdist5","Barrel: 3D Distance(track-cluster) vs. cluster eta; cluster eta; | distance | (cm)",100,-1.,1.,50,0,25);
00280 hA[200]=h=new TH2F("muWClustET_eta","W: 2x2 cluster ET vs. cluster eta; cluster detector eta; 2x2 cluster ET",100,-1.,1.,60,0,60);
00281 hA[201]=h=new TH2F("muWTrPT_eta","W: Track PT vs. cluster eta; cluster detector eta; Track PT",100,-1.,1.,100,0,100);
00282 hA[202]=h=new TH2F("muBTrPT_eta","Background: Track PT vs. cluster eta; cluster detector eta; Track PT",100,-1.,1.,100,0,100);
00283 hA[203]=h=new TH2F("muWE2P_eta","W: E/P vs. cluster eta; cluster detector eta; E/P",100,-1.,1.,15,0,3);
00284 hA[204]=h=new TH2F("muBE2P_eta","Background: E/P vs. cluster eta; cluster detector eta; E/P",100,-1.,1.,15,0,3);
00285 hA[205]=h=new TH2F("muWTr2D1","lastHit on track (W candidate); detector eta ; detector phi (rad)",100,-1.1,1.1,240,-PI,PI);
00286 hA[206]=h=new TH2F("muWClustET_eta_clgt14","2x2 cluster ET vs. cluster eta (Cluster ET > 14); cluster detector eta; 2x2 cluster ET",100,-1.,1.,60,0,60);
00287 hA[207]=h=new TH2F("muWClustET_eta_towerIso","2x2 cluster ET vs. cluster eta (2x2/4x4 frac); cluster detector eta; 2x2 cluster ET",100,-1.,1.,60,0,60);
00288 hA[208]=h=new TH2F("muWClustET_eta_delR","2x2 cluster ET vs. cluster eta (track-clust match); cluster detector eta; 2x2 cluster ET",100,-1.,1.,60,0,60);
00289 hA[209]=h=new TH2F("muWClustET_eta_nearIso","2x2 cluster ET vs. cluster eta (2x2/near frac); cluster detector eta; 2x2 cluster ET",100,-1.,1.,60,0,60);
00290
00291
00292 hA[210]=h=new TH2F("muETOWnearET_eta_tr2cl","ETOW near cone ET vs. cluster eta; cluster detector eta; ETOW near ET",100,-1.,1.,20,0,20);
00293 hA[211]=h=new TH2F("muETOWnearET_phi_etaGT0_tr2cl","ETOW near cone ET vs. cluster phi (cluster eta > 0); cluster detector phi; ETOW near ET",240,-PI,PI,20,0,20);
00294 hA[212]=h=new TH2F("muETOWnearET_phi_etaLT0_tr2cl","ETOW near cone ET vs. cluster phi (cluster eta < 0); cluster detector phi; ETOW near ET",240,-PI,PI,20,0,20);
00295
00296
00297 string bx7name[4]={"[0,30]","[30,40]","[40,110]","[110,120]"};
00298 for(int i=0; i<4; i++){
00299 hA[215+i]=h=new TH2F(Form("muBTOW_adcGT10goodVer_7bx%d",i),Form("BTOW tower ADC > 10 (after good vertex found) bx7=%s; detector eta; detector phi",bx7name[i].data()),40,-1,1,120,-PI,PI);
00300 hA[219+i]=h=new TH2F(Form("muBTOW_etGT2goodVer_7bx%d",i),Form("BTOW tower ET > 2.0 GeV (after good vertex found) bx7=%s; detector eta; detector phi",bx7name[i].data()),40,-1,1,120,-PI,PI);
00301 hA[223+i]=h=new TH2F(Form("muBTOW_adcGT10noVer_7bx%d",i),Form("BTOW tower ADC > 10 (no rank>0 vertex found) bx7=%s; detector eta; detector phi",bx7name[i].data()),40,-1,1,120,-PI,PI);
00302 hA[227+i]=h=new TH2F(Form("muETOW_adcGT10goodVer_7bx%d",i),Form("ETOW tower ADC > 10 (after good vertex found) bx7=%s; detector eta bin; detector phi bin",bx7name[i].data()),12,0,12,60,0,60);
00303 hA[231+i]=h=new TH2F(Form("muETOW_etGT2goodVer_7bx%d",i),Form("ETOW tower ET > 2.0 GeV (after good vertex found) bx7=%s; detector eta bin; detector phi bin",bx7name[i].data()),12,0,12,60,0,60);
00304 hA[235+i]=h=new TH2F(Form("muETOW_adcGT10noVer_7bx%d",i),Form("ETOW tower ADC > 10 (no rank>0 vertex found) bx7=%s; detector eta bin; detector phi bin",bx7name[i].data()),12,0,12,60,0,60);
00305 }
00306
00307 hA[240]=h=new TH2F("muclustEt_etaFinal_noE","Barrel: 2x2 Cluster ET vs. lepton eta, final selection (no EEMC in veto); lepton eta in LAB; lepton 2x2 Cluster ET (GeV)",32,-2.,2.,60,0.,60.);
00308 hA[241]=h=new TH2F("muclustEt_etaFinal","Barrel: 2x2 Cluster ET vs. lepton eta, final selection; lepton eta in LAB; lepton 2x2 Cluster ET (GeV)",32,-2.,2.,60,0.,60.);
00309
00310
00311 for(int i=0;i<mxHA;i++) {
00312 if( hA[i]==0) continue;
00313 HList->Add( hA[i]);
00314 }
00315
00316
00317 LOG_INFO<<Form("%s::initHistos done",GetName())<<endm;
00318
00319 }
00320
00321
00322
00323
00324
00325
00326
00327
