[u/mrichter/AliRoot.git] / PWG1 / TRD / AliTRDv0Monitor.cxx

/**************************************************************************
* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
*                                                                        *
* Author: The ALICE Off-line Project.                                    *
* Contributors are mentioned in the code where appropriate.              *
*                                                                        *
* Permission to use, copy, modify and distribute this software and its   *
* documentation strictly for non-commercialf purposes is hereby granted   *
* without fee, provided that the above copyright notice appears in all   *
* copies and that both the copyright notice and this permission notice   *
* appear in the supporting documentation. The authors make no claims     *
* about the suitability of this software for any purpose. It is          *
* provided "as is" without express or implied warranty.                  *
**************************************************************************/

/////////////////////////////////////////////////////
//
//  Monitor V0 for TRD
//
//  Authors:                                          
//  Markus Heide <mheide@uni-muenster.de> 
//////////////////////////////////////////////////////
    
//#include "TPDGCode.h"
/*#include "TFile.h"
#include "TTree.h"
#include "TEventList.h"*/
#include "TObjArray.h"
#include "TH2.h"
#include "TH2F.h"
#include "TH1I.h"
#include "TCanvas.h"
// #include "TPad.h"
// #include "TLegend.h"

#include "AliLog.h"
#include "AliESDtrack.h"
#include "AliPID.h"
#include "AliTRDv0Monitor.h"
#include "info/AliTRDv0Info.h"


ClassImp(AliTRDv0Monitor)

//________________________________________________________________________
AliTRDv0Monitor::AliTRDv0Monitor() 
  :AliTRDrecoTask()
  ,fhQualityReductions(NULL)
  ,fV0s(NULL)
  ,fData(NULL)
  ,fInfo(NULL)
  ,fP(-1.) 
{
  //
  // Default constructor
  //
  SetNameTitle("v0Monitor", "Reference V0 Monitor");
}

//________________________________________________________________________
AliTRDv0Monitor::AliTRDv0Monitor(const char *name, const char *title) 
  :AliTRDrecoTask(name, title)
  ,fhQualityReductions(NULL)
  ,fV0s(NULL)
  ,fData(NULL)
  ,fInfo(NULL)
  ,fP(-1.)
{
  //
  // Default constructor
  //
  DefineInput(2, TObjArray::Class()); // v0 list
  DefineInput(3, TObjArray::Class()); // pid info list 
}


//____________________________________________________________________
void AliTRDv0Monitor::MakeSummary(){
  TCanvas *cOut = new TCanvas("v0MonitorSummary1", "Summary 1 for task V0Monitor", 1024, 768);
  cOut->cd();
  //GetRefFigure(4);
  cOut->SaveAs("V0MonitorSummary.gif");

  cOut = new TCanvas("v0MonitorSummary2","Summary 2 for task V0Monitor", 1024, 768);
  cOut->cd();
  //GetRefFigure(5);
  cOut->SaveAs("V0MonitorSummary2.gif");
}

//________________________________________________________________________
Bool_t AliTRDv0Monitor::GetRefFigure(Int_t /*ifig*/)
{

  AliInfo("Implementation on going ...");
  return kTRUE;
}

//________________________________________________________________________
void AliTRDv0Monitor::UserCreateOutputObjects()
{
  // Create histograms
  // Called once


fContainer = new TObjArray(kNPlots);
fContainer->SetName("V0Monitoring");


const char *samplename[AliPID::kSPECIES] = {"electrons","muons","pions","kaons","protons"};
const char *decayname[AliTRDv0Info::kNDecays] = {"gamma","K0s","Lambda","AntiLambda"};
const char *detectorname[kNDets] = {"ITS","TPC","TOF"};

fhQualityReductions = new TH1I(Form("fhQualityReductions"),Form("Number of tracks cut out by different quality cut steps"),11,-9,2);
fContainer->Add(fhQualityReductions);

for(Int_t ipart = 0;ipart < AliPID::kSPECIES; ipart++){
  fhCutReductions[ipart] = new TH1I(Form("fhCutReductions_%s",samplename[ipart]),Form("Number of tracks cut out by different cut steps for %s",samplename[ipart]),19,-17,2);
  fContainer->Add(fhCutReductions[ipart]);
  for(Int_t idetector = 0; idetector < kNDets; idetector++){
    fhDetPID[idetector][ipart] = new TH2F(Form("fhDetector_%s_%s",detectorname[idetector],samplename[ipart]),Form("%s Likelihood for %s vs. momentum",detectorname[idetector], samplename[ipart]),100,0.,10.,100, 0., 1.);

    fContainer->Add(fhDetPID[idetector][ipart]);
  }
  fhComPID[ipart] = new TH2F(Form("fhComPID_%s",samplename[ipart]),Form("Combined TPC/TOF PID: Likelihood for %s",samplename[ipart]),100,0.,10.,100,0.,1.);

    fContainer->Add(fhComPID[ipart]);

    for(Int_t cutstep = 0; cutstep < kNCutSteps; cutstep++){
      fhTPCdEdx[ipart][cutstep] = new TH2F(Form("fhTPCdEdx_%s_[%d]",samplename[ipart],cutstep),Form("TPC dE/dx for %s [%d]",samplename[ipart],cutstep),100,0.,10.,300,0.,300.);
      
      fContainer->Add(fhTPCdEdx[ipart][cutstep]);
    }
}

for(Int_t iDecay = 0; iDecay < AliTRDv0Info::kNDecays; iDecay++){   
  for(Int_t cutstep =0; cutstep < kNCutSteps; cutstep++){
    fhV0Chi2ndf[iDecay][cutstep] =  new TH2F(Form("fhV0Chi2ndf_%s_[%d]",decayname[iDecay],cutstep),Form("Chi2/NDF vs. momentum"),100,0.,10.,500, 0., 500.);
    
    fContainer->Add(fhV0Chi2ndf[iDecay][cutstep]);
    
    fhPsiPair[iDecay][cutstep] =  new TH2F(Form("fhV0PsiPair_%s_[%d]",decayname[iDecay],cutstep),Form("Psi_pair vs. momentum"),100,0.,10.,200, 0., 1.6);
    
    fContainer->Add(fhPsiPair[iDecay][cutstep]);

    fhPointAngle[iDecay][cutstep] =  new TH2F(Form("fhPointAngle_%s_[%d]",decayname[iDecay],cutstep),Form("Pointing Angle vs. momentum"),100,0.,10.,500, 0., 1.6);     
    fContainer->Add(fhPointAngle[iDecay][cutstep]);

    fhDCA[iDecay][cutstep] =  new TH2F(Form("fhDCA_%s_[%d]",decayname[iDecay],cutstep),Form("V0 Daughter DCA vs. momentum"),100,0.,10.,500, 0., 1.);
    
    fContainer->Add(fhDCA[iDecay][cutstep]);

    fhOpenAngle[iDecay][cutstep] =  new TH2F(Form("fhOpenAngle_%s_[%d]",decayname[iDecay],cutstep),Form("Opening Angle vs. momentum"),100,0.,10.,500, 0., 1.6);
    
    fContainer->Add(fhOpenAngle[iDecay][cutstep]);

    fhRadius[iDecay][cutstep] =  new TH2F(Form("fhRadius_%s_[%d]",decayname[iDecay],cutstep),Form("V0 Generation Radius vs. momentum"),100,0.,10.,500, 0., 150.);
    
    fContainer->Add(fhRadius[iDecay][cutstep]);
  }

  fhInvMass[iDecay] =  new TH2F(Form("fhInvMass_%s",decayname[iDecay]),Form("Invariant Mass vs. momentum"),100,0.,10.,500, 0., 2.);
    
  fContainer->Add(fhInvMass[iDecay]); 

} 

/*TH1F *hV0mcPID[AliPID::kSPECIES][AliPID::kSPECIES];
  Int_t nPBins = 200;
  
  
  for(Int_t iSpecies = 0; iSpecies < AliPID::kSPECIES; iSpecies++){
  for(Int_t iSample = 0; iSample < AliPID::kSPECIES; iSample++){
  
  fhV0mcPID[iSample][iSpecies] = new TH1F(Form("fhV0mcPID_%s_is_%s",name[iSample],name[iSpecies]),Form("%s contained in %s sample",name[iSpecies],name[iSample]), nPBins, 0.2, 13.);
  }
  }*/
}
//________________________________________________________________________
void AliTRDv0Monitor::UserExec(Option_t *) 
{
  // Main loop
  // Called for each event
  
  if(!(fTracks = dynamic_cast<TObjArray*>(GetInputData(1)))) return;
  if(!(fV0s    = dynamic_cast<TObjArray*>(GetInputData(2)))) return;
  if(!(fInfo   = dynamic_cast<TObjArray*>(GetInputData(3)))) return;
  
  
  AliTRDtrackInfo     *track = NULL;
  AliTRDv0Info *v0(NULL);


  for(Int_t itrk=0; itrk<fTracks->GetEntriesFast(); itrk++){
    track = (AliTRDtrackInfo*)fTracks->UncheckedAt(itrk);
    for(Int_t iv(0); iv<fV0s->GetEntriesFast(); iv++){
      if(!(v0 = (AliTRDv0Info*)fV0s->At(iv))) continue;
      if(!v0->HasTrack(track)) continue;
      ULong_t status = track->GetStatus();
      if(!(status&AliESDtrack::kTRDpid)) continue;
      
      fhQualityReductions->Fill(v0->fQuality);//fills integer codes for tracks cut out by track/V0 quality cuts
      
      if(!(v0->fQuality == 1))continue;

      for(Int_t part = 0; part < AliPID::kSPECIES; part++){
  fhCutReductions[part]->Fill(v0->GetPID(part,track));//fill in numbers of tracks eliminated by different PID cuts
      }

    
      for(Int_t idecay(0), part(-1); idecay <  AliTRDv0Info::kNDecays; idecay++){//loop over decay types considered for reference data
  
  if(idecay ==  AliTRDv0Info::kLambda){ //protons and pions from Lambda
    part = AliPID::kProton;
  } else if(idecay == AliTRDv0Info::kAntiLambda) { //antiprotons and pions from Anti-Lambda     
    part = AliPID::kProton;
  } else if(idecay ==   AliTRDv0Info::kK0s) {//pions from K0s
    part = AliPID::kPion;
  } else if(idecay ==  AliTRDv0Info::kGamma) {//electrons from conversions
    part = AliPID::kElectron;
  } 
  
  //fill histograms with track/V0 quality cuts only
  fhPsiPair[idecay][0]->Fill(v0->fV0Momentum,v0->fPsiPair);//Angle between daughter momentum plane and plane perpendicular to magnetic field
  fhInvMass[idecay]->Fill(v0->fV0Momentum,v0->fInvMass[idecay]);//Invariant mass
  fhPointAngle[idecay][0]->Fill(v0->fV0Momentum,v0->fPointingAngle);// = TMath::ACos(esdv0->GetV0CosineOfPointingAngle()); // Cosine of pointing angle
  fhOpenAngle[idecay][0]->Fill(v0->fV0Momentum,v0->fOpenAngle);// opening angle between daughters
  fhDCA[idecay][0]->Fill(v0->fV0Momentum,v0->fDCA);// Distance of closest approach of daughter tracks	
  fhV0Chi2ndf[idecay][0]->Fill(v0->fV0Momentum,v0->fChi2ndf[idecay]);//Kalman Filter Chi2/NDF
  fhRadius[idecay][0]->Fill(v0->fV0Momentum,v0->fRadius);//distance of decay/conversion from primary vertex in x-y plane

  if(v0->HasTrack(track) == -1){	  
    fhTPCdEdx[part][0]->Fill(v0->fTrackN->P(),v0->fTPCdEdx[AliTRDv0Info::kNeg]);//TPC dE/dx for negative track
  } else if(v0->HasTrack(track) == 1){
    fhTPCdEdx[part][0]->Fill(v0->fTrackP->P(),v0->fTPCdEdx[AliTRDv0Info::kPos]);//TPC dE/dx for positive track
  }

  //fill histograms after invariant mass cuts
  if((v0->fInvMass[idecay] < v0->fUpInvMass[idecay][0])&&(v0->fInvMass[idecay]> v0->fDownInvMass[idecay])){
    fhV0Chi2ndf[idecay][1]->Fill(v0->fV0Momentum,v0->fChi2ndf[idecay]);
    fhPsiPair[idecay][1]->Fill(v0->fV0Momentum,v0->fPsiPair);
    fhPointAngle[idecay][1]->Fill(v0->fV0Momentum,v0->fPointingAngle);
    fhOpenAngle[idecay][1]->Fill(v0->fV0Momentum,v0->fOpenAngle);
    fhDCA[idecay][1]->Fill(v0->fV0Momentum,v0->fDCA);
    fhRadius[idecay][1]->Fill(v0->fV0Momentum,v0->fRadius);
    if(v0->HasTrack(track) == -1)
      fhTPCdEdx[part][1]->Fill(v0->fTrackN->P(),v0->fTPCdEdx[AliTRDv0Info::kNeg]);
    else if(v0->HasTrack(track) == 1)
      fhTPCdEdx[part][1]->Fill(v0->fTrackP->P(),v0->fTPCdEdx[AliTRDv0Info::kPos]);

  }

  //fill histograms after all reference selection cuts
  if(v0->GetPID(part,track)==1){
    fhV0Chi2ndf[idecay][2]->Fill(v0->fV0Momentum,v0->fChi2ndf[idecay]);
    fhPsiPair[idecay][2]->Fill(v0->fV0Momentum,v0->fPsiPair);
    fhPointAngle[idecay][2]->Fill(v0->fV0Momentum,v0->fPointingAngle);
    fhOpenAngle[idecay][2]->Fill(v0->fV0Momentum,v0->fOpenAngle);
    fhDCA[idecay][2]->Fill(v0->fV0Momentum,v0->fDCA);
    fhRadius[idecay][2]->Fill(v0->fV0Momentum,v0->fRadius);
    if(v0->HasTrack(track) == -1)
      fhTPCdEdx[part][2]->Fill(v0->fTrackN->P(),v0->fTPCdEdx[AliTRDv0Info::kNeg]);
    else if(v0->HasTrack(track) == 1)
      fhTPCdEdx[part][2]->Fill(v0->fTrackP->P(),v0->fTPCdEdx[AliTRDv0Info::kPos]);

  }
      }
    }
  }
  PostData(1, fContainer);
}
Commit	Line	Data
0ed6f095	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercialf purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/////////////////////////////////////////////////////
	17	//
	18	// Monitor V0 for TRD
	19	//
	20	// Authors:
	21	// Markus Heide <mheide@uni-muenster.de>
	22	//////////////////////////////////////////////////////
dea861a3	23
	24	//#include "TPDGCode.h"
	25	/*#include "TFile.h"
0ed6f095	26	#include "TTree.h"
dea861a3	27	#include "TEventList.h"*/
0ed6f095	28	#include "TObjArray.h"
	29	#include "TH2.h"
	30	#include "TH2F.h"
	31	#include "TH1I.h"
fa6441b5	32	#include "TCanvas.h"
dea861a3	33	// #include "TPad.h"
dea861a3	34	// #include "TLegend.h"
0ed6f095	35
	36	#include "AliLog.h"
	37	#include "AliESDtrack.h"
	38	#include "AliPID.h"
	39	#include "AliTRDv0Monitor.h"
	40	#include "info/AliTRDv0Info.h"
	41
	42
	43	ClassImp(AliTRDv0Monitor)
	44
	45	//________________________________________________________________________
	46	AliTRDv0Monitor::AliTRDv0Monitor()
	47	:AliTRDrecoTask()
dea861a3	48	,fhQualityReductions(NULL)
	49	,fV0s(NULL)
	50	,fData(NULL)
	51	,fInfo(NULL)
	52	,fP(-1.)
0ed6f095	53	{
	54	//
	55	// Default constructor
	56	//
	57	SetNameTitle("v0Monitor", "Reference V0 Monitor");
	58	}
	59
	60	//________________________________________________________________________
	61	AliTRDv0Monitor::AliTRDv0Monitor(const char name, const char title)
	62	:AliTRDrecoTask(name, title)
dea861a3	63	,fhQualityReductions(NULL)
	64	,fV0s(NULL)
	65	,fData(NULL)
	66	,fInfo(NULL)
	67	,fP(-1.)
0ed6f095	68	{
	69	//
	70	// Default constructor
	71	//
	72	DefineInput(2, TObjArray::Class()); // v0 list
	73	DefineInput(3, TObjArray::Class()); // pid info list
0ed6f095	74	}
	75
	76
dea861a3	77
fa6441b5	78	//____________________________________________________________________
	79	void AliTRDv0Monitor::MakeSummary(){
	80	TCanvas *cOut = new TCanvas("v0MonitorSummary1", "Summary 1 for task V0Monitor", 1024, 768);
	81	cOut->cd();
	82	//GetRefFigure(4);
	83	cOut->SaveAs("V0MonitorSummary.gif");
	84
	85	cOut = new TCanvas("v0MonitorSummary2","Summary 2 for task V0Monitor", 1024, 768);
	86	cOut->cd();
	87	//GetRefFigure(5);
	88	cOut->SaveAs("V0MonitorSummary2.gif");
	89	}
129a3fe2	90
	91	//________________________________________________________________________
	92	Bool_t AliTRDv0Monitor::GetRefFigure(Int_t /ifig/)
	93	{
dea861a3	94
129a3fe2	95	AliInfo("Implementation on going ...");
	96	return kTRUE;
	97	}
	98
