StRoot: StEEmcPool/StEEsoloPi0/StEEsoloPi0Maker.cxx Source File

00001 // *-- Author : Jan Balewski
00002 // 
00003 // $Id: StEEsoloPi0Maker.cxx,v 1.14 2011/04/11 19:35:42 fisyak Exp $
00004 
00005 #include <TFile.h>
00006 
00007 #include "StEEsoloPi0Maker.h"
00008 
00009 #include "TChain.h"
00010 #include "TClonesArray.h"
00011 #include "StL0Trigger.h"
00012 #include "StEventInfo.h"
00013 #include "StEventSummary.h"
00014 #include "StarClassLibrary/StThreeVectorF.hh"
00015 #include "StMuDSTMaker/COMMON/StMuEvent.h"
00016 #include "StMuDSTMaker/COMMON/StMuTrack.h"
00017 #include "StMuDSTMaker/COMMON/StMuDstMaker.h"
00018 #include "StMuDSTMaker/COMMON/StMuTriggerIdCollection.h"
00019 #include "StEvent/StTriggerId.h"
00020 
00021 
00022 #include "StEEmcUtil/database/EEmcDbItem.h"
00023 #include "StEEmcUtil/database/StEEmcDb.h"
00024 
00025 #include <StMessMgr.h>
00026 
00027 ClassImp(StEEsoloPi0Maker)
00028 
00029 //________________________________________________
00030 //________________________________________________
00031 StEEsoloPi0Maker::StEEsoloPi0Maker(const char* self ,const char* muDstMakerName) : StMaker(self){
00032   mMuDstMaker = (StMuDstMaker*)GetMaker(muDstMakerName);
00033   assert(mMuDstMaker);
00034   MCflag=0;// default off
00035 
00036 }
00037 
00038 
00039 //___________________ _____________________________
00040 //________________________________________________
00041 StEEsoloPi0Maker::~StEEsoloPi0Maker(){
00042  // Save all objects in this file
00043   // hfile->Write();
00044 }
00045 
00046 //________________________________________________
00047 //________________________________________________
00048 Int_t StEEsoloPi0Maker::InitRun(int runNo){
00049   // call initRun() only once for M-C
00050 
00051   if( !MCflag) initRun(runNo);
00052   static int first=1;
00053   if(first) return kStOK;
00054   initRun(runNo);
00055   first=0;
00056   return kStOK;
00057 }
00058  
00059 //________________________________________________
00060 //________________________________________________
00061 Int_t StEEsoloPi0Maker::Init(){
00062   eeDb = (StEEmcDb*)this->GetDataSet("StEEmcDb");
00063   EEsoloPi0::init();
00064   LOG_INFO << "has MCflag="<< MCflag<<endm;
00065   return StMaker::Init();
00066 }
00067 
00068 //________________________________________________
00069 //________________________________________________
00070 Int_t StEEsoloPi0Maker::Finish(){
00071   finish();
00072   return kStOK;
00073 }
00074  
00075  
00076  
00077 
00078 //________________________________________________
00079 //________________________________________________
00080 Int_t StEEsoloPi0Maker::Make(){
00081   clear();  
00082   static int n0=0,n1=0,n2=0,n3=0;
00083   // printf("%s::Make() is called ..........n0,1,2,3= %d %d %d %d \n",StMaker::GetName(),n0,n1,n2,n3);
00084   n0++;
00085   //............. trigger sort
00086   if( !MCflag&& !unpackMuTrig()) return kStOK;
00087   n1++;
00088 
00089   //................. CTB sort
00090   float sum=getCtbSum();  sum=sum;
00091   //tmp  if(sum<75 || sum > 800) return kStOK; 
00092   n2++;
00093     
00094   // ............. acquire EEMC data 
00095   if(!unpackMuEemc())   return kStOK;
00096   n3++;
00097 
00098   findTowerClust();
00099   findTowerPi0();
00100 
00101   return kStOK;
00102 }
00103 
00104 
00105 
00106 //________________________________________________
00107 //________________________________________________
00108 bool StEEsoloPi0Maker::unpackMuTrig(){ // trigger filter, has wrongID
00109   
00110   //  printf("%s::unpackMuTrig() is called ..........\n",StMaker::GetName());
00111  
00112   // Access to muDst .......................
00113   StMuEvent* muEve = mMuDstMaker->muDst()->event();
00114   
00115   StMuTriggerIdCollection& trgIdColl=muEve->triggerIdCollection();
00116 
00117   const StTriggerId& oflTrgId=trgIdColl.nominal();
00118   vector<unsigned int> trgId=oflTrgId.triggerIds();
00119 
00120 #if 0
00121   StEventInfo &info=muEve->eventInfo();
00122   int nPrim = mMuDstMaker->muDst()->primaryTracks()->GetEntries(); 
00123   printf("\n\n ==================== processing eventID %d nPrim=%d nTrig=%d==============\n", info.id(),nPrim, trgId.size());
00124 #endif
00125 
00126   bool isGood=false;
00127   UInt_t i;
00128   for(i = 0; i < trgId.size() ; i++){
00129     //    printf("i=%d trgId=%d\n",i,trgId[i]);
00130 #if 0
00131     //.......... minB trig in pp200 in 2004
00132     if(trgId[i]==10) isGood=true;
00133     if(trgId[i]==45010) isGood=true;
00134     if(trgId[i]==45020) isGood=true;
00135 #endif 
00136     //.......... some trigs in ppTrans in 2006 
00137     if(trgId[i]==127641) isGood=true; //e-http-l2gam
00138     // if(trgId[i]==127652) isGood=true; //e-jp0-l2jet
00139     // if(trgId[i]==127551) isGood=true; //e-jp0
00140   }
00141   
00142 #if 0 // TPC vertex, not used (yet)
00143   StEventSummary &smry=muEve->eventSummary();
00144   StThreeVectorF ver=smry.primaryVertexPosition();
00145 
00146   b.zTpc=ver.z();
00147   if( fabs(ver.x())<0.000001 &&fabs(ver.y())<0.000001 &&fabs(ver.z())<0.000001 )  b.zTpc=999;
00148 #endif
00149 
00150   return isGood;
00151 }
00152 
00153 //________________________________________________
00154 //________________________________________________
00155 bool StEEsoloPi0Maker::unpackMuEemc(){
00156   
00157   nInpEve++;
00158   gMessMgr->Debug() <<GetName()<<"::unpackMuDst() is called "<<endm;
00159 
00160   // Access to muDst .......................
00161   StMuEmcCollection* emc = mMuDstMaker->muDst()->muEmcCollection();
00162   if (!emc) {
00163     gMessMgr->Warning() <<"No EMC data for this event"<<endm;    return false;
00164   }
00165 
00166   int i, n1=0;
00167   //printf("aaa %d %d \n",eeDb->mfirstSecID,eeDb->mlastSecID);
00168   //.........................  T O W E R S .....................
00169   for (i=0; i< emc->getNEndcapTowerADC(); i++) {
00170     int sec,eta,sub,val;
00171     //muDst  ranges:sec:1-12, sub:1-5, eta:1-12
00172     emc->getEndcapTowerADC(i,val,sec,sub,eta);
00173     assert(sec>0 && sec<=MaxSectors);// total corruption of muDst
00174 
00175     //tmp, for fasted analysis use only hits from sectors init in DB
00176     if(sec<eeDb->getFirstSector() || sec>eeDb->getLastSector()) continue;
00177  
00178     const EEmcDbItem *x=eeDb->getTile(sec,'A'+sub-1,eta,'T');
00179     assert(x); // it should never happened for muDst
00180 
00181     if(x->fail ) continue; // drop broken channels
00182 
00183     int iphi=(x->sec-1)*MaxSubSec+(x->sub-'A');
00184     int ieta=x->eta-1;
00185     assert(iphi>=0 && iphi<MaxPhiBins);
00186     assert(ieta>=0 && ieta<MaxEtaBins);
00187     int irad=iphi*MaxEtaBins+ieta; // unified spiral index
00188     assert(irad>=0 && irad<EEsoloPi0::MxTw);
00189 
00190     float  ene=-102;
00191  
00192     float  rawAdc=val;
00193     if(rawAdc<x->thr) continue;
00194     float  adc=rawAdc-x->ped; 
00195     if(x->gain) ene=adc/x->gain;
00196 
00197     if(MCflag) ene/=0.8; //fug factor for sampling fraction in M-C being ~4% insted of assumed 5%
00198 
00199     int ii=(x->sec-1)*60+(x->sub-'A')*12+x->eta-1;
00200     assert(ii==irad);
00201     //    if(adc>0)
00202     // printf("adc=%f ene/GeV=%f rawAdc=%f idG=%f,hSec=%d %s\n",adc,ene,rawAdc,mfixEmTgain[ieta],x->sec,x->name); 
00203     
00204     //aa int iT=0;// for towers   
00205     //aa tileAdc[iT][ieta][iphi]=adc;// store towers
00206     //aa tileThr[iT][ieta][iphi]=rawAdc>x->thr;
00207 
00208     n1++;
00209     //aa tileEne[iT][ieta][iphi]=ene;
00210     float recoEner=ene/scaleFactor; // ideal
00211     // printf("new ii=%d ene=%f del=%f mcf=%d\n",idar,recoEner,recoEner-soloMip[irad].e,MCflag);
00212     soloMip[irad].e= recoEner;
00213     
00214   }// end of loop over towers
00215 
00216   gMessMgr->Debug() <<GetName()<<"::unpackMuDst(),  found total" <<n1<<" towers with ADC>thr"<<endm;
00217 
00218   return true;
00219 }
00220 
00221 //________________________________________________
00222 //________________________________________________
00223 float StEEsoloPi0Maker::getCtbSum(){
00224   
00225   // printf("%s::getCtbSum() is called ..........\n",StMaker::GetName());
00226  
00227   // Access to muDst .......................
00228   StMuEvent* muEve = mMuDstMaker->muDst()->event();
00229   StCtbTriggerDetector* ctbDet = &(muEve->ctbTriggerDetector());
00230   
00231   assert(ctbDet);
00232   float ctbSum = 0;
00233   int nHit=0;
00234   for (UInt_t slat = 0; slat < ctbDet->numberOfSlats(); slat++) {
00235     for (UInt_t tray = 0; tray < ctbDet->numberOfTrays(); tray++) {
00236       float  adc = ctbDet->mips(tray,slat,0);
00237       ctbSum += adc;      
00238       if(adc > 5) nHit++;
00239     }
00240   }
00241   // printf("CTB %d hit ADC>5  sumADC=%f (all) \n",nHit, ctbSum);
00242   hA[7]->Fill(ctbSum); 
00243   return ctbSum;
00244 }
00245 
00246 
00247 
00248 
00249 // $Log: StEEsoloPi0Maker.cxx,v $
00250 // Revision 1.14  2011/04/11 19:35:42  fisyak
00251 // Replace uint by UInt_t, use TMath
00252 //
00253 // Revision 1.13  2009/02/04 20:33:21  ogrebeny
00254 // Moved the EEMC database functionality from StEEmcDbMaker to StEEmcUtil/database. See ticket http://www.star.bnl.gov/rt2/Ticket/Display.html?id=1388
00255 //
00256 // Revision 1.12  2007/10/19 23:18:42  balewski
00257 // 2008 cleanup, now works only w/ regular muDst
00258 //
00259 // Revision 1.11  2004/10/27 18:07:50  balewski
00260 // practical use of 'sim' flavor for M-C
00261 //
00262 // Revision 1.10  2004/10/21 13:31:25  balewski
00263 // to match new name of emcCollection in muDst
00264 //
00265 // Revision 1.9  2004/09/29 18:04:44  balewski
00266 // now it runs on M-C as well
00267 //
00268 // Revision 1.8  2004/09/03 04:50:52  balewski
00269 // big clenup
00270 //
00271 // Revision 1.7  2004/08/26 04:39:40  balewski
00272 // towards pi0
00273 //
00274 // Revision 1.6  2004/08/17 15:46:56  balewski
00275 // clenup for pp200 in 2004
00276 //
00277 // Revision 1.5  2004/08/09 20:28:31  balewski
00278 // add trig selection
00279 //
00280 // Revision 1.4  2004/05/07 21:38:38  balewski
00281 // gamma finder with SMD
00282 //
00283 // Revision 1.3  2004/04/29 04:16:43  perev
00284 // typo fixed. getTile not getTail
00285 //
00286 // Revision 1.2  2004/04/14 19:34:01  balewski
00287 // access to trigger data
00288 //
00289 // Revision 1.1  2004/04/14 17:09:09  balewski
00290 // new copy of pi0finder with towers only, should work on ezTree as well (after small cleanup)
00291 //