[u/mrichter/AliRoot.git] / TPC / CalibMacros / CalibTime.C

/// \file CalibTime.C
///
/// Macro to visualize and analyze time dependent variables
///
/// ~~~
/// .L $ALICE_ROOT/TPC/CalibMacros/CalibTime.C
///
/// // 1. Init - Load libraries tools 
/// Init()
///
/// // 2. Load debug streamers 
/// InitChains()
/// ~~~

gSystem->Load("libANALYSIS");
gSystem->Load("libSTAT");
gSystem->Load("libTPCcalib");

TFile f("CalibObjects.root");
AliTPCcalibTime *calibTime = (AliTPCcalibTime *)f->Get("TPCCalib")->FindObject("calibTime");

TPair * addPair= calibTime->GetMapDz()->FindObject(" D0SCO ");
THnSparse* addHist=dynamic_cast<THnSparseF*>(addPair->Value());

TGraph * gr = AliTPCcalibBase::FitSlices(addHist,2,0,100,100);
gr->SetMarkerColor(2);
gr->Draw("same*");
//Make Fit

AliSplineFit fit;
fit.SetGraph(gr)
fit->SetMinPoints(gr->GetN()+1);
fit->InitKnots(gr,2,0,0.001)
fit.SplineFit(0)
TGraph * grfit = fit.MakeGraph(gr->GetX()[0],gr->GetX()[gr->GetN()-1],50000,0);
gr->SetMarkerStyle(25);
//gr->Draw("alp");
grfit->SetLineColor(2);
grfit->Draw("lu");

//
// Chain Based analysis
//


TChain * chainLaser=0, *chainDz=0, *chaindEdx=0; 

void Init(){  
  /// Load neccesary libraries

  gSystem->Load("libANALYSIS");
  gSystem->Load("libSTAT");
  gSystem->Load("libTPCcalib");
  gROOT->Macro("~/NimStyle.C");
  gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros");
  gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+");
}

void InitChains(){  
  AliXRDPROOFtoolkit tool;   
  //load debug streams
  chainLaser = tool.MakeChain("time.txt","laserInfo",0,10200);
  chainLaser->Lookup();
  // cosmic dZ - drfit velocity part 
  chainDz = tool.MakeChain("time.txt","cosmicDz",0,10200);
  chainDz->Lookup();
  //
  // cosmic dE - drift velocity part 
  chaindEdx = tool.MakeChain("time.txt","cosmicdEdx",0,10200);
  chaindEdx->Lookup();
  //
  // Set Alias
  //
  chaindEdx->SetAlias("side","(-1+(p0.fP[1]>0)*2)");
  chaindEdx->SetAlias("der","side*(dedx0-dedx1)/(dedx0+dedx1)");
  chaindEdx->SetAlias("derIO","side*(dedx0Out-dedx0In)/(dedx0In+dedx0Out)");
  chaindEdx->SetAlias("dr","(1-abs(p0.fP[1]/250))");
  chaindEdx->SetAlias("isOK","dedx0In>0&&dedx0Out>0&&dedx1In>0&&dedx1Out>0");
  chaindEdx->SetAlias("dedxM","(dedx0+dedx1)*0.5");

}


void MakeTglFitCosmic(){
  /// Fit the z correction factor

  TStatToolkit toolkit;
  Double_t chi2=0;
  Int_t    npoints=0;
  TVectorD fitParamC,fitParamIO;
  TMatrixD covMatrix;  
  //
  // dedx delta
  // 
  TString fstring="";
  fstring+="side++";
  fstring+="p0.fP[3]++";
  //fstring+="p0.fP[3]*dr++";
  TString *strTheta = toolkit.FitPlane(chaindEdx,"der",fstring->Data(), "isOK", chi2,npoints,fitParamC,covMatrix);
  chaindEdx->SetAlias("derc",strTheta.Data());
  strTheta->Tokenize("+")->Print(); 
  //
  // dedx delta
  // 
  fstring="";
  fstring+="side++";
  fstring+="p0.fP[3]++";
  //fstring+="p0.fP[3]*dr++";
  TString *strTheta = toolkit.FitPlane(chaindEdx,"derIO",fstring->Data(), "isOK", chi2,npoints,fitParamIO,covMatrix);
  chaindEdx->SetAlias("derIOc",strTheta.Data());
  strTheta->Tokenize("+")->Print(); 
  //
  // Make Plot
  //
  chaindEdx->Draw("der:p0.fP[3]>>his(20,-0.5,0.5)","","prof");
  chaindEdx->Draw("der-derc:p0.fP[3]>>hisC(20,-0.5,0.5)","","prof");
  his->SetXTitle("tan(#theta)");
  his->SetYTitle("(Q_{u}-Q_{d})/(Q_{u}+Q_{d})");
  his->Draw("");
  hisC->Draw("same");
}


