void TrsRead(const Int_t Nevents=1)
{
if (gClassTable->GetID("StChain") < 0) Load();
chain = new StChain("trs");
chain->SetDebug();
// Create the makers to be called by the current chain
// ROOT Db
const char* mainDB = "$STAR/StDb";
St_db_Maker *dbMk = new St_db_Maker("db",mainDB);
dbMk->SetDebug();
chain->SetInput("params","db:StDb/params");
const char* calibDB = "$STAR_ROOT/calib";
St_db_Maker *calibMk = new St_db_Maker("calib",calibDB);
chain->SetInput("calib","calib:calib");
calibMk->SetDebug();
StTrsMaker *tpc_raw = new StTrsMaker("Trs");
// MySQL DB
// const char *mainDB = "MySQL:Geometry";
// St_db_Maker *dbMk = new St_db_Maker("Geometry",mainDB);
// dbMk->SetDebug();
// dbMk->Init();
// dbMk->GetDataBase("Geometry/tpc");
// const char *calibDB = "MySQL:Calib";
// St_db_Maker *calibMk = new St_db_Maker("Calib",calibDB);
// calibMk->SetDebug();
// calibMk->Init();
// calibMk->GetDataBase("Calibrations/tpc");
// StTpcDbMaker *tpcDbMk = new StTpcDbMaker("tpcDb");
// tpcDbMk->Init();
// tpcDbMk->Make();
// cout << "The Db: " << gStTpcDb << endl;
// Tell TRS to read the file.
tpc_raw->readFile("test.trs");
// Init the main chain and all its makers
int iInit = chain->Init();
if (iInit) chain->Fatal(iInit,"on init");
gBenchmark->Start("trs");
for (Int_t i =1; i <= Nevents; i++) {
if (chain->Make(i)) break;
if (i != Nevents) chain->Clear();
printf ("===========================================\n");
printf ("=========================================== Done with Event no. %d\n",i);
printf ("===========================================\n");
}
if (Nevents > 1) {
chain->Finish();
gBenchmark->Print("trs");
}
else b = new TBrowser;
}
void runMuDstAna(TString inputFileList, int nFiles, const char* outName) {
gROOT->Macro("loadMuDst.C");
gSystem->Load("StTpcDb.so");
gSystem->Load("StDetectorDbMaker.so");
gSystem->Load("StDbUtilities");
gSystem->Load("StDbLib.so");
gSystem->Load("StDbBroker.so");
gSystem->Load("St_db_Maker");
gSystem->Load("StTriggerUtilities");
gSystem->Load("MuDstAna.so");
// List of member links in the chain
StChain* chain = new StChain;
StMuDstMaker* muDstMaker = new StMuDstMaker(0, 0, "", inputFileList.Data(), "MuDst", nFiles);
St_db_Maker* dbMk = 0;
dbMk = new St_db_Maker("StarDb", "MySQL:StarDb", "MySQL:StarDb", "$STAR/StarDb");
dbMk->SetDebug(0);
StDetectorDbTriggerID* detdbid = StDetectorDbTriggerID::instance();
// Turn off everything but Primary tracks in order to speed up the analysis and eliminate IO
muDstMaker-> SetStatus("*", 0); // Turn off all branches
muDstMaker-> SetStatus("PrimaryVertices", 1);
muDstMaker-> SetStatus("PrimaryTracks", 1);
muDstMaker-> SetStatus("MuEvent", 1); //.. turn on the useful branches
muDstMaker-> SetStatus("DetectorStates", 1);
cout << "input list " << inputFileList << endl;
MuDstAna* analysisMaker = new MuDstAna(muDstMaker, outName);
// Loop over the links in the chain
chain->Init();
chain->EventLoop(1,1000000);
chain->Finish();
// Cleanup
delete chain;
}
/**
* This is an alternative to bfcMixer_pythia() allowing the user to explicitly
* specify BFC options for all three chains.
*
* Sample input values:
*
* daqFileList = "@run10148002.list"
* fzdFile = "/path/to/my_pythia.fzd"
* flag = "W"
*
*/
void EmbeddingReco(const int nEvents, const std::string daqFileList, const std::string fzdFile,
const std::string bfcOptionsChain1, const std::string bfcOptionsChain2, const std::string bfcOptionsChain3,
const std::string flag)
{
// Dynamically link some shared libs
gROOT->LoadMacro("bfc.C");
if (gClassTable->GetID("StBFChain") < 0) Load();
//______________Create the main chain object______________________________________
Chain = new StChain("Embedding");
//________________________________________________________________________________
bfc(-1, bfcOptionsChain1.c_str(), daqFileList.c_str());
chain1 = chain;
chain1->SetName("One");
chain1->SetAttr(".call", "SetActive(0)", "St_db_Maker::"); // Use DB cache to reduce overhead
Chain->cd();
//________________________________________________________________________________
bfc(-1, bfcOptionsChain2.c_str(), fzdFile.c_str());
chain2 = chain;
chain2->SetName("Two");
Chain->cd();
if (chain2->GetOption("TRS")) {
StTrsMaker *trsMk = (StTrsMaker *) chain2->GetMaker("Trs");
if (!trsMk) {
cout << "Cannot find Trs in chain2" << endl;
return;
}
trsMk->setNormalFactor(1.32);
trsMk->SetMode(0);
}
//________________________________________________________________________________
// gSystem->Load("StFtpcMixerMaker");
// StFtpcMixerMaker *ftpcmixer = new StFtpcMixerMaker("FtpcMixer","daq","trs");
//________________________________________________________________________________
TString OutputFileName(gSystem->BaseName(fzdFile.c_str()));
OutputFileName.ReplaceAll("*", "");
OutputFileName.ReplaceAll(".fzd", "");
// OutputFileName.Append("_emb.root");
OutputFileName.Append(".root");
bfc(-1, bfcOptionsChain3.c_str(), 0, OutputFileName);
chain3 = chain;
chain3->SetName("Three");
Chain->cd();
//________________________________________________________________________________
StTpcMixerMaker *mixer = (StTpcMixerMaker *) chain3->Maker("TpcMixer");
mixer->SetInput("Input1", "TpxRaw/.data/Event");
if (bfcOptionsChain2.find("TpcRS") != std::string::npos) {
mixer->SetInput("Input2", "TpcRS/Event");
}
else {
mixer->SetInput("Input2", "Trs/.const/Event");
}
Chain->cd();
if (flag == "W") {
// blacklist some detectors to save DB load
St_db_Maker *dbMk = (St_db_Maker *)chain3->GetMaker("db");
dbMk->SetAttr("blacklist", "svt");
dbMk->SetAttr("blacklist", "ssd");
// dbMk->SetAttr("blacklist", "ftpc"); // S.F.
}
//------------------------------------ EMC MIXERS ------------------------------------
// Add BEMC mixer to chain3
StEmcRawMaker *emcRaw = (StEmcRawMaker *)chain3->GetMaker("emcRaw");
emcRaw->getBemcRaw()->saveAllStEvent(true); // use all 4800 BEMC towers
gSystem->Load("StEmcMixerMaker");
StEmcMixerMaker *bemcMixer = new StEmcMixerMaker;
chain3->AddAfter("EmcSimulator", bemcMixer);
// Set EEMC fast and slow simulator in embedding mode
StEEmcFastMaker *eefs = (StEEmcFastMaker *)chain3->GetMaker("eefs");
eefs->SetEmbeddingMode(); // Use local StEmcCollection
eefs->UseFullTower(true); // Use full ETOW detector
StEEmcSlowMaker *eess = new StEEmcSlowMaker;
eess->setEmbeddingMode(true);
// Add EEMC mixer to chain3
StEEmcMixerMaker *eemcMixer = new StEEmcMixerMaker;
//------------------------------------------------------------------------------------
//----------------------------- TRIGGER FILTER -----------------------------
// We want to achieve the following ordering for makers:
// 1. BBC simulator
// 2. BEMC simulator
// 3. BEMC mixer
// 4. EEMC fast simulator
// 5. EEMC slow simulator
// 6. EEMC mixer
// 7. Pythia event maker
// 8. Trigger simulator
// 9. Trigger filter
// 10. TPC maker
// Place TPC chain after EMC makers
chain3->AddAfter("eefs", chain3->GetMaker("tpcChain"));
chain3->AddAfter("eefs", eemcMixer);
chain3->AddAfter("eefs", eess);
if (flag == "Jet") {
// Place Pythia maker after GEANT maker
// and trigger filter after EMC makers
gSystem->Load("StJetSkimEvent");
gSystem->Load("StMCAsymMaker");
gSystem->Load("StBfcTriggerFilterMaker");
StPythiaEventMaker *pythia = new StPythiaEventMaker;
TString pyfile = gSystem->BaseName(fzdFile.c_str());
pyfile.ReplaceAll(".fzd", ".pythia.root");
pythia->SetPythiaFile(pyfile);
chain3->AddAfter("geant", pythia);
// Place trigger simulator after EMC makers
gSystem->Load("StTriggerUtilities");
StTriggerSimuMaker *trgsim = new StTriggerSimuMaker;
trgsim->setMC(1);
// BBC was not used in Run 9 jet triggers
//trgsim->useBbc();
//trgsim->bbc->setSource("StEvent");
trgsim->useBemc();
trgsim->bemc->setConfig(StBemcTriggerSimu::kOnline);
trgsim->useEemc();
trgsim->eemc->setSource("StEvent");
/*
trgsim->bemc->setBarrelJetPatchTh0(32);
trgsim->bemc->setBarrelJetPatchTh1(43);
trgsim->bemc->setBarrelJetPatchTh2(64);
*/
trgsim->setBarrelJetPatchTh(0, 32); // set for Run-11 based on trigger versioning page 20130728
trgsim->setBarrelJetPatchTh(1, 43);
trgsim->setBarrelJetPatchTh(2, 64);
trgsim->setOverlapJetPatchTh(0, 32);
trgsim->setOverlapJetPatchTh(1, 43);
trgsim->setOverlapJetPatchTh(2, 64);
trgsim->setEndcapJetPatchTh(0, 32);
trgsim->setEndcapJetPatchTh(1, 43);
trgsim->setEndcapJetPatchTh(2, 64);
trgsim->setBarrelHighTowerTh(0, 11);
trgsim->setBarrelHighTowerTh(1, 18);
trgsim->setBarrelHighTowerTh(2, 25);
trgsim->setBarrelHighTowerTh(3, 31);
trgsim->setEndcapHighTowerTh(0, 25);
trgsim->setEndcapHighTowerTh(1, 31);
// StBfcTriggerFilterMaker* trgfilt = new StBfcTriggerFilterMaker;
// no trigger filter for Run-11 VPDMB
// trgfilt->SetOkAllEvents(1);
// The BFC trigger filter will select only JP1, AJP and BHT3 events
// trgfilt->SetJP1();
// trgfilt->SetAJP();
// trgfilt->SetBHT3();
// Lower all jet patch thresholds by one unit from
// their values obtained from the database using
// the current timestamp.
//trgfilt->changeJPThresh(-1);
chain3->AddBefore("tpcChain", trgsim);
// chain3->AddBefore("tpcChain",trgfilt); // no trigger filter for Run-11 VPDMB
// Move these makers after trigger decision
// *** VERY IMPORTANT ***
// The order of TpxRaw and TpcRS *must* be preserved
// or the embedding will *not* work. [RT# 2299]
// http://www.star.bnl.gov/rt2/Ticket/Display.html?id=2299
StTpcRSMaker *TpcRS = (StTpcRSMaker *)chain2->GetMaker("TpcRS");
StTpcHitMaker *TpxRaw = (StTpcHitMaker *)chain1->GetMaker("TpxRaw");
chain3->AddBefore("TpcMixer", TpxRaw);
chain3->AddBefore("TpcMixer", TpcRS);
#if 0
// Turn on debugging of DB maker
St_db_Maker *db = (St_db_Maker *)chain1->GetMaker("db");
db->SetDebug(2);
#endif
//--------------------------------------------------------------------------
TString trgfile = gSystem->BaseName(fzdFile.c_str());
trgfile.ReplaceAll(".fzd", ".trig.root");
TFile *ofile = TFile::Open(trgfile, "recreate");
assert(ofile);
TH2F *hBarrelHighTowerSimu = new TH2F("hBarrelHighTowerSimu", "BEMC high tower simu;trigger patch;high tower", 300, 0, 300, 64, 0, 64);
TH2F *hBarrelPatchSumSimu = new TH2F("hBarrelPatchSumSimu", "BEMC patch sum simu;trigger patch;patch sum", 300, 0, 300, 64, 0, 64);
TH2F *hEndcapHighTowerSimu = new TH2F("hEndcapHighTowerSimu", "EEMC high tower simu;trigger patch;high tower", 90, 0, 90, 64, 0, 64);
TH2F *hEndcapPatchSumSimu = new TH2F("hEndcapPatchSumSimu", "EEMC patch sum simu;trigger patch;patch sum", 90, 0, 90, 64, 0, 64);
TH2F *hBarrelJetPatchSimu = new TH2F("hBarrelJetPatchSimu", "BEMC jet patch;jet patch;adc", 18, 0, 18, 160, 0, 160);
TH2F *hEndcapJetPatchSimu = new TH2F("hEndcapJetPatchSimu", "EEMC jet patch;jet patch;adc", 6, 0, 6, 160, 0, 160);
TH2F *hOverlapJetPatchSimu = new TH2F("hOverlapJetPatchSimu", "BEMC-EEMC-overlap;jet patch;adc", 6, 0, 6, 160, 0, 160);
}
//--------------------------------------------------------------------------
// Initialize chain
Chain->Init();
PrintTimer(Chain);
puts("Order of makers in BFCMIXER:");
StMaker::lsMakers(Chain);
// Event loop
int mNTotal = 0;
int mNFailed = 0;
TBenchmark evnt;
StIOMaker *inputStream = (StIOMaker *)chain1->GetMaker("inputStream");
for (int iEvent = 1; iEvent <= nEvents; ++iEvent) {
evnt.Reset();
evnt.Start("QAInfo:");
Chain->Clear();
int iMake = Chain->Make(iEvent);
if (iMake == kStErr) ++mNFailed;
if (inputStream->GetMakeReturn() % 10 == kStEOF) {
inputStream->Rewind();
--iEvent;
continue;
}
//if (iMake == kStSkip) continue;
if (iMake % 10 == kStEOF || iMake % 10 == kStFatal) break;
++mNTotal;
PrintTimer(Chain);
//--------------------------------------------------------------------------
if (flag == "Jet") {
// BEMC high towers and trigger patches
for (int triggerpatch = 0; triggerpatch < 300; ++triggerpatch) {
hBarrelHighTowerSimu->Fill(triggerpatch, trgsim->bemc->getBEMC_FEE_HT_ADC()[triggerpatch]);
hBarrelPatchSumSimu->Fill(triggerpatch, trgsim->bemc->getBEMC_FEE_TP_ADC()[triggerpatch]);
} // for triggerpatch
// BEMC jet patches
for (int jetpatch = 0; jetpatch < 18; ++jetpatch) {
hBarrelJetPatchSimu->Fill(jetpatch, trgsim->bemc->barrelJetPatchAdc(jetpatch));
} // for jetpatch
// EEMC high towers and trigger patches
for (int triggerpatch = 0; triggerpatch < 90; ++triggerpatch) {
hEndcapHighTowerSimu->Fill(triggerpatch, trgsim->eemc->getEndcapHighTower(triggerpatch));
hEndcapPatchSumSimu->Fill(triggerpatch, trgsim->eemc->getEndcapPatchSum(triggerpatch));
} // for triggerpatch
// EEMC jet patches
for (int jetpatch = 0; jetpatch < 6; ++jetpatch) {
hEndcapJetPatchSimu->Fill(jetpatch, trgsim->eemc->endcapJetPatchAdc(jetpatch));
} // for jetpatch
// BEMC-EEMC-overlap jet patches
for (int i = 0; i < 2; ++i) {
int jetpatch, adc;
trgsim->emc->getOverlapJetPatchAdc(i, jetpatch, adc);
hOverlapJetPatchSimu->Fill(jetpatch, adc);
} // for i
} //end of if("jet")
//--------------------------------------------------------------------------
evnt.Stop("QAInfo:");
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", iEvent, Chain->GetRunNumber(), Chain->GetEventNumber(), Chain->GetDate(), Chain->GetTime(), chain3->GetMakeReturn(), evnt.GetRealTime("QAInfo:"), evnt.GetCpuTime("QAInfo:"));
} // End event loop
printf("QAInfo:EventLoop completed code %d\n", iMake);
gSystem->Exec("date");
TDatime t;
printf("QAInfo:Run is finished at Date/Time %i/%i; Total events processed: %i and not completed: %i\n", t.GetDate(), t.GetTime(), mNTotal, mNFailed);
//--------------------------------------------------------------------------
if (flag == "Jet") {
ofile->Write();
ofile->Close();
}
}
int fgtRawSpectraQA( const Char_t *filenameIn = "testfile.daq",
const Char_t *filebaseOut = "testfile.rawQA",
Int_t date = 20120115,
Int_t time = 0,
Float_t pedRelThres = 0, // 0 = no cut, otherwise: # of sigma above pedestal for cut
Int_t nevents = 2000,
Int_t timebin = 2,
UInt_t statusMask,
Bool_t cutShortEvents = 1 ){
LoadLibs();
Int_t ierr = 0;
//
// START CONSTRUCTING THE CHAIN
//
cout << "Constructing the chain" << endl;
analysisChain = new StChain("eemcAnalysisChain");
//
// THE DATABASE
//
// note: DB is used to convert elec coords into geoIds in
// StEvent/StFgtStrip, and do determine the disc and quad given
// an rdo/arm/apv combination.
cout << "Constructing St_db_Maker" << endl;
TString dir0 = "MySQL:StarDb";
TString dir1 = "$STAR/StarDb";
St_db_Maker *dbMkr = new St_db_Maker( "dbMkr", dir0, dir1 );
dbMkr->SetDateTime(date,time);
cout << "Constructing StFgtDbMaker" << endl;
fgtDbMkr = new StFgtDbMaker( "fgtDbMkr" );
//
// NOW THE OTHER READERS AND MAKERS
//
cout << "Constructing the daq reader" << endl;
daqRdr = new StFgtRawDaqReader( "daqReader", filenameIn, "" );
daqRdr->setIsCosmic( 0 );
daqRdr->cutShortEvents( cutShortEvents );
// a2cMkr = new StFgtA2CMaker( "a2cMkr" );
// a2cMkr->setStatusMask( statusMask );
// a2cMkr->setAbsThres( -10000 ); // set to below -4096 to skip cut
// a2cMkr->setRelThres( pedRelThres ); // set to zero to skip cut
cout << "Constructing the QA Makers" << endl;
qaMkr = new StFgtQaRawOctAdc( "qaMkr", 2 );
qaMkr->setTimeBin( 2 );
// debug
analysisChain->ls(4);
cout << "Initializing" << endl;
ierr = analysisChain->Init();
if( ierr ){
cout << "Error initializing" << endl;
return;
};
if( nevents < 0 )
nevents = 1<<30; // a big number
cout << "max nevents = " << nevents << endl;
for( int i=0; i<nevents && !ierr; ++i ){
if( i+1 % 100 == 0 )
cout << "\ton event number " << i << endl;
//cout << "clear" << endl;
analysisChain->Clear();
//cout << "make" << endl;
ierr = analysisChain->Make();
};
//
// Calls the ::Finish() method on all makers
//
cout << "finish" << endl;
analysisChain->Finish();
cout << "Making plots" << endl;
gROOT->SetStyle("Plain");
gStyle->SetOptStat(0);
gStyle->SetPalette(1);
can = new TCanvas( "fgtRawQA", "fgtRawQA", 850, 1100);
can->Divide( 1, 4 );
// open output file
can->Print( (std::string(filebaseOut) + ".ps[").data() );
//std::vector< TH2F* >& hist = qaMkr->getHistVec();
cout << "First hist at " << qaMkr->getHist(0) << endl;
cout << "Finding max" << endl;
// get max
Float_t max = 0;
for( Int_t idx = 0; idx < 48; ++idx ){
if( qaMkr->getHist(idx) )
if( qaMkr->getHist(idx)->GetMaximum() > max )
max = qaMkr->getHist(idx)->GetMaximum();
};
max = 0.9*max;
cout << "Max is " << max << endl;
TH2F *dummy[4];
const Char_t dummyNames[4][20] = { "dummy1", "dummy2", "dummy3", "dummy4" };
Float_t xMin = qaMkr->getHist(0)->GetXaxis()->GetXmin();
Float_t xMax = qaMkr->getHist(0)->GetXaxis()->GetXmax();
Float_t yMin = qaMkr->getHist(0)->GetYaxis()->GetXmin();
Float_t yMax = qaMkr->getHist(0)->GetYaxis()->GetXmax();
for( Int_t i=0; i<4; ++i )
dummy[i] = new TH2F( dummyNames[i], "", 5, xMin, xMax, 1, yMin, yMax );
std::stringstream ss;
Int_t subpad = 1;
idx = 0;
for( Int_t rdo = 1; rdo < 3; ++rdo ){
for( Int_t arm = 0; arm < 6; ++arm ){
for( Int_t startIdx = 0; startIdx < 4; ++startIdx, ++idx ){
for( Int_t oct = 0; oct < 1; ++oct ){
cout << "rdo/arm/disc/oct = " << rdo << '/' << arm << '/' << startIdx << '/' << oct << endl;
if( qaMkr->getHist(idx)->GetEntries() > 0 ){
can->cd(subpad);
gPad->SetLeftMargin( 0.06 );
gPad->SetRightMargin( 0.05 );
gPad->SetBottomMargin( 0.11 );
qaMkr->getHist(idx)->SetMaximum( max );
dummy[subpad-1]->GetXaxis()->SetTitleOffset(0.9);
dummy[subpad-1]->GetXaxis()->SetTitleSize(0.06);
dummy[subpad-1]->GetXaxis()->SetLabelSize(0.06);
dummy[subpad-1]->GetYaxis()->SetTitleOffset(0.55);
dummy[subpad-1]->GetYaxis()->SetTitleSize(0.06);
dummy[subpad-1]->GetYaxis()->SetLabelSize(0.06);
gPad->SetGridx(0);
gPad->SetGridy(0);
for( Int_t i = 0; i<5; ++i ){
ss.str("");
ss.clear();
ss << "channels in APV "; // startArray[startIdx]+i;
dummy[subpad-1]->GetXaxis()->SetBinLabel( i+1, ss.str().data() );
};
dummy[subpad-1]->GetXaxis()->SetNdivisions(222,0);
dummy[subpad-1]->SetMinimum( 0 );
dummy[subpad-1]->SetMaximum( max );
dummy[subpad-1]->SetTitle( qaMkr->getHist(idx)->GetTitle() );
qaMkr->getHist(idx)->SetTitle("");
dummy[subpad-1]->GetYaxis()->SetTitle( qaMkr->getHist(idx)->GetYaxis()->GetTitle() );
qaMkr->getHist(idx)->GetYaxis()->SetTitle("");
dummy[subpad-1]->Draw("COLZ");
qaMkr->getHist(idx)->Draw("COLZ SAME");
gPad->Update();
gPad->Modified();
TPaveText *title = (TPaveText*)(gPad->GetPrimitive("title"));
if( title ){
title->SetX1NDC( 0.045 );
title->SetX2NDC( 0.55 );
title->SetY1NDC( 0.91 ) ;
title->SetY2NDC( 0.999 );
title->SetBorderSize(0);
title->SetTextAlign( 12 );
title->SetTextColor(kBlue);
title->Draw();
};
TPave *palette = (TPave*)(gPad->GetPrimitive("palette"));
if( palette ){
palette->SetX1NDC( 0.955 );
palette->SetX2NDC( 0.985 );
palette->Draw();
};
++subpad;
if( subpad == 5 ){
subpad = 1;
can->Print( (std::string(filebaseOut) + ".ps").data() );
};
};
};
};
};
};
if( subpad != 1 )
can->Print( (std::string(filebaseOut) + ".ps").data() );
can->Print( (std::string(filebaseOut) + ".ps]").data() );
gSystem->Exec(( std::string("ps2pdf -dAutoRotatePages=/None ") + filebaseOut + ".ps" ).data());
cerr << "\tall done" << endl;
return;
};
void rdMu2TrigSimu(char *file="/star/data47/reco/pp200/pythia6_410/9_11gev/cdf_a/y2006c/gheisha_on/p07ic/rcf1309_*_2000evts.MuDst.root"){
int nevents = 200;
int flagMC=1; // 0/1 == Using Real/Simulation data files
int useEemc=1; // 0/1 == Exclude/Include EEMC in Trigger Decisions
int useBemc=1; // 0/1 == Exclude/Include BEMC in Trigger Decisions
int useL2=1; // 0/1 == Exclude/Include L2 in Trigger Decisions
int L2ConfigYear=2006; // possible: 2006, 2008
int bemcConfig=2; // Online==1, Offline==2, Expert==3
int playConfig=0; // jan:100_199
int emcEveDump=0; // extrating raw EMC data in a custom format
int outputL2Histo=0;//output L2 histos to directory outL2
TString outDir="./outL2/";
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
assert( !gSystem->Load("StDetectorDbMaker"));
assert( !gSystem->Load("StDbUtilities"));
assert( !gSystem->Load("StDbBroker"));
assert( !gSystem->Load("St_db_Maker"));
assert( !gSystem->Load("StEEmcUtil")); // needed by eemcDb
assert( !gSystem->Load("StEEmcDbMaker"));
assert( !gSystem->Load("StDaqLib")); // needed by bemcDb
assert( !gSystem->Load("StEmcRawMaker"));
assert( !gSystem->Load("StEmcADCtoEMaker"));
if (flagMC) {
assert( !gSystem->Load("StMcEvent"));
assert( !gSystem->Load("StMcEventMaker"));
assert( !gSystem->Load("StEmcSimulatorMaker"));
assert( !gSystem->Load("StEEmcSimulatorMaker"));
assert( !gSystem->Load("StEpcMaker"));
}
assert( !gSystem->Load("StTriggerUtilities"));
gROOT->Macro("LoadLogger.C");
cout << " loading done " << endl;
chain= new StChain("StChain");
if (flagMC){
TString geantFile;
geantFile += file;
geantFile.ReplaceAll("MuDst.root", "geant.root");
printf("geantFile=%s\n", geantFile.Data());
StIOMaker* ioMaker = new StIOMaker();
ioMaker->SetFile(geantFile);
ioMaker->SetIOMode("r");
ioMaker->SetBranch("*",0,"0"); //deactivate all branches
ioMaker->SetBranch("geantBranch",0,"r"); //activate geant Branch
StMcEventMaker *evtMaker = new StMcEventMaker();
}
//Need MuDstMaker to get data
printf(" Analyzing file=%s\n",file);
StMuDstMaker* muDstMaker =new StMuDstMaker(0,0,"",file,"",1000);
//Database -- get a real calibration from the database
St_db_Maker* dbMk =0;
if(useEemc || useL2) // full DB access
dbMk = new St_db_Maker("StarDb","MySQL:StarDb","MySQL:StarDb","$STAR/StarDb");
else // only Barrel is uploaded, is faster
dbMk = new St_db_Maker("Calibrations","MySQL:Calibrations_emc");
//If MC then must set database time and date
//If Endcap fast simu is used tower gains in DB do not matter,JB
if (flagMC) dbMk->SetDateTime(20060522, 55000);//timestamp R7142018
//Collect all output histograms
TObjArray* HList=new TObjArray;
//Endcap DB
if(useEemc || useL2) new StEEmcDbMaker("eemcDb");
//Get BEMC adc values
if (flagMC && useBemc) {
StEmcSimulatorMaker* emcSim = new StEmcSimulatorMaker(); //use this instead to "redo" converstion from geant->adc
if (bemcConfig == 1) {
emcSim->setCheckStatus(kBarrelEmcTowerId,false); //this returns hits regardless of offline tower status
}
emcSim->setCalibSpread(kBarrelEmcTowerId,0.15);//spread gains by 15%
}
if (flagMC==0 && useBemc){
StEmcADCtoEMaker *bemcAdc = new StEmcADCtoEMaker();//for real data this sets calibration and status
if (bemcConfig == 1) {
bemcAdc->setCheckStatus(kBarrelEmcTowerId,false);
}
}
//must use slow simulator to get pedestals correct for L2
if (flagMC==1 && useEemc){
StEEmcSlowMaker *slowSim = new StEEmcSlowMaker("slowSim");
slowSim->setSamplingFraction(0.0384); // effectively scales all Tower energies with a factor of 1.3 (added by: Ilya Selyuzhenkov; April 11, 2008)
slowSim->setAddPed(true);
slowSim->setSmearPed(true);
}
//Get TriggerMaker
StTriggerSimuMaker *simuTrig = new StTriggerSimuMaker("StarTrigSimu");
simuTrig->setHList(HList);
simuTrig->setMC(flagMC); // must be before individual detectors, to be passed
simuTrig->useBbc();
if(useEemc) {
simuTrig->useEemc(0);//default=0:just process ADC, 1,2:comp w/trgData,see .
simuTrig->eemc->setSetupPath(eemcSetupPath);
}
if(useBemc){
simuTrig->useBemc();
simuTrig->bemc->setConfig(bemcConfig);
}
if(flagMC && useEemc==2){
// pass one argument to M-C as generic switch
// Endcap specific params -- ok Jan you need to change this to a default "online" setup
int eemcDsmSetup[20]; // see StEemcTriggerSimu::initRun() for definition
memset(eemcDsmSetup, 0,sizeof(eemcDsmSetup));// clear all, may be a bad default
eemcDsmSetup[0]=3; // HTthr0
eemcDsmSetup[1]=12; // HTthr1
eemcDsmSetup[2]=22; // HTthr2
eemcDsmSetup[3]=1; // TPthr0
eemcDsmSetup[4]=17; // TPthr1
eemcDsmSetup[5]=31; // TPthr2
eemcDsmSetup[10]=2; //HTTPthrSelc, 2=use_thres_#1
simuTrig->eemc->setDsmSetup(eemcDsmSetup);
}
if(useL2) {
/*
reads all input/setup files from L2setup-yyyymmdd/
writes all output files to L2out-yyyymmdd
depending on the DB time stamp
both dierectiorie MUST exist, setup must be reasonable
*/
StGenericL2Emulator* simL2Mk=0;
if(L2ConfigYear==2006) simL2Mk= new StL2_2006EmulatorMaker;
else if(L2ConfigYear==2008) simL2Mk= new StL2_2008EmulatorMaker;
assert(simL2Mk);
simL2Mk->setSetupPath(eemcSetupPath);
simL2Mk->setOutPath(outDir.Data());
if (flagMC) simL2Mk->setMC();
simuTrig->useL2(simL2Mk);
}
//if(emcEveDump) new StJanEventMaker;
StTriggerSimuPlayMaker *playMk= new StTriggerSimuPlayMaker; // to develope user analysis of trigQA
playMk->setConfig(playConfig);
playMk->setHList(HList);
chain->ls(3);
chain->Init();
for (Int_t iev=0;iev<nevents; iev++) {
cout << "\n****************************************** " << endl;
cout << "Working on eventNumber:\t" << iev <<"\tof:\t"<<nevents<<endl;
cout << "****************************************** " << endl;
chain->Clear();
int iret = chain->Make(iev);
total++;
if (iret % 10 == kStEOF || iret % 10 == kStFatal) {
cout << "Bad return code!" << endl;
break;
}
int trigID[3]={127213,117211,137611};
StMuDst *muDst = muDstMaker->muDst();
StMuEvent *muEvent = muDst->event();
StMuTriggerIdCollection trig = muEvent -> triggerIdCollection();
StTriggerId l1trig = trig.nominal();
if( l1trig.isTrigger(trigID[0])) {
cout<<" SimuTrigger 127213 ="<<simuTrig->isTrigger(trigID[0])<<" BEMC="<<simuTrig->bemc->triggerDecision(trigID[0])<<" L2="<<simuTrig->lTwo->triggerDecision(trigID[0])<<endl;
}
if( l1trig.isTrigger(trigID[1])) {
cout<<" SimuTrigger 117211 ="<<simuTrig->isTrigger(trigID[1])<<" BEMC="<<simuTrig->bemc->triggerDecision(trigID[1])<<" L2="<<simuTrig->lTwo->triggerDecision(trigID[1])<<endl;
}
if( l1trig.isTrigger(trigID[2])) {
cout<<" SimuTrigger 137611 ="<<simuTrig->isTrigger(trigID[2])<<" BEMC="<<simuTrig->bemc->triggerDecision(trigID[2])<<" L2="<<simuTrig->lTwo->triggerDecision(trigID[2])<<endl;
}
StTriggerSimuResult trigResult = simuTrig->detailedResult(trigID[2]);
if (trigResult.bemcDecision()==1){
vector<short> towerId = trigResult.highTowerIds();
for (unsigned i=0; i<towerId.size(); i++) {
cout<<" LO Trigger BEMC Tower="<<towerId[i]<<" adc="<<trigResult.highTowerAdc(towerId[i])<<endl;
}
}
if (trigResult.l2Decision()==1){
vector<short> towerId = trigResult.highTowerIds();
for (unsigned i=0; i<towerId.size(); i++) {
cout<<" L2 Trigger BEMC Tower="<<towerId[i]<<" adc="<<trigResult.highTowerAdc(towerId[i])<<endl;
}
}
}
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << total << endl;
cout << "****************************************** " << endl;
if (outputL2Histo==1) {
TString fileMu=file;
printf("=%s=\n",fileMu.Data());
if(fileMu.Contains(".lis")) fileMu.ReplaceAll(".lis",".trgSim");
if(fileMu.Contains(".MuDst.root")) fileMu.ReplaceAll(".MuDst.root",".trgSim");
TString outF=outDir+fileMu;
outF+=".hist.root";
printf("=%s=\n",outF.Data());
hf=new TFile(outF,"recreate");
if(hf->IsOpen()) {
HList->Write();
printf("\n Histo saved -->%s<\n",outF.Data());
}
else {
printf("\n Failed to open Histo-file -->%s<, continue\n",outF.Data());
}
}
//cout <<Form("sorting done %d of nEve=%d, CPU rate=%.1f Hz, total time %.1f minute(s) \n\n",total,nEntries,rate,tMnt)<<endl;
}
void
show
(
char* inpDir = "", // MuDST directory
char* inpFile = "show.lis", // MuDST file(s);
char* outFile = "show.root",// output tree file
Int_t nFiles = 50, // # of MuDST file(s)
Int_t nEvents = 100 // # of events
)
// remeber to adjust dbase timestamp below !!!!
// what a ... design
{
//gErrorIgnoreLevel=1999;
// load root/root4star libraries
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
// load more libraries :)
gSystem->Load("libmysqlclient");
gSystem->Load("libGeom");
gSystem->Load("StDbLib");
gSystem->Load("StDbBroker");
gSystem->Load("St_db_Maker");
// load even more libraries (EEMC stuff)
gSystem->Load("StEEmcUtil");
gSystem->Load("StEEmcDbMaker");
gSystem->Load("StEEmcPoolTTM");
// create the chain
chain = new StChain("StChain");
//
now = new TDatime;
// for display
TCanvas *c1 = new TCanvas("eemc","eemc",10,10,1000,1000);
TPaveLabel *tlab = new TPaveLabel(-0.99,+0.99,+0.99,+0.90,
"EEMC TOWERS & TPC TRACKS Piotr A Zolnierczuk (IU)");
eventInfo = new TPaveText (-0.99,-0.99,+0.0 ,-0.75);
dateInfo = new TPaveLabel(+0.60,-0.99,+0.99,-0.95,now->AsString());
TGeoManager *gm = new TGeoManager("eemc", "eemc tower display");
TGeoVolume *top = gm->MakeBox("star",0, 200., 200., 350.);
TGeoVolume *smbox = gm->MakeBox("smbox1",0, 2., 2., 2.);
smbox->SetLineColor(kRed);
// eemc
eemc = new EEmcTTDisplay();
eemc->SetMagneticField(0.5); // in Tesla
eemc->SetShowExtrapolatedTracks(true);
TGeoTranslation *etra = new TGeoTranslation(0.0,0.0,0.5*(eemc->getZ1()+eemc->getZ2()));
top->AddNode(smbox, 1,NULL);
top->AddNode(eemc(),1,etra);
gm->SetTopVolume(top);
gm->CloseGeometry();
gm->SetVisLevel(4);
gm->SetVisOption(0);
c1->SetTheta(90);
c1->SetPhi(0);
top->Draw();
tlab->Draw();
gPad->Update();
// now we add Makers to the chain... some of that is black magic :)
muDstMk = new StMuDstMaker(0,0,inpDir,inpFile,"",nFiles); // muDST main chain
StMuDbReader *db = StMuDbReader::instance(); // need the database
St_db_Maker *dbMk = new St_db_Maker("StarDb", "MySQL:StarDb"); // need another db(?)
new StEEmcDbMaker("eemcDb"); // need EEMC database
// now comment in/out/change the below if you want it your way
dbMk->setTimeStampDay(20040331); // format: yyyymmdd
// finally after so many lines we arrive at the good stuff
ttm = new EEmcTTMMaker ("TTM",muDstMk,eemcDbMk);
ttm->Summary(cout); //
StMuDebug::setLevel(0);
chain->Init();
StEEmcDb *eemcDb = (StEEmcDb*)chain->GetDataSet("StEEmcDb");
eemcDb->setSectors(1,12); // request EEMC DB for sectors you need (dafault:1-12)
eemcDb->setPreferedFlavor("onlped","eemcPMTped"); // request alternative flavor
chain->ls(3);
//---------------------------------------------------
next();
}
void TestIOMaker(Int_t nevents=3,
//This file points to 30 events of 10 neg muons w/ pt=.9
const char* MainFile="/star/data22/ITTF/data/simple_geant/DEV_10_8_01/*.event.root")
//This file points to 5 muons /event
//const char* MainFile="/star/data22/ITTF/EvalData/MCNtuple/muon_100_neg.event.root")
//This file points to 110 events from mevsim (homebrew had. cocktail)
//const char* MainFile="/star/data22/ITTF/data/mevsim/10_9_01/*.event.root")
//This file points to a nightly low density hadronic cocktail reconstruction.
//const char* MainFile="/star/rcf/test/dev/trs_redhat61/Tue/year_2001/hc_lowdensity/*.event.root")
{
// Dynamically link needed shared libs
cout <<"Loading St_base"<<endl;
gSystem->Load("St_base");
cout <<"Loading StChain"<<endl;
gSystem->Load("StChain");
cout <<"Loading St_Tables"<<endl;
gSystem->Load("St_Tables");
cout <<"Loading StUtilities"<<endl;
gSystem->Load("StUtilities");
cout <<"Loading StIOMaker"<<endl;
gSystem->Load("StIOMaker");
cout <<"Loading StarClassLibrary"<<endl;
gSystem->Load("StarClassLibrary");
cout <<"Loading DataBase"<<endl;
gSystem->Load("StDbUtilities");
gSystem->Load("StDbLib");
gSystem->Load("StDbBroker");
gSystem->Load("St_db_Maker");
gSystem->Load("StTpcDb");
cout <<"Loading StEvent"<<endl;
gSystem->Load("StEvent");
cout <<"Loading StEventMaker"<<endl;
gSystem->Load("StEventMaker");
cout <<"Loading StEmcUtil"<<endl;
gSystem->Load("StEmcUtil");
cout <<"Loading StMcEvent"<<endl;
gSystem->Load("StMcEvent");
cout <<"Loading StMcEventMaker"<<endl;
gSystem->Load("StMcEventMaker");
cout <<"Loading AssociationMaker"<<endl;
gSystem->Load("StAssociationMaker");
// create a new instance of the chain
chain = new StChain("StChain");
chain->SetDebug();
// add makers to the chain
StIOMaker* ioMaker = new StIOMaker("IO","r",MainFile,"bfcTree");
ioMaker->SetDebug();
ioMaker->SetIOMode("r");
ioMaker->SetBranch("*",0,"0"); //deactivate all branches
ioMaker->SetBranch("geantBranch",0,"r"); //activate geant Branch
ioMaker->SetBranch("dstBranch",0,"r"); //activate Event Branch
ioMaker->SetBranch("runcoBranch",0,"r"); //activate runco Branch
//Calibration Maker (StarDB,not a real Database!)
const char* calibDB = "MySQL:StarDb";
const char* paramsDB = "$STAR/StarDb";
St_db_Maker* calibMk = new St_db_Maker("StarDb",calibDB,paramsDB);
calibMk->SetDateTime("year_2b");
//calibMk->SetDateTime("year_1h");
calibMk->SetDebug();
//Read Tpc Database access
StTpcDbMaker *tpcDbMk = new StTpcDbMaker("tpcDb");
//StEventMaker
StEventMaker* eventReader = new StEventMaker("events","title");
eventReader->doPrintEventInfo = 0;
//StMcEventMaker
StMcEventMaker* mcEventReader = new StMcEventMaker();
//Association
StAssociationMaker* assocMaker = new StAssociationMaker();
cout <<"Calling Init() Methods "<<endl;
chain->Init();
cout <<"Starting Event Loop"<<endl;
int istat=0,iev=1;
EventLoop: if (iev<=nevents && !istat) {
chain->Clear();
cout << "---------------------- Processing Event : " << iev << endl;
istat = chain->Make(iev);
if (istat) {
cout << "Last Event Processed. Status = " << istat << endl;
}
iev++; goto EventLoop;
}
return;
}
void fPostStatus_db(
const char* opt = "",
const char* year = "17sim",
const char* input = "fPostStatus.txt",
int readdate = 0,
int readtime = 0)
{
// storeTime is beginning time for validity range in case of WRITING DB
TString option(opt), yr(year), storeTime;
int date, time; // time for READING DB
std::cout <<"year = " <<year <<std::endl;
if (yr.Contains("17sim"))
{
storeTime = "2016-12-10 00:00:00";
date = 20161210;
time = 0;
}
else if (yr.Contains("17ofl"))
{
storeTime = "2016-12-20 00:00:00";
date = 20161220;
time = 0;
}
else { std::cout << "Please specify valid year tag\n"; exit; }
if (readdate > 0) date = readdate;
if (readtime > 0) time = readtime;
std::cout << "Opt =" << opt << "\n";
std::cout << "write = " << option.Contains("writedb") << "\n";
std::cout << "storetime = " << storeTime << "\n";
std::cout << "date, time = " << date <<" "<< time << "\n";
gROOT->Macro("./loadlib.C");
//-------------------------------------------
const Int_t MAX_DB_INDEX = 241;
fpostStatus_st in[MAX_DB_INDEX];
FILE *FP = fopen(input, "r");
if (!FP) { printf("Could not open %s\n", input); exit; }
printf("\nReading %s\n", input);
char line[1000];
int n = 0;
while (fgets(line, 1000, FP) != NULL)
{
sscanf(line,"%hu %hu", &in[n].slatid, &in[n].status);
printf("slatId=%3d status=%1d\n", in[n].slatid, in[n].status);
n++;
}
printf("Found %d entries\n", n);
fclose(FP);
printf("File Closed\n", n);
//-------------------------------------------
#if 1
if (option.Contains("writedb"))
{
gSystem->Setenv("DB_ACCESS_MODE", "write");
StDbManager* mgr = StDbManager::Instance();
StDbConfigNode* node = mgr->initConfig("Calibrations_fps");
StDbTable* dbtable = node->addDbTable("fpostStatus");
mgr->setStoreTime(storeTime.Data());
dbtable->SetTable((char*)&in, MAX_DB_INDEX);
if (yr.Contains("sim")) dbtable->setFlavor("sim");
mgr->storeDbTable(dbtable);
std::cout << "INFO: table saved to database" << std::endl;
}
//-------------------------------------------
if (option.Contains("readdb"))
{
std::cout << "INFO: Reading database" << std::endl;
gSystem->Unsetenv("DB_ACCESS_MODE");
//gSystem->Unsetenv("DB_SERVER_LOCAL_CONFIG");
St_db_Maker* dbMk = new St_db_Maker("db", "MySQL:StarDb", "$STAR/StarDb");
dbMk->SetDebug();
dbMk->SetDateTime(date, time); // event or run start time, set to your liking
if (yr.Contains("ofl")) { dbMk->SetFlavor("ofl"); }
else if (yr.Contains("sim")) { dbMk->SetFlavor("sim"); }
dbMk->Init();
dbMk->Make();
TDataSet *DB = 0;
DB = dbMk->GetDataBase("Calibrations/fps/fpostStatus");
if (!DB) std::cout << "ERROR: no table found in db, or malformed local db config" << std::endl;
St_fpostStatus *dataset = 0;
dataset = (St_fpostStatus*) DB->Find("fpostStatus");
if (!dataset) { std::cout << "ERROR: dataset does not contain requested table" << std::endl; return; }
Int_t rows = dataset->GetNRows();
if (rows > 1) std::cout << "INFO: found INDEXED table with " << rows << " rows" << std::endl;
TDatime val[2];
dbMk->GetValidity((TTable*)dataset, val);
std::cout << "Dataset validity range: [ " << val[0].GetDate() << "." << val[0].GetTime() << " - "
<< val[1].GetDate() << "." << val[1].GetTime() << " ] " << std::endl;
fpostStatus_st *table = (fpostStatus_st*)dataset->GetTable();
for (Int_t i = 0; i < rows; i++)
{
std::cout << Form("Row=%3d slatid=%3d status=%1d\n", i, table[i].slatid, table[i].status);
}
}
#endif
return;
}//Main
void runTrigTreeMaker( const Char_t *file="rcf10063_2_6000evts.MuDst.root", const Char_t* fileOut = "rcf10063_2_6000evts.TrigTree-A.root", Char_t trigVer = 'g', Int_t nevents = 10000000,
Int_t flagMC = 1, const Char_t *eemcSetupPath="/star/u/sgliske/Share/StarTrigSimuSetup/" ){
//int flagMC=0; // 0/1 == Using Real/Simulation data files
int useEemc=2; // 0/1 == Exclude/Include EEMC in Trigger Decisions
int useBemc=0; // 0/1 == Exclude/Include BEMC in Trigger Decisions
int useL2=1; // 0/1 == Exclude/Include L2 in Trigger Decisions
int L2ConfigYear=2006; // possible: 2006, 2008
int bemcConfig=1; // Online==1, Offline==2, Expert==3
int playConfig=0; // jan:100_199
int emcEveDump=0; // extrating raw EMC data in a custom format
int outputL2Histo=1;//output L2 histos to directory outL2
TString outDir="./outL2/";
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
assert( !gSystem->Load("StDetectorDbMaker"));
assert( !gSystem->Load("StTpcDb"));
assert( !gSystem->Load("StDbUtilities"));
assert( !gSystem->Load("StDbBroker"));
assert( !gSystem->Load("St_db_Maker"));
assert( !gSystem->Load("StEEmcUtil")); // needed by eemcDb
assert( !gSystem->Load("StEEmcDbMaker"));
assert( !gSystem->Load("StDaqLib")); // needed by bemcDb
assert( !gSystem->Load("StEmcRawMaker"));
assert( !gSystem->Load("StEmcADCtoEMaker"));
if (flagMC) {
assert( !gSystem->Load("StMcEvent"));
assert( !gSystem->Load("StMcEventMaker"));
assert( !gSystem->Load("StEmcSimulatorMaker"));
assert( !gSystem->Load("StEEmcSimulatorMaker"));
assert( !gSystem->Load("StEpcMaker"));
}
assert( !gSystem->Load("StTriggerUtilities"));
assert( !gSystem->Load("StSpinDbMaker") );
assert( !gSystem->Load("StEEmcPoolEEmcTreeContainers") );
assert( !gSystem->Load("StEEmcTreeMaker") );
gROOT->Macro("LoadLogger.C");
cout << " loading done " << endl;
chain= new StChain("StChain");
if (flagMC){
TString geantFile;
geantFile += file;
geantFile.ReplaceAll("MuDst.root", "geant.root");
printf("geantFile=%s\n", geantFile.Data());
StIOMaker* ioMaker = new StIOMaker();
ioMaker->SetFile(geantFile);
ioMaker->SetIOMode("r");
ioMaker->SetBranch("*",0,"0"); //deactivate all branches
ioMaker->SetBranch("geantBranch",0,"r"); //activate geant Branch
StMcEventMaker *evtMaker = new StMcEventMaker();
}
//Need MuDstMaker to get data
printf(" Analyzing file=%s\n",file);
StMuDstMaker* muDstMaker =new StMuDstMaker(0,0,"",file,"",1000);
// just to make sure StEvent is found for the trig simu
StMuDst2StEventMaker *muDst2StEvent = new StMuDst2StEventMaker();
//Database -- get a real calibration from the database
St_db_Maker* dbMk =0;
if(useEemc || useL2) // full DB access
dbMk = new St_db_Maker("StarDb","MySQL:StarDb","MySQL:StarDb","$STAR/StarDb");
else // only Barrel is uploaded, is faster
dbMk = new St_db_Maker("Calibrations","MySQL:Calibrations_emc");
dbMk->SetAttr("blacklist", "fgt");
dbMk->SetAttr("blacklist", "svt");
dbMk->SetAttr("blacklist", "tpc");
dbMk->SetAttr("blacklist", "ftpc");
// CORRECTION: do not always need to set date and time for MC
//If MC then must set database time and date
//If Endcap fast simu is used tower gains in DB do not matter,JB
//if (flagMC) dbMk->SetDateTime(20060522, 55000);//timestamp R7142018
//Collect all output histograms
TObjArray* HList=new TObjArray;
//Endcap DB
if(useEemc || useL2) new StEEmcDbMaker("eemcDb");
//Get BEMC adc values
if (flagMC && useBemc) {
StEmcSimulatorMaker* emcSim = new StEmcSimulatorMaker(); //use this instead to "redo" converstion from geant->adc
if (bemcConfig == 1) {
emcSim->setCheckStatus(kBarrelEmcTowerId,false); //this returns hits regardless of offline tower status
}
emcSim->setCalibSpread(kBarrelEmcTowerId,0.15);//spread gains by 15%
}
if (flagMC==0 && useBemc){
StEmcADCtoEMaker *bemcAdc = new StEmcADCtoEMaker();//for real data this sets calibration and status
if (bemcConfig == 1) {
bemcAdc->setCheckStatus(kBarrelEmcTowerId,false);
}
}
// original comment: must use slow simulator to get pedestals
// correct for L2. New comment: I think this is needed if slow simulator was
// used in the MC production.
if (0 && flagMC==1 && useEemc){
StEEmcSlowMaker *slowSim = new StEEmcSlowMaker("slowSim");
// note: changing the sampling fraction does not effect the
// L0 trigger decision. Also, it is probably better to use the value
// in the slow simulator, so that it matches data. (S. Gliske,
// Sept 2012)
// This 1.3 comes from the iron/air MC, where the sampling
// fraction was believed to be near 4%, but slow/fast simulator
// used 5%.
// Now that the slow/fast simulator uses 0.048, it is a 1.25 scaling factor
// slowSim->setSamplingFraction(0.0384); // effectively scale all Tower energies with a factor of 1.3 (added by: Ilya Selyuzhenkov; April 11, 2008)
slowSim->setAddPed(true);
slowSim->setSmearPed(true);
}
//Get TriggerMaker
StTriggerSimuMaker *simuTrig = new StTriggerSimuMaker("StarTrigSimu");
simuTrig->setHList(HList);
simuTrig->setMC(flagMC); // must be before individual detectors, to be passed
simuTrig->useBbc();
if(useEemc) {
simuTrig->useEemc(0);//default=0:just process ADC, 1,2:comp w/trgData,see .
simuTrig->eemc->setSetupPath(eemcSetupPath);
}
if(useBemc){
simuTrig->useBemc();
simuTrig->bemc->setConfig(bemcConfig);
}
if(flagMC && useEemc==2){
// pass one argument to M-C as generic switch
// Endcap specific params -- ok Jan you need to change this to a default "online" setup
int eemcDsmSetup[20]; // see StEemcTriggerSimu::initRun() for definition
memset(eemcDsmSetup, 0,sizeof(eemcDsmSetup));// clear all, may be a bad default
if( trigVer == 'd' ){
// valid runs 7130037-7132029
eemcDsmSetup[0]=6; // HTthr0
eemcDsmSetup[1]=12; // HTthr1
eemcDsmSetup[2]=22; // HTthr2
eemcDsmSetup[3]=1; // TPthr0
eemcDsmSetup[4]=20; // TPthr1
eemcDsmSetup[5]=31; // TPthr2
} else if ( trigVer == 'e' ){
// valid runs 7132045-7133051
eemcDsmSetup[0]=6; // HTthr0
eemcDsmSetup[1]=17; // HTthr1
eemcDsmSetup[2]=22; // HTthr2
eemcDsmSetup[3]=1; // TPthr0
eemcDsmSetup[4]=20; // TPthr1
eemcDsmSetup[5]=31; // TPthr2
} else if ( trigVer == 'f' ){
// valid runs 7133052-7156040
eemcDsmSetup[0]=6; // HTthr0
eemcDsmSetup[1]=16; // HTthr1
eemcDsmSetup[2]=22; // HTthr2
eemcDsmSetup[3]=1; // TPthr0
eemcDsmSetup[4]=20; // TPthr1
eemcDsmSetup[5]=31; // TPthr2
} else if ( trigVer == 'g' ){
// test trigger, ht1 & tp2 are 1 DSM (roughly 5%) higher than 'e'
eemcDsmSetup[0]=6; // HTthr0
eemcDsmSetup[1]=18; // HTthr1
eemcDsmSetup[2]=22; // HTthr2
eemcDsmSetup[3]=1; // TPthr0
eemcDsmSetup[4]=21; // TPthr1
eemcDsmSetup[5]=31; // TPthr2
} else if ( trigVer == 'h' || trigVer == 'i' ){
// test trigger, ht1 & tp2 are 2 DSM (roughly 10%) higher than 'e'
eemcDsmSetup[0]=6; // HTthr0
eemcDsmSetup[1]=19; // HTthr1
eemcDsmSetup[2]=22; // HTthr2
eemcDsmSetup[3]=1; // TPthr0
eemcDsmSetup[4]=22; // TPthr1
eemcDsmSetup[5]=31; // TPthr2
} else {
cerr << "Invalid trigger version" << endl;
return;
};
eemcDsmSetup[10]=2; //HTTPthrSelc, 2=use_thres_#1
simuTrig->eemc->setDsmSetup(eemcDsmSetup);
}
StGenericL2Emulator* simL2Mk=0;
if(useL2) {
/*
NOTE: the following comment was in the rdMu2TrigSimu.C macro,
but does not seem applicable. Just use the given
StarTrigSimuSetup directory, and make sure the given outDir
exists, all will be well.
reads all input/setup files from L2setup-yyyymmdd/
writes all output files to L2out-yyyymmdd
depending on the DB time stamp
both dierectiorie MUST exist, setup must be reasonable
*/
if(L2ConfigYear==2006) simL2Mk= new StL2_2006EmulatorMaker;
else if(L2ConfigYear==2008) simL2Mk= new StL2_2008EmulatorMaker;
assert(simL2Mk);
simL2Mk->setSetupPath(eemcSetupPath);
simL2Mk->setOutPath(outDir.Data());
if (flagMC) simL2Mk->setMC();
simuTrig->useL2(simL2Mk);
}
// now the maker to write the tree
StTrigTreeMaker_t *treeMkr = new StTrigTreeMaker_t( "trigTreeMaker", fileOut, 1 );
treeMkr->addTrigger( -999 ); // i.e. no trigger cut now
chain->ls(3);
cout << "******* Initialize *******" << endl ;
chain->Init();
cout << "******* done with initializing *******" << endl;
for (Int_t iev=0;iev<nevents; iev++) {
if( iev % 100 == 1 )
cout << "Working on eventNumber:\t" << iev <<"\tof:\t"<<nevents<<endl;
chain->Clear();
int iret = chain->Make(iev);
total++;
if (iret % 10 == kStEOF || iret % 10 == kStFatal) {
cout << iev << " Bad return code!" << endl;
break;
}
//struct L2gammaResult2006* result = dynamic_cast< L2gammaResult2006* >( simL2Mk->result() );
//cout << "ZZZ " << result->trigger << endl;
if( simL2Mk ){
// const unsigned int* result = simL2Mk->result();
// cout << "ZZZ " << std::hex << result[0] << result[1] << result[2] << result[3] << result[4] << result[5] << result[6] << result[7] << std::dec
// << " " << simL2Mk->mAcceptTriggerList.size() << ' ' << simL2Mk->mVetoTriggerList.size() << endl;
// cout << "ZZZ accept: ";
// for( Int_t i=0; i<simL2Mk->mAcceptTriggerList.size(); ++i )
// cout << simL2Mk->mAcceptTriggerList.at(i) << ' ';
// cout << endl;
// cout << "ZZZ veto: ";
// for( Int_t i=0; i<simL2Mk->mVetoTriggerList.size(); ++i )
// cout << simL2Mk->mVetoTriggerList.at(i) << ' ';
// cout << endl;
// cout << "ZZZ " << simuTrig->isTrigger( 137641 ) << endl;
}
}
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << total << endl;
cout << "****************************************** " << endl;
if (outputL2Histo==1) {
TString fileMu=file;
printf("=%s=\n",fileMu.Data());
if(fileMu.Contains(".lis")) fileMu.ReplaceAll(".lis",".trgSim");
if(fileMu.Contains(".MuDst.root")) fileMu.ReplaceAll(".MuDst.root",".trgSim");
TString outF=outDir+fileMu;
outF+=".hist.root";
printf("=%s=\n",outF.Data());
hf=new TFile(outF,"recreate");
if(hf->IsOpen()) {
cout << "HList not saved, or seg faults" << endl;
//HList->Write();
//printf("\n Histo saved -->%s<\n",outF.Data());
}
else {
printf("\n Failed to open Histo-file -->%s<, continue\n",outF.Data());
}
}
}
void Run2006SimuTriggerMaker(const char *dir ="",
const char* file="/star/data32/reco/pp200/pythia6_205/above_35gev/cdf_a/y2004y/gheisha_on/p05ih/rcf1230_10_4000evts.MuDst.root",
const char *fname="/star/data32/reco/pp200/pythia6_205/above_35gev/cdf_a/y2004y/gheisha_on/p05ih/rcf1230_10_4000evts.event.root",
const char *filter = "")
{
int nevents = 10;
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
gSystem->Load("StTpcDb");
gSystem->Load("StDetectorDbMaker");
gSystem->Load("StDbUtilities");
gSystem->Load("StMcEvent");
gSystem->Load("StMcEventMaker");
gSystem->Load("StDaqLib");
gSystem->Load("StEmcRawMaker");
gSystem->Load("StEmcADCtoEMaker");
gSystem->Load("StPreEclMaker");
gSystem->Load("StEpcMaker");
gSystem->Load("StEmcSimulatorMaker");
gSystem->Load("StEmcUtil");
gSystem->Load("StDbBroker");
gSystem->Load("St_db_Maker");
gSystem->Load("StEEmcUtil");
gSystem->Load("StEEmcDbMaker");
gSystem->Load("StEmcTriggerMaker");
double pi = atan(1.0)*4.0;
cout << " loading done " << endl;
chain= new StChain("StChain");
chain->SetDebug(1);
gMessMgr->SwitchOff("D");
gMessMgr->SwitchOff("I");
StIOMaker* ioMaker = new StIOMaker();
ioMaker->SetFile(fname);
ioMaker->SetIOMode("r");
ioMaker->SetBranch("*",0,"0"); //deactivate all branches
ioMaker->SetBranch("geantBranch",0,"r"); //activate geant Branch
class StMcEventMaker *mcEventMaker = new StMcEventMaker();
mcEventMaker->doPrintEventInfo = false;
mcEventMaker->doPrintMemoryInfo = false;
StMuDebug::setLevel(1);
StMuDstMaker* muDstMaker = new StMuDstMaker(0,0,dir,file,filter,1e6,"MuDst");
//Database -- get a real calibration from the database
St_db_Maker* dbMk = new St_db_Maker("StarDb","MySQL:StarDb","MySQL:StarDb","$STAR/StarDb");
dbMk->SetDateTime(20060506,214129 ); // for simulation
//to use ideal ideal gains/status/calib uncommment correct line below
//dbMk->SetFlavor("sim","bemcPed");
//dbMk->SetFlavor("sim","bemcStatus");
//dbMk->SetFlavor("sim","bemcCalib");
//dbMk->SetFlavor("sim","bemcGain");
//dbMk->SetFlavor("sim","eemcPMTcal");
//dbMk->SetFlavor("sim","eemcPIXcal");
chain->ls(3);
//Database interface
StDetectorDbMaker* detDbMk = new StDetectorDbMaker();
//Endcap DB
StEEmcDbMaker* eemcb = new StEEmcDbMaker("eemcDb");
//StEmcADCtoEMaker *adc = new StEmcADCtoEMaker(); // this will just convert what's in MuDst to ADC, use for data only!
StEmcSimulatorMaker* emcSim = new StEmcSimulatorMaker(); //use this instead to "redo" converstion from geant->adc
emcSim->setCalibSpread(kBarrelEmcTowerId, 0.15);
StPreEclMaker* preEcl = new StPreEclMaker(); //need this to fill new StEvent information
//StEmcTrigger
StEmcTriggerMaker *emcTrig = new StEmcTriggerMaker("bemctrigger");
chain->Init();
chain->PrintInfo();
chain->ls(3);
TChain* fileChain = muDstMaker->chain();
int BL1_ADC[6];
int hold=-1;
for (Int_t iev=0;iev<nevents; iev++) {
cout << "****************************************** " << endl;
cout << "Working on eventNumber:\t" << iev <<"\tof:\t"<<nevents<<endl;
cout << "*************************1***************** " << endl;
chain->Clear();
int iret = chain->Make(iev);
total++;
if (iret) {
cout << "Bad return code!" << endl;
break;
}
cout<<" "<<endl;
map<int,int>::iterator iter;
map<int,int> tower;
map<int,int> tpatch;
map<int,int> jpatch;
cout<<endl;
cout<<endl;
cout<<" 1101 ht1 = "<<emcTrig->isTrigger(1101)<<" = "<<emcTrig->is2003HT1()<<endl;
tower=emcTrig->barrelTowersAboveThreshold(1101);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
cout<<"old tower id="<<emcTrig->get2003HT1_ID()<<" adc="<<emcTrig->get2003HT1_ADC()<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 2101 ht1 = "<<emcTrig->isTrigger(2201)<<" = "<<emcTrig->is2003HT1()<<endl;
tower=emcTrig->barrelTowersAboveThreshold(2201);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
cout<<"old tower id="<<emcTrig->get2003HT1_ID()<<" adc="<<emcTrig->get2003HT1_ADC()<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 2202 ht2 = "<<emcTrig->isTrigger(2202)<<" = "<<emcTrig->is2003HT2()<<endl;
tower=emcTrig->barrelTowersAboveThreshold(2202);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
cout<<"old tower id="<<emcTrig->get2003HT2_ID()<<" adc="<<emcTrig->get2003HT2_ADC()<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 45201 ht1 = "<<emcTrig->isTrigger(45201)<<" = "<<emcTrig->is2004HT1()<<endl;
tower=emcTrig->barrelTowersAboveThreshold(45201);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
cout<<"old tower id="<<emcTrig->get2004HT1_ID()<<" adc="<<emcTrig->get2004HT1_ADC()<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 45202 ht1 = "<<emcTrig->isTrigger(45202)<<" = "<<emcTrig->is2004HT2()<<endl;
tower=emcTrig->barrelTowersAboveThreshold(45202);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
cout<<"old tower id="<<emcTrig->get2004HT2_ID()<<" adc="<<emcTrig->get2004HT2_ADC()<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 45206 jp1 = "<<emcTrig->isTrigger(45206)<<" = "<<emcTrig->is2004JP1()<<endl;
jpatch=emcTrig->barrelJetPatchesAboveThreshold(45206);
cout<<"Total #'s of jpatches="<<jpatch.size()<<endl;
for ( iter=jpatch.begin();iter !=jpatch.end(); iter++){
cout<<"jpatch id="<<iter->first<<" adc="<<iter->second<<endl;
cout<<"old patch id="<<emcTrig->get2004JP1_ID()<<" adc="<<emcTrig->get2004JP1_ADC()<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 45207 jp2 = "<<emcTrig->isTrigger(45207)<<" = "<<emcTrig->is2004JP2()<<endl;
jpatch=emcTrig->barrelJetPatchesAboveThreshold(45207);
cout<<"Total #'s of jpatches="<<jpatch.size()<<endl;
for ( iter=jpatch.begin();iter !=jpatch.end(); iter++){
cout<<"jpatch id="<<iter->first<<" adc="<<iter->second<<endl;
cout<<"old patch id="<<emcTrig->get2004JP2_ID()<<" adc="<<emcTrig->get2004JP2_ADC()<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 127212 ht2 (matrix0)="<<emcTrig->isTrigger(127212)<<endl;
tower=emcTrig->barrelTowersAboveThreshold(127212);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 127213 ht2 (matrix1)="<<emcTrig->isTrigger(127213)<<endl;
tower=emcTrig->barrelTowersAboveThreshold(127213);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 137213 ht2 (matrix3) ="<<emcTrig->isTrigger(127213)<<endl;
tower=emcTrig->barrelTowersAboveThreshold(127213);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 127501 jp0 (matrix1)="<<emcTrig->isTrigger(127501)<<endl;
jpatch=emcTrig->barrelJetPatchesAboveThreshold(127501);
cout<<"Total #'s of jpatches="<<jpatch.size()<<endl;
for ( iter=jpatch.begin();iter !=jpatch.end(); iter++){
cout<<"jpatch id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 137501 jp0 (matrix3) ="<<emcTrig->isTrigger(137501)<<endl;
jpatch=emcTrig->barrelJetPatchesAboveThreshold(137501);
cout<<"Total #'s of jpatches="<<jpatch.size()<<endl;
for ( iter=jpatch.begin();iter !=jpatch.end(); iter++){
cout<<"jpatch id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 127221 jp1 (matrix1)="<<emcTrig->isTrigger(127221)<<endl;
jpatch=emcTrig->barrelJetPatchesAboveThreshold(127221);
cout<<"Total #'s of jpatches="<<jpatch.size()<<endl;
for ( iter=jpatch.begin();iter !=jpatch.end(); iter++){
cout<<"jpatch id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 137221 jp1 (matrix2)="<<emcTrig->isTrigger(137221)<<endl;
jpatch=emcTrig->barrelJetPatchesAboveThreshold(137221);
cout<<"Total #'s of jpatches="<<jpatch.size()<<endl;
for ( iter=jpatch.begin();iter !=jpatch.end(); iter++){
cout<<"jpatch id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 137222 jp1 (matrix3)="<<emcTrig->isTrigger(137222)<<endl;
jpatch=emcTrig->barrelJetPatchesAboveThreshold(127501);
cout<<"Total #'s of jpatches="<<jpatch.size()<<endl;
for ( iter=jpatch.begin();iter !=jpatch.end(); iter++){
cout<<"jpatch id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 127821 http (matrix1)="<<emcTrig->isTrigger(127821)<<endl;
tpatch=emcTrig->barrelTriggerPatchesAboveThreshold(127821);
cout<<"Total #'s of tpatches="<<tpatch.size()<<endl;
for ( iter=tpatch.begin();iter !=tpatch.end(); iter++){
cout<<"tpatch id="<<iter->first<<" adc="<<iter->second<<endl;
}
tower=emcTrig->barrelTowersAboveThreshold(127821);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
int triggerPatch = emcTrig->barrelTriggerPatchForTower(iter->first);
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<" tp (id,adc)=("<< triggerPatch << "," << tpatch[triggerPatch] << ")" <<endl;
}
cout<<endl;
cout<<endl;
cout<<" 137821 http (matrix3) ="<<emcTrig->isTrigger(137821)<<endl;
tpatch=emcTrig->barrelTriggerPatchesAboveThreshold(137821);
cout<<"Total #'s of tpatches="<<tpatch.size()<<endl;
for ( iter=tpatch.begin();iter !=tpatch.end(); iter++){
cout<<"tpatch id="<<iter->first<<" adc="<<iter->second<<endl;
}
tower=emcTrig->barrelTowersAboveThreshold(137821);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
int triggerPatch = emcTrig->barrelTriggerPatchForTower(iter->first);
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<" tp (id,adc)=("<< triggerPatch << "," << tpatch[triggerPatch] << ")" <<endl;
}
cout<<endl;
cout<<endl;
cout<<" 137822 http (matrix4)="<<emcTrig->isTrigger(137822)<<endl;
tpatch=emcTrig->barrelTriggerPatchesAboveThreshold(137822);
cout<<"Total #'s of tpatches="<<tpatch.size()<<endl;
for ( iter=tpatch.begin();iter !=tpatch.end(); iter++){
cout<<"tpatch id="<<iter->first<<" adc="<<iter->second<<endl;
}
tower=emcTrig->barrelTowersAboveThreshold(137822);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
int triggerPatch = emcTrig->barrelTriggerPatchForTower(iter->first);
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<" tp (id,adc)=("<< triggerPatch << "," << tpatch[triggerPatch] << ")" <<endl;
}
cout<<endl;
cout<<endl;
cout<<" 117705 jpsi (matrix1)="<<emcTrig->isTrigger(117705)<<endl;
tower=emcTrig->barrelTowersAboveThreshold(117705);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 117602 upsilon (matrix1)="<<emcTrig->isTrigger(117602)<<endl;
tower=emcTrig->barrelTowersAboveThreshold(117602);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 137602 upsilon (matrix3)="<<emcTrig->isTrigger(137602)<<endl;
tower=emcTrig->barrelTowersAboveThreshold(137602);
cout<<"Total #'s of towers="<<tower.size()<<endl;
for ( iter=tower.begin();iter !=tower.end(); iter++){
cout<<"tower id="<<iter->first<<" adc="<<iter->second<<endl;
}
cout<<endl;
cout<<endl;
cout<<" 127611 http-L2="<<emcTrig->isTrigger(127611)<<endl;
cout<<endl;
cout<<endl;
cout<<" 127622 jp0-etot-L2="<<emcTrig->isTrigger(127622)<<endl;
cout<<endl;
cout<<endl;
cout<<" 137622 jp0-etot-L2="<<emcTrig->isTrigger(137622)<<endl;
for (int j=0;j<6;j++){
BL1_ADC[j]=0.0;
emcTrig->get2006BL1_ADC(j,&hold);
BL1_ADC[j]=hold;
cout<<" 2x1 Patch "<<j<<" = "<<BL1_ADC[j]<<endl;
}
}
chain->Finish();
cout << "****************************************** " << endl;
cout << "total number of events " << total << endl;
cout << "****************************************** " << endl;
}
int rdEztPanitkin( int trigID=0,//96294,
int nEve=100,
Int_t nFiles =20,
char* file="runList/R7101015.lis", //R6121034.lis,
char* inDir = "./"
){
//inDir="/star/data05/scratch/eemcdb/muDst/2005/048/";
// file="st_physics_6048024_raw_2010006.MuDst.root";
inDir=" myMuDst/";
file="st_physics_7166007_raw_1040002.MuDst.root";
gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
loadSharedLibraries();
cout << " loading done " << endl;
assert( !gSystem->Load("StEEmcPoolmuEztPanitkin"));
assert( !gSystem->Load("StEEmcUtil"));
// libraries below are needed for DB interface
assert( !gSystem->Load("StDbBroker"));
assert( !gSystem->Load("St_db_Maker"));
assert( !gSystem->Load("StEEmcDbMaker"));
// create chain
chain = new StChain("StChain");
printf("adding muDst from '%s' ....\n",file);
// Now we add Makers to the chain...
muMk = new StMuDstMaker(0,0,inDir,file,"MuDst.root",nFiles);
TChain* tree=muMk->chain(); assert(tree);
int nEntries=(int) tree->GetEntries();
printf("total eve in chain =%d\n",nEntries);
printf("in=%s%s=\n",inDir,file);
//return;
St_db_Maker *stDb = new St_db_Maker("StarDb", "MySQL:StarDb");
stDb->SetFlavor("onlPed","eemcPMTped");
new StEEmcDbMaker("EEmcDbMaker");
HList=new TObjArray;
myMk3=new MuEzPanitkinMaker("myPanitkin","MuDst");
myMk3->SetHList(HList);
myMk3->SetTrigIdFilter(trigID);
//myMk3->SetTrigIdFilter(66300); //zeroB
//if(trigID==66007)
// myMk3->SetHistoPixels();
// myMk3->SetSpy();
{
// Endcap specific DSM params
int eemcDsmSetup[20]; // see StEemcTriggerSimu::initRun() for definition
memset(eemcDsmSetup, 0,sizeof(eemcDsmSetup));// clear all, may be a bad default
eemcDsmSetup[0]=3; // HTthr0
eemcDsmSetup[1]=12; // HTthr1
eemcDsmSetup[2]=22; // HTthr2
eemcDsmSetup[3]=1; // TPthr0
eemcDsmSetup[4]=17; // TPthr1
eemcDsmSetup[5]=31; // TPthr2
eemcDsmSetup[10]=2; //HTTPthrSelc, 2=use_thres_#1
// myMk3->SetDsmSetup(eemcDsmSetup);
}
gMessMgr->SwitchOff("D");
gMessMgr->SwitchOn("I");
muMk->SetStatus("*",0);
muMk->SetStatus("EztAll",1);
muMk->SetStatus("MuEvent",1);
chain->Init();
chain->ls(3);
// muMk->printArrays();
printf("All Ezt-branches set\n");
int eventCounter=0;
int stat=0;
int t1=time(0);
StMuTimer timer;
timer.start();
//---------------------------------------------------
while ( 1) {// loop over events
if(eventCounter>=nEve) break;
eventCounter++;
chain->Clear();
stat = chain->Make();
if(stat) break;
if(eventCounter%200!=0)continue;
printf("\n\n ====================%d processing ==============\n", eventCounter);
}
printf("sorting done, nEve=%d of %d\n",nEve, nEntries);
int t2=time(0);
if(t1==t2) t2++;
float rate=1.*eventCounter/(t2-t1);
float nMnts=(t2-t1)/60.;
printf("sorting done %d of nEve=%d, elapsed rate=%.1f Hz, tot %.1f minutes\n",eventCounter,nEntries,rate,nMnts);
if (eventCounter) {
cout << "CPU time/event= " << timer.elapsedTime()/eventCounter << " sec "
<< " rate= " << eventCounter/timer.elapsedTime() << " Hz" << endl;
}
myMk3->saveHisto("eemcQA");
// TString txt="R6035106-"; txt+=trigID; myMk3->saveHisto(txt);
}
int fgtSingleEvents( const Char_t *filenameIn = "testfile.daq",
const Char_t *filenameOut = "testfile.root",
Int_t nevents = 200,
Bool_t isCosmic = 0,
Bool_t cutShortEvents = 0 ){
LoadLibs();
Int_t ierr = 0;
if( isCosmic )
cout << "Is Cosmic" << endl;
else
cout << "Is not cosmic" << endl;
//
// START CONSTRUCTING THE CHAIN
//
cout << "Constructing the chain" << endl;
analysisChain = new StChain("eemcAnalysisChain");
std::string fgtDbMkrName = "";
if( !isCosmic ){
// always cut short events if it is cosmic data
cutShortEvents = 1;
cout << "Loading St_db_Maker" << endl;
gSystem->Load("libStDb_Tables.so");
gSystem->Load("StDbLib.so");
gSystem->Load("St_db_Maker");
gSystem->Load("StDbBroker");
TString dir0 = "MySQL:StarDb";
TString dir1 = "$STAR/StarDb";
St_db_Maker *dbMkr = new St_db_Maker( "dbMkr", dir0, dir1 );
dbMkr->SetDateTime(20120115,1);
cout << "Loading StFgtDbMaker" << endl;
gSystem->Load("StFgtDbMaker");
cout << "Constructing StFgtDbMaker" << endl;
fgtDbMkr = new StFgtDbMaker( "fgtDbMkr" );
//fgtDbMkr->SetFlavor("ideal",""); // mapping is wrong, but at least the code runs...
fgtDbMkrName = fgtDbMkr->GetName();
};
//
// NOW THE OTHER READERS AND MAKERS
//
cout << "Constructing the daq reader" << endl;
daqRdr = new StFgtRawDaqReader( "daqReader", filenameIn, fgtDbMkrName.data() );
daqRdr->setIsCosmic( isCosmic );
daqRdr->cutShortEvents( cutShortEvents );
cout << "Constructing the A2C converter" << endl;
a2cMkr = new StFgtA2CMaker( "a2cMaker" );
//a2cMkr->setTimeBinMask( 0xFF );
a2cMkr->setAbsThres( 100 ); // set to below -4096 to skip cut
a2cMkr->setRelThres( 0 ); // set to zero to skip cut
if( !isCosmic )
a2cMkr->setFgtDb( fgtDbMkr->getDbTables() );
//a2cMkr->doRemoveOtherTimeBins( 0 );
//a2cMkr->doCutBadStatus( 0 );
cout << "Constructing the QA Maker" << endl;
qaMkr = new StFgtSingleEventQA( "fgtSingleEventQA" );
qaMkr->setFilename( filenameOut );
// debug
// analysisChain->ls(4);
cout << "Initializing" << endl;
ierr = analysisChain->Init();
if( ierr ){
cout << "Error initializing" << endl;
return;
};
if( nevents < 0 )
nevents = 1<<30; // a big number
cout << "max nevents = " << nevents << endl;
for( int i=0; i<nevents && !ierr; ++i ){
if( i+1 % 100 == 0 )
cout << "\ton event number " << i << endl;
//cout << "clear" << endl;
analysisChain->Clear();
//cout << "make" << endl;
ierr = analysisChain->Make();
};
//
// Calls the ::Finish() method on all makers
//
cout << "finish" << endl;
analysisChain->Finish();
cerr << "\tall done" << endl;
return;
};
void lasertest(TString filename,int ftpc, int lsec, int straight, int gfit,char* t0, char* gas,float gastemp,float mbfield)
{
Int_t minz, maxz;
Int_t minrad = 0;
Int_t maxrad = 30;
cout<<"Starting lasertest.C:"<<endl;
cout<<" filename = "<<filename<<".root"<<endl;
cout<<" ftpc = "<<ftpc;
if ( ftpc == 1 ) cout<<" FTPC West"<<endl;
if ( ftpc == 2 ) cout<<" FTPC East"<<endl;
cout<<" lsec = "<<lsec<<endl;
cout<<" straight = "<<straight<<endl;
cout<<" gfit = "<<gfit<<endl;
// for FTPC West
if ( ftpc == 1 ) {
minz = 0;
maxz = 300;
}
// for FTPC East
if ( ftpc == 2 ) {
minz = -300;
maxz = 0;
}
cout<<" minz = "<<minz<<endl;
cout<<" maxz = "<<maxz<<endl;
cout<<" minrad = "<<minrad<<endl;
cout<<" maxrad = "<<maxrad<<endl;
cout<<" t0 = "<<t0<<endl;
cout<<" gas = "<<gas<<endl;
cout<<" gastemp = "<<gastemp<<endl;
cout<<" mbfield = "<<mbfield<<endl;
cout<<endl;
if (gClassTable->GetID("TTable") < 0) gSystem->Load("libStar");
gSystem->Load("libMinuit.so");
gSystem->Load("libSpectrum.so");
gSystem->Load("St_base");
gSystem->Load("StChain");
gSystem->Load("libtpc_Tables");
gSystem->Load("StUtilities");
gSystem->Load("StarClassLibrary");
gSystem->Load("StEvent");
gSystem->Load("StarMagField");
gSystem->Load("StMagF");
gSystem->Load("libStDb_Tables.so");
gSystem->Load("StDetectorDbMaker.so");
gSystem->Load("StTpcDb");
gSystem->Load("StDbUtilities");
gSystem->Load("StDbLib.so");
gSystem->Load("StDbBroker.so");
gSystem->Load("libStDb_Tables.so");
gSystem->Load("St_db_Maker.so");
gSystem->Load("StFtpcCalibMaker");
gSystem->Load("libftpc_Tables");
gSystem->Load("StFtpcClusterMaker");
gSystem->Load("StFtpcTrackMaker");
// Create the makers to be called by the current chain
const char *mysqlDB = "MySQL:StarDb";
const char *paramsDB = "$STAR/StarDb";
//const char *paramsDB = "$PWD/StarDb";
StChain *chain = new StChain();
StFtpcCalibMaker *laser=new StFtpcCalibMaker();
laser->GetRunInfo(filename);
cout<<" date = "<<laser->Date()<<" time = "<<laser->Time()<<endl;
St_db_Maker *dbMk = new St_db_Maker("db",mysqlDB,paramsDB);
dbMk->SetDateTime(laser->Date(),laser->Time());
dbMk->Init();
dbMk->Make();
cout<<"dbDate = "<<dbMk->GetDateTime().GetDate()<<endl;
cout<<"After Database init !!!"<<endl;
cout<<endl;
if (laser->DbInit(mbfield) == kStWarn) {
delete laser;
break;
}
// Interation over tzero
for (int i=-5;i<7;i++)
{
float step=i/10.0;
// Interation over gas composition
for (int k=-5;k<8;k++)
{
float step2=k/10.0;
char t[3];char g[3];
sprintf(t,"%.2f",step);
sprintf(g,"%.2f",step2);
cout<<endl;
cout<<"laser->DoLaserCalib: deltat0 = "<<t<<" und deltaGas = "<<g<<endl;
if (k==0 && i==0)
cout<<"Kommt zum Schluss !!!"<<endl;
else
laser->DoLaserCalib(filename,ftpc,lsec,straight,gfit,minz,maxz,minrad,maxrad,t,g,gastemp,mbfield);
}
}
laser->DoLaserCalib(filename,ftpc,lsec,straight,gfit,minz,maxz,minrad,maxrad,"0","0",gastemp,mbfield);
delete laser;
}
