StRoot: St_ctf_Maker/St_ctf_Maker.cxx Source File

00001 // $Id: St_ctf_Maker.cxx,v 1.17 2009/11/10 21:09:59 fisyak Exp $
00002 // $Log: St_ctf_Maker.cxx,v $
00003 // Revision 1.17  2009/11/10 21:09:59  fisyak
00004 // pams clean up
00005 //
00006 // Revision 1.16  2007/04/28 17:55:47  perev
00007 // Redundant StChain.h removed
00008 //
00009 // Revision 1.15  2003/04/30 20:39:14  perev
00010 // Warnings cleanup. Modified lines marked VP
00011 //
00012 // Revision 1.14  2000/06/26 22:13:06  fisyak
00013 // remove params
00014 //
00015 // Revision 1.13  1999/07/15 13:57:48  perev
00016 // cleanup
00017 //
00018 // Revision 1.12  1999/03/11 03:55:06  perev
00019 // new schema
00020 //
00021 // Revision 1.11  1999/03/02 00:11:49  llope
00022 // reduced table maxlen's for mdc2
00023 //
00024 // Revision 1.10  1999/02/26 14:49:54  kathy
00025 // fixed histogram limits as per Bill's request
00026 //
00027 // Revision 1.9  1999/02/25 21:11:46  kathy
00028 // fix histograms
00029 //
00030 // Revision 1.8  1999/02/25 19:38:50  kathy
00031 // fix up histograms
00032 //
00033 // Revision 1.7  1999/02/23 21:25:42  llope
00034 // fixed histograms, added 1/beta vs p
00035 //
00036 // Revision 1.6  1999/02/23 01:09:47  fisyak
00037 // Add Bill Llope dst tables
00038 //
00039 // Revision 1.5  1999/02/06 00:15:46  fisyak
00040 // Add adc/tdc histograms
00041 //
00042 // Revision 1.4  1999/01/25 23:39:12  fisyak
00043 // Add tof
00044 //
00045 // Revision 1.3  1999/01/21 00:52:31  fisyak
00046 // Cleanup
00047 //
00048 // Revision 1.2  1999/01/02 19:08:14  fisyak
00049 // Add ctf
00050 //
00051 // Revision 1.1  1999/01/01 02:39:38  fisyak
00052 // Add ctf Maker
00053 //
00054 // Revision 1.7  1998/10/31 00:25:45  fisyak
00055 // Makers take care about branches
00056 //
00057 // Revision 1.6  1998/10/06 18:00:29  perev
00058 // cleanup
00059 //
00060 // Revision 1.5  1998/10/02 13:46:08  fine
00061 // DataSet->DataSetIter
00062 //
00063 // Revision 1.4  1998/08/14 15:25:58  fisyak
00064 // add options
00065 //
00066 // Revision 1.3  1998/08/10 02:32:07  fisyak
00067 // Clean up
00068 //
00069 // Revision 1.2  1998/07/20 15:08:15  fisyak
00070 // Add tcl and tpt
00071 //
00073 //                                                                      //
00074 // St_ctf_Maker class for Makers                                        //
00075 //                                                                      //
00077 
00078 #include "St_ctf_Maker.h"
00079 #include "St_DataSetIter.h"
00080 #include "ctf/St_ctg_Module.h"
00081 #include "ctf/St_ctg_Module.h"
00082 #include "ctf/St_cts_Module.h"
00083 #include "ctf/St_ctu_Module.h"
00084 #ifdef __TPT__
00085 #include "ctf/St_fill_dst_tof_Module.h"
00086 #endif
00087 #include "TH1.h"
00088 #include "TH2.h"
00089 
00090 ClassImp(St_ctf_Maker)
00091 
00092 //_____________________________________________________________________________
00093 St_ctf_Maker::St_ctf_Maker(const char *name):StMaker(name),
00094 m_ctb(0),
00095 m_ctb_slat_phi(0),
00096 m_ctb_slat_eta(0),
00097 m_ctb_slat(0),
00098 m_cts_ctb(0),
00099 m_tof(0),
00100 m_tof_slat_phi(0),
00101 m_tof_slat_eta(0),
00102 m_tof_slat(0),
00103 m_cts_tof(0)
00104 {
00105    drawinit=kFALSE;
00106 }
00107 //_____________________________________________________________________________
00108 St_ctf_Maker::~St_ctf_Maker(){
00109 }
00110 //_____________________________________________________________________________
00111 Int_t St_ctf_Maker::Init(){
00112   int iInit=0;
00113 
00114 // Create tables
00115   St_DataSet *ctfg = GetInputDB("ctf/ctg");
00116   assert (ctfg);
00117   St_DataSetIter       gime(ctfg);
00118   m_ctb          = (St_ctg_geo      *) gime("ctb");
00119   m_ctb_slat_phi = (St_ctg_slat_phi *) gime("ctb_slat_phi");
00120   m_ctb_slat_eta = (St_ctg_slat_eta *) gime("ctb_slat_eta");
00121   m_ctb_slat     = (St_ctg_slat     *) gime("ctb_slat");
00122 
00123   Int_t Res_ctg_ctb  =  ctg (m_ctb,m_ctb_slat_phi,m_ctb_slat_eta,m_ctb_slat);
00124   if (Res_ctg_ctb!=kSTAFCV_OK) iInit = kStWarn;
00125   m_tof          = (St_ctg_geo      *) gime("tof");
00126   m_tof_slat_phi = (St_ctg_slat_phi *) gime("tof_slat_phi");
00127   m_tof_slat_eta = (St_ctg_slat_eta *) gime("tof_slat_eta");
00128   m_tof_slat     = (St_ctg_slat     *) gime("tof_slat");
00129 
00130   Int_t Res_ctg_tof  =  ctg (m_tof,m_tof_slat_phi,m_tof_slat_eta,m_tof_slat); if(Res_ctg_tof){};
00131   // Special treatment for double names
00132   //  m_cts          = (St_cts_mpara    *) params("ctf/cts")->GetList()->FindObject("cts");
00133   St_DataSet *ctfs = GetInputDB("ctf/cts");
00134   assert (ctfs);
00135   gime.Reset(ctfs);
00136   m_cts_ctb          = (St_cts_mpara    *) gime("cts_ctb");
00137   m_cts_tof          = (St_cts_mpara    *) gime("cts_tof");
00138   // Create Histograms  
00139   m_adcc  = new TH1F("CtfCtbrawAdc","CTB ADCs",128,0.,1024.);
00140   m_adct  = new TH1F("CtfTofrawAdc","TOF ADCs",128,0.,2048.);
00141   m_tsvsp = new TH2F("CtfDsttoftrkTsvsp","TOF 1/beta vs. Ptrk",80,0.1,4.1,100,0.,10.);
00142   m_tsvsp->SetXTitle("TPC track momentum (GeV/c)");
00143   m_tsvsp->SetYTitle("1/beta from TOFp");
00144   m_tsvsp1= new TH2F("CtfDsttoftrkTsvsp1","TOF 1/beta vs. Ptrk",80,0.1,4.1,100,0.,10.);
00145   m_tsvsp1->SetXTitle("TPC track momentum (GeV/c)");
00146   m_tsvsp1->SetYTitle("1/beta from TOFp, slat Nhits=1");
00147   return StMaker::Init();
00148 }
00149 //_____________________________________________________________________________
00150 Int_t St_ctf_Maker::Make(){
00151 //  PrintInfo();
00152 
00153   int iMake=kStOK;
00154   
00155   St_DataSet *gea = GetInputDS("geant");
00156   if (gea) {//          Geant input exists
00157     St_DataSetIter geant(gea);
00158     St_g2t_track   *g2t_track   = (St_g2t_track *)   geant("g2t_track");
00159     St_g2t_ctf_hit *g2t_ctb_hit = (St_g2t_ctf_hit *) geant("g2t_ctb_hit");
00160     if (g2t_ctb_hit) {
00161       St_cts_mslat *ctb_mslat = new  St_cts_mslat("ctb_mslat", 240); m_DataSet->Add(ctb_mslat);
00162       St_cts_event *ctb_event = new  St_cts_event("ctb_event",   1); m_DataSet->Add(ctb_event);
00163       St_ctu_raw   *ctb_raw   = new  St_ctu_raw("ctb_raw",     240); m_DataSet->Add(ctb_raw);
00164       St_ctu_cor   *ctb_cor   = new  St_ctu_cor("ctb_cor",     240); m_DataSet->Add(ctb_cor);
00165       Int_t Res_cts_ctb = cts(g2t_ctb_hit, g2t_track,
00166                               m_ctb,  m_ctb_slat, m_ctb_slat_phi, m_ctb_slat_eta, m_cts_ctb,
00167                               ctb_event, ctb_mslat, ctb_raw); if(Res_cts_ctb){};
00168       Int_t Res_ctu_ctb = ctu(m_ctb,  m_ctb_slat, ctb_raw, ctb_cor);if(Res_ctu_ctb){};
00169       ctu_raw_st *raw   = ctb_raw->GetTable();
00170       for (Int_t i=0; i<ctb_raw->GetNRows();i++,raw++){
00171         m_adcc->Fill((Float_t) raw->adc);
00172       }
00173     }
00174 
00175     St_g2t_ctf_hit *g2t_tof_hit = (St_g2t_ctf_hit *) geant("g2t_tof_hit");
00176     if (g2t_tof_hit) {
00177       St_cts_mslat *tof_mslat = new  St_cts_mslat("tof_mslat",500); m_DataSet->Add(tof_mslat); //safe factor of 5 for mdc2
00178       St_cts_event *tof_event = new  St_cts_event("tof_event",  1); m_DataSet->Add(tof_event);
00179       St_ctu_raw   *tof_raw   = new  St_ctu_raw("tof_raw",     50); m_DataSet->Add(tof_raw);
00180       St_ctu_cor   *tof_cor   = new  St_ctu_cor("tof_cor",     50); m_DataSet->Add(tof_cor); 
00181       Int_t Res_cts_tof = cts(g2t_tof_hit, g2t_track,
00182                               m_tof,  m_tof_slat, m_tof_slat_phi, m_tof_slat_eta, m_cts_tof,
00183                               tof_event, tof_mslat, tof_raw);if(Res_cts_tof){}
00184 #ifdef __TPT__
00185       St_DataSet *tpc_tracks = GetDataSet("tpc_tracks");
00186       St_tpt_track  *tptrack = 0;
00187       St_tte_mctrk  *mctrk   = 0;
00188       if (tpc_tracks) {
00189          St_DataSetIter tpcI(tpc_tracks);
00190          tptrack = (St_tpt_track *) tpcI["tptrack"];
00191          mctrk   = (St_tte_mctrk *) tpcI["mctrk"];
00192       }
00193       St_DataSet *global = GetDataSet("global");
00194       St_dst_vertex     *vertex      = 0;
00195       if (global) {
00196          St_DataSetIter globalI(global);
00197          vertex  = (St_dst_vertex *) globalI["dst/vertex"];
00198       }
00199       if (g2t_track && tptrack && mctrk && vertex ) {
00200          St_dst_tof_trk *dst_tof_trk = new St_dst_tof_trk("dst_tof_trk",250); //safe factor of 5 for mdc2
00201          m_DataSet->Add(dst_tof_trk);
00202          St_dst_tof_evt *dst_tof_evt = new St_dst_tof_evt("dst_tof_evt",1);
00203          m_DataSet->Add(dst_tof_evt);
00204          Int_t Res_fill_dst_tof = fill_dst_tof(g2t_tof_hit,g2t_track,
00205                                               tptrack,mctrk,vertex,
00206                                               m_tof,m_tof_slat,
00207                                               m_tof_slat_phi,m_tof_slat_eta,
00208                                               m_cts_tof,tof_mslat,
00209                                               dst_tof_trk,dst_tof_evt); if(Res_fill_dst_tof){};
00210          dst_tof_trk_st *dst = dst_tof_trk->GetTable();
00211          for (Int_t i=0; i<dst_tof_trk->GetNRows();i++,dst++){
00212            m_tsvsp->Fill(dst->ptot,dst->ts_mtime);
00213            if (dst->n_hits == 1) {
00214             m_tsvsp1->Fill(dst->ptot,dst->ts_mtime);
00215            }
00216          }
00217       }  
00218       Int_t Res_ctu_tof = ctu(m_tof,  m_tof_slat, tof_raw, tof_cor);if(Res_ctu_tof){}
00219       ctu_raw_st *raw   = tof_raw->GetTable();
00220       for (Int_t i=0; i<tof_raw->GetNRows();i++,raw++){
00221         m_adct->Fill((Float_t) raw->adc);
00222       }
00223 #endif
00224     }
00225   }
00226   return iMake;
00227   
00228 }
00229 //_____________________________________________________________________________
00230