Commit ec86ba3e authored by Jérémie Dudouet's avatar Jérémie Dudouet Committed by Jérémie Dudouet
Browse files

add FOM, second interaction position point and interaction type (photo,compton,pair) in the Tree

parent 3ca1bd03
......@@ -477,6 +477,11 @@ void TB_AGATA_Tracking::InitTree()
fTree->Branch("trackX1", trackX1, "trackX1[nbTrack]/F");
fTree->Branch("trackY1", trackY1, "trackY1[nbTrack]/F");
fTree->Branch("trackZ1", trackZ1, "trackZ1[nbTrack]/F");
fTree->Branch("trackX2", trackX2, "trackX2[nbTrack]/F");
fTree->Branch("trackY2", trackY2, "trackY2[nbTrack]/F");
fTree->Branch("trackZ2", trackZ2, "trackZ2[nbTrack]/F");
fTree->Branch("trackType", trackType, "trackType[nbTrack]/I");
fTree->Branch("trackFOM", trackFOM, "trackFOM[nbTrack]/F");
fTree->Branch("trackT", trackT, "trackT[nbTrack]/F");
fTree->Branch("trackCrystalID", trackCrystalID, "trackCrystalID[nbTrack]/I");
fTree->Branch("TStrack", &TSTrack, "TStrack/l");
......@@ -525,6 +530,9 @@ void TB_AGATA_Tracking::Process(Int_t idet)
trackE[i] = gamma1->GetE();
trackType[i] = gamma1->GetType();
trackFOM[i] = gamma1->GetFOM();
trackT[i] = hit1->GetT();
trackCrystalID[i] = hit1->GetID(1);
......@@ -532,6 +540,19 @@ void TB_AGATA_Tracking::Process(Int_t idet)
trackY1[i] = gamma1->GetY();
trackZ1[i] = gamma1->GetZ();
//if tracked gamma is photo or pair, the second int point is a copy of the first one (because not useful)
if(trackType[i]==2){
trackX2[i] = gamma1->GetX();
trackY2[i] = gamma1->GetY();
trackZ2[i] = gamma1->GetZ();
}
else{
trackX2[i] = trackX1[i];
trackY2[i] = trackY1[i];
trackZ2[i] = trackZ1[i];
}
if(trackE[i]==0)
nbTrack--;
}
......
......@@ -115,6 +115,12 @@ protected:
Float_t trackX1[MaxTrackedGamma]; // X position of the first interaction point
Float_t trackY1[MaxTrackedGamma]; // Y position of the first interaction point
Float_t trackZ1[MaxTrackedGamma]; // Z position of the first interaction point
Float_t trackX2[MaxTrackedGamma]; // X position of the Second interaction point
Float_t trackY2[MaxTrackedGamma]; // Y position of the Second interaction point
Float_t trackZ2[MaxTrackedGamma]; // Z position of the Second interaction point
Int_t trackType[MaxTrackedGamma]; // Interaction type (1-photo 2-compt 3-pair)
Float_t trackFOM[MaxTrackedGamma]; // figure of merit
Float_t trackT[MaxTrackedGamma]; // time of the track gamma
Int_t trackCrystalID[MaxTrackedGamma]; // CrystalID of the first interaction point of the track
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment