for(Int_t la=0;la<4;la++){//loop on layers
Double_t parp[3];Double_t park[3];Double_t parpi[3];
sp->GetParFitLayer(la,fMom,parp,park,parpi);
- // cout<<"lay = "<<la<<endl;
- //for(Int_t i=0;i<3;i++)cout<<parp[i]<<" "<<park[i]<<" "<<parpi[i]<<endl;
+
Double_t range[6];
range[0]=0.3*parp[1];
range[1]=2.*parp[1];
}
- cout<<"PROTONS"<<endl;
fPBayesp=CookCombinedBayes(condFun,prior,0);
- for(Int_t la=0;la<4;la++)cout<<" lay "<<la<<"->prob "<<condFun[la][0];cout<<endl;cout<<"bayes = "<<fPBayesp<<endl;
- cout<<"KAONS"<<endl;
fPBayesk=CookCombinedBayes(condFun,prior,1);
- for(Int_t la=0;la<4;la++)cout<<" lay "<<la<<"->prob "<<condFun[la][1];cout<<endl;cout<<"bayes = "<<fPBayesk<<endl;
- cout<<"PIONS"<<endl;
fPBayespi=CookCombinedBayes(condFun,prior,2);
- for(Int_t la=0;la<4;la++)cout<<" lay "<<la<<"->prob "<<condFun[la][2];cout<<endl;cout<<"bayes = "<<fPBayespi<<endl;
-
-
fdEdx=dEdx;
}
}
- cout<<"PROTONS"<<endl;
fPBayesp=CookCombinedBayes(condFun,prior,0);
- for(Int_t la=0;la<4;la++)cout<<" lay "<<la<<"->prob "<<condFun[la][0];cout<<endl;cout<<"bayes = "<<fPBayesp<<endl;
- cout<<"KAONS"<<endl;
fPBayesk=CookCombinedBayes(condFun,prior,1);
- for(Int_t la=0;la<4;la++)cout<<" lay "<<la<<"->prob "<<condFun[la][1];cout<<endl;cout<<"bayes = "<<fPBayesk<<endl;
- cout<<"PIONS"<<endl;
fPBayespi=CookCombinedBayes(condFun,prior,2);
- for(Int_t la=0;la<4;la++)cout<<" lay "<<la<<"->prob "<<condFun[la][2];cout<<endl;cout<<"bayes = "<<fPBayespi<<endl;
fdEdx=trackITS->GetdEdx();
}
//______________________________________________________________
void AliITSSteerPid::InitLayer(TString fileITS,TString fileFitPar){
// it opens the files useful for the PID
- printf("££££ InitLayer.....\n");
TFile *fClonarr2=new TFile (fileITS,"r");
fVect2=(TClonesArray*)fClonarr2->Get("vectfitits_0");//truncated mean
fVect2lay1=(TClonesArray*)fClonarr2->Get("vectfitits_1");//lay 1
TFile *fFitPar=new TFile (fileFitPar);
fFitTree=(TTree*)fFitPar->Get("tree");
- printf("££££ InitLayer.....DONE \n");
}
//______________________________________________________________
//______________________________________________________________
void AliITSSteerPid::GetParFitLayer(Int_t nolay,Float_t mom,Double_t *parp,Double_t *park,Double_t *parpi){
- //it gives the parameters of the convoluted functions (WL, MP, WG) for protons, kaons and pions for a given momentum and ITS layer
+ //it gives the parameters of the convoluted functions (WL, MP, WG) for
+ //protons, kaons and pions for a given momentum and ITS layer
Double_t parfit0pro[3]={0,0,0};
Double_t parfit1pro[3]={0,0,0};
}//______________________________________________________________
void AliITSSteerPid::GetLangausProPars(Float_t mom,Double_t *parfit0,Double_t *parfit1,Double_t *parfit3,Double_t *par){
- //It finds the parameters of the convoluted Landau-Gaussian response function for protons (Width Landau, Most Probable, Width Gaussian)
+ //It finds the parameters of the convoluted Landau-Gaussian response
+ //function for protons (Width Landau, Most Probable, Width Gaussian)
par[0]=parfit0[0]+parfit0[1]/mom;
par[1]=parfit1[0]/(mom*mom)+parfit1[1]/(mom*mom)*TMath::Log(mom*mom)+parfit1[2];
par[2]=parfit3[0]/(mom*mom)+parfit3[1]/(mom*mom)*TMath::Log(mom*mom)+parfit3[2];
}
//______________________________________________________________
void AliITSSteerPid::GetLangausKaoPars(Float_t mom,Double_t *parfit0,Double_t *parfit1,Double_t *parfit3,Double_t *par){
- // It finds the parameters of the convoluted Landau-Gaussian response function for kaons (Width Landau, Most Probable, Width Gaussian)
+ // It finds the parameters of the convoluted Landau-Gaussian response
+ //function for kaons (Width Landau, Most Probable, Width Gaussian)
par[0]=parfit0[0]+parfit0[1]/(mom*mom);
par[1]=parfit1[0]/(mom*mom)+parfit1[1]/(mom*mom)*TMath::Log(mom*mom)+parfit1[2];
//______________________________________________________________
void AliITSSteerPid::GetLangausPiPars(Float_t mom,Double_t *parfit0,Double_t *parfit1,Double_t *parfit3,Double_t *par){
- //It finds the parameters of the convoluted Landau-Gaussian response function for pions (Width Landau, Most Probable, Width Gaussian)
+ //It finds the parameters of the convoluted Landau-Gaussian response
+ //function for pions (Width Landau, Most Probable, Width Gaussian)
par[0]=parfit0[0]/(mom*mom)+parfit0[1]/(mom*mom)*TMath::Log(mom*mom)+parfit0[2];
par[1]=parfit1[0]/(mom)+parfit1[1]/(mom)*TMath::Log(mom*mom)+parfit1[2];