void KVReconstructedEvent::IdentifyEvent()
{
//All particles which have not been previously identified (IsIdentified=kFALSE), and which
//may be identified independently of all other particles in their group according to the 1st
//order coherency analysis (KVReconstructedNucleus::GetStatus=0), will be identified.
//Particles stopping in first member of a telescope (KVReconstructedNucleus::GetStatus=3) will
//have their Z estimated from the energy loss in the detector (if calibrated).
KVReconstructedNucleus *d;
while ((d = GetNextParticle())) {
if (!d->IsIdentified()){
if(d->GetStatus() == KVReconstructedNucleus::kStatusOK){
// identifiable particles
d->Identify();
}
else if(d->GetStatus() == KVReconstructedNucleus::kStatusStopFirstStage) {
// particles stopped in first member of a telescope
// estimation of Z (minimum) from energy loss (if detector is calibrated)
UInt_t zmin = d->GetStoppingDetector()->FindZmin(-1., d->GetMassFormula());
if( zmin ){
d->SetZ( zmin );
d->SetIsIdentified();
// "Identifying" telescope is taken from list of ID telescopes
// to which stopping detector belongs
d->SetIdentifyingTelescope( (KVIDTelescope*)d->GetStoppingDetector()->GetIDTelescopes()->At(0) );
}
}
}
}
}
bool KVINDRABackwardGroupReconstructor::DoCoherencyAnalysis(KVReconstructedNucleus& PART)
{
// Coherency analysis for backward rings 10-17 of INDRA
PART.SetParameter("UseFullChIoEnergyForCalib", !(theChio && theChio->GetNHits() > 1));
bool ok = false;
if (PART.GetStoppingDetector()->IsType("CSI")) {
// particles stopping in CsI detectors
// check coherency of CsI-R/L and Si-CsI identifications
ok = CoherencyChIoCsI(PART);
}
else {
// particle stopped in ChIo (=> Zmin)
ok = PART.IsIdentified() && identifying_telescope;
}
return ok;
}
Bool_t KVINDRAEtalonGroupReconstructor::CoherencyEtalons(KVReconstructedNucleus& PART)
{
// Called by Identify() for particles stopping in etalon modules of Rings 10-17.
KVIdentificationResult* IDcsi, *IDsilicsi, *IDsi75sili, *IDcisi75, *IDcicsi;
IDcsi = IDsilicsi = IDsi75sili = IDcisi75 = IDcicsi = nullptr;
if (PART.GetStoppingDetector()->IsType("CSI")) {
IDcsi = PART.GetIdentificationResult("CSI_R_L");
IDcicsi = PART.GetIdentificationResult("CI_CSI");
IDsilicsi = PART.GetIdentificationResult("SILI_CSI");
}
IDsi75sili = PART.GetIdentificationResult("SI75_SILI");
IDcisi75 = PART.GetIdentificationResult("CI_SI75");
PART.SetParameter("PileupChIo", kFALSE);
PART.SetParameter("IncludeEtalonsInCalibration", kTRUE);
Bool_t haveCsI = IDcsi && IDcsi->IDOK;
Bool_t haveSiLiCsI = IDsilicsi && IDsilicsi->IDOK;
Bool_t haveSi75SiLi = IDsi75sili && IDsi75sili->IDOK;
Bool_t haveChIoSi75 = IDcisi75 && IDcisi75->IDOK;
KVIDTelescope* idt_csi, *idt_silicsi, *idt_si75sili, *idt_cisi75, *idt_cicsi;
idt_csi = idt_silicsi = idt_si75sili = idt_cisi75 = idt_cicsi = nullptr;
if (PART.GetStoppingDetector()->IsType("CSI")) {
idt_csi =
(KVIDTelescope*)PART.GetReconstructionTrajectory()->GetIDTelescopes()->FindObjectByType(IDcsi->GetType());
idt_silicsi =
(KVIDTelescope*)PART.GetReconstructionTrajectory()->GetIDTelescopes()->FindObjectByType(IDsilicsi->GetType());
idt_cicsi =
(KVIDTelescope*)PART.GetReconstructionTrajectory()->GetIDTelescopes()->FindObjectByType(IDcicsi->GetType());
}
idt_si75sili = (KVIDTelescope*)PART.GetReconstructionTrajectory()->GetIDTelescopes()->FindObjectByType(IDsi75sili->GetType());
idt_cisi75 = (KVIDTelescope*)PART.GetReconstructionTrajectory()->GetIDTelescopes()->FindObjectByType(IDcisi75->GetType());
// Treat cases where particle hit etalon telescope
if (idt_csi) {
if (haveCsI) {
partID = *IDcsi;
identifying_telescope = idt_csi;
if (haveSi75SiLi) {
// check for heavy fragment in Si75-SiLi
if (IDsi75sili->Z > partID.Z) {
if (haveChIoSi75) {
// check we don't have a better identification in ChIo-Si75
if (IDcisi75->IDquality < IDsi75sili->IDquality && IDcisi75->Z > partID.Z) {
PART.SetParameter("PileupSi75", kTRUE);
IDcisi75->SetComment("CsI identification with another particle stopped in Si75");
PART.SetParameter("UseFullChIoEnergyForCalib", kFALSE); // calculate ChIo energy for this particle
}
else {
PART.SetParameter("PileupSiLi", kTRUE);
IDsi75sili->SetComment("CsI identification with another particle stopped in SiLi");
PART.SetParameter("UseFullChIoEnergyForCalib", kFALSE); // calculate ChIo energy for this particle
}
}
}
}
else if (haveChIoSi75) {
// check for heavy fragment in ChIo-Si75
if (IDcisi75->Z > partID.Z) {
PART.SetParameter("PileupSi75", kTRUE);
IDcisi75->SetComment("CsI identification with another particle stopped in Si75");
PART.SetParameter("UseFullChIoEnergyForCalib", kFALSE); // calculate ChIo energy for this particle
}
}
return kTRUE;
}
else if (haveSiLiCsI) {
partID = *IDsilicsi;
identifying_telescope = idt_silicsi;
if (haveChIoSi75) {
// check for heavy fragment in ChIo-Si75
if (IDcisi75->Z > partID.Z) {
PART.SetParameter("PileupSi75", kTRUE);
IDcisi75->SetComment("CsI identification with another particle stopped in Si75");
PART.SetParameter("UseFullChIoEnergyForCalib", kFALSE); // calculate ChIo energy for this particle
}
}
return kTRUE;
}
}
else if (PART.GetStoppingDetector()->IsType("SILI")) {
if (haveSi75SiLi) {
partID = *IDsi75sili;
identifying_telescope = idt_si75sili;
// check ChIo-Si75 id is coherent (either no id or Z<=this one)
if (haveChIoSi75) {
if (IDcisi75->Z > partID.Z) PART.SetParameter("PileupChIo", kTRUE);
}
return kTRUE;
}
else if (haveChIoSi75) {
partID = *IDcisi75;
identifying_telescope = idt_cisi75;
return kTRUE;
}
}
else if (PART.GetStoppingDetector()->IsType("SI75")) { // MFR march 2014 SiLi is not the real stopping detector
if (haveChIoSi75) {
partID = *IDcisi75;
identifying_telescope = idt_cisi75;
return kTRUE;
} // end MFR march 2014
}
return kFALSE;
}