ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitAna/DataTree/interface/SuperCluster.h

Comparing UserCode/MitAna/DataTree/interface/SuperCluster.h (file contents):
Revision 1.12 by loizides, Wed Apr 8 10:24:33 2009 UTC vs.
Revision 1.17 by sixie, Fri Mar 26 15:09:11 2010 UTC

#	Line 26 | Line 26 | namespace mithep
26		SuperCluster() : fEnergy(0), fEtaWidth(0), fPreshowerEnergy(0),
27		fPhiWidth(0), fRawEnergy(0) {}
28
29	<	void                    AddCluster(const BasicCluster *c)      { fClusters.Add(c);           }
30	<	const BasicCluster     *Cluster(UInt_t i)       const { return fClusters.At(i);              }
31	<	UInt_t                  ClusterSize()           const { return fClusters.Entries();          }
32	<	Double_t                Energy()                const { return fEnergy;                      }
33	<	Double_t                Eta()                   const { return fPoint.Eta();                 }
34	<	Double_t                EtaWidth()              const { return fEtaWidth;                    }
35	<	EObjType                ObjType()               const { return kSuperCluster;                }
36	<	Double_t                Phi()                   const { return fPoint.Phi();                 }
37	<	Double_t                PhiWidth()              const { return fPhiWidth;                    }
38	<	ThreeVectorC            Point()                 const { return fPoint.V();                   }
39	<	void                    Print(Option_t *opt="") const;
40	<	Double_t                PreshowerEnergy()       const { return fPreshowerEnergy;             }
41	<	Double_t                RawEnergy()             const { return fRawEnergy;                   }
42	<	Double_t                Rho()                   const { return fPoint.Rho();                 }
43	<	const BasicCluster     *Seed()                  const { return fSeedRef.Obj();               }
44	<	void                    SetEnergy(Double_t energy)                 { fEnergy = energy;       }
45	<	void                    SetEtaWidth(Double_t etaWidth)             { fEtaWidth = etaWidth;   }
46	<	void                    SetPhiWidth(Double_t phiWidth)             { fPhiWidth = phiWidth;   }
47	<	void                    SetPreshowerEnergy(Double_t e)             { fPreshowerEnergy = e;   }
48	<	void                    SetRawEnergy(Double_t rawEnergy)           { fRawEnergy = rawEnergy; }
49	<	void                    SetSeed(const BasicCluster *s)             { fSeedRef = s;           }
50	<	void                    SetXYZ(Double_t x, Double_t y, Double_t z) { fPoint.SetXYZ(x,y,z);   }
29	>	void                   AddCluster(const BasicCluster *c)          { fClusters.Add(c);        }
30	>	const BasicCluster    *Cluster(UInt_t i)       const { return fClusters.At(i);               }
31	>	UInt_t                 ClusterSize()           const { return fClusters.Entries();           }
32	>	Int_t                  Compare(const TObject *o) const;
33	>	Double_t               Energy()                const { return fEnergy;                       }
34	>	Double_t               Et()                    const;
35	>	Double_t               Eta()                   const { return fPoint.Eta();                  }
36	>	Double_t               EtaWidth()              const { return fEtaWidth;                     }
37	>	Double_t               HcalDepth1Energy()      const { return fHcalDepth1Energy;             }
38	>	Double_t               HcalDepth2Energy()      const { return fHcalDepth2Energy;             }
39	>	Double_t               HadDepth1OverEm()       const { return fHcalDepth1Energy/fEnergy;     }
40	>	Double_t               HadDepth2OverEm()       const { return fHcalDepth2Energy/fEnergy;     }
41	>	Double_t               HadOverEm()             const { return (fHcalDepth1Energy+
42	>	fHcalDepth2Energy)/fEnergy;   }
43	>	Bool_t                 IsSortable()            const { return kTRUE;                         }
44	>	EObjType               ObjType()               const { return kSuperCluster;                 }
45	>	Double_t               Phi()                   const { return fPoint.Phi();                  }
46	>	Double_t               PhiWidth()              const { return fPhiWidth;                     }
47	>	ThreeVectorC           Point()                 const { return fPoint.V();                    }
48	>	void                   Print(Option_t *opt="") const;
49	>	Double_t               PreshowerEnergy()       const { return fPreshowerEnergy;              }
50	>	Double_t               RawEnergy()             const { return fRawEnergy;                    }
51	>	Double_t               Rho()                   const { return fPoint.Rho();                  }
52	>	const BasicCluster    *Seed()                  const { return fSeedRef.Obj();                }
53	>	void                   SetEnergy(Double_t energy)                 { fEnergy = energy;        }
54	>	void                   SetEtaWidth(Double_t etaWidth)             { fEtaWidth = etaWidth;    }
55	>	void                   SetPhiWidth(Double_t phiWidth)             { fPhiWidth = phiWidth;    }
56	>	void                   SetPreshowerEnergy(Double_t e)             { fPreshowerEnergy = e;    }
57	>	void                   SetRawEnergy(Double_t rawEnergy)           { fRawEnergy = rawEnergy;  }
58	>	void                   SetHcalDepth1Energy(Double_t x)            { fHcalDepth1Energy = x;   }
59	>	void                   SetHcalDepth2Energy(Double_t x)            { fHcalDepth2Energy = x;   }
60	>	void                   SetSeed(const BasicCluster *s)             { fSeedRef = s;            }
61	>	void                   SetXYZ(Double_t x, Double_t y, Double_t z) { fPoint.SetXYZ(x,y,z);    }
62
63		protected:
64	<	Vect3C                  fPoint;           //centroid Position
65	<	Double32_t              fEnergy;          //[0,0,14]super cluster energy
66	<	Double32_t              fEtaWidth;        //[0,0,14]width in Phi
67	<	Double32_t              fPreshowerEnergy; //[0,0,14]energy in the preshower
68	<	Double32_t              fPhiWidth;        //[0,0,14]width in Phi
69	<	Double32_t              fRawEnergy;       //[0,0,14]super cluster raw energy
70	<	RefArray<BasicCluster>  fClusters;        //assigned basic clusters
71	<	Ref<BasicCluster>       fSeedRef;         //seed cluster
64	>	Vect3C                  fPoint;               //centroid Position
65	>	Double32_t              fEnergy;              //[0,0,14]super cluster energy
66	>	Double32_t              fEtaWidth;            //[0,0,14]width in Phi
67	>	Double32_t              fPreshowerEnergy;     //[0,0,14]energy in the preshower
68	>	Double32_t              fPhiWidth;            //[0,0,14]width in Phi
69	>	Double32_t              fRawEnergy;           //[0,0,14]super cluster raw energy
70	>	Double32_t              fHcalDepth1Energy;    //[0,0,14] hcal depth1 over ECAL energy
71	>	Double32_t              fHcalDepth2Energy;    //[0,0,14] hcal depth2 over ECAL energy
72	>	RefArray<BasicCluster>  fClusters;            //assigned basic clusters
73	>	Ref<BasicCluster>       fSeedRef;             //seed cluster
74
75	<	ClassDef(SuperCluster, 1) // Super cluster class
75	>	ClassDef(SuperCluster, 2) // Super cluster class
76		};
77		}
78	+
79	+	//--------------------------------------------------------------------------------------------------
80	+	inline Double_t mithep::SuperCluster::Et() const
81	+	{
82	+	// Return transverse energy.
83	+
84	+	return fEnergy*fPoint.Rho()/fPoint.V().R();
85	+	}
86	+
87	+	//--------------------------------------------------------------------------------------------------
88	+	inline Int_t mithep::SuperCluster::Compare(const TObject *o) const
89	+	{
90	+	// Default compare function for sorting according to transverse momentum.
91	+	// Returns -1 if this object is smaller than given object, 0 if objects are
92	+	// equal and 1 if this is larger than given object.
93	+
94	+	const mithep::SuperCluster *s = dynamic_cast<const mithep::SuperCluster *>(o);
95	+	if (!s)
96	+	return 1;
97	+
98	+	Double_t mye = Energy();
99	+	Double_t e   = s->Energy();
100	+	if (mye>e)
101	+	return -1;
102	+	else if (e>mye)
103	+	return +1;
104	+	return 0;
105	+	}
106		#endif

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines