rdMu2spinTest.C

class StChain;
class StMuEmcCollection;

class StEEmcDbMaker;
class StMuDstMaker;
class TChain;

StEEmcDbMaker  *myDb;
StMuDstMaker* muMk;
StChain *chain=0;


int rdMu2spinTest( char* Rrun    ="R6172038",
                   Int_t nFiles  = 20,
                   char* inDir   = "./lis/",
                   int nEve=500)
{

    #define  USE_DB  
  char *outPath="iterX/";


  gROOT->LoadMacro("$STAR/StRoot/StMuDSTMaker/COMMON/macros/loadSharedLibraries.C");
  loadSharedLibraries();
  cout << " loading done " << endl;

  TString fileL=Rrun; fileL+=".lis";
  int runNo=atoi(Rrun+1);

  // tmp : 2006
   inDir   = "/star/data05/scratch/balewski/2006-prod-muDst/7097090/";
   fileL="st_physics_7097090_raw_1010001.MuDst.root";

  // tmp : 2005
  // inDir   = "/star/data05/scratch/balewski/2005-prod-muDst/6174052/";
  // fileL="st_physics_6174068_raw_2030018.MuDst.root";
 

  // runNo=999; Rrun="R999";

#ifdef USE_DB
  gSystem->Load("StDbLib");
  gSystem->Load("StDbBroker");
  gSystem->Load("St_db_Maker");
  gSystem->Load("StSpinDbMaker");
#endif
    
  // create chain    
  chain = new StChain("StChain"); 
  // inDir   ="/star/u/jwebb/work/2005/TEST/";
  // fileL="test.lis";

  printf("adding muDst from '%s' ....\n",fileL.Data());
  // Now we add Makers to the chain...   
  muMk = new StMuDstMaker(0,0,inDir,fileL,"MuDst.root",nFiles);
  TChain* tree=muMk->chain(); assert(tree); 
  int nEntries=(int)tree->GetEntries();
  printf("total eve in chain =%d\n",nEntries);

#ifdef USE_DB
  St_db_Maker *stDb = new St_db_Maker("StarDb", "MySQL:StarDb");
  spDb=new StSpinDbMaker("spinDb");
  
  gMessMgr->SwitchOff("D");
  gMessMgr->SwitchOn("I");
  hbxI=new TH1F("bXI","Intended fill pattern vs. STAR bXing; bXing at STAR IP",120,-0.5,119.5);
#endif

  hbx48=new TH1F("bX48","Rate vs. raw bx48; bXing= raw bx48",120,-0.5,119.5);
  hbx7=new TH1F("bX7","Rate vs. raw bx7; bXing= raw bx7",120,-0.5,119.5);
  hbx748=new TH2F("bX748","BX ID correlation ;  raw bx7;  raw bx48",120,-0.5,119.5,120,-0.5,119.5);

  hbx48c=new TH1F("bX48c","Rate vs. yellBX using bx48 ; bXing= bx48+offset",120,-0.5,119.5);
  hbx7c=new TH1F("bX7c","Rate vs. yellBX using bx7; bXing= bx7+offset",120,-0.5,119.5);

  hbx48cm=new TH1F("bX48cm","Masking ON, Rate vs. STAR IP bXing(bx48) ; bXing= bx48+offset",120,-0.5,119.5);

  int t1=time(0);
  chain->Init();
  int t2=time(0);
  printf("%d seconds used to Init() the chain\n",t2-t1);

  chain->ls(3);
  //  return;
  int eventCounter=0;
  int stat=0;
  t1=time(0);

  //---------------------------------------------------
  while ( 1 ) {// loop over events
    if(eventCounter>=nEve) break;
    chain->Clear();
    stat = chain->Make();
    if(stat) break; // end of events
    eventCounter++;
    // dump spinDB only once:
    // if(eventCounter==1)  spDb->print(0);  // 1-large printout, 0-smaler
  
   // Access to muDst .......................
    StMuEvent* muEve = muMk->muDst()->event();
    int nPrim = muMk->muDst()->primaryTracks()->GetEntries();  // get number of primary tracks
    StEventInfo &info=muEve->eventInfo();

    StL0Trigger *trig=&(muEve->l0Trigger());
    int bx48=trig->bunchCrossingId();
    int bx7=trig->bunchCrossingId7bit(muEve->runId());
    hbx48->Fill(bx48);
    hbx7->Fill(bx7);
    hbx748->Fill(bx7,bx48);
    //  printf("w/o DB bx7=%d  bx48=%d    del=%d \n", bx7, bx48,(bx7-bx48+120)%120);

#ifdef USE_DB 
    int bxStar48=spDb->BXyellowUsingBX48(bx48);
    int bxStar7=spDb->BXyellowUsingBX7(bx7);

    //printf("eveID=%d bx7=%d bx48=%d --> yellBX:=%d ? %d --> del=%d\n",info.id(),bx7,bx48, bxStar7,bxStar48, (bxStar7-bxStar48+120)%120);
    assert(spDb->offsetBX48minusBX7(bx48,bx7)==0);
    hbx48c->Fill(bxStar48);
    hbx7c->Fill(bxStar7);
    if(!spDb->isMaskedUsingBX48(bx48) ) hbx48cm->Fill(bxStar48);
#endif

   if(eventCounter%500==0) printf("\n\n ====================%d  processing eventID %d nPrim=%d ==============\n", eventCounter++,info.id(),nPrim);

  }

  if(eventCounter>0) printf("sorting done, nEve=%d of :%d :%d :is missing\n",eventCounter, nEntries,-eventCounter+nEntries);
  t2=time(0);
  if(t2==t1) t2=t1+1;
  float rate=1.*eventCounter/(t2-t1);
  float nMnts=(t2-t1)/60.;  printf("sorting done %d of   nEve=%d, CPU rate=%.1f Hz tot %.1f minutes\n",eventCounter,nEntries,rate,nMnts);
  
  TString fileH=outPath; fileH+=Rrun; fileH+=".hist.root";
  printf(" saving -->%s\n",fileH.Data());
  new TFile(fileH,"recreate");
  hbx48->Write();
  hbx7->Write();
  hbx748->Write();
  hbx48c->Write();
  hbx7c->Write();
  hbx48cm->Write();
  
#ifdef USE_DB
  //play with spinDb information
  // spDb->print(0); // 0=short, 1=huge
  const int * spin8bits=spDb->getSpin8bits();
  for(int bx=0;bx<120;bx++){
    bool isFilled=(spin8bits[bx] & 0x11)==0x11;
    if(isFilled) hbxI->Fill(bx);
    assert(isFilled==spDb->isBXfilledUsingBXyellow(bx));
    // spDb->isBXmaskedUsingBXyellow(bx);
    //    if(spDb->isBXfilledUsingBX48(bx48)) hbxI->Fill(bxStar);
    //printf("bxStar=%3d   bx48=%3d   fill=%d\n",bxStar,bx48,spDb->isBXfilledUsingBX48(bx48));
  }
  hbxI->Write();
#endif


  chain->Finish();
  
  
  
}