StRoot: macros/embedding/bfcMixer

00001 
00002 //
00003 // $Id: bfcMixer_DaqFzd.C,v 1.2 2012/02/17 20:47:36 fisyak Exp $
00004 //
00005 // $Log: bfcMixer_DaqFzd.C,v $
00006 // Revision 1.2  2012/02/17 20:47:36  fisyak
00007 // Remove nodefault option from chain3
00008 //
00009 // Revision 1.1  2011/10/14 20:05:23  jeromel
00010 // Saved version ... floating aorund and being used or tested
00011 //
00012 //
00013 // JET EMBEDDING MACRO
00014 //
00015 // Pibero Djawotho <pibero@tamu.edu>
00016 // Texas A&M University
00017 // 27 July 2011
00018 //
00020 
00021 class StChain;
00022 StChain* Chain = 0;
00023 
00024 class StBFChain;
00025 StBFChain* chain1 = 0;
00026 StBFChain* chain2 = 0;
00027 StBFChain* chain3 = 0;
00028 
00029 //_____________________________________________________________________
00030 void bfcMixer_Jet2(const Int_t Nevents = 1000,
00031                   const Char_t* daqfile = "/star/data03/daq/2009/WembAll/st_zerobias_adc_10085134_raw_4450001.daq",
00032                   const Char_t* fzdfile = "/star/institutions/uky/gdwebb/pp500Embedding2009/output/pTbin_4_5/st_zerobias_adc_10085134_raw_4450001.fzd",
00033                   const Char_t* prodName = "P09igpp500")
00034 {
00035   // Production chains for P08ic - p+p, Au+Au 9 GeV and d+Au
00036   TString prodP08iepp("DbV20081117 B2008a ITTF IAna ppOpt l3onl emcDY2 fpd ftpc trgd ZDCvtx NosvtIT NossdIT Corr4 OSpaceZ2 OGridLeak3D VFMCE -hitfilt");
00037   TString prodP08iedAu("DbV20090213 P2008 ITTF OSpaceZ2 OGridLeak3D beamLine VFMCE TpxClu -VFMinuit -hitfilt");
00038   TString prodP10iapp("DbV20091001 pp2009c TpcRS ITTF OSpaceZ2 OGridLeak3D beamLine, VFMCE TpcRS -VFMinuit -hitfilt");
00039 
00040   // BES Run10 chains
00041   TString prodP10ihAuAu39("DbV20100909 P2010a,btof,BEmcChkStat,Corr4,OSpaceZ2,OGridLeak3D,VFMCE TpxClu -VFMinuit -hitfilt");
00042   TString prodP10ihAuAu11("DbV20100821 P2010a,btof,BEmcChkStat,Corr4,OSpaceZ2,OGridLeak3D,VFMCE TpxClu -VFMinuit -hitfilt");
00043   TString prodP10ihAuAu7("DbV20100821 P2010a,btof,BEmcChkStat,Corr4,OSpaceZ2,OGridLeak3D,VFMCE TpxClu -VFMinuit -hitfilt");
00044 
00045   // Run 9 500 GeV chain for W from Jan Balewski
00046   TString prodP09igpp500("DbV20091225 pp2009c ITTF VFPPVnoCTB  beamLine BEmcChkStat btof Corr4 OSpaceZ2 OGridLeak3D -dstout -evout");
00047 
00048   // Run 9 200 GeV chain for jet embedding
00049   TString prodP10icpp200("DbV20100130 pp2009c ITTF VFPPVnoCTB beamLine BEmcChkStat btof fmsdat Corr4 OSpaceZ2 OGridLeak3D");
00050 
00051   // Additional tags needed for embedding
00052   prodP09igpp500 += " TpxClu -VFMinuit VFPPVnoCTB beamLine -hitfilt";
00053   prodP10icpp200 += " TpxClu -VFMinuit VFPPVnoCTB beamLine -hitfilt";
00054 
00055   TString geomP08ic("ry2008");
00056   TString geomP10ih("ry2010");
00057   TString geomP09ig("ry2009a");
00058   TString geomP10ic("ry2009a");
00059   TString chain1Opt("in,magF,tpcDb,NoDefault,TpxRaw,-ittf,NoOutput");
00060   TString chain2Opt("fzin,gen_T,geomT,sim_T,TpcRS,-ittf,-tpc_daq,nodefault");
00061   chain2Opt += " ";
00062 
00063   TString chain3Opt;
00064   if (prodName == "P08icpp")           { chain3Opt = prodP08icpp;       chain2Opt += geomP08ic; }
00065   else if (prodName == "P08iepp")      { chain3Opt = prodP08iepp;       chain2Opt += geomP08ic; }
00066   else if (prodName == "P08icAuAu9")   { chain3Opt = prodP08icAuAu9;    chain2Opt += geomP08ic; }
00067   else if (prodName == "P08icdAu")     { chain3Opt = prodP08icdAu;      chain2Opt += geomP08ic; }
00068   else if (prodName == "P08iedAu")     { chain3Opt = prodP08iedAu;      chain2Opt += geomP08ic; }
00069   else if (prodName == "P08icAuAu200") { chain3Opt = prodP08icAuAu200;  chain2Opt += geomP08ic; }
00070   else if (prodName == "P10iapp")      { chain3Opt = prodP10iapp;       chain2Opt += geomP10ih; }
00071   else if (prodName == "P10ihAuAu39")  { chain3Opt = prodP10ihAuAu39;   chain2Opt += geomP10ih; }
00072   else if (prodName == "P10ihAuAu11")  { chain3Opt = prodP10ihAuAu11;   chain2Opt += geomP10ih; }
00073   else if (prodName == "P10ihAuAu7")   { chain3Opt = prodP10ihAuAu7;    chain2Opt += geomP10ih; }
00074   else if (prodName == "P09igpp500")   { chain3Opt = prodP09igpp500;    chain2Opt += geomP09ig; }
00075   else if (prodName == "P10icpp200")   { chain3Opt = prodP10icpp200;    chain2Opt += geomP10ic; }
00076   else {
00077     cout << "Choice prodName " << prodName << " does not correspond to known chain. Processing impossible." << endl;
00078     return;
00079   }
00080   chain3Opt += ",Embedding,TpcMixer,GeantOut,MiniMcMk,McAna,-in,NoInput,useInTracker"; 
00081   chain3Opt += ",";
00082 
00083   if (prodName == "P08icpp")           { chain3Opt += geomP08ic; }
00084   else if (prodName == "P08iepp")      { chain3Opt += geomP08ic; }
00085   else if (prodName == "P08icAuAu9")   { chain3Opt += geomP08ic; }
00086   else if (prodName == "P08icdAu")     { chain3Opt += geomP08ic; }
00087   else if (prodName == "P08iedAu")     { chain3Opt += geomP08ic; }
00088   else if (prodName == "P08icAuAu200") { chain3Opt += geomP08ic; }
00089   else if (prodName == "P10iapp")      { chain3Opt += geomP10ih; }
00090   else if (prodName == "P10ihAuAu39")  { chain3Opt += geomP10ih; }
00091   else if (prodName == "P10ihAuAu11")  { chain3Opt += geomP10ih; }
00092   else if (prodName == "P10ihAuAu7")   { chain3Opt += geomP10ih; }
00093   else if (prodName == "P09igpp500")   { chain3Opt += geomP09ig; }
00094   else if (prodName == "P10icpp200")   { chain3Opt += geomP10ic; }
00095   else {
00096     cout << "Choice prodName " << prodName << " does not correspond to known chain. Processing impossible. " << endl;
00097     return;
00098   }
00099 
00100   // Add BEMC simulators to chain
00101   chain3Opt += ",emcSim";
00102 
00103   // Add EEMC fast simulator to chain
00104   chain3Opt += ",EEfs";
00105 
00106   // Dynamically link some shared libs
00107   gROOT->LoadMacro("bfc.C");
00108   if (gClassTable->GetID("StBFChain") < 0) Load();
00109   //______________Create the main chain object______________________________________
00110   Chain = new StChain("Embedding");
00111   //________________________________________________________________________________
00112   bfc(-1,chain1Opt,daqfile);
00113   chain1 = chain;
00114   chain1->SetName("One"); 
00115   Chain->cd();
00116   //________________________________________________________________________________  
00117   bfc(-1,chain2Opt,fzdfile);
00118   chain2 = chain;
00119   chain2->SetName("Two");
00120   Chain->cd();
00121   if (chain2->GetOption("TRS")){
00122     StTrsMaker *trsMk = (StTrsMaker *) chain2->GetMaker("Trs");
00123     if (!trsMk) {
00124       cout << "Cannot find Trs in chain2" << endl;
00125       return;
00126     }
00127     trsMk->setNormalFactor(1.32);
00128     trsMk->SetMode(0);
00129   }
00130   //________________________________________________________________________________
00131   //  gSystem->Load("StFtpcMixerMaker");
00132   //  StFtpcMixerMaker  *ftpcmixer = new StFtpcMixerMaker("FtpcMixer","daq","trs");
00133   //________________________________________________________________________________
00134   TString OutputFileName(gSystem->BaseName(fzdfile));
00135   OutputFileName.ReplaceAll("*","");
00136   OutputFileName.ReplaceAll(".fzd","");
00137   //  OutputFileName.Append("_emb.root");
00138   OutputFileName.Append(".root");
00139   bfc(-1,chain3Opt,0,OutputFileName);
00140   chain3 = chain;
00141   chain3->SetName("Three"); 
00142   Chain->cd();
00143   //________________________________________________________________________________
00144   // Read raw TPX hits and generate simulated TPX hits AFTER
00145   // the event pass the trigger filter
00146   StTpcHitMaker* TpxRaw = (StTpcHitMaker*)chain1->GetMaker("TpxRaw");
00147   StTpcRSMaker* TpcRS = (StTpcRSMaker*)chain2->GetMaker("TpcRS");
00148   chain3->AddBefore("TpcMixer",TpcRS);
00149   chain3->AddBefore("TpcMixer",TpxRaw);
00150   //________________________________________________________________________________
00151   StTpcMixerMaker  *mixer = (StTpcMixerMaker *) chain3->Maker("TpcMixer");
00152   if( prodName == "P08icAuAu200")
00153         {
00154           mixer->SetInput("Input1","MixerEvent");
00155         }
00156   else
00157         {
00158           mixer->SetInput("Input1","TpxRaw/.data/Event");
00159         }
00160 
00161   if (chain2Opt.Contains("TpcRS",TString::kIgnoreCase)) {
00162     mixer->SetInput("Input2","TpcRS/Event");
00163   } else {
00164     mixer->SetInput("Input2","Trs/.const/Event");
00165   }
00166   Chain->cd();
00167 
00168   //------------------------------------ EMC MIXERS ------------------------------------
00169   // Add BEMC mixer to chain3
00170   StEmcRawMaker* emcRaw = (StEmcRawMaker*)chain3->GetMaker("emcRaw");
00171   emcRaw->getBemcRaw()->saveAllStEvent(true); // use all 4800 BEMC towers
00172   gSystem->Load("StEmcMixerMaker");
00173   StEmcMixerMaker* bemcMixer = new StEmcMixerMaker;
00174   chain3->AddAfter("EmcSimulator",bemcMixer);
00175   // Set EEMC fast and slow simulator in embedding mode
00176   StEEmcFastMaker* eefs = (StEEmcFastMaker*)chain3->GetMaker("eefs");
00177   eefs->SetEmcCollectionLocal(false); // don't create local StEmcCollection
00178   eefs->UseFullTower(true); // Use full ETOW detector
00179   StEEmcSlowMaker* eess = new StEEmcSlowMaker;
00180   eess->setEmbeddingMode(false); // don't use local StEmcCollection
00181   eess->setAddPed(true);    // add pedestals                 
00182   eess->setSmearPed(true);  // smear pedestals                 
00183   eess->setDropBad(true);   // force bad channels to be dropped in db
00184   eess->setOverwrite(true); // overwrite ADC values
00185   eess->setSource("StEvent");             
00186   // Add EEMC mixer to chain3
00187   //StEEmcMixerMaker* eemcMixer = new StEEmcMixerMaker;
00188   //------------------------------------------------------------------------------------
00189 
00190   //----------------------------- TRIGGER FILTER -----------------------------
00191   // We want to achieve the following ordering for makers:
00192   // 1. BBC simulator
00193   // 2. BEMC simulator
00194   // 3. BEMC mixer
00195   // 4. EEMC fast simulator
00196   // 5. EEMC slow simulator
00197   // 6. EEMC mixer
00198   // 7. Pythia event maker
00199   // 8. Trigger simulator
00200   // 9. Trigger filter
00201   // 10. TPC maker
00202 
00203   // Place TPC chain after EMC makers
00204   chain3->AddAfter("eefs",chain3->GetMaker("tpcChain"));
00205   //chain3->AddAfter("eefs",eemcMixer);
00206   chain3->AddAfter("eefs",eess);
00207 
00208   // Place Pythia maker after GEANT maker
00209   // and trigger filter after EMC makers
00210   gSystem->Load("StJetSkimEvent");
00211   gSystem->Load("StBfcTriggerFilterMaker");
00212 
00213   StPythiaEventMaker* pythia = new StPythiaEventMaker;
00214   TString pyfile = gSystem->BaseName(fzdfile);
00215   pyfile.ReplaceAll(".fzd",".pythia.root");
00216   pythia->SetPythiaFile(pyfile);
00217   chain3->AddAfter("geant",pythia);
00218 
00219   // Place trigger simulator after EMC makers
00220   gSystem->Load("StTriggerUtilities");
00221   StTriggerSimuMaker* trgsim = new StTriggerSimuMaker;
00222   trgsim->setMC(1);
00223   // BBC was not used in Run 9 jet triggers
00224   //trgsim->useBbc();
00225   //trgsim->bbc->setSource("StEvent");
00226   trgsim->useBemc();
00227   trgsim->bemc->setConfig(StBemcTriggerSimu::kOnline);
00228   trgsim->useEemc();
00229   trgsim->eemc->setSource("StEvent");
00230 
00231   trgsim->bemc->setBarrelJetPatchTh0(19);
00232   trgsim->bemc->setBarrelJetPatchTh1(26);
00233   trgsim->bemc->setBarrelJetPatchTh2(34);
00234 
00235   trgsim->emc->setOverlapJetPatchTh0(18);
00236   trgsim->emc->setOverlapJetPatchTh1(25);
00237   trgsim->emc->setOverlapJetPatchTh2(33);
00238 
00239   trgsim->eemc->setEndcapJetPatchTh0(17);
00240   trgsim->eemc->setEndcapJetPatchTh1(24);
00241   trgsim->eemc->setEndcapJetPatchTh2(31);
00242 
00243   trgsim->bemc->setBarrelHighTowerTh0(11);
00244   trgsim->bemc->setBarrelHighTowerTh1(15);
00245   trgsim->bemc->setBarrelHighTowerTh2(18);
00246   trgsim->bemc->setBarrelHighTowerTh3(24);
00247 
00248   trgsim->eemc->setEndcapHighTowerTh0(16);
00249   trgsim->eemc->setEndcapHighTowerTh1(25);
00250 
00251   StBfcTriggerFilterMaker* trgfilt = new StBfcTriggerFilterMaker;
00252   // The BFC trigger filter will select only JP1, AJP and BHT3 events
00253   trgfilt->SetJP1();
00254   trgfilt->SetAJP();
00255   trgfilt->SetBHT3();
00256   // Lower all jet patch thresholds by one unit from
00257   // their values obtained from the database using
00258   // the current timestamp.
00259   //trgfilt->changeJPThresh(-1);
00260   chain3->AddBefore("tpcChain",trgsim);
00261   chain3->AddBefore("tpcChain",trgfilt);
00262 #if 0
00263   //--------------------------------------------------------------------------
00264   {
00265     TDatime t;
00266     gMessMgr->QAInfo() << Form("Run is started at Date/Time %i/%i",t.GetDate(),t.GetTime()) << endm;
00267   }
00268   gMessMgr->QAInfo() << Form("Run on %s in %s",gSystem->HostName(),gSystem->WorkingDirectory()) << endm;
00269   gMessMgr->QAInfo() << Form("with %s", Chain->GetCVS()) << endm;
00270 
00271   TAttr::SetDebug(0);
00272   Chain->SetAttr(".Privilege",0,"*"                );   //All  makers are NOT priviliged
00273   Chain->SetAttr(".Privilege",1,"StBFChain::*" );       //StBFChain is priviliged
00274   Chain->SetAttr(".Privilege",1,"StIOInterFace::*" );   //All IO makers are priviliged
00275   Chain->SetAttr(".Privilege",1,"St_geant_Maker::*");   //It is also IO maker
00276   //Chain->SetAttr(".Privilege",1,"StPrepEmbedMaker::*"); //It is also IO maker
00277   //  Chain->SetDEBUG(0);
00278   if (Nevents < 0) return;
00279   Int_t iInit = Chain->Init();
00280   if (iInit >=  kStEOF) {Chain->FatalErr(iInit,"on init"); return;}
00281   StMaker *treeMk = Chain->GetMaker("outputStream");
00282   cout << "Order of makers in BFCMIXER:" << endl;
00283   StMaker::lsMakers(Chain);
00284   Chain->EventLoop(Nevents,treeMk);
00285   gMessMgr->QAInfo() << "Run completed " << endm;
00286   gSystem->Exec("date");
00287 #endif
00288 #if 0
00289   // Run chain
00290   int istat = chain->Init();
00291   if (istat) {
00292     cout << "Chain initialization failed" << endl;
00293     chain->Fatal(istat,"during Init()");
00294   }
00295   // Detailed listing of makers in the chain
00296   //chain->ls(0);
00297   cout << "Order of makers in BFC:" << endl;
00298   StMaker::lsMakers(chain);
00299   chain->EventLoop(Nevents);
00300 #endif
00301   TString trgfile = gSystem->BaseName(fzdfile);
00302   trgfile.ReplaceAll(".fzd",".trig.root");
00303   TFile* ofile = TFile::Open(trgfile,"recreate");
00304   assert(ofile);
00305   TH2F* hBarrelHighTowerSimu = new TH2F("hBarrelHighTowerSimu","BEMC high tower simu;trigger patch;high tower",300,0,300,64,0,64);
00306   TH2F* hBarrelPatchSumSimu = new TH2F("hBarrelPatchSumSimu","BEMC patch sum simu;trigger patch;patch sum",300,0,300,64,0,64);
00307   TH2F* hEndcapHighTowerSimu = new TH2F("hEndcapHighTowerSimu","EEMC high tower simu;trigger patch;high tower",90,0,90,64,0,64);
00308   TH2F* hEndcapPatchSumSimu = new TH2F("hEndcapPatchSumSimu","EEMC patch sum simu;trigger patch;patch sum",90,0,90,64,0,64);
00309   TH2F* hBarrelJetPatchSimu = new TH2F("hBarrelJetPatchSimu",";jet patch;adc",18,0,18,160,0,160);
00310   TH2F* hEndcapJetPatchSimu = new TH2F("hEndcapJetPatchSimu",";jet patch;adc",6,0,6,160,0,160);
00311   TH2F* hOverlapJetPatchSimu = new TH2F("hOverlapJetPatchSimu",";jet patch;adc",6,0,6,160,0,160);
00312   // Initialize chain
00313   Chain->Init();
00314   cout << "Order of makers in BFCMIXER:" << endl;
00315   StMaker::lsMakers(Chain);
00316   // Event loop
00317   for (int iEvent = 1; iEvent <= Nevents; ++iEvent) {
00318     Chain->Clear();
00319     int status = Chain->Make(iEvent);
00320     if (status == kStSkip) continue;
00321     if (status % 10 == kStEOF || status % 10 == kStFatal) break;
00322     // BEMC high towers and trigger patches
00323     for (int triggerpatch = 0; triggerpatch < 300; ++triggerpatch) {
00324       hBarrelHighTowerSimu->Fill(triggerpatch,trgsim->bemc->getBEMC_FEE_HT_ADC()[triggerpatch]);
00325       hBarrelPatchSumSimu->Fill(triggerpatch,trgsim->bemc->getBEMC_FEE_TP_ADC()[triggerpatch]);
00326     } // for triggerpatch
00327     // BEMC jet patches
00328     for (int jetpatch = 0; jetpatch < 18; ++jetpatch) {
00329       hBarrelJetPatchSimu->Fill(jetpatch,trgsim->bemc->barrelJetPatchAdc(jetpatch));
00330     } // for jetpatch
00331     // EEMC high towers and trigger patches
00332     for (int triggerpatch = 0; triggerpatch < 90; ++triggerpatch) {
00333       hEndcapHighTowerSimu->Fill(triggerpatch,trgsim->eemc->getOutHT(triggerpatch));
00334       hEndcapPatchSumSimu->Fill(triggerpatch,trgsim->eemc->getOutTPsum(triggerpatch));
00335     } // for triggerpatch
00336     // EEMC jet patches
00337     for (int jetpatch = 0; jetpatch < 6; ++jetpatch) {
00338       hEndcapJetPatchSimu->Fill(jetpatch,trgsim->eemc->endcapJetPatchAdc(jetpatch));
00339     } // for jetpatch
00340     // BEMC-EEMC-overlap jet patches
00341     for (int i = 0; i < 2; ++i) {
00342       int jetpatch, adc;
00343       trgsim->emc->getOverlapJetPatchAdc(i,jetpatch,adc);
00344       hOverlapJetPatchSimu->Fill(jetpatch,adc);
00345     } // for i
00346   } // for iEvent
00347   ofile->Write();
00348   ofile->Close();
00349 }
macros/embedding/bfcMixer_DaqFzd.C
