theAnalysis->FillList(mcpiplus,"PionLoosePlus","PidAlgoIdealCharged");

pimatchplus.Cleanup();
pimatchminus.Cleanup();

		theAnalysis->FillList(mcpiplus,"PionLoosePlus","PidAlgoIdealCharged");
		theAnalysis->FillList(mcpiminus,"PionLooseMinus","PidAlgoIdealCharged");
theAnalysis->FillList(piall,"PionAll");


for ( int ii=0 ; ii<piall.GetLength() ; ++ii ){
if ( theAnalysis->McTruthMatch(piall[ii]) ){
if ( 211 == piall[ii]->PdgCode()  ){
pimatchplus.Add( piall[ii] );
}
if ( -211 == piall[ii]->PdgCode()  ){
pimatchminus.Add( piall[ii] );
}

}
}

// consistency check (show suspicious events)
if ( ( pimatchplus.GetLength()!= mcpiplus.GetLength() ) || ( pimatchminus.GetLength() != mcpiminus.GetLength() ) ){
cout << "evt " << i << endl;
cout << "pimatchplus = " << pimatchplus.GetLength() << "     piplus = " << mcpiplus.GetLength() << endl;
cout << "pimatchminus = " << pimatchminus.GetLength() << "     piminus = " << mcpiminus.GetLength() << endl;
			++nsuspiciousevents;
}

		npimatchplus += pimatchplus.GetLength();
		npimatchminus += pimatchminus.GetLength();

		nmcpiplus += mcpiplus.GetLength();
		nmcpiminus += mcpiminus.GetLength();