StRoot: StFtpcDriftMapMaker/macros/FtpcDriftMapMaker.C Source File

00001 // $Id: FtpcDriftMapMaker.C,v 1.5 2009/11/10 11:07:09 jcs Exp $
00002 // $Log: FtpcDriftMapMaker.C,v $
00003 // Revision 1.5  2009/11/10 11:07:09  jcs
00004 // Now necessary to also load StDetectorDbMaker library
00005 //
00006 // Revision 1.4  2009/11/10 10:59:26  jcs
00007 // change map to Map to avoid conflict with cint
00008 //
00009 // Revision 1.3  2006/08/02 13:59:16  jcs
00010 // add deltaAr argument to allow user to change gas compostion (default: deltaAr=0)
00011 //
00012 // Revision 1.2  2005/12/12 09:30:52  jcs
00013 // load StarMagField library
00014 //
00015 // Revision 1.1  2005/12/12 09:09:11  jcs
00016 // Macro used to run StFtpcDriftMapMakerZ
00017 //
00018 //
00019 //======================================================================
00020 // owner:  Janet Seyboth  (jcs@mppmu.mpg.de)
00021 // what it does: compute drift map for FTPCs
00022 //
00023 //    const Int_t      Map     = 2           use mapped field values
00024 //                             = 1           use constant field values
00025 //    const Float_t   |factor| > 0.8         scale from full field
00026 //                    factor   > 0           normal field
00027 //                    factor   < 0           reversed field
00028 //                                           factor is the field scaling factor
00029 //    const Float_t   deltaAr                change gas composition: 
00030 //                                                   % Ar + deltaAr
00031 //                                                   % CO2 - deltaAr
00032 //
00033 //    Default:  const Int_t   Map    = 2
00034 //              const Float_t factor = 1.0
00035 //              const Float_t deltaAr = 0.0
00036 //   
00037 //======================================================================
00038 
00039 void FtpcDriftMapMaker(const Int_t Map=2, const Float_t factor=1.0, const Float_t deltaAr = 0.0)
00040 {
00041     cout<<"FtpcDriftMapMaker.C called with:"<<endl;
00042     cout<<"                                  Map     = "<<Map<<endl;
00043     cout<<"                                  factor  = "<<factor<<endl;
00044     cout<<"                                  deltaAr = "<<deltaAr<<endl;
00045     if (gClassTable->GetID("TTable") < 0) gSystem->Load("libTable");
00046     gSystem->Load("St_base");
00047     gSystem->Load("StChain");
00048     gSystem->Load("St_Tables");
00049     gSystem->Load("StUtilities");
00050     gSystem->Load("StarClassLibrary");
00051     gSystem->Load("StDetectorDbMaker");
00052     gSystem->Load("StDbUtilities");
00053     gSystem->Load("StDbLib");
00054     gSystem->Load("StDbBroker");
00055     gSystem->Load("St_db_Maker");
00056     gSystem->Load("StarMagField");
00057     gSystem->Load("StMagF");
00058     gSystem->Load("libftpc_Tables");
00059     gSystem->Load("StFtpcDriftMapMaker");
00060     gSystem->Load("StFtpcClusterMaker");
00061 
00062     //  Create the makers to be called by the current chain
00063     const char *mysqlDB =  "MySQL:StarDb";
00064     const char *paramsDB = "$STAR/StarDb";
00065     //const char *paramsDB = "$PWD/StarDb";
00066     StChain *chain =  new StChain();
00067     
00068     St_db_Maker *dbMk = new St_db_Maker("db",mysqlDB,paramsDB);
00069     dbMk->SetDateTime(20010501,00000);
00070 
00071   // Full Field Positive ?
00072   if ( factor > 0.8 ) {
00073      dbMk->SetFlavor("ffp10kv","ftpcVDrift");
00074      dbMk->SetFlavor("ffp10kv","ftpcdVDriftdP");
00075      dbMk->SetFlavor("ffp10kv","ftpcDeflection");
00076      dbMk->SetFlavor("ffp10kv","ftpcdDeflectiondP");
00077      gMessMgr->Info() << "StFtpcDriftMapMaker::InitRun: flavor set to ffp10kv"<<endm;
00078   }
00079   else if ( factor > 0.2 ) {
00080      dbMk->SetFlavor("hfp10kv","ftpcVDrift");
00081      dbMk->SetFlavor("hfp10kv","ftpcdVDriftdP");
00082      dbMk->SetFlavor("hfp10kv","ftpcDeflection");
00083      dbMk->SetFlavor("hfp10kv","ftpcdDeflectiondP");
00084      gMessMgr->Info() << "StFtpcDriftMapMaker::InitRun: flavor set to hfp10kv"<<endm;
00085   }
00086   else if ( factor > -0.2 ) {
00087      dbMk->SetFlavor("zf10kv","ftpcVDrift");
00088      dbMk->SetFlavor("zf10kv","ftpcdVDriftdP");
00089      dbMk->SetFlavor("zf10kv","ftpcDeflection");
00090      dbMk->SetFlavor("zf10kv","ftpcdDeflectiondP");
00091      gMessMgr->Info() << "StFtpcDriftMapMaker::InitRun: flavor set to zf10kv"<<endm;
00092   }
00093   else if ( factor > -0.8 ) {
00094      dbMk->SetFlavor("hfn10kv","ftpcVDrift");
00095      dbMk->SetFlavor("hfn10kv","ftpcdVDriftdP");
00096      dbMk->SetFlavor("hfn10kv","ftpcDeflection");
00097      dbMk->SetFlavor("hfn10kv","ftpcdDeflectiondP");
00098      gMessMgr->Info() << "StFtpcDriftMapMaker::InitRun: flavor set to hfn10kv"<<endm;
00099   }
00100   else {
00101      dbMk->SetFlavor("ffn10kv","ftpcVDrift");
00102      dbMk->SetFlavor("ffn10kv","ftpcdVDriftdP");
00103      dbMk->SetFlavor("ffn10kv","ftpcDeflection");
00104      dbMk->SetFlavor("ffn10kv","ftpcdDeflectiondP");
00105      gMessMgr->Info() << "StFtpcDriftMapMaker::InitRun: flavor set to ffn10kv"<<endm;
00106   }
00107     dbMk->Init();
00108     dbMk->Make();
00109     
00110     //b=new TBrowser();
00111     
00112     StFtpcDriftMapMaker *ftpcDriftMapMk = new StFtpcDriftMapMaker(Map,factor,deltaAr);
00113 }