UtilitiesAMSEvent.h

Go to the documentation of this file.

#pragma once
 
#include <root.h>
 
 
namespace AMSEventUtilities
{
    //______________________________________________________________________________________
    double FindMaxEkinTrdMCClusters(AMSEventR * pev, int ntrdmccluster, int mc0_PID)
    {
        
        double max_clusterZ = 0.0;      
        double max_clusterZEkin = -1;  
 
        for (int iTRDMCCluster = 0; iTRDMCCluster < ntrdmccluster; iTRDMCCluster++)
        {
            TrdMCClusterR ctr = pev->TrdMCCluster(iTRDMCCluster);
            
            if (ctr.Xgl[2] > max_clusterZ && ctr.GtrkID == 1 && ctr.ParticleNo == mc0_PID)
            {
                max_clusterZ = ctr.Xgl[2];
                max_clusterZEkin = ctr.Ekin;
            }
        }
 
        return max_clusterZEkin;
    }
}