StRoot: macros/examples/bfcMixer_v4

00001 
00002 //
00003 // Macro for running chain with different inputs
00004 //
00005 // Owner:  Yuri Fisyak
00006 //
00007 // $Id: bfcMixer_v4_noFTPC.C,v 1.2 2006/08/23 22:30:06 fisyak Exp $
00008 //
00010 
00011 //TBrowser *b = 0;
00012 class StBFChain;
00013 StBFChain  *chain=0;
00014 class StMaker;
00015 StMaker    *treeMk=0;
00016 StBFChain *chain1, *chain2, *chain3;
00017 class StEvent;
00018 StEvent *Event;
00019 class St_geant_Maker;
00020 class StIOMaker;
00021 class StEventDisplayMaker; StEventDisplayMaker *dsMk = 0;
00022 class StEventMaker; StEventMaker *evMk = 0;
00023 class StMixerMaker;
00024 //class StFtpcMixerMaker;
00025 class StEvtHddr;
00026 //_____________________________________________________________________
00027 void Load(){
00028   gSystem->Load("St_base");
00029   gSystem->Load("StChain");
00030   gSystem->Load("StUtilities");
00031   gSystem->Load("StBFChain");
00032   //Extra things to load for the acceptance filter
00033   //  gSystem->Load("StarClassLibrary");
00034   gSystem->Load("StAnalysisUtilities");
00035   //  gSystem->Load("StV0AccMaker.so");
00036 
00037   if (chain) delete chain;
00038 }
00039 //_____________________________________________________________________
00040 void bfcMixer_v4_noFTPC(const Int_t Nevents=10,
00041                         //             const Char_t *file1="/auto/pdsfdv08/starprod/daq/2004/production62GeV/ReversedFullField/st_physics_adc_5090009_raw_2060002.daq",
00042                         const Char_t *file1="/star/data03/daq/2004/090/st_physics_adc_5090009_raw_2060002.daq",
00043                         //           const Char_t *file2="/home/starofl/embedding/GSTAR/gtest.fz",
00044                         const Char_t *file2="./gtest.fz",
00045                         //             const Char_t *file3="/home/starofl/embedding/GSTAR/st_physics_2270008_raw_0030.vertices.dat",
00046                         //             const Char_t *file3="/home/starofl/embedding/GSTAR/st_physics_2270008_raw_0030.vertices.dat",
00047                         const Char_t *file3="./st_physics_adc_5090009_raw_2060002.vertices.dat",
00048              const Float_t zvertex_low=-175.0,
00049              const Float_t zvertex_high=175.0,
00050              const Char_t *mode="strange",
00051              const Char_t *acc_mode="off" )
00052 {
00053   // Dynamically link some shared libs
00054   if (gClassTable->GetID("StBFChain") < 0) Load();
00055 
00056   // Create the main chain object
00057   chain = new StBFChain("Embedding");
00058   
00059   StMaker *saveMk = 0;
00060   chain->SetFlags("-ittf,NoDefault");
00061   // Create chain1 object
00062   chain1 = new StBFChain("One");
00063   
00064   saveMk = chain1->cd();
00065   //  chain1->SetFlags("in NoDefault");
00066   //  chain1->SetFlags("in alltrigger NoDefault");
00067   //  chain1->SetFlags("in Physics DbV20020226 NoDefault");
00068   chain1->SetFlags("in Physics DbV20060421 NoDefault -ittf");
00069   chain1->Set_IO_Files(file1);
00070   chain1->Load();
00071   chain1->Instantiate();
00072 
00073   saveMk->cd();
00074   
00075   // Create chain2 object
00076   chain2 = new StBFChain("Two");
00077   saveMk = chain2->cd();
00078   //  chain2->SetFlags("fzin DbV20020226 gen_T geomT sim_T tpc trs -tcl -tpt -PreVtx -tpc_daq");   // 
00079   chain2->SetFlags("fzin DbV20060421 gen_T geomT sim_T tpc trs -tcl -tpt -PreVtx -tpc_daq -ittf");   // 
00080   chain2->Set_IO_Files(file2);
00081   chain2->Load();
00082   chain2->Instantiate();
00083   St_geant_Maker *geantMk = chain2->GetMaker("geant");
00084   if (geantMk) geantMk->SetMode(1);   // Mixer mode - do not modify EvtHddr
00085   
00086   if (chain2->GetOption("TRS")){
00087     StTrsMaker *trsMk = (StTrsMaker *) chain2->GetMaker("Trs");
00088     trsMk->setNormalFactor(1.15);
00089   }
00090 
00091   // Add the acceptance filter maker before TRS  
00092   if (!strcmp(mode,"strange")){
00093     if (!strcmp(acc_mode,"on")){
00094       
00095       Char_t *extraMaker = "StV0AccMaker";
00096       if (gClassTable->GetID(extraMaker) < 0) gSystem->Load(extraMaker);
00097       StMaker *extraMk = (StMaker *)chain1->GetMaker(extraMaker);
00098       if(extraMk) delete extraMk;
00099       extraMk = chain->New(extraMaker,"before");
00100       if (extraMk) {
00101         Char_t *before = "Trs";
00102         StMaker *trsmk = chain1->GetMaker(before);
00103         if (trsmk) chain1->AddBefore(before,extraMk);
00104         StV0AccCuts *cuts = ((StV0AccMaker *)extraMk)->GetCutsPtr();
00105         cuts->SetFilter();
00106         cuts->SetV0MinDecayLen(0.);
00107         cuts->SetV0DaughterMinImpact(0);
00108         cuts->SetV0DaughterMinHit(10.);
00109         cuts->SetXiV0MaxImpact(5);
00110         cuts->SetXiMinDecayLen(2.);
00111         cuts->SetXiV0PiMinImpact(0.);
00112         cuts->SetXiDaughterMinHit(10.);
00113         cuts->SetKinkMinDecayRad(128.);
00114         cuts->SetKinkMaxDecayRad(184.);
00115       }
00116     }
00117   }
00118   // end additional maker code
00119 
00120   saveMk->cd();
00121 
00122   // Mixer
00123   gSystem->Load("StMixerMaker");
00124   StMixerMaker  *mixer = new StMixerMaker("Mixer","daq","trs");
00125   chain1->SetInput("Input1","StDAQReader");
00126   chain2->SetInput("Input2","Event");
00127   mixer->writeFile("mixer.trs",Nevents);
00128 
00129   //  gSystem->Load("StFtpcMixerMaker");
00130   //  StFtpcMixerMaker  *ftpcmixer = new StFtpcMixerMaker("FtpcMixer","daq","trs");
00131 
00132   // Create chain3 object
00133   chain3 = new StBFChain("Three");
00134   saveMk = chain3->cd();
00135 
00136   //  Base CuCu chain, fails. ELH: standard P06ib cucu options:
00137   //chain3->SetFlags("Simu NoDefault NoInput onlraw -onlcl P2005 DbV20060421 useCDV ITTF tofDat -SvtIT SCEbyE OGridLeak OShortR OSpaceZ2");
00138 
00139 // use Simu NoDefault NoInput onlraw -onlcl and standard chain options
00140 //  expand P2005, get rid of corrections and ftpc:
00141 //chain3->SetFlags("Simu NoDefault NoInput onlraw -onlcl DbV20050515 ry2005b tpc_daq tpc emcDY2 global dst Kalman event evout QA Tree GeantOut fcf ctf -Prevtx -nohits CMuDST ZDCvtx tofDat Xi2 Kink2 EST ToF svt_daq SvtD svtdEdx xiSvt l3onl fpd eemcD pmdRaw"); 
00142 //  chain3->SetFlags("Simu NoDefault NoInput onlraw -onlcl B2005 l3onl fcf emcDY2 fpd trgd ZDCvtx DbV20060421 useCDV ITTF tofDat -SvtIT SCEbyE");
00143 
00144   chain3->SetFlags("Simu NoDefault NoInput onlraw -onlcl ry2005b,tpc_daq,tpcI,svt_daq,SvtD,Physics,Idst,l0,tags,Tree,evout l3onl fcf emcDY2 fpd trgd ZDCvtx DbV20060421 useCDV ITTF tofDat -SvtIT MuDST -trg VFMinuit");
00145   //  chain3->SetFlags("Simu NoDefault NoInput onlraw -onlcl ry2005b,tpc_daq,tpcI,svt_daq,SvtD,Physics,Idst,l0,Tree,evout l3onl fcf emcDY2 fpd trgd ZDCvtx DbV20060421 useCDV ITTF tofDat -SvtIT MuDST -trg VFMinuit");
00146 
00147 
00148 
00149   //  StRTSClientFCF *fcfMk = (StRTSClientFCF *) chain3->GetMaker("");
00150   //  fcfMk->SetMode("0x1");
00151 
00152 
00153   TString OutputFileName(gSystem->BaseName(file1));
00154   OutputFileName.ReplaceAll("*","");
00155   OutputFileName.ReplaceAll(".daq","");
00156   OutputFileName.Append(".root");
00157   cout <<"AAR - Setting file output to: " <<OutputFileName.Data()<<endl;
00158   
00159   chain3->Set_IO_Files(0,OutputFileName.Data());
00160   chain3->Load();
00161   chain3->Instantiate();
00162   St_geant_Maker *geantMk = (St_geant_Maker *) chain->GetMaker("geant");
00163   geantMk->SetActive(kTRUE);
00164   StMaker *tpcdaqMk = chain3->GetMaker("tpc_raw");
00165   if(!tpcdaqMk )
00166     {
00167       cout <<" Error: no tpc daq maker. End. "<<endl;
00168       return;
00169     }
00170   tpcdaqMk->SetMode(1);   // Trs
00171   tpcdaqMk->SetInput("Event","MixerEvent");
00172 
00173   //  StMaker *ftpccluMk = chain3->GetMaker("ftpc_hits");
00174   //  ftpccluMk->SetInput("ftpc_raw","FtpcMixer");
00175   cout <<"AAR - going to file"<<endl;
00176 
00177   saveMk->cd();
00178   {
00179     TDatime t;
00180     printf ("QAInfo:Run is started at Date/Time%i/%i\n",t.GetDate(),t.GetTime());
00181   }
00182   printf ("QAInfo:Run on %s in %s\n",
00183         gSystem->HostName(),
00184         gSystem->WorkingDirectory());
00185   printf ("QAInfo: with %s\n", chain->GetCVS());
00186 
00187   // Init the chain and all its makers
00188   cout <<"AAR - initializing the chain"<<endl;
00189 
00190   chain->SetDEBUG(0);
00191   if (Nevents >= 0) {
00192     cout << "AAR - in init if"<<endl;
00193     Int_t iInit = chain->Init();
00194     cout <<"AAR - init value"<<iInit<<endl;
00195   }
00196 
00197  cout <<"AAR - Getting Tree"<<endl;
00198   // chain->SetDEBUG();
00199   treeMk = chain->GetMaker("tree");
00200   TBenchmark evnt;
00201   Int_t iMake = 0, i = 1, iBad = 0;
00202 
00203  cout <<"AAR - Setting IOMaker"<<endl;
00204 
00205   StIOMaker *inpMk = (StIOMaker *)chain1->GetMaker("inputStream");
00206   FILE *fp = fopen(file3,"r");
00207   Float_t x, y, z;
00208   Int_t ncols, eventnumber, mult, skip=0, oldskip = 0, skiptest=0;
00209   printf("zvertex_low = %f zvertex_high = %f\n",zvertex_low, zvertex_high);
00210 
00211   // vtxMk = (StVertexMaker*) chain3->GetMaker("vertex");
00212 
00213 
00214   cout <<"AAR - Entering Event Loop"<<endl;
00215 EventLoop: if (i <= Nevents && iMake != kStEOF && iMake != kStFatal) {
00216    evnt.Reset();
00217    evnt.Start("QAInfo:");
00218   printf ("ELH3\n");
00219    chain->Clear();
00220   printf ("ELH4\n");
00221    
00222    ncols = fscanf(fp,"%d %d %d %f %f %f",&eventnumber,&skip,&mult,&x,&y,&z);
00223    if(ncols<0) break;
00224 
00225    printf("\nEvent number: %d Multiplicity = %d\n",eventnumber,mult);
00226    printf("vertex in bfcMixer: %16f %16f %16f\n",x,y,z);
00227 
00228    // skip events in chain1 (daq file) only.
00229    // fz file does not have these events.
00230    // strange initialization due to details of inpMk->Skip
00231    // ...not well understood.  elh
00232 
00233    printf ("bfcMixer: i = %i\n",i);
00234    printf ("bfcMixer: skip = %i\n",skip);
00235    printf ("bfcMixer: oldskip = %i\n",oldskip);
00236    if (inpMk && skip>0) {
00237      if(i == 1) {skip++;}
00238      skiptest = inpMk->Skip(skip-oldskip-1);
00239      printf("bfcMixer: skiptest1 = %i\n",skiptest);
00240      skiptest = inpMk->Skip();
00241      printf("bfcMixer: skiptest2 = %i\n",skiptest);
00242      if(i == 1) {skip--;}
00243    }
00244    oldskip = skip;
00245 
00246    use this to force the vertex position to match:
00247    // vtxMk->FixVertex(x,y,z);
00248 
00249    iMake = chain->Make(i);
00250    if (treeMk && iMake == kStErr) {treeMk->Make(i); iBad++;}
00251    StEvtHddr *fEvtHddr = (StEvtHddr*)chain->GetDataSet("EvtHddr");
00252    StEvtHddr *fEvtHddrDaq = (StEvtHddr*)chain1->GetDataSet("EvtHddr");
00253    *fEvtHddr = *fEvtHddrDaq;
00254    // gSystem->Exec("ps ux");
00255    evnt.Stop("QAInfo:");
00256    // evnt->Show("QAInfo:");
00257    printf ("QAInfo: Done with Event [no. %d/run %d/evt. %d/Date.Time%d.%d/sta %d] Real Time = %10.2f seconds Cpu Time =  %10.2f seconds \n", i,chain->GetRunNumber(),chain->GetEventNumber(),chain->GetDate(), chain->GetTime(),
00258           iMake,evnt.GetRealTime("QAInfo:"),evnt.GetCpuTime("QAInfo:"));
00259 
00260 
00261    // Be sure to get the same event from daq and fz files...  elh
00262    if(eventnumber == chain->GetEventNumber()){
00263      printf("bfcMixer:  Event number %d found in both .fz file and .daq file\n",eventnumber);
00264    }else{
00265      printf("bfcMixer: Error! Event number from .fz file: %d\n",eventnumber);
00266      printf("bfcMixer: Error! Event number from .daq file: %d\n",chain->GetEventNumber());
00267      printf("bfcMixer: Error! Events are out-of-sync, exiting!\n");
00268      break;
00269    }
00270    i++;
00271    goto EventLoop;
00272  }
00273 
00274   fflush(stdout);
00275   printf ("QAInfo:Run completed ");
00276   gSystem->Exec("date");
00277   if (evMk) Event = (StEvent *) chain->GetInputDS("StEvent");
00278   {
00279     TDatime t;
00280     printf ("\nQAInfo:Run is finished at Date/Time%i/%i\n",t.GetDate(),t.GetTime());
00281   }
00282 }
00283
macros/examples/bfcMixer_v4_noFTPC.C
