Char_t name[50];
if(!fInit){
- Warning("Exec","Initilization not succesfull. Aborting.");
+ Warning("Exec","(opt=%s) Initilization not succesfull. Aborting.\n",opt);
return;
} // end if !fInit
//-----------------------------------------------------------
Int_t AliITSPid::GetPcode(TClonesArray* rps,Float_t pm)
{
+ Info("GetPcode","method not implemented - Inputs TClonesArray *%x , Float_t %f",rps,pm);
return 0;
}
//-----------------------------------------------------------
}
// -----------------------------------------------------------------------
-void AliITSRiemannFit::InitPoints(Int_t evnt,Int_t ntracks,AliITS *ITS,
+void AliITSRiemannFit::InitPoints(Int_t ntracks,AliITS *ITS,
TTree *TR,Int_t nparticles){
//////////////////////////////////////////////////////////////////////
// Fill the class member fPointRecs with the reconstructed points
}
//-------------------------------------------------------------------
-Int_t AliITSRiemannFit::FitHelix(Int_t tracknumber,Int_t charge,Double_t Px,Double_t Py,Double_t Pz,Double_t& fd0,
+Int_t AliITSRiemannFit::FitHelix(Int_t tracknumber,Double_t Px,Double_t Py,Double_t Pz,Double_t& fd0,
Double_t& fphi,Double_t& u0, Double_t& v0, Double_t& rho,Double_t& omega, Double_t& z0,
Double_t& vpar,Double_t& chisql, Double_t& fCorrLin,Double_t& fFit,
Int_t first,Int_t second,Int_t third,Int_t fourth,Int_t fifth,Int_t sixth) {
Float_t GetdY(Int_t i) const {return this->fPointRecs[i]->fdy;}
Float_t GetdZ(Int_t i) const {return this->fPointRecs[i]->fdz;}
- void InitPoints(Int_t evnt,Int_t ntracks,AliITS *ITS,TTree *TR,Int_t nparticles);
+ void InitPoints(Int_t ntracks,AliITS *ITS,TTree *TR,Int_t nparticles);
void WritePoints(void);
void ReadPoints(void);
static Int_t SolveCubic(Double_t a,Double_t b,Double_t c,Double_t&,Double_t&,Double_t&);
- Int_t FitHelix(Int_t tracknumber,Int_t charge,Double_t Px,Double_t Py,Double_t Pz,
+ Int_t FitHelix(Int_t tracknumber,Double_t Px,Double_t Py,Double_t Pz,
Double_t& fd0,Double_t& fphi,Double_t& u0, Double_t& v0, Double_t& rho,
Double_t& omega, Double_t& z0,
Double_t& vpar,Double_t& chisql,Double_t& fCorrLin,Double_t& fFit,
}
-void AliITSTrackV1::AddMS(AliITSRad *rl, Double_t mass) {
+void AliITSTrackV1::AddMS(Double_t mass) {
//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
////////// Modification of the covariance matrix to take into account multiple scattering ///////////
fC44+=theta2*q30*q30;
}
-void AliITSTrackV1::PrimaryTrack(AliITSRad *rl) {
+void AliITSTrackV1::PrimaryTrack() {
//Origin A. Badala' and G.S. Pappalardo: e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
// calculation of part of covariance matrix for vertex constraint
fd2(i)=fC22;
ftgl2(i)=fC33;
fdtgl(i)=fC32;
- AddMS(rl);
+ AddMS();
AddEL(-1,0);
}
}
Double_t GetSigmaphi() const{return fC00;} // gets the phi variance
Double_t GetSigmaZ() const{return fC11;} // gets the Z variance
void AddEL(Double_t signdE, Bool_t flagtot, Double_t mass=0.1396); // adds the energy loss
- void AddMS(AliITSRad *rl,Double_t mass=0.1396); // modify the covariance matrix to take into account the multiple scattering
+ void AddMS(Double_t mass=0.1396); // modify the covariance matrix to take into account the multiple scattering
void Correct(Double_t rk); // correct the track to take into account the real detector geometry
void SetDv(Double_t x) {fDv=x;} // sets the radial impact parameter for vertex constraint
void SetZv(Double_t x) {fZv=x;} // sets longitudinal impact parameter for vertex constraint
void Setfnoclust() {fnoclust++;} //modify fnoclust
Double_t GetsigmaDv() const {return fsigmaDv;} // gets sigma for Dv extraction
Double_t GetsigmaZv() const {return fsigmaZv;} // gets sigma for Zv extraction
- void PrimaryTrack(AliITSRad *rl); // calculation of part of covariance matrix for vertex constraint
+ void PrimaryTrack(); // calculation of part of covariance matrix for vertex constraint
void Setd2(TVector &x) {for(Int_t i=0; i<6; i++){fd2(i)=x(i);}} // sets the vector fd2
void Settgl2(TVector &x) {for(Int_t i=0; i<6; i++){ftgl2(i)=x(i);}} // sets the vector ftgl2
void Setdtgl(TVector &x) {for(Int_t i=0; i<6; i++){fdtgl(i)=x(i);}} // sets the vector fdtgl
primaryTrack.SetZv(zv);
primaryTrack.SetsigmaDv(sigmaDv);
primaryTrack.SetsigmaZv(sigmaZv);
- primaryTrack.PrimaryTrack(frl);
+ primaryTrack.PrimaryTrack();
TVector d2=primaryTrack.Getd2();
TVector tgl2=primaryTrack.Gettgl2();
TVector dtgl=primaryTrack.Getdtgl();
// cout<<" chi2pred = "<<chi2pred<<"\n";
// if(chi2pred>fChi2max) continue; //aggiunto il 30-7-2001
if(iriv == 0) flaghit=1;
- (*newTrack).AddMS(frl); // add the multiple scattering
+ (*newTrack).AddMS(); // add the multiple scattering
//matrix to the covariance matrix
(*newTrack).AddEL(1.,0);
AliITSTrackV1 *newTrack = new AliITSTrackV1(*trackITS);
(*newTrack).Setfnoclust();
//(*newTrack).SetLayer((*trackITS).GetLayer()-1);
- (*newTrack).AddMS(frl); // add the multiple scattering matrix
+ (*newTrack).AddMS(); // add the multiple scattering matrix
// to the covariance matrix
(*newTrack).AddEL(1.,0);
listoftrack.AddLast(newTrack);
delete hITSYv;
Char_t name[30];
// sprintf(name,"Vertex_%d",evnumber);
+ if(fDebug>0)Info("FindVertexForCurrentEvent","Vertex found for event %d",evnumber);
sprintf(name,"Vertex");
fCurrentVertex = new AliITSVertex(Position,Resolution,SNR,name);
return fCurrentVertex;
TString treeName = "TreeT_ITS_";
treeName += evnumb;
// TTree *trkTree=(TTree*)fInFile->Get(treeName.Data()); masera
- TTree *trkTree;
+ TTree *trkTree=0;
if(!trkTree) return fCurrentVertex;
//--------------------------------------------------------------------
//This functions propagates reconstructed ITS tracks back
//--------------------------------------------------------------------
- Error("PropagateBack","This method is not converted to NewIO yet\n");
+ Error("PropagateBack","This method is not converted to NewIO yet (Args: %x , %x )\n",inp,out);
return 1;
/*
TFile *in=(TFile*)inp;
// This functions refits ITS tracks using the
// "inward propagated" TPC tracks
//--------------------------------------------------------------------
- Error("RefitInward","This method is not converted to NewIO yet\n");
+ Error("RefitInward","This method is not converted to NewIO yet (Args: %x , %x )\n",in,out);
return 1;
/*
TFile *in=(TFile*)inp;