[u/mrichter/AliRoot.git] / RICH / RichMenu.C

AliRun *a;    AliStack *s;  AliRunLoader *al; //globals for easy manual manipulations
AliRICH   *r; AliLoader    *rl;

//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Bool_t AliceRead()
{
  Info("ReadAlice","Tring to read ALICE from SIMULATED FILE...");
  if(gAlice){
    delete gAlice->GetRunLoader();
    delete gAlice;
  }      
  
  if(gSystem->Exec("ls galice.root>/dev/null")==256){//there is no galice.root in current directory
    AliceNew();
    RichGet();
    return kFALSE; //new session
  }else{
    if(!(al=AliRunLoader::Open())){//if not possible to read from galice.root, then remove grabage and reinvoke AliceRead()
      gSystem->Exec("rm -rf *.root *.dat");
      AliceRead();
    }
    al->LoadgAlice();//before this gAlice is 0;
    if(!gAlice) Fatal("menu.C::ReadAlice","No gAlice in file");
    a=al->GetAliRun();//provides pointer to AliRun object
    Info("AliceRead","Run contains %i event(s)",a->GetEventsPerRun());      
    RichGet();
    return kTRUE;   //old session opened from file
  }        
}//AliceRead()
//__________________________________________________________________________________________________
void AliceNew()
{
  Info("AliceNew","Init new session");
  new AliRun("gAlice","Alice experiment system");  gAlice->Init(); a=gAlice; al=gAlice->GetRunLoader();
}//AliceNew()    
//__________________________________________________________________________________________________
void RichGet()
{
  if(!(r=(AliRICH*)al->GetAliRun()->GetDetector("RICH")))     Warning("RICH/menu.C::ReadAlice","No RICH in file");
  if(!(rl=al->GetDetectorLoader("RICH")))           Warning("RICH/menu.C::ReadAlice","No RICH loader in file");        
}

//__________________________________________________________________________________________________
void RichMenu()
{   
  TControlBar *pMenu = new TControlBar("vertical","MAIN");
       
  if(AliceRead()){//it's from file, show some info
    pMenu->AddButton("Display single chambers"         ,"r->Display();"  , "Display Fast");
    pMenu->AddButton("Display ALL chambers"            ,"r->DisplayEvent(0,0);"  , "Display Fast");
    pMenu->AddButton("Recon with stack"                ,"AliRICHReconstructor::CheckPR(        )","Create RSR.root with ntuple hn");    
    pMenu->AddButton("RichAna no Recon"                ,"AliRICHReconstructor::RichAna(0,0,kFALSE)","Create RichAna.root with ntuple hn without PatRec");    
    pMenu->AddButton("RichAna with Recon"              ,"AliRICHReconstructor::RichAna(0,0,kTRUE )","Create RichAna.root with ntuple hn with PatRec");    
    pMenu->AddButton("HITS Print"                      ,"h();"      ,"To print hits: h()");
    pMenu->AddButton("HITS QA"                         ,"hqa()"     ,"QA plots for hits: hqa()");
    pMenu->AddButton("Print sdigits"                   ,"s();"      ,"To print sdigits: s()");
    pMenu->AddButton("Print digits"                    ,"d();"      ,"To print digits: d()");
    pMenu->AddButton("Clusters print"                  ,"c();"      ,"To print clusters: c()");  
    pMenu->AddButton("Clusters QA"                     ,"cqa();"    ,"Clusters QA: cqa() or AliRICHReconstructor::CluQA(al)");  
    pMenu->AddButton("Print ESD"                       ,"e();"      ,"To print ESD status");  
    pMenu->AddButton("Print Matrix"                    ,"m();"      ,"To print probability matrix");  
    pMenu->AddButton("Print occupancy"                 ,"r->OccupancyPrint(-1);" ,"To print occupancy");  
    pMenu->AddButton("Print event summary  "           ,"r->SummaryOfEvent();"   ,"To print a summary of the event");  
  }else{//it's aliroot, simulate
    pMenu->AddButton("Debug ON",     "DebugON();",   "Switch debug on-off");   
    pMenu->AddButton("Debug OFF",    "DebugOFF();",  "Switch debug on-off");   
    pMenu->AddButton("Run",          "a->Run(1)",  "Process!");
  }
  pMenu->AddButton("Test segmentation"  ,"rp->TestSeg()"  ,"Test AliRICHParam segmentation methods"     );
  pMenu->AddButton("Test response"      ,"rp->TestResp()" ,"Test AliRICHParam response methods"         );
  pMenu->AddButton("Test transformation","rp->TestTrans()","Test AliRICHParam transformation methods"   );
  pMenu->AddButton("Geo GUI"            ,"geo();"                                                          ,"Shows geometry"                             ); 
  pMenu->AddButton("Debug ON"           ,"AliLog::SetGlobalDebugLevel(AliLog::kDebug);"                    ,"Switch debug on"                            );   
  pMenu->AddButton("Debug OFF"          ,"AliLog::SetGlobalDebugLevel(0);"                                 ,"Switch debug off"                           );   
  pMenu->AddButton("Browser"            ,"new TBrowser;"                                                   ,"Start ROOT TBrowser"                        );
  pMenu->AddButton("Quit"               ,".q"                                                              ,"Close session"                              );
  pMenu->Show();
}//menu()
//__________________________________________________________________________________________________
void DebugOFF(){  Info("DebugOFF","");  AliLog::SetGlobalDebugLevel(0);}
void DebugON() {  Info("DebugON","");   AliLog::SetGlobalDebugLevel(AliLog::kDebug);}
//void geo (            )   { gGeoManager->SetVisOption(0);gGeoManager->GetTopVolume()->Draw(); AliRICHParam::DrawAxis();}
void geo()                  {  if(!gGeoManager) TGeoManager::Import("geometry.root");gGeoManager->GetTopVolume()->Draw();AliRICHParam::DrawAxis();}
  
