Adding invert function -

change the curvature to opposite

(Marian)

diff --git a/TPC/AliTPCcalibCosmic.cxx b/TPC/AliTPCcalibCosmic.cxx
index 42b0f5a5b420b013529897ffa486a35df22d3226..5cc31d2490d9ac55413da1cb1cf245cf20a1d8e3 100644 (file)
--- a/TPC/AliTPCcalibCosmic.cxx
+++ b/TPC/AliTPCcalibCosmic.cxx
@@ -477,7 +477,19 @@ void AliTPCcalibCosmic::BinLogX(TH1 *h) {
   
 }
 
-
+AliExternalTrackParam *AliTPCcalibCosmic::Invert(AliExternalTrackParam *input)
+{
+  //
+  // Invert paramerameter  - not covariance yet
+  //
+  AliExternalTrackParam *output = new AliExternalTrackParam(*input);
+  Double_t * param = (Double_t*)output->GetParameter();
+  param[0]*=-1;
+  param[3]*=-1;
+  param[4]*=-1;
+  //
+  return output;
+}
 
 
 /*
@@ -537,3 +549,58 @@ TString * thetaParam = TStatToolkit::FitPlane(chain,"Tr0.fP[3]+Tr1.fP[3]", strFi
 */
 
 
+/*
+gSystem->Load("libANALYSIS");
+gSystem->Load("libSTAT");
+gSystem->Load("libTPCcalib");
+
+TStatToolkit toolkit;
+Double_t chi2;
+TVectorD fitParam;
+TMatrixD covMatrix;
+Int_t npoints;
+//
+TCut cutT("cutT","abs(Tr1.fP[3]+Tr0.fP[3])<0.03");  // OK
+TCut cutD("cutD","abs(Tr0.fP[0]+Tr1.fP[0])<5");     // OK
+TCut cutPt("cutPt","abs(Tr1.fP[4]+Tr0.fP[4])<0.2&&abs(Tr0.fP[4])+abs(Tr1.fP[4])<10");
+TCut cutN("cutN","min(Orig0.fTPCncls,Orig1.fTPCncls)>110");
+TCut cutA="abs(norm-1)<0.3"+cutT+cutD+cutPt+cutN;
+
+
+
+TTree * chain = Track0;
+
+
+chain->SetAlias("norm","signalTot0.fElements[3]/signalTot1.fElements[3]");
+//
+chain->SetAlias("dr1","(signalTot0.fElements[1]/signalTot0.fElements[3])");
+chain->SetAlias("dr2","(signalTot0.fElements[2]/signalTot0.fElements[3])");
+chain->SetAlias("dr4","(signalTot0.fElements[4]/signalTot0.fElements[3])");
+chain->SetAlias("dr5","(signalTot0.fElements[5]/signalTot0.fElements[3])");
+
+TString fstring="";  
+fstring+="dr1++";
+fstring+="dr2++";
+fstring+="dr4++";
+fstring+="dr5++";
+//
+fstring+="dr1*dr2++";
+fstring+="dr1*dr4++";
+fstring+="dr1*dr5++";
+fstring+="dr2*dr4++";
+fstring+="dr2*dr5++";
+fstring+="dr4*dr5++";
+
+
+
+TString *strqdedx = toolkit.FitPlane(chain,"norm",fstring->Data(), cutA, chi2,npoints,fitParam,covMatrix,-1,0,200000);
+  
+chain->SetAlias("corQT",strqdedx->Data());
+
+*/
+
+
+
+
+
+

diff --git a/TPC/AliTPCcalibCosmic.h b/TPC/AliTPCcalibCosmic.h
index fe1664825db68e2e7a367afdc55cae5c6170edc4..a9b058c16b64be1108d1dbea22071a63a798fc14 100644 (file)
--- a/TPC/AliTPCcalibCosmic.h
+++ b/TPC/AliTPCcalibCosmic.h
@@ -30,6 +30,7 @@ public:
   void              SetGainMap(AliTPCCalPad *GainMap){fGainMap = GainMap;};
   static void       CalculateBetheParams(TH2F *hist, Double_t * initialParam);
   static Double_t   CalculateMIPvalue(TH1F * hist);
+  AliExternalTrackParam *Invert(AliExternalTrackParam *input);
   //
   TH1F   *          GetHistNTracks(){return fHistNTracks;};
   TH1F   *          GetHistClusters(){return fClusters;};