21 #include "Utilities/func/MathUtil.h" 30 #include "NovaDAQConventions/DAQConventions.h" 106 RC_Tree = tfs->
make<TTree>(
"RC_Tree",
"Information from RC track info");
143 if (p->
PdgCode()!=42 )
continue;
145 unsigned ntruehits_Tot = flshits.size();
147 if ( ntruehits_Tot==0 )
continue;
148 double momentum2 = (p->
P())*(p->
P());
149 double mass2 = (p->
Mass())*(p->
Mass());
153 for (
unsigned h = 0;
h!= flshits.size(); ++
h) {
164 const int _xMax = 383;
165 const int _yMax = 383;
166 const int _zMax = 895;
167 const int _xdelt = 35;
168 const int _ydelt = 35;
169 const int _zdelt = 15;
171 int PX_max, PY_max, PX_min, PY_min;
172 double TX_max, TY_max, TX_min, TY_min;
193 for (
unsigned sliceIdx = 0; sliceIdx != slices->size(); ++sliceIdx) {
213 bool Intrusion =
false;
215 for (
unsigned i = 0;
i != slice.
NCell(); ++
i) {
217 if (chit->
ADC() < 501)
continue;
220 float tns = chit->
TNS();
229 if (PX_max < chit->Plane() ) PX_max = chit->
Plane();
230 if (PX_min > chit->
Plane() ) PX_min = chit->
Plane();
231 if (TX_max < tns ) TX_max = tns;
232 if (TX_min > tns ) TX_min = tns;
235 if (chit->
Cell() < _xdelt){
239 else if (chit->
Cell() > _xMax - _xdelt){
246 else if (chit->
Cell() > _xMax-_xdelt && entry_id!=2)
_penetrated =
true;
252 if (PY_max < chit->Plane() ) PY_max = chit->
Plane();
253 if (PY_min > chit->
Plane() ) PY_min = chit->
Plane();
254 if (TY_max < chit->TNS() ) TY_max = chit->
TNS();
255 if (TY_min > chit->
TNS() ) TY_min = chit->
TNS();
258 if (chit->
Cell() < _ydelt) {
262 else if (chit->
Cell() > _yMax - _ydelt) {
270 else if (chit->
Cell() > _yMax-_ydelt && entry_id!=4)
_penetrated =
true;
275 if (chit->
Plane() < _zdelt) {
279 else if (chit->
Plane() > _zMax - _zdelt) {
T max(const caf::Proxy< T > &a, T b)
SubRunNumber_t subRun() const
back track the reconstruction to the simulation
unsigned int NCell(geo::View_t view) const
Number of cells in view view.
const sim::ParticleNavigator & ParticleNavigator() const
Get a reference to the ParticleNavigator.
unsigned short Plane() const
std::vector< sim::FLSHit > ParticleToFLSHit(const int &trackID) const
All the FLSHits that were created by the track id trackID, sorted from most to least light...
list_type::const_iterator const_iterator
Vertical planes which measure X.
A collection of associated CellHits.
DEFINE_ART_MODULE(TestTMapFile)
unsigned short Cell() const
double P(const int i=0) const
std::string fClusterInput
Input folder from cluster reco.
Collect Geo headers and supply basic geometry functions.
EDAnalyzer(Table< Config > const &config)
FmmTriggerAna(fhicl::ParameterSet const &pset)
EventNumber_t event() const
art::Ptr< rb::CellHit > Cell(geo::View_t view, unsigned int viewIdx) const
Get the ith cell from view view.
T * make(ARGS...args) const
int16_t ADC(uint32_t i) const
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
void analyze(const art::Event &evt)
T min(const caf::Proxy< T > &a, T b)
Encapsulate the geometry of one entire detector (near, far, ndos)