void MakeTglFitCosmic(){
  ///

  TStatToolkit toolkit;
  Double_t chi2=0;
  Int_t    npoints=0;
  TVectorD fitParamC,fitParamIO;
  TMatrixD covMatrix;  
  TString fstring="";
  fstring="";
  fstring+="ppit++";
  fstring+="(trigger==1)++";
  fstring+="(trigger==2)++";
  fstring+="(trigger==4)++";
  fstring+="(trigger==8)++";
  TString *strPress = toolkit.FitPlane(chaindEdx,"(dedx0+dedx1)*0.5",fstring->Data(), "isOK", chi2,npoints,fitParamIO,covMatrix);
  //
  //
  chaindEdx->SetAlias("de",strTheta.Data());
  strTheta->Tokenize("+")->Print(); 


}
Commit	Line	Data
35d81915	1	/// \file CalibTime.C
	2	///
	3	/// Macro to visualize and analyze time dependent variables
	4	///
	5	/// ~~~
	6	/// .L $ALICE_ROOT/TPC/CalibMacros/CalibTime.C
	7	///
	8	/// // 1. Init - Load libraries tools
	9	/// Init()
	10	///
	11	/// // 2. Load debug streamers
	12	/// InitChains()
	13	/// ~~~
2be25cc9	14
	15	gSystem->Load("libANALYSIS");
	16	gSystem->Load("libSTAT");
	17	gSystem->Load("libTPCcalib");
	18
2be25cc9	19	TFile f("CalibObjects.root");
	20	AliTPCcalibTime calibTime = (AliTPCcalibTime )f->Get("TPCCalib")->FindObject("calibTime");
	21
2be25cc9	22	TPair * addPair= calibTime->GetMapDz()->FindObject(" D0SCO ");
	23	THnSparse* addHist=dynamic_cast<THnSparseF*>(addPair->Value());
	24
2be25cc9	25	TGraph * gr = AliTPCcalibBase::FitSlices(addHist,2,0,100,100);
	26	gr->SetMarkerColor(2);
	27	gr->Draw("same*");
	28	//Make Fit
	29
	30	AliSplineFit fit;
	31	fit.SetGraph(gr)
	32	fit->SetMinPoints(gr->GetN()+1);
	33	fit->InitKnots(gr,2,0,0.001)
	34	fit.SplineFit(0)
	35	TGraph * grfit = fit.MakeGraph(gr->GetX()[0],gr->GetX()[gr->GetN()-1],50000,0);
	36	gr->SetMarkerStyle(25);
	37	//gr->Draw("alp");
	38	grfit->SetLineColor(2);
	39	grfit->Draw("lu");
	40
