StRoot: StSpinPool/StWalgoB2009/St2009pubSpinMaker.cxx Source File

00001 // $Id: St2009pubSpinMaker.cxx,v 1.12 2010/05/01 01:31:45 balewski Exp $
00002 //
00003 //*-- Author : Jan Balewski, MIT
00004 // 
00005 #include <StMuDSTMaker/COMMON/StMuDstMaker.h>
00006 #include <StMuDSTMaker/COMMON/StMuDst.h>
00007 #include <StMuDSTMaker/COMMON/StMuEvent.h>
00008 #include <StEvent/StEventInfo.h> // just to get time
00009 #include "StSpinPool/StSpinDbMaker/StSpinDbMaker.h"
00010 
00011 #include "St2009WMaker.h"
00012 
00013 #include "St2009pubSpinMaker.h"
00014 
00015 ClassImp(St2009pubSpinMaker)
00016 
00017 //_____________________________________________________________________________
00018 //
00019 St2009pubSpinMaker::St2009pubSpinMaker(const char *name):StMaker(name){
00020   wMK=0;HList=0;
00021   core=name;
00022   par_QPTlow=0.010;
00023 
00024   par_QPThighET0=25; 
00025   par_QPThighET1=50; 
00026   par_QPThighA=0.08; 
00027   par_QPThighB=0.0013; 
00028   par_leptonEta1=-1.; par_leptonEta2=1.;
00029   par_useNoEEMC=0;
00030  }
00031 
00032 
00033 //_____________________________________________________________________________
00034 //
00035 Int_t 
00036 St2009pubSpinMaker::Init(){
00037   assert(wMK);
00038   assert(HList);
00039   initHistos();
00040   nRun=0;
00041   hbxIdeal=0;
00042   return StMaker::Init();
00043 }
00044 
00045 
00046 //_____________________________________________________________________________
00047 //
00048 Int_t 
00049 St2009pubSpinMaker::FinishRun  (int runNo){
00050   char txt[1000];
00051   sprintf(txt,"events T= %d %d",Tfirst,Tlast);
00052   printf("Finish run=%d , events time range %s\n",runNo,txt);
00053   hbxIdeal->GetYaxis()->SetTitle(txt);
00054   return kStOK;
00055 }
00056 
00057 //_____________________________________________________________________________
00058 //
00059 Int_t 
00060 St2009pubSpinMaker::InitRun  (int runNo){
00061   assert(runNo>= 10081007); // F10407, first pp500 long fill with defined pol pattern
00062   assert(runNo<=10103046); // F10536, last pp500 fill in run9
00063 
00064   char txt[1000],txt0[100];
00065   sprintf(txt0,"bxIdeal%d",nRun);
00066   sprintf(txt,"intended fill pattern  R%d-%d vs. bXing; %s", runNo,nRun,spinDb->getV124comment());
00067   nRun++;
00068   Tfirst=int(2e9); Tlast=-Tfirst;
00069   hbxIdeal=new TH1F(core+txt0,txt,128,-0.5,127.5);
00070   hbxIdeal->SetFillColor(kYellow);
00071   HList->Add(hbxIdeal);
00072   
00073   spinDb->print(0); // 0=short, 1=huge
00074   for(int bx=0;bx<120;bx++){
00075     if(spinDb->isBXfilledUsingInternalBX(bx))  hbxIdeal->Fill(bx);   
00076   }
00077 
00078   sprintf(txt,"bXing= bx48+off=%d",spinDb->BX48offset());
00079   hA[3]->GetXaxis()->SetTitle(txt);
00080 
00081   sprintf(txt,"bXing= bx7+off=%d",spinDb->BX7offset());
00082   hA[4]->GetXaxis()->SetTitle(txt);
00083 
00084   LOG_INFO<<Form("::InitRun(%d) done, W-spin sorting  params: exclude |Q/PT| < %.2f OR |Q/PT| above line %.3f*(ET-%.1f)-%6e if ET<%.1f, for AL use leptonEta in[%.1f,%.1f] useNoEEMC=%d", runNo,
00085                  par_QPTlow,par_QPThighET0, par_QPThighA ,par_QPThighB,par_QPThighET1,par_leptonEta1, par_leptonEta2,par_useNoEEMC
00086                  )<<endm;        
00087   return kStOK;
00088 }
00089 
00090 //_____________________________________________________________________________
00091 //
00092 Int_t 
00093 St2009pubSpinMaker::Make(){
00094   int T=wMK->mMuDstMaker->muDst()->event()->eventInfo().time();
00095   if(Tlast<T) Tlast=T;
00096   if(Tfirst>T) Tfirst=T;
00097 
00098   bXingSort();
00099   return kStOK;
00100 }
00101 
00102 //_____________________________________________________________________________
00103 //
00104 void 
00105 St2009pubSpinMaker::bXingSort(){
00106   //has access to whole W-algo-maker data via pointer 'wMK'
00107   
00108   hA[0]->Fill("inp",1.);
00109 
00110   assert(spinDb->isValid());  // all 3 DB records exist 
00111   assert(spinDb->isPolDirLong());  // you do not want mix Long & Trans by accident
00112   
00113   if(wMK->wEve.vertex.size()<=0) return; 
00114   //......... require: L2W-trig (ET or rnd) & vertex is reasonable .......
00115   
00116   int bx48=wMK->wEve.bx48;
00117   int bx7=wMK->wEve.bx7;
00118   if(spinDb->offsetBX48minusBX7(bx48,bx7)) {
00119    printf("BAD bx7=%d bx48=%d del=%d\n",bx7,bx48,spinDb->offsetBX48minusBX7(bx48,bx7));
00120    hA[0]->Fill("badBx48",1.);
00121    return; // both counters must be in sync
00122   }
00123 
00124   //remove events tagged as Zs
00125   if(wMK->wEve.zTag) return;
00126   hA[0]->Fill("noZ",1.);
00127  
00128 
00129   hA[1]->Fill(bx48);
00130   hA[2]->Fill(bx7);
00131 
00132   int bxStar48= spinDb->BXstarUsingBX48(bx48);
00133   int bxStar7=spinDb->BXstarUsingBX7(bx7);
00134   hA[3]->Fill(bxStar48);
00135   hA[4]->Fill(bxStar7);
00136 
00137   int spin4=spinDb->spin4usingBX48(bx48); 
00138   hA[5]->Fill(bxStar7,spin4);
00139 
00140   float par_maxDsmThr=58;
00141   float par_myET=25; // monitoring cut
00142   if( wMK->wEve.l2bitRnd) { // lumi monitor BHT3-random
00143     StMuEvent* muEve = wMK->mMuDstMaker->muDst()->event();
00144     int max=0;
00145     for (int m=0;m<300;m++)  {    // access L0-HT data
00146       int val=muEve->emcTriggerDetector().highTower(m);
00147       if(max<val) max=val;
00148     }
00149     // avoid too much energy - can be W-events (1/milion :)
00150     if(max<par_maxDsmThr)  {  hA[6]->Fill(spin4);  hA[0]->Fill("BG1",1.);}
00151     return;
00152   }
00153   
00154   if( wMK->wEve.l2bitET==0) return; 
00155   //..... it is guaranteed ..... L2W-ET>13 did fired  ......
00156   
00157   
00158   // search for  Ws ............
00159   for(uint iv=0;iv<wMK->wEve.vertex.size();iv++) {
00160     WeveVertex &V=wMK->wEve.vertex[iv];
00161     for(uint it=0;it<V.eleTrack.size();it++) {
00162       WeveEleTrack &T=V.eleTrack[it];
00163       if(T.pointTower.id==0) continue;
00164 
00165       /* Collect QCD background for lumi monitors */
00166       float frac24=T.cluster.ET/(T.cl4x4.ET);
00167       if(iv==0 && it==0 && frac24<wMK->par_clustFrac24) {
00168         hA[31]->Fill(T.cluster.ET);
00169         if( T.cluster.ET <20. ) { hA[7]->Fill(spin4);  hA[0]->Fill("BG2",1.);}
00170       }
00171 
00172       if(T.isMatch2Cl==false) continue;
00173       assert(T.cluster.nTower>0); // internal logical error
00174       assert(T.nearTotET>0); // internal logical error
00175 
00176       int iQ=0; // plus
00177       float p_Q=T.prMuTrack->charge();
00178       if( p_Q<0 ) iQ=1;// minus
00179       float ET=T.cluster.ET;
00180       
00181 
00182       //put final W cut here
00183       bool isW= T.cluster.ET /T.nearTotET> wMK->par_nearTotEtFrac;  // near cone
00184       if(par_useNoEEMC) 
00185         isW=isW && T.sPtBalance_noEEMC>wMK->par_ptBalance; // awayET
00186       else
00187         isW=isW && T.sPtBalance>wMK->par_ptBalance; // awayET
00188     
00189       if(!isW) { // AL(QCD)
00190         if(ET>15 &&ET<20 ) hA[16+iQ]->Fill(spin4);
00191         continue;
00192       }
00193 
00194       hA[0]->Fill("Wcut",1.);
00195 
00196       hA[30]->Fill(T.prMuTrack->eta());
00197       if(fabs(T.primP.Eta()) > wMK->par_leptonEta) continue;     
00198       // allows further cuts on eta
00199       if(T.prMuTrack->eta()<par_leptonEta1) continue;
00200       if(T.prMuTrack->eta()>par_leptonEta2) continue;
00201       hA[0]->Fill("eta",1.);
00202 
00203       //::::::::::::::::::::::::::::::::::::::::::::::::
00204       //:::::accepted W events for x-section :::::::::::
00205       //::::::::::::::::::::::::::::::::::::::::::::::::
00206 
00207       if(ET>par_myET) hA[0]->Fill("W25",1.);
00208       float q2pt=T.prMuTrack->charge()/T.prMuTrack->pt();
00209       if(ET>par_myET) hA[8]->Fill(q2pt);
00210       hA[9]->Fill(ET,q2pt);
00211       
00212       // apply cut on reco charge
00213       if( fabs(q2pt)< par_QPTlow) continue;
00214       if(ET>par_myET) hA[0]->Fill("Qlow",1.);
00215 
00216       if(par_QPTlow>0) { // abaility to skip all Q/PT cuts
00217         if( fabs(q2pt)< par_QPTlow) continue;
00218         float highCut=par_QPThighA - (ET-par_QPThighET0)*par_QPThighB;
00219         // printf("fff ET=%f q2pr=%f highCut=%f passed=%d\n",ET, q2pt,highCut,fabs(q2pt)<highCut);
00220         if( ET>par_myET && ET<par_QPThighET1 && fabs(q2pt)>highCut) continue;
00221       }
00222 
00223       if(ET>par_myET) {
00224         hA[0]->Fill("Qhigh",1.);
00225         if(p_Q>0) hA[0]->Fill("Q +",1.);
00226         else  hA[0]->Fill("Q -",1.);
00227       }
00228 
00229      
00230       hA[10+iQ]->Fill(ET);
00231       if(ET>25 &&ET<50 ) hA[12+iQ]->Fill(spin4);
00232       if(ET>32 &&ET<44 ) hA[14+iQ]->Fill(spin4);
00233      
00234       hA[18+iQ]->Fill(spin4,ET);         
00235      
00236     } // end of loop over tracks
00237   }// end of loop ove vertices
00238 
00239 }
00240 
00241 
00242 // $Log: St2009pubSpinMaker.cxx,v $
00243 // Revision 1.12  2010/05/01 01:31:45  balewski
00244 // added W->JJ code & JES calibration
00245 //
00246 // Revision 1.11  2010/04/14 20:00:08  balewski
00247 // added AL w/o endcap
00248 //
00249 // Revision 1.10  2010/03/22 17:18:32  balewski
00250 // now the < > sign is right
00251 //
00252 // Revision 1.9  2010/03/22 16:11:42  balewski
00253 // better computation of AL(QCD)
00254 //
00255 // Revision 1.8  2010/03/20 19:19:05  balewski
00256 // added ability to drop Q/PT cut for spin analysis
00257 //
00258 // Revision 1.7  2010/03/15 17:05:46  balewski
00259 // cleanup, used for W AL sort March 15, 2010
00260 //
00261 // Revision 1.6  2010/03/14 22:50:31  balewski
00262 // *** empty log message ***
00263 //
00264 // Revision 1.5  2010/02/04 03:48:12  balewski
00265 // add ET for lumi monitor
00266 //
00267 // Revision 1.4  2010/01/28 20:10:05  balewski
00268 // added eta dependent spin sorting
00269 //
00270 // Revision 1.3  2010/01/28 03:42:55  balewski
00271 // cleanup
00272 //
00273 // Revision 1.2  2010/01/27 22:12:25  balewski
00274 // spin code matched to x-section code
00275 //
00276 // Revision 1.1  2009/11/23 23:00:18  balewski
00277 // code moved spin-pool
00278 //
00279 // Revision 1.1  2009/11/23 21:11:18  balewski
00280 // start
00281 //