StRoot: St_QA_Maker/StQAMakerBase.cxx Source File

00001 // $Id: StQAMakerBase.cxx,v 2.38 2012/03/05 03:42:32 genevb Exp $ 
00002 // $Log: StQAMakerBase.cxx,v $
00003 // Revision 2.38  2012/03/05 03:42:32  genevb
00004 // Remove TPC XY dist, add TPC RPhi charge
00005 //
00006 // Revision 2.37  2012/02/08 22:10:35  genevb
00007 // Updates for Run 12
00008 //
00009 // Revision 2.36  2011/05/26 19:59:38  genevb
00010 // Cleanup in destructors
00011 //
00012 // Revision 2.35  2009/11/19 20:34:38  genevb
00013 // Remove Event Summary (using defunct old software monitors)
00014 //
00015 // Revision 2.34  2009/03/27 21:18:36  genevb
00016 // Add Jet Patch trigger histograms
00017 //
00018 // Revision 2.33  2007/12/12 19:50:55  genevb
00019 // Update for trigger words
00020 //
00021 // Revision 2.32  2007/11/30 05:38:50  genevb
00022 // Changes for Run8: mostly silicon removal, TOF addition
00023 //
00024 // Revision 2.31  2007/04/24 00:33:38  genevb
00025 // Always do PMD hists now
00026 //
00027 // Revision 2.30  2005/02/22 19:38:39  genevb
00028 // Do PMD hists only for year 4 and later (real data)
00029 //
00030 // Revision 2.29  2005/02/08 17:22:46  genevb
00031 // PMD histo changes, handle estGlobal/ITTF tracks
00032 //
00033 // Revision 2.28  2004/12/13 15:52:37  genevb
00034 // Numerous updates: PMD, primtrk, FPD, QAShift lists
00035 //
00036 // Revision 2.27  2004/02/12 16:54:24  genevb
00037 // Separate MinBias histos
00038 //
00039 // Revision 2.26  2004/02/12 05:03:14  genevb
00040 // Year 4 AuAu changes. New SVT histos.
00041 //
00042 // Revision 2.25  2004/01/10 01:10:18  genevb
00043 // Preparations for Year 5, added some svt plots
00044 //
00045 // Revision 2.24  2003/04/19 00:17:50  genevb
00046 // Updated for dAu/pp running
00047 //
00048 // Revision 2.23  2003/02/28 16:01:09  genevb
00049 // Further improvements for previous check-in
00050 //
00051 // Revision 2.22  2003/02/28 06:17:56  genevb
00052 // Allow StQAMakerBase::Make to be called for all events
00053 //
00054 // Revision 2.21  2003/02/20 20:09:54  genevb
00055 // Several changes for new trigger scheme, dAu data
00056 //
00057 // Revision 2.20  2003/02/19 06:38:29  genevb
00058 // Rework trigger and mult/event class sections
00059 //
00060 // Revision 2.19  2003/02/15 22:00:52  genevb
00061 // Add tpcSectors, fix ftpc east/west charge
00062 //
00063 // Revision 2.18  2002/04/23 01:59:56  genevb
00064 // Addition of BBC/FPD histos
00065 //
00066 // Revision 2.17  2002/03/11 20:59:58  genevb
00067 // Fixed bug with placement of trigger hists
00068 //
00069 // Revision 2.16  2002/03/01 22:51:20  genevb
00070 // Small set-to-zero possible bug fix
00071 //
00072 // Revision 2.15  2002/02/12 18:42:00  genevb
00073 // Additional FTPC histograms
00074 //
00075 // Revision 2.14  2002/01/26 03:04:07  genevb
00076 // Fixed some problems with fcl histos
00077 //
00078 // Revision 2.13  2002/01/21 22:09:24  genevb
00079 // Include some ftpc histograms from StFtpcClusterMaker
00080 //
00081 // Revision 2.12  2001/12/28 09:19:13  genevb
00082 // Adjustments for pp running
00083 //
00084 // Revision 2.11  2001/12/20 03:11:08  genevb
00085 // pp trigger words 0x2XXX
00086 //
00087 // Revision 2.10  2001/09/10 18:00:13  genevb
00088 // Another trigger word
00089 //
00090 // Revision 2.9  2001/08/29 20:45:15  genevb
00091 // Trigger word histos
00092 //
00093 // Revision 2.8  2001/08/07 07:51:28  lansdell
00094 // primvtx check for different multiplicities crashed for MC data, now fixed
00095 //
00096 // Revision 2.7  2001/05/23 00:14:53  lansdell
00097 // more changes for qa_shift histograms
00098 //
00099 // Revision 2.6  2001/05/16 20:57:03  lansdell
00100 // new histograms added for qa_shift printlist; some histogram ranges changed; StMcEvent now used in StEventQA
00101 //
00102 // Revision 2.5  2001/04/28 22:05:13  genevb
00103 // Added EMC histograms
00104 //
00105 // Revision 2.4  2001/04/26 16:34:40  genevb
00106 // Fixed some histogram ranges
00107 //
00108 // Revision 2.3  2000/08/25 22:06:15  genevb
00109 // Raised high mult bin to 2500 global tracks
00110 //
00111 // Revision 2.2  2000/08/25 20:29:34  lansdell
00112 // year1 print list changed slightly; cosmetic improvement of some event summary histograms
00113 //
00114 // Revision 2.1  2000/08/25 16:04:10  genevb
00115 // Introduction of files
00116 //
00117 // 
00119 //                                                                       //
00120 //  StQAMakerBase base class for QA Histogram Makers                     //
00121 //                                                                       //
00123 
00124 #include "StMessMgr.h"
00125 #include "StQAMakerBase.h"
00126 #include "StQABookHist.h"
00127 #include "QAH.h"
00128 #include "TList.h"
00129 
00130 ClassImp(StQAMakerBase)
00131 
00132 //_____________________________________________________________________________
00133 StQAMakerBase::StQAMakerBase(const char *name, const char *title, const char* type) :
00134  StMaker(name,title), QAMakerType(type) {
00135 
00136   hists = 0;
00137   histsList.Expand(32);
00138   histsList.SetOwner();
00139   Int_t i;
00140   for (i=0;i<32;i++) histsList.AddAt(0,i);
00141   histsSet = StQA_Undef;
00142   eventClass = 3;
00143   ITTF = kFALSE;
00144   EST = -1; // -1 = unknown
00145   allTrigs = kFALSE;
00146 
00147 //  - Set all the histogram booking constants
00148 
00149   ntrk   = 50;
00150   nmnpt  = 50;
00151   nmneta = 50;
00152   nxyz   = 50;
00153 
00154 //  - Zero all histogram pointers
00155   mNullPrimVtx = 0; // histogram for number of events without primary vertex
00156   mMultClass = 0;   // histogram for number of events in mult classes
00157   mTrigWord = 0;    // histogram for event trigger words
00158   mTrigBits = 0;    // histogram for event trigger bits
00159   for (i=0; i<24; i++) mTpcSectorPlot[i] = 0;
00160   m_pnt_rpTQW = 0;
00161   m_pnt_rpTQE = 0;
00162 
00163 // FTPC histograms
00164   m_ftpc_chargestepW=0; 
00165   m_ftpc_chargestepE=0; 
00166   m_ftpc_fcl_radialW=0;  
00167   m_ftpc_fcl_radialE=0;  
00168 
00169 }
00170 //_____________________________________________________________________________
00171 StQAMakerBase::~StQAMakerBase() {
00172   if (mNullPrimVtx) delete mNullPrimVtx;
00173   if (mMultClass) delete mMultClass;
00174   if (mTrigWord) delete mTrigWord;
00175   if (mTrigBits) delete mTrigBits;
00176   for (Int_t i=0; i<24; i++) { if (mTpcSectorPlot[i]) delete mTpcSectorPlot[i]; };
00177   if (m_pnt_rpTQW) delete m_pnt_rpTQW;
00178   if (m_pnt_rpTQE) delete m_pnt_rpTQE;
00179 }
00180 //_____________________________________________________________________________
00181 Int_t StQAMakerBase::Init() {
00182 // Histogram booking must wait until first event Make() to determine event type
00183   eventCount = 0;
00184   firstEventClass = kTRUE;
00185   return StMaker::Init();
00186 }
00187 //_____________________________________________________________________________
00188 void StQAMakerBase::Clear(Option_t* opt) {
00189   firstEventClass = kTRUE;
00190   StMaker::Clear(opt);
00191 }
00192 //_____________________________________________________________________________
00193 Int_t StQAMakerBase::Make() {
00194 
00195   if (Debug())
00196     gMessMgr->Info(" In StQAMakerBase::Make()");
00197 
00198   if (!mNullPrimVtx) BookHist();
00199   // See BookHist() for definitions of histsSet values,
00200   // which should be set during Make() of the derived QA Maker class
00201   // event class also decided in derived Make()
00202   switch (histsSet) {
00203     case (StQA_AuAuOld) :
00204       mMultClass->Fill((float) eventClass);
00205     default : {
00206       if (!eventClass) { hists=0; return kStOk; }
00207       hists = (StQABookHist*) histsList.At((--eventClass));
00208     }
00209   }
00210 
00211   if (!hists) NewQABookHist();
00212     
00213 
00214   if (!fillHists) return kStOk;
00215   // Call methods to fill histograms
00216 
00217 
00218   // Those not divided by event class:
00219   if (firstEventClass) {
00220     firstEventClass = kFALSE;
00221   }
00222 
00223   // Those divided by event class:
00224   // histograms from table globtrk
00225   MakeHistGlob();
00226   // histograms from table primtrk - must be done after global tracks
00227   MakeHistPrim();
00228   // histograms from table primtrk & dst_dedx
00229   MakeHistPID();
00230   // histograms from table dst_dedx
00231   MakeHistDE();
00232   // histograms from table point
00233   MakeHistPoint();
00234   // histograms from table dst_vertex,dst_v0_vertex,dst_xi_vertex,dst_kinkVertex
00235   MakeHistVertex();
00236   // histograms from EMC in StEvent
00237   MakeHistEMC();
00238   // histograms from geant and reco tables 
00239   if (histsSet==StQA_MC) MakeHistEval();
00240   // histograms from BBC in StEvent
00241   if (histsSet==StQA_AuAuOld) MakeHistBBC();
00242   // histograms from FPD in StEvent
00243   if (histsSet==StQA_AuAuOld) MakeHistFPD();
00244   // histograms from PMD in StEvent
00245   if (!(histsSet==StQA_run12all ||
00246         histsSet==StQA_run12      )) MakeHistPMD();
00247   // histograms from TOF in StEvent
00248   if (histsSet==StQA_run8     ||
00249       histsSet==StQA_run12all ||
00250       histsSet==StQA_run12      ) MakeHistTOF();
00251 
00252   eventCount++;
00253   return kStOk;
00254 }
00255 //_____________________________________________________________________________
00256 void StQAMakerBase::NewQABookHist() {
00257   const char* pre = prefix[eventClass].Data();
00258   if (Debug())
00259     gMessMgr->Info() <<
00260       "StQAMakerBase: booking histograms with prefix: " <<
00261       pre << endm;
00262   QAH::maker = (StMaker*) (this);
00263   hists = new StQABookHist(pre);
00264   histsList.AddAt(hists,eventClass);
00265   hists->BookHist(histsSet);
00266 }
00267 //_____________________________________________________________________________
00268 TH2F* StQAMakerBase::MH1F(const Text_t* name, const Text_t* title,
00269                           Int_t nbinsx, Axis_t xlow, Axis_t xup) {
00270   TH2F* h = QAH::MH1F(name,title,nbinsx,xlow,xup,eventClass);
00271   if (eventClass>1) {
00272     h->Rebin(0,"low mult");
00273     h->Rebin(1,"mid mult");
00274     h->Rebin(2,"high mult");
00275   }
00276   return h;
00277 }
00278 //_____________________________________________________________________________
00279 void StQAMakerBase::BookHist() {  
00280 // book histograms
00281 
00282   Int_t tempClass = eventClass;
00283 
00284   if (allTrigs) {
00285     prefix[0] = QAMakerType;
00286     eventClass = 1;
00287   } else {
00288   switch (histsSet) {
00289 
00290   // Real data with three multiplicity classes (low, medium, high)
00291 
00292     case (StQA_AuAuOld) : {
00293       (prefix[0] = QAMakerType) += "LM";
00294       (prefix[1] = QAMakerType) += "MM";
00295       (prefix[2] = QAMakerType) += "HM";
00296       eventClass = 3;
00297       break; }
00298 
00299   // Real data with event classes for different triggers
00300 
00301     case (StQA_run12) :
00302     case (StQA_run8) :
00303     case (StQA_AuAu) : {
00304       (prefix[0] = QAMakerType) += "MB";  // Minbias
00305       (prefix[1] = QAMakerType) += "CL";  // Central
00306       (prefix[2] = QAMakerType) += "HT";  // HighTower
00307       (prefix[3] = QAMakerType) += "XX";  // OtherPhysics
00308       (prefix[4] = QAMakerType) += "JP";  // JetPatch
00309       eventClass = 5;
00310       break; }
00311 
00312     case (StQA_dAu) : {
00313       prefix[0] = QAMakerType;  // Minbias
00314       (prefix[1] = QAMakerType) += "HP";
00315       (prefix[2] = QAMakerType) += "XX";
00316       eventClass = 3;
00317       break; }
00318 
00319   // the following data sets use the defaults
00320  
00321    // Generic data (e.g. Monte Carlo) with just one event class
00322     case (StQA_MC) :
00323 
00324    // pp data with just one event class
00325     case (StQA_pp) :
00326     case (StQA_run12all) :
00327 
00328     default  : {
00329       prefix[0] = QAMakerType;
00330       eventClass = 1;
00331     }
00332   }
00333   }
00334   
00335   QAH::maker = (StMaker*) (this);
00336   QAH::preString = QAMakerType;
00337 
00338   BookHistTrigger();
00339   BookHistGeneral();
00340   BookHistFcl();
00341 
00342   Int_t tempClass2 = eventClass;
00343   // Must book the histograms with no special prefix now
00344   for (eventClass=0; eventClass<tempClass2; eventClass++) {
00345     if (!(QAMakerType.CompareTo(prefix[eventClass]))) {
00346       NewQABookHist();
00347       break;
00348     }
00349   }
00350   eventClass = tempClass;
00351   
00352 }
00353 //_____________________________________________________________________________
00354 void StQAMakerBase::BookHistGeneral(){  
00355 
00356   mNullPrimVtx = QAH::H1F("QaNullPrimVtx","event primary vertex check",40,-2,2);
00357   mNullPrimVtx->SetXTitle("has primary vertex? (yes = 1, no = -1)");
00358   mNullPrimVtx->SetYTitle("# of events");
00359 
00360   if (histsSet == StQA_AuAuOld) {
00361     mMultClass = QAH::H1F("QaMultClass","event multiplicity class",5,-0.5,4.5);
00362     mMultClass->SetXTitle("mult class (0=?/MC, 1=LM, 2=MM, 3=HM)");
00363     mMultClass->SetYTitle("# of events");
00364   }
00365 
00366   char namebuf[32];
00367   char titlebuf[64];
00368   for (Int_t i=0; i<24; i++) {
00369     sprintf(namebuf ,"QaTpcSector%d",i+1);
00370     sprintf(titlebuf,"Hits in TPC Sector %d",i+1);
00371     //mTpcSectorPlot[i] = QAH::H2F(namebuf,titlebuf,104,-52.,52.,150,50.,200.);
00372     mTpcSectorPlot[i] = QAH::H2F(namebuf,titlebuf,104,-52.,52.,47,-0.5,46.5);
00373     mTpcSectorPlot[i]->SetXTitle("along padrows [cm]");
00374     //mTpcSectorPlot[i]->SetYTitle("across padrows [cm]");
00375     mTpcSectorPlot[i]->SetYTitle("padrow");
00376   }
00377   m_pnt_rpTQW    = QAH::H2F("QaPointRPTpcQW","point: r-phi distribution of charge, tpcW",20,58.75,196.75,72,0,TMath::TwoPi());
00378   m_pnt_rpTQE    = QAH::H2F("QaPointRPTpcQE","point: r-phi distribution of charge, tpcE",20,58.75,196.75,72,0,TMath::TwoPi());
00379 }
00380 //_____________________________________________________________________________
00381 void StQAMakerBase::BookHistTrigger(){  
00382 
00383   QAH::maker = (StMaker*) (this);
00384   QAH::preString = QAMakerType;
00385   if (mTrigWord) return;
00386   mTrigWord = QAH::H1F("QaTrigWord","trigger word",8,0.5,8.5);
00387   mTrigWord->SetXTitle("1:MinBias 2:Central 3:HiPt 4:Jet 5:HiTower 6:OtherPhys");
00388   mTrigBits = QAH::H1F("QaTrigBits","trigger bits",32,-0.5,31.5);
00389 }
00390 //_____________________________________________________________________________
00391 void StQAMakerBase::BookHistFcl(){
00392 
00393   // Get fcl histograms from FTPC makers
00394   if (!(m_ftpc_chargestepW)) {
00395     // First try to get histograms from StFtpcClusterMaker named "ftpc_hits"
00396     StMaker* fhMaker = GetMaker("ftpc_hits");
00397     if (fhMaker) {
00398       m_ftpc_chargestepW = (TH1F*) (fhMaker->GetHist("fcl_chargestepW"));
00399       AddHist(m_ftpc_chargestepW);
00400       m_ftpc_chargestepE = (TH1F*) (fhMaker->GetHist("fcl_chargestepE"));
00401       AddHist(m_ftpc_chargestepE);
00402       m_ftpc_fcl_radialW = (TH1F*) (fhMaker->GetHist("fcl_radialW"));
00403       AddHist(m_ftpc_fcl_radialW);
00404       m_ftpc_fcl_radialE = (TH1F*) (fhMaker->GetHist("fcl_radialE"));
00405       AddHist(m_ftpc_fcl_radialE);
00406     } else {
00407       // "ftpc_hits" maker doesn't exist, so look in hist branch
00408       // *** Currently isn't working for bfcread_event_QAhist.C ***
00409       St_DataSet* hDS = GetDataSet("histBranch");
00410       if (hDS) {
00411         // hDS->ls(9);
00412         St_DataSet* fhDS = hDS->Find("ftpc_hitsHist");
00413         if (fhDS) {
00414           m_ftpc_chargestepW =
00415             (TH1F*) (fhDS->FindObject("fcl_chargestepW"));
00416           AddHist(m_ftpc_chargestepW);
00417           m_ftpc_chargestepE =
00418             (TH1F*) (fhDS->FindObject("fcl_chargestepE"));
00419           AddHist(m_ftpc_chargestepE);
00420           m_ftpc_fcl_radialW =
00421             (TH1F*) (fhDS->FindObject("fcl_radialW"));
00422           AddHist(m_ftpc_fcl_radialW);
00423           m_ftpc_fcl_radialE =
00424             (TH1F*) (fhDS->FindObject("fcl_radialE"));
00425           AddHist(m_ftpc_fcl_radialE);
00426         }
00427       }
00428     }
00429   }
00430 }
00431 //_____________________________________________________________________________