2be25cc9	41	//
	42	// Chain Based analysis
	43	//
	44
	45
	46	TChain * chainLaser=0, chainDz=0, chaindEdx=0;
	47
	48	void Init(){
35d81915	49	/// Load neccesary libraries
35d81915	50
2be25cc9	51	gSystem->Load("libANALYSIS");
	52	gSystem->Load("libSTAT");
	53	gSystem->Load("libTPCcalib");
	54	gROOT->Macro("~/NimStyle.C");
	55	gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros");
	56	gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+");
	57	}
	58
	59	void InitChains(){
	60	AliXRDPROOFtoolkit tool;
	61	//load debug streams
	62	chainLaser = tool.MakeChain("time.txt","laserInfo",0,10200);
	63	chainLaser->Lookup();
	64	// cosmic dZ - drfit velocity part
	65	chainDz = tool.MakeChain("time.txt","cosmicDz",0,10200);
	66	chainDz->Lookup();
	67	//
	68	// cosmic dE - drift velocity part
	69	chaindEdx = tool.MakeChain("time.txt","cosmicdEdx",0,10200);
	70	chaindEdx->Lookup();
	71	//
	72	// Set Alias
	73	//
	74	chaindEdx->SetAlias("side","(-1+(p0.fP[1]>0)*2)");
	75	chaindEdx->SetAlias("der","side*(dedx0-dedx1)/(dedx0+dedx1)");
	76	chaindEdx->SetAlias("derIO","side*(dedx0Out-dedx0In)/(dedx0In+dedx0Out)");
	77	chaindEdx->SetAlias("dr","(1-abs(p0.fP[1]/250))");
	78	chaindEdx->SetAlias("isOK","dedx0In>0&&dedx0Out>0&&dedx1In>0&&dedx1Out>0");
	79	chaindEdx->SetAlias("dedxM","(dedx0+dedx1)*0.5");
	80
	81	}
	82
	83
	84	void MakeTglFitCosmic(){
35d81915	85	/// Fit the z correction factor
35d81915	86
2be25cc9	87	TStatToolkit toolkit;
	88	Double_t chi2=0;
	89	Int_t npoints=0;
	90	TVectorD fitParamC,fitParamIO;
	91	TMatrixD covMatrix;
	92	//
	93	// dedx delta
	94	//
	95	TString fstring="";
	96	fstring+="side++";
	97	fstring+="p0.fP[3]++";
	98	//fstring+="p0.fP[3]*dr++";
	99	TString *strTheta = toolkit.FitPlane(chaindEdx,"der",fstring->Data(), "isOK", chi2,npoints,fitParamC,covMatrix);
	100	chaindEdx->SetAlias("derc",strTheta.Data());
	101	strTheta->Tokenize("+")->Print();
	102	//
	103	// dedx delta
	104	//
	105	fstring="";
	106	fstring+="side++";
	107	fstring+="p0.fP[3]++";
	108	//fstring+="p0.fP[3]*dr++";
	109	TString *strTheta = toolkit.FitPlane(chaindEdx,"derIO",fstring->Data(), "isOK", chi2,npoints,fitParamIO,covMatrix);
	110	chaindEdx->SetAlias("derIOc",strTheta.Data());
	111	strTheta->Tokenize("+")->Print();
	112	//
	113	// Make Plot
	114	//
	115	chaindEdx->Draw("der:p0.fP[3]>>his(20,-0.5,0.5)","","prof");
	116	chaindEdx->Draw("der-derc:p0.fP[3]>>hisC(20,-0.5,0.5)","","prof");
	117	his->SetXTitle("tan(#theta)");
	118	his->SetYTitle("(Q_{u}-Q_{d})/(Q_{u}+Q_{d})");
	119	his->Draw("");
	120	hisC->Draw("same");
	121	}
	122
	123
	124	void MakeTglFitCosmic(){
35d81915	125	///
35d81915	126
2be25cc9	127	TStatToolkit toolkit;
	128	Double_t chi2=0;
	129	Int_t npoints=0;
	130	TVectorD fitParamC,fitParamIO;
	131	TMatrixD covMatrix;
	132	TString fstring="";
	133	fstring="";
	134	fstring+="ppit++";
	135	fstring+="(trigger==1)++";
	136	fstring+="(trigger==2)++";
	137	fstring+="(trigger==4)++";
	138	fstring+="(trigger==8)++";
	139	TString strPress = toolkit.FitPlane(chaindEdx,"(dedx0+dedx1)0.5",fstring->Data(), "isOK", chi2,npoints,fitParamIO,covMatrix);
	140	//
	141	//
	142	chaindEdx->SetAlias("de",strTheta.Data());
	143	strTheta->Tokenize("+")->Print();
	144
	145
	146	}