StRoot: StEmcPool/StEmcOfflineCalibrationMaker/macros/bemcCalibMacro.C Source File

00001 void bemcCalibMacro(const char* dir="./",
00002                            const char* name  = "test.root",
00003                            const char* filelist = "test.list",
00004                            int nFiles = 1,
00005                            int nEvents = 1000,
00006                            const char* outPath = "./")
00007 {
00008         gROOT->Macro("LoadLogger.C");
00009         gROOT->Macro("loadMuDst.C");
00010     gSystem->Load("StarMagField.so");
00011     gSystem->Load("StMagF");
00012     gSystem->Load("StDetectorDbMaker");
00013     gSystem->Load("StTpcDb");
00014     gSystem->Load("St_db_Maker");
00015     gSystem->Load("StDbUtilities");
00016     gSystem->Load("StMcEvent");
00017     gSystem->Load("StMcEventMaker");
00018     gSystem->Load("StDaqLib");
00019     gSystem->Load("StEmcRawMaker");
00020     gSystem->Load("StEmcADCtoEMaker");
00021     gSystem->Load("StEpcMaker");
00022     gSystem->Load("StTriggerUtilities");
00023     //gSystem->Load("StEmcTriggerMaker");
00024     gSystem->Load("StDbBroker");
00025     gSystem->Load("libgeometry_Tables");
00026     gSystem->Load("StEEmcUtil");
00027     gSystem->Load("StEEmcDbMaker");
00028     gSystem->Load("StEmcOfflineCalibrationMaker");
00029         
00030     StChain* chain = new StChain("StChain");
00031         //chain->SetDebug();
00032         
00033     StMuDstMaker* muDstMaker = new StMuDstMaker(0,0,"",filelist,"",nFiles);
00034     St_db_Maker *dbMaker = new St_db_Maker("StarDb","MySQL:StarDb");
00035     StEEmcDbMaker* eemcb = new StEEmcDbMaker("eemcDb");
00036         StEmcADCtoEMaker *adc = new StEmcADCtoEMaker();
00037         /*
00038         StDetectorId towerid = static_cast<StDetectorId>(kBarrelEmcTowerId);
00039         StDetectorId prsid = static_cast<StDetectorId>(kBarrelEmcTowerId+1);
00040         StDetectorId smdeid = static_cast<StDetectorId>(kBarrelEmcTowerId+2);
00041         StDetectorId smdpid = static_cast<StDetectorId>(kBarrelEmcTowerId+3);
00042         adc->setDoZeroSuppression(towerid,0);
00043         adc->setDoZeroSuppression(prsid,0);
00044         adc->setDoZeroSuppression(smdeid,0);
00045         adc->setDoZeroSuppression(smdpid,0);
00046         */      
00047         //get control table so we can turn off BPRS zero-suppression and save hits from "bad" caps
00048 
00049         controlADCtoE_st* control_table = adc->getControlTable();
00050         control_table->CutOff[1] = -1;
00051         control_table->CutOffType[1] = 0;
00052         control_table->DeductPedestal[1] = 2;
00053         adc->saveAllStEvent(kTRUE);
00054         //StEmcTriggerMaker *emcTrig = new StEmcTriggerMaker("bemctrigger");
00055         //emcTrig->setDbMaker(dbMaker);
00056 
00057     StTriggerSimuMaker* trigsim = new StTriggerSimuMaker();
00058     //trigsim->useBbc();
00059     trigsim->useBemc();
00060     trigsim->useEemc();
00061     trigsim->bemc->setConfig(StBemcTriggerSimu::kOffline);
00062     //StGenericL2Emulator* simL2Mk = new StL2_2008EmulatorMaker;
00063     //assert(simL2Mk);
00064     //simL2Mk->setSetupPath("/afs/rhic.bnl.gov/star/users/kocolosk/public/StarTrigSimuSetup/");
00065     //simL2Mk->setOutPath(outPath);
00066     //trigsim->useL2(simL2Mk);
00067 
00068         
00069         TString outfile(dir);
00070         outfile += "/";
00071         outfile += name;
00072         StEmcOfflineCalibrationMaker* bemcCalibMaker = new StEmcOfflineCalibrationMaker("bemcCalibMaker",outfile.Data());
00073         //bemcCalibMaker->subtractPedestals = true;
00074 
00075         /*
00076         //2006 Triggers 
00077         bemcCalibMaker->addMinBiasTrigger(200000);
00078         bemcCalibMaker->addMinBiasTrigger(200002);
00079         bemcCalibMaker->addMinBiasTrigger(200012);
00080         
00081         bemcCalibMaker->addHighTowerTrigger(117211);
00082         bemcCalibMaker->addHighTowerTrigger(117212);
00083         bemcCalibMaker->addHighTowerTrigger(127212);
00084         bemcCalibMaker->addHighTowerTrigger(127213);
00085         bemcCalibMaker->addHighTowerTrigger(137213);
00086         */
00087 
00088         //2008 Triggers
00089         /*
00090         bemcCalibMaker->addMinBiasTrigger(220000);
00091         bemcCalibMaker->addHighTowerTrigger(220500);//bht0
00092         bemcCalibMaker->addHighTowerTrigger(220510);//bht1
00093         bemcCalibMaker->addHighTowerTrigger(220520);//bht2
00094         bemcCalibMaker->addHighTowerTrigger(1);//bht2-mon
00095 
00096         bemcCalibMaker->addFastTrigger(220900);//fms-fast
00097         bemcCalibMaker->addFastTrigger(220901);//fms-fast
00098         bemcCalibMaker->addFastTrigger(220910);//fms-fast tpc line
00099         bemcCalibMaker->addFastTrigger(220920);//fpd e fast
00100         bemcCalibMaker->addFastTrigger(220710);//tof reading only tpx and tof
00101         bemcCalibMaker->addFastTrigger(19);//bbc-fast
00102         */
00103 
00104         //2009 500 GeV triggers
00105         /*
00106         bemcCalibMaker->addMinBiasTrigger(230010);//BBCMB
00107         bemcCalibMaker->addMinBiasTrigger(230011);//BBCMB-Cat0
00108         bemcCalibMaker->addMinBiasTrigger(230012);//Cat1
00109         bemcCalibMaker->addMinBiasTrigger(230013);//Cat2
00110         bemcCalibMaker->addMinBiasTrigger(230014);//Cat3
00111         bemcCalibMaker->addMinBiasTrigger(230015);//Cat4
00112 
00113         bemcCalibMaker->addMinBiasTrigger(230020);//VPBMB didn't work
00114 
00115         bemcCalibMaker->addHighTowerTrigger(230531);//bht3
00116         bemcCalibMaker->addHighTowerTrigger(230601);//l2w
00117         */
00118 
00119         //2009 200 GeV triggers
00120         bemcCalibMaker->addMinBiasTrigger(240010);//BBCMB
00121         bemcCalibMaker->addMinBiasTrigger(240110);//BBCMB
00122         bemcCalibMaker->addMinBiasTrigger(240120);//BBCMB
00123         bemcCalibMaker->addMinBiasTrigger(240220);//BBCMB
00124 
00125         bemcCalibMaker->addMinBiasTrigger(240011);//BBCMB-Cat0
00126         bemcCalibMaker->addMinBiasTrigger(240111);//BBCMB-Cat0
00127         bemcCalibMaker->addMinBiasTrigger(240121);//BBCMB-Cat0
00128         bemcCalibMaker->addMinBiasTrigger(240221);//BBCMB-Cat0
00129 
00130         bemcCalibMaker->addMinBiasTrigger(240012);//BBCMB-Cat1
00131         bemcCalibMaker->addMinBiasTrigger(240112);//BBCMB-Cat1
00132         bemcCalibMaker->addMinBiasTrigger(240122);//BBCMB-Cat1
00133         bemcCalibMaker->addMinBiasTrigger(240222);//BBCMB-Cat1
00134 
00135         bemcCalibMaker->addMinBiasTrigger(240013);//BBCMB-Cat2
00136         bemcCalibMaker->addMinBiasTrigger(240113);//BBCMB-Cat2
00137         bemcCalibMaker->addMinBiasTrigger(240123);//BBCMB-Cat2
00138         bemcCalibMaker->addMinBiasTrigger(240223);//BBCMB-Cat2
00139 
00140         bemcCalibMaker->addMinBiasTrigger(240014);//BBCMB-Cat3
00141         bemcCalibMaker->addMinBiasTrigger(240114);//BBCMB-Cat3
00142         bemcCalibMaker->addMinBiasTrigger(240124);//BBCMB-Cat3
00143         bemcCalibMaker->addMinBiasTrigger(240224);//BBCMB-Cat3
00144 
00145         bemcCalibMaker->addMinBiasTrigger(240015);//BBCMB-Cat4
00146         bemcCalibMaker->addMinBiasTrigger(240115);//BBCMB-Cat4
00147         bemcCalibMaker->addMinBiasTrigger(240125);//BBCMB-Cat4
00148         bemcCalibMaker->addMinBiasTrigger(240225);//BBCMB-Cat4
00149 
00150         bemcCalibMaker->addMinBiasTrigger(240020);//VPDMB no vertex
00151         bemcCalibMaker->addMinBiasTrigger(240025);//VPDMB vertex < 30
00152 
00153         bemcCalibMaker->addMinBiasTrigger(240130);//ZDCMB
00154 
00155         bemcCalibMaker->addHighTowerTrigger(240550);//bht0*VPDMB*!BHT1
00156         bemcCalibMaker->addHighTowerTrigger(240570);//bht0*VPDMB*!BHT2
00157         bemcCalibMaker->addHighTowerTrigger(240560);//bht1*VPDMB
00158         bemcCalibMaker->addHighTowerTrigger(240540);//bht2*VPDMB
00159         bemcCalibMaker->addHighTowerTrigger(240530);//bht3
00160 
00161         bemcCalibMaker->addFastTrigger(240801);//fms-fast
00162         bemcCalibMaker->addFastTrigger(240802);//fms-fast
00163 
00164         bemcCalibMaker->addFastTrigger(240640);//upsilon broken so skip
00165 
00166 
00167 
00168         StMemStat memory;
00169         memory.PrintMem(NULL);
00170         
00171         chain->Init();
00172         cout<<"chain initialized"<<endl;
00173         
00174         TStopwatch total;
00175         TStopwatch timer;
00176         
00177         int i=0;
00178         while(i<nEvents && chain->Make()==kStOk)
00179         {
00180                 if(i % 500 == 0){
00181                         cout<<"done with event "<<i;
00182                         cout<<"\tcpu: "<<timer.CpuTime()<<"\treal: "<<timer.RealTime()<<"\tratio: "<<timer.CpuTime()/timer.RealTime();//<<endl;
00183                                 timer.Start();
00184                                 memory.PrintMem(NULL);
00185                 }
00186                 i++;
00187                 chain->Clear();
00188         }
00189         
00190         chain->ls(3);
00191         chain->Finish();
00192         printf("my macro processed %i events in %s",i,name);
00193         cout<<"\tcpu: "<<total.CpuTime()<<"\treal: "<<total.RealTime()<<"\tratio: "<<total.CpuTime()/total.RealTime()<<endl;
00194 }