0ed6f095	99	//________________________________________________________________________
	100	void AliTRDv0Monitor::UserCreateOutputObjects()
	101	{
	102	// Create histograms
	103	// Called once
	104
0ed6f095	105
dea861a3	106	fContainer = new TObjArray(kNPlots);
	107	fContainer->SetName("V0Monitoring");
	108
	109
	110	const char *samplename[AliPID::kSPECIES] = {"electrons","muons","pions","kaons","protons"};
	111	const char *decayname[AliTRDv0Info::kNDecays] = {"gamma","K0s","Lambda","AntiLambda"};
	112	const char *detectorname[kNDets] = {"ITS","TPC","TOF"};
	113
	114	fhQualityReductions = new TH1I(Form("fhQualityReductions"),Form("Number of tracks cut out by different quality cut steps"),11,-9,2);
	115	fContainer->Add(fhQualityReductions);
	116
	117	for(Int_t ipart = 0;ipart < AliPID::kSPECIES; ipart++){
	118	fhCutReductions[ipart] = new TH1I(Form("fhCutReductions_%s",samplename[ipart]),Form("Number of tracks cut out by different cut steps for %s",samplename[ipart]),19,-17,2);
	119	fContainer->Add(fhCutReductions[ipart]);
	120	for(Int_t idetector = 0; idetector < kNDets; idetector++){
	121	fhDetPID[idetector][ipart] = new TH2F(Form("fhDetector_%s_%s",detectorname[idetector],samplename[ipart]),Form("%s Likelihood for %s vs. momentum",detectorname[idetector], samplename[ipart]),100,0.,10.,100, 0., 1.);
	122
	123	fContainer->Add(fhDetPID[idetector][ipart]);
	124	}
	125	fhComPID[ipart] = new TH2F(Form("fhComPID_%s",samplename[ipart]),Form("Combined TPC/TOF PID: Likelihood for %s",samplename[ipart]),100,0.,10.,100,0.,1.);
	126
	127	fContainer->Add(fhComPID[ipart]);
	128
	129	for(Int_t cutstep = 0; cutstep < kNCutSteps; cutstep++){
	130	fhTPCdEdx[ipart][cutstep] = new TH2F(Form("fhTPCdEdx_%s_[%d]",samplename[ipart],cutstep),Form("TPC dE/dx for %s [%d]",samplename[ipart],cutstep),100,0.,10.,300,0.,300.);
	131
	132	fContainer->Add(fhTPCdEdx[ipart][cutstep]);
	133	}
	134	}
	135
	136	for(Int_t iDecay = 0; iDecay < AliTRDv0Info::kNDecays; iDecay++){
	137	for(Int_t cutstep =0; cutstep < kNCutSteps; cutstep++){
	138	fhV0Chi2ndf[iDecay][cutstep] = new TH2F(Form("fhV0Chi2ndf_%s_[%d]",decayname[iDecay],cutstep),Form("Chi2/NDF vs. momentum"),100,0.,10.,500, 0., 500.);
	139
	140	fContainer->Add(fhV0Chi2ndf[iDecay][cutstep]);
	141
	142	fhPsiPair[iDecay][cutstep] = new TH2F(Form("fhV0PsiPair_%s_[%d]",decayname[iDecay],cutstep),Form("Psi_pair vs. momentum"),100,0.,10.,200, 0., 1.6);
	143
	144	fContainer->Add(fhPsiPair[iDecay][cutstep]);
	145
	146	fhPointAngle[iDecay][cutstep] = new TH2F(Form("fhPointAngle_%s_[%d]",decayname[iDecay],cutstep),Form("Pointing Angle vs. momentum"),100,0.,10.,500, 0., 1.6);
	147	fContainer->Add(fhPointAngle[iDecay][cutstep]);
	148
	149	fhDCA[iDecay][cutstep] = new TH2F(Form("fhDCA_%s_[%d]",decayname[iDecay],cutstep),Form("V0 Daughter DCA vs. momentum"),100,0.,10.,500, 0., 1.);
	150
	151	fContainer->Add(fhDCA[iDecay][cutstep]);
	152
	153	fhOpenAngle[iDecay][cutstep] = new TH2F(Form("fhOpenAngle_%s_[%d]",decayname[iDecay],cutstep),Form("Opening Angle vs. momentum"),100,0.,10.,500, 0., 1.6);
0ed6f095	154
dea861a3	155	fContainer->Add(fhOpenAngle[iDecay][cutstep]);
0ed6f095	156
dea861a3	157	fhRadius[iDecay][cutstep] = new TH2F(Form("fhRadius_%s_[%d]",decayname[iDecay],cutstep),Form("V0 Generation Radius vs. momentum"),100,0.,10.,500, 0., 150.);
	158
	159	fContainer->Add(fhRadius[iDecay][cutstep]);
	160	}
	161
	162	fhInvMass[iDecay] = new TH2F(Form("fhInvMass_%s",decayname[iDecay]),Form("Invariant Mass vs. momentum"),100,0.,10.,500, 0., 2.);
	163
	164	fContainer->Add(fhInvMass[iDecay]);
	165
	166	}
	167
	168	/TH1F hV0mcPID[AliPID::kSPECIES][AliPID::kSPECIES];
	169	Int_t nPBins = 200;
	170
	171
	172
	173	for(Int_t iSpecies = 0; iSpecies < AliPID::kSPECIES; iSpecies++){
	174	for(Int_t iSample = 0; iSample < AliPID::kSPECIES; iSample++){
	175
	176	fhV0mcPID[iSample][iSpecies] = new TH1F(Form("fhV0mcPID_%s_is_%s",name[iSample],name[iSpecies]),Form("%s contained in %s sample",name[iSpecies],name[iSample]), nPBins, 0.2, 13.);
	177	}
	178	}*/
0ed6f095	179	}
	180	//________________________________________________________________________
	181	void AliTRDv0Monitor::UserExec(Option_t *)
	182	{
	183	// Main loop
	184	// Called for each event
	185
	186	if(!(fTracks = dynamic_cast<TObjArray*>(GetInputData(1)))) return;
	187	if(!(fV0s = dynamic_cast<TObjArray*>(GetInputData(2)))) return;
	188	if(!(fInfo = dynamic_cast<TObjArray*>(GetInputData(3)))) return;
	189
	190
	191	AliTRDtrackInfo *track = NULL;
	192	AliTRDv0Info *v0(NULL);
	193
dea861a3	194
0ed6f095	195
	196	for(Int_t itrk=0; itrk<fTracks->GetEntriesFast(); itrk++){
	197	track = (AliTRDtrackInfo*)fTracks->UncheckedAt(itrk);
	198	for(Int_t iv(0); iv<fV0s->GetEntriesFast(); iv++){
	199	if(!(v0 = (AliTRDv0Info*)fV0s->At(iv))) continue;
	200	if(!v0->HasTrack(track)) continue;
	201	ULong_t status = track->GetStatus();
	202	if(!(status&AliESDtrack::kTRDpid)) continue;
	203
fa6441b5	204	fhQualityReductions->Fill(v0->fQuality);//fills integer codes for tracks cut out by track/V0 quality cuts
0ed6f095	205
	206	if(!(v0->fQuality == 1))continue;
	207
	208	for(Int_t part = 0; part < AliPID::kSPECIES; part++){
dea861a3	209	fhCutReductions[part]->Fill(v0->GetPID(part,track));//fill in numbers of tracks eliminated by different PID cuts
0ed6f095	210	}
0ed6f095	211
dea861a3	212
0ed6f095	213	for(Int_t idecay(0), part(-1); idecay < AliTRDv0Info::kNDecays; idecay++){//loop over decay types considered for reference data
dea861a3	214
	215	if(idecay == AliTRDv0Info::kLambda){ //protons and pions from Lambda
	216	part = AliPID::kProton;
	217	} else if(idecay == AliTRDv0Info::kAntiLambda) { //antiprotons and pions from Anti-Lambda
	218	part = AliPID::kProton;
	219	} else if(idecay == AliTRDv0Info::kK0s) {//pions from K0s
	220	part = AliPID::kPion;
	221	} else if(idecay == AliTRDv0Info::kGamma) {//electrons from conversions
	222	part = AliPID::kElectron;
	223	}
	224
	225	//fill histograms with track/V0 quality cuts only
	226	fhPsiPair[idecay][0]->Fill(v0->fV0Momentum,v0->fPsiPair);//Angle between daughter momentum plane and plane perpendicular to magnetic field
	227	fhInvMass[idecay]->Fill(v0->fV0Momentum,v0->fInvMass[idecay]);//Invariant mass
	228	fhPointAngle[idecay][0]->Fill(v0->fV0Momentum,v0->fPointingAngle);// = TMath::ACos(esdv0->GetV0CosineOfPointingAngle()); // Cosine of pointing angle
	229	fhOpenAngle[idecay][0]->Fill(v0->fV0Momentum,v0->fOpenAngle);// opening angle between daughters
	230	fhDCA[idecay][0]->Fill(v0->fV0Momentum,v0->fDCA);// Distance of closest approach of daughter tracks
	231	fhV0Chi2ndf[idecay][0]->Fill(v0->fV0Momentum,v0->fChi2ndf[idecay]);//Kalman Filter Chi2/NDF
	232	fhRadius[idecay][0]->Fill(v0->fV0Momentum,v0->fRadius);//distance of decay/conversion from primary vertex in x-y plane
	233
	234	if(v0->HasTrack(track) == -1){
	235	fhTPCdEdx[part][0]->Fill(v0->fTrackN->P(),v0->fTPCdEdx[AliTRDv0Info::kNeg]);//TPC dE/dx for negative track
	236	} else if(v0->HasTrack(track) == 1){
	237	fhTPCdEdx[part][0]->Fill(v0->fTrackP->P(),v0->fTPCdEdx[AliTRDv0Info::kPos]);//TPC dE/dx for positive track
0ed6f095	238	}
0ed6f095	239
dea861a3	240	//fill histograms after invariant mass cuts
	241	if((v0->fInvMass[idecay] < v0->fUpInvMass[idecay][0])&&(v0->fInvMass[idecay]> v0->fDownInvMass[idecay])){
	242	fhV0Chi2ndf[idecay][1]->Fill(v0->fV0Momentum,v0->fChi2ndf[idecay]);
	243	fhPsiPair[idecay][1]->Fill(v0->fV0Momentum,v0->fPsiPair);
	244	fhPointAngle[idecay][1]->Fill(v0->fV0Momentum,v0->fPointingAngle);
	245	fhOpenAngle[idecay][1]->Fill(v0->fV0Momentum,v0->fOpenAngle);
	246	fhDCA[idecay][1]->Fill(v0->fV0Momentum,v0->fDCA);
	247	fhRadius[idecay][1]->Fill(v0->fV0Momentum,v0->fRadius);
	248	if(v0->HasTrack(track) == -1)
	249	fhTPCdEdx[part][1]->Fill(v0->fTrackN->P(),v0->fTPCdEdx[AliTRDv0Info::kNeg]);
	250	else if(v0->HasTrack(track) == 1)
	251	fhTPCdEdx[part][1]->Fill(v0->fTrackP->P(),v0->fTPCdEdx[AliTRDv0Info::kPos]);
0ed6f095	252
dea861a3	253	}
0ed6f095	254
dea861a3	255	//fill histograms after all reference selection cuts
	256	if(v0->GetPID(part,track)==1){
	257	fhV0Chi2ndf[idecay][2]->Fill(v0->fV0Momentum,v0->fChi2ndf[idecay]);
	258	fhPsiPair[idecay][2]->Fill(v0->fV0Momentum,v0->fPsiPair);
	259	fhPointAngle[idecay][2]->Fill(v0->fV0Momentum,v0->fPointingAngle);
	260	fhOpenAngle[idecay][2]->Fill(v0->fV0Momentum,v0->fOpenAngle);
	261	fhDCA[idecay][2]->Fill(v0->fV0Momentum,v0->fDCA);
	262	fhRadius[idecay][2]->Fill(v0->fV0Momentum,v0->fRadius);
	263	if(v0->HasTrack(track) == -1)
	264	fhTPCdEdx[part][2]->Fill(v0->fTrackN->P(),v0->fTPCdEdx[AliTRDv0Info::kNeg]);
	265	else if(v0->HasTrack(track) == 1)
	266	fhTPCdEdx[part][2]->Fill(v0->fTrackP->P(),v0->fTPCdEdx[AliTRDv0Info::kPos]);
0ed6f095	267
dea861a3	268	}
	269	}
	270	}
	271	}
	272	PostData(1, fContainer);
	273	}