void dis (Int_t evt=-1)   {r->Display      (evt);}                //utility display 
void dum (            )   {r->Dump         (   );}                //utility display 

void h   (Int_t evt=0 )   {r->HitPrint  (evt);}   //print hits for requested event
void hqa (            )   {r->HitQA     (   );}   //hits QA plots for all events 

void s   (Int_t evt=0 )   {r->SDigPrint (evt);}   //print sdigits for requested event

void d   (Int_t evt=0 )   {r->DigPrint  (evt);}   //print digits for requested event
void dqa (            )   {AliRICHReconstructor::DigQA     (al );}   //digits QA plots for all events


void c   (Int_t evt=0 )   {r->CluPrint  (evt);}   //print clusters for requested event
void cqa (            )   {AliRICHReconstructor::CluQA     (al );}   //clusters QA plots for all events
void ct  (            )   {AliRICHReconstructor::Test      (   );}   //test clusters by predifined list of digits

void t   (Int_t evt=0          )   {AliRICHParam::Stack(evt);}    
void tid (Int_t tid,Int_t evt=0)   {AliRICHParam::Stack(evt,tid);} 
void e   (                     )   {AliRICHTracker::EsdPrint();}                   
void m   (Double_t probCut=0.7 )   {AliRICHTracker::MatrixPrint(probCut);}                   

Int_t nem (Int_t evt=0) {AliRICHParam::StackCount(kElectron  ,evt);} //utility number of electrons
Int_t nep (Int_t evt=0) {AliRICHParam::StackCount(kPositron  ,evt);} //utility number of positrons
Int_t nmup(Int_t evt=0) {AliRICHParam::StackCount(kMuonPlus  ,evt);} //utility number of positive muons
Int_t nmum(Int_t evt=0) {AliRICHParam::StackCount(kMuonMinus ,evt);} //utility number of negative muons
Int_t npi0(Int_t evt=0) {AliRICHParam::StackCount(kPi0       ,evt);} //utility number of neutral pions 
Int_t npip(Int_t evt=0) {AliRICHParam::StackCount(kPiPlus    ,evt);} //utility number of positive pions
Int_t npim(Int_t evt=0) {AliRICHParam::StackCount(kPiMinus   ,evt);} //utility number of negative pions
Int_t nk0 (Int_t evt=0) {AliRICHParam::StackCount(kK0        ,evt);} //utility number of neutral kaons
Int_t nkp (Int_t evt=0) {AliRICHParam::StackCount(kKPlus     ,evt);} //utility number of positive kaons
Int_t nkm (Int_t evt=0) {AliRICHParam::StackCount(kKMinus    ,evt);} //utility number of negative kaons
Int_t npp (Int_t evt=0) {AliRICHParam::StackCount(kProton    ,evt);} //utility number of protons
Int_t npm (Int_t evt=0) {AliRICHParam::StackCount(kProtonBar ,evt);} //utility number of antiprotons
Commit	Line	Data
	1	AliRun a; AliStack s; AliRunLoader *al; //globals for easy manual manipulations
	2	AliRICH r; AliLoader rl;
	3
	4	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	5	Bool_t AliceRead()
	6	{
	7	Info("ReadAlice","Tring to read ALICE from SIMULATED FILE...");
	8	if(gAlice){
	9	delete gAlice->GetRunLoader();
	10	delete gAlice;
	11	}
	12
	13	if(gSystem->Exec("ls galice.root>/dev/null")==256){//there is no galice.root in current directory
	14	AliceNew();
	15	RichGet();
	16	return kFALSE; //new session
	17	}else{
	18	if(!(al=AliRunLoader::Open())){//if not possible to read from galice.root, then remove grabage and reinvoke AliceRead()
	19	gSystem->Exec("rm -rf .root .dat");
	20	AliceRead();
	21	}
	22	al->LoadgAlice();//before this gAlice is 0;
	23	if(!gAlice) Fatal("menu.C::ReadAlice","No gAlice in file");
	24	a=al->GetAliRun();//provides pointer to AliRun object
	25	Info("AliceRead","Run contains %i event(s)",a->GetEventsPerRun());
	26	RichGet();
	27	return kTRUE; //old session opened from file
	28	}
	29	}//AliceRead()
	30	//__________________________________________________________________________________________________
	31	void AliceNew()
	32	{
	33	Info("AliceNew","Init new session");
	34	new AliRun("gAlice","Alice experiment system"); gAlice->Init(); a=gAlice; al=gAlice->GetRunLoader();
	35	}//AliceNew()
	36	//__________________________________________________________________________________________________
	37	void RichGet()
	38	{
	39	if(!(r=(AliRICH*)al->GetAliRun()->GetDetector("RICH"))) Warning("RICH/menu.C::ReadAlice","No RICH in file");
	40	if(!(rl=al->GetDetectorLoader("RICH"))) Warning("RICH/menu.C::ReadAlice","No RICH loader in file");
	41	}
	42
	43	//__________________________________________________________________________________________________
	44	void RichMenu()
	45	{
	46	TControlBar *pMenu = new TControlBar("vertical","MAIN");
	47
	48	if(AliceRead()){//it's from file, show some info
	49	pMenu->AddButton("Display single chambers" ,"r->Display();" , "Display Fast");
	50	pMenu->AddButton("Display ALL chambers" ,"r->DisplayEvent(0,0);" , "Display Fast");
	51	pMenu->AddButton("Recon with stack" ,"AliRICHReconstructor::CheckPR( )","Create RSR.root with ntuple hn");
	52	pMenu->AddButton("RichAna no Recon" ,"AliRICHReconstructor::RichAna(0,0,kFALSE)","Create RichAna.root with ntuple hn without PatRec");
	53	pMenu->AddButton("RichAna with Recon" ,"AliRICHReconstructor::RichAna(0,0,kTRUE )","Create RichAna.root with ntuple hn with PatRec");
	54	pMenu->AddButton("HITS Print" ,"h();" ,"To print hits: h()");
	55	pMenu->AddButton("HITS QA" ,"hqa()" ,"QA plots for hits: hqa()");
	56	pMenu->AddButton("Print sdigits" ,"s();" ,"To print sdigits: s()");
	57	pMenu->AddButton("Print digits" ,"d();" ,"To print digits: d()");
	58	pMenu->AddButton("Clusters print" ,"c();" ,"To print clusters: c()");
	59	pMenu->AddButton("Clusters QA" ,"cqa();" ,"Clusters QA: cqa() or AliRICHReconstructor::CluQA(al)");
	60	pMenu->AddButton("Print ESD" ,"e();" ,"To print ESD status");
	61	pMenu->AddButton("Print Matrix" ,"m();" ,"To print probability matrix");
	62	pMenu->AddButton("Print occupancy" ,"r->OccupancyPrint(-1);" ,"To print occupancy");
	63	pMenu->AddButton("Print event summary " ,"r->SummaryOfEvent();" ,"To print a summary of the event");
	64	}else{//it's aliroot, simulate
	65	pMenu->AddButton("Debug ON", "DebugON();", "Switch debug on-off");
	66	pMenu->AddButton("Debug OFF", "DebugOFF();", "Switch debug on-off");
	67	pMenu->AddButton("Run", "a->Run(1)", "Process!");
	68	}
	69	pMenu->AddButton("Test segmentation" ,"rp->TestSeg()" ,"Test AliRICHParam segmentation methods" );
	70	pMenu->AddButton("Test response" ,"rp->TestResp()" ,"Test AliRICHParam response methods" );
	71	pMenu->AddButton("Test transformation","rp->TestTrans()","Test AliRICHParam transformation methods" );
	72	pMenu->AddButton("Geo GUI" ,"geo();" ,"Shows geometry" );
	73	pMenu->AddButton("Debug ON" ,"AliLog::SetGlobalDebugLevel(AliLog::kDebug);" ,"Switch debug on" );
	74	pMenu->AddButton("Debug OFF" ,"AliLog::SetGlobalDebugLevel(0);" ,"Switch debug off" );
	75	pMenu->AddButton("Browser" ,"new TBrowser;" ,"Start ROOT TBrowser" );
	76	pMenu->AddButton("Quit" ,".q" ,"Close session" );
	77	pMenu->Show();
	78	}//menu()
	79	//__________________________________________________________________________________________________
	80	void DebugOFF(){ Info("DebugOFF",""); AliLog::SetGlobalDebugLevel(0);}
	81	void DebugON() { Info("DebugON",""); AliLog::SetGlobalDebugLevel(AliLog::kDebug);}
	82	//void geo ( ) { gGeoManager->SetVisOption(0);gGeoManager->GetTopVolume()->Draw(); AliRICHParam::DrawAxis();}
	83	void geo() { if(!gGeoManager) TGeoManager::Import("geometry.root");gGeoManager->GetTopVolume()->Draw();AliRICHParam::DrawAxis();}
	84
	85	void dis (Int_t evt=-1) {r->Display (evt);} //utility display
	86	void dum ( ) {r->Dump ( );} //utility display
	87
	88	void h (Int_t evt=0 ) {r->HitPrint (evt);} //print hits for requested event
	89	void hqa ( ) {r->HitQA ( );} //hits QA plots for all events
	90
	91	void s (Int_t evt=0 ) {r->SDigPrint (evt);} //print sdigits for requested event
	92
	93	void d (Int_t evt=0 ) {r->DigPrint (evt);} //print digits for requested event
	94	void dqa ( ) {AliRICHReconstructor::DigQA (al );} //digits QA plots for all events
	95
	96
	97	void c (Int_t evt=0 ) {r->CluPrint (evt);} //print clusters for requested event
	98	void cqa ( ) {AliRICHReconstructor::CluQA (al );} //clusters QA plots for all events
	99	void ct ( ) {AliRICHReconstructor::Test ( );} //test clusters by predifined list of digits
	100
	101	void t (Int_t evt=0 ) {AliRICHParam::Stack(evt);}
	102	void tid (Int_t tid,Int_t evt=0) {AliRICHParam::Stack(evt,tid);}
	103	void e ( ) {AliRICHTracker::EsdPrint();}
	104	void m (Double_t probCut=0.7 ) {AliRICHTracker::MatrixPrint(probCut);}
	105
	106	Int_t nem (Int_t evt=0) {AliRICHParam::StackCount(kElectron ,evt);} //utility number of electrons
	107	Int_t nep (Int_t evt=0) {AliRICHParam::StackCount(kPositron ,evt);} //utility number of positrons
	108	Int_t nmup(Int_t evt=0) {AliRICHParam::StackCount(kMuonPlus ,evt);} //utility number of positive muons
	109	Int_t nmum(Int_t evt=0) {AliRICHParam::StackCount(kMuonMinus ,evt);} //utility number of negative muons
	110	Int_t npi0(Int_t evt=0) {AliRICHParam::StackCount(kPi0 ,evt);} //utility number of neutral pions
	111	Int_t npip(Int_t evt=0) {AliRICHParam::StackCount(kPiPlus ,evt);} //utility number of positive pions
	112	Int_t npim(Int_t evt=0) {AliRICHParam::StackCount(kPiMinus ,evt);} //utility number of negative pions
	113	Int_t nk0 (Int_t evt=0) {AliRICHParam::StackCount(kK0 ,evt);} //utility number of neutral kaons
	114	Int_t nkp (Int_t evt=0) {AliRICHParam::StackCount(kKPlus ,evt);} //utility number of positive kaons
	115	Int_t nkm (Int_t evt=0) {AliRICHParam::StackCount(kKMinus ,evt);} //utility number of negative kaons
	116	Int_t npp (Int_t evt=0) {AliRICHParam::StackCount(kProton ,evt);} //utility number of protons
	117	Int_t npm (Int_t evt=0) {AliRICHParam::StackCount(kProtonBar ,evt);} //utility number of antiprotons