#include <iostream>#include <cassert>#include "TParticle.h"#include "DataUtil/TruthHelper.h"#include "MessageService/MsgService.h"#include "JobControl/JobCModuleRegistry.h"#include "JobControl/JobCommand.h"#include "MinosObjectMap/MomNavigator.h"#include "TClonesArray.h"#include "MCNtuple/Module/NtpMCModule.h"#include "MCNtuple/NtpMCRecord.h"#include "MCNtuple/NtpMCTruth.h"#include "MCNtuple/NtpMCStdHep.h"#include "MCNtuple/NtpMCDigiScintHit.h"#include "MCNtuple/NtpMCDetSimResult.h"#include "MCNtuple/NtpMCPhotonResult.h"#include "MCNtuple/NtpMCFluxInfo.h"#include "MCNtuple/NtpMCFluxWgt.h"#include "MCNtuple/NtpMCStdHepHit.h"#include "StandardNtuple/NtpStRecord.h"#include "Record/SimSnarlRecord.h"#include "Record/SimSnarlHeader.h"#include "Plex/PlexPlaneId.h"#include "Plex/PlexStripEndId.h"#include "REROOT_Classes/REROOT_NeuKin.h"#include "REROOT_Classes/REROOT_FluxInfo.h"#include "REROOT_Classes/REROOT_FluxWgt.h"#include "Conventions/Munits.h"#include "Conventions/Mphysical.h"#include "Util/UtilMath.h"#include "UgliGeometry/UgliGeomHandle.h"#include "UgliGeometry/UgliStripHandle.h"#include "Digitization/DigiScintHit.h"#include "RawData/RawRecord.h"#include "RawData/RawDigitDataBlock.h"#include "RawData/RawDataBlock.h"#include "RawData/RawDigit.h"#include "DataUtil/Truthifier.h"#include "PhotonTransport/PhotonEventResult.h"#include "DetSim/SimEventResult.h"Go to the source code of this file.
Functions | |
| CVSID ("$Id: NtpMCModule.cxx,v 1.41 2009/05/21 18:29:36 rhatcher Exp $") | |
| JOBMODULE (NtpMCModule,"NtpMCModule","A module for filling MC ntuple records.") | |
| bool | compareP4 (Float_t *p4neu, TParticle *part) |
|
||||||||||||
|
Definition at line 67 of file NtpMCModule.cxx. References UtilMath::sameValue(). Referenced by NtpMCModule::FillNtpMCTruth(). 00068 {
00069 bool result = true;
00070 Float_t part_p4[4] = { part->Px(), part->Py(), part->Pz(), part->Energy() };
00071 for (int k=0; k<4; ++k) {
00072 // use float strictness scaled by 100 (vs default 10)
00073 // because for fava values are exchanged as ASCII messages
00074 // and more precision gets lost
00075 result &= UtilMath::sameValue(p4neu[k],part_p4[k],false,100.);
00076 }
00077 return result;
00078 }
|
|
||||||||||||
|
|
|
||||||||||||||||
|
|
1.3.9.1