V0Validator Class Reference

#include <Validation/RecoVertex/src/V0Validator.cc>

Inheritance diagram for V0Validator:

Public Member Functions
	V0Validator (const edm::ParameterSet &)
	~V0Validator ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginRun (const edm::Run &, const edm::EventSetup &)
virtual void	endRun (const edm::Run &, const edm::EventSetup &)

Private Attributes
std::string	dirName
edm::EDGetTokenT < edm::SimTrackContainer >	edmSimTrackContainerToken_
edm::EDGetTokenT < edm::SimVertexContainer >	edmSimVertexContainerToken_
edm::EDGetTokenT< edm::View < reco::Track > >	edmView_recoTrack_Token_
MonitorElement *	fakeKsMass
MonitorElement *	fakeLamMass
int	genK0s
int	genLam
MonitorElement *	goodKsMass
MonitorElement *	goodLamMass
float	K0sCandEta
int	k0sCandFound
float	K0sCandpT
float	K0sCandR
unsigned int	K0sCandStatus
float	K0sGenEta
float	K0sGenpT
float	K0sGenR
unsigned int	K0sGenStatus
unsigned int	K0sPiCandStatus [2]
unsigned int	K0sPiEff [2]
int	k0sTracksFound
MonitorElement *	ksAbsoluteDistResolution
MonitorElement *	ksCandStatus
MonitorElement *	ksEffVsEta
MonitorElement *	ksEffVsEta_denom
MonitorElement *	ksEffVsEta_num
MonitorElement *	ksEffVsPt
MonitorElement *	ksEffVsPt_denom
MonitorElement *	ksEffVsPt_num
MonitorElement *	ksEffVsR
MonitorElement *	ksEffVsR_denom
MonitorElement *	ksEffVsR_num
MonitorElement *	ksFakeDauRadDist
MonitorElement *	ksFakeVsEta
MonitorElement *	ksFakeVsEta_denom
MonitorElement *	ksFakeVsEta_num
MonitorElement *	ksFakeVsPt
MonitorElement *	ksFakeVsPt_denom
MonitorElement *	ksFakeVsPt_num
MonitorElement *	ksFakeVsR
MonitorElement *	ksFakeVsR_denom
MonitorElement *	ksFakeVsR_num
float	KsGenX
float	KsGenY
float	KsGenZ
MonitorElement *	ksMassAll
MonitorElement *	ksTkEffVsEta
MonitorElement *	ksTkEffVsEta_num
MonitorElement *	ksTkEffVsPt
MonitorElement *	ksTkEffVsPt_num
MonitorElement *	ksTkEffVsR
MonitorElement *	ksTkEffVsR_num
MonitorElement *	ksTkFakeVsEta
MonitorElement *	ksTkFakeVsEta_num
MonitorElement *	ksTkFakeVsPt
MonitorElement *	ksTkFakeVsPt_num
MonitorElement *	ksTkFakeVsR
MonitorElement *	ksTkFakeVsR_num
MonitorElement *	ksXResolution
MonitorElement *	ksYResolution
MonitorElement *	ksZResolution
MonitorElement *	lamAbsoluteDistResolution
float	LamCandEta
int	lamCandFound
float	LamCandpT
float	LamCandR
unsigned int	LamCandStatus
MonitorElement *	lamCandStatus
MonitorElement *	lamEffVsEta
MonitorElement *	lamEffVsEta_denom
MonitorElement *	lamEffVsEta_num
MonitorElement *	lamEffVsPt
MonitorElement *	lamEffVsPt_denom
MonitorElement *	lamEffVsPt_num
MonitorElement *	lamEffVsR
MonitorElement *	lamEffVsR_denom
MonitorElement *	lamEffVsR_num
MonitorElement *	lamFakeDauRadDist
MonitorElement *	lamFakeVsEta
MonitorElement *	lamFakeVsEta_denom
MonitorElement *	lamFakeVsEta_num
MonitorElement *	lamFakeVsPt
MonitorElement *	lamFakeVsPt_denom
MonitorElement *	lamFakeVsPt_num
MonitorElement *	lamFakeVsR
MonitorElement *	lamFakeVsR_denom
MonitorElement *	lamFakeVsR_num
float	LamGenEta
float	LamGenpT
float	LamGenR
unsigned int	LamGenStatus
float	LamGenX
float	LamGenY
float	LamGenZ
MonitorElement *	lamMassAll
unsigned int	LamPiCandStatus [2]
unsigned int	LamPiEff [2]
MonitorElement *	lamTkEffVsEta
MonitorElement *	lamTkEffVsEta_num
MonitorElement *	lamTkEffVsPt
MonitorElement *	lamTkEffVsPt_num
MonitorElement *	lamTkEffVsR
MonitorElement *	lamTkEffVsR_num
MonitorElement *	lamTkFakeVsEta
MonitorElement *	lamTkFakeVsEta_num
MonitorElement *	lamTkFakeVsPt
MonitorElement *	lamTkFakeVsPt_num
MonitorElement *	lamTkFakeVsR
MonitorElement *	lamTkFakeVsR_num
int	lamTracksFound
MonitorElement *	lamXResolution
MonitorElement *	lamYResolution
MonitorElement *	lamZResolution
MonitorElement *	nKs
MonitorElement *	nLam
int	noTPforK0sCand
int	noTPforLamCand
int	realK0sFound
int	realK0sFoundEff
int	realLamFound
int	realLamFoundEff
edm::EDGetTokenT < reco::RecoToSimCollection >	recoRecoToSimCollectionToken_
edm::EDGetTokenT < reco::SimToRecoCollection >	recoSimToRecoCollectionToken_
edm::EDGetTokenT < reco::VertexCompositeCandidateCollection >	recoVertexCompositeCandidateCollection_k0s_Token_
edm::EDGetTokenT < reco::VertexCompositeCandidateCollection >	recoVertexCompositeCandidateCollection_lambda_Token_
std::string	theDQMRootFileName
DQMStore *	theDQMstore
edm::EDGetTokenT < TrackingParticleCollection >	trackingParticleCollection_Eff_Token_
edm::EDGetTokenT < TrackingParticleCollection >	trackingParticleCollectionToken_
edm::EDGetTokenT< std::vector < reco::Vertex > >	vec_recoVertex_Token_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: Creates validation histograms for RecoVertex/V0Producer

Implementation: <Notes on="" implementation>="">

Definition at line 87 of file V0Validator.h.

Constructor & Destructor Documentation

V0Validator::V0Validator ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 31 of file V0Validator.cc.

References genK0s, genLam, k0sCandFound, lamCandFound, noTPforK0sCand, noTPforLamCand, cppFunctionSkipper::operator, realK0sFound, realK0sFoundEff, realLamFound, realLamFoundEff, and theDQMstore.

  : theDQMRootFileName(iConfig.getParameter<std::string>("DQMRootFileName"))
  , dirName(iConfig.getParameter<std::string>("dirName"))
  , recoRecoToSimCollectionToken_( consumes<reco::RecoToSimCollection>( edm::InputTag( std::string( "trackingParticleRecoTrackAsssociation" ) ) ) )
  , recoSimToRecoCollectionToken_( consumes<reco::SimToRecoCollection>( edm::InputTag( std::string( "trackingParticleRecoTrackAsssociation" ) ) ) )
  , trackingParticleCollection_Eff_Token_( consumes<TrackingParticleCollection>( edm::InputTag( std::string( "mix" )
                                                  , std::string( "MergedTrackTruth" )
                                                )
                                         )
                       )
  , trackingParticleCollectionToken_( consumes<TrackingParticleCollection>( edm::InputTag( std::string( "mix" )
                                             , std::string( "MergedTrackTruth" )
                                               )
                                        )
                      )
  , edmView_recoTrack_Token_( consumes< edm::View<reco::Track> >( edm::InputTag( std::string( "generalTracks" ) ) ) )
  , edmSimTrackContainerToken_( consumes<edm::SimTrackContainer>( edm::InputTag( std::string( "g4SimHits" ) ) ) )
  , edmSimVertexContainerToken_( consumes<edm::SimVertexContainer>( edm::InputTag( std::string( "g4SimHits" ) ) ) )
  , vec_recoVertex_Token_( consumes< std::vector<reco::Vertex> >( edm::InputTag( std::string( "offlinePrimaryVertices" ) ) ) )
  , recoVertexCompositeCandidateCollection_k0s_Token_( consumes<reco::VertexCompositeCandidateCollection>( iConfig.getParameter<edm::InputTag>( "kShortCollection" ) ) )
  , recoVertexCompositeCandidateCollection_lambda_Token_( consumes<reco::VertexCompositeCandidateCollection>( iConfig.getParameter<edm::InputTag>( "lambdaCollection" ) ) ) {
  genLam = genK0s = realLamFoundEff = realK0sFoundEff = lamCandFound = 
    k0sCandFound = noTPforK0sCand = noTPforLamCand = realK0sFound = realLamFound = 0;
  theDQMstore = edm::Service<DQMStore>().operator->();
}

V0Validator::~V0Validator

(

)

Definition at line 58 of file V0Validator.cc.

                          {

}

Member Function Documentation

void V0Validator::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 462 of file V0Validator.cc.

References funct::abs(), gather_cfg::cout, edmSimTrackContainerToken_, edmSimVertexContainerToken_, edmView_recoTrack_Token_, fakeKsMass, fakeLamMass, MonitorElement::Fill(), genK0s, genLam, edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), goodKsMass, goodLamMass, i, edm::Ref< C, T, F >::isNonnull(), j, K0sCandEta, K0sCandpT, K0sCandR, K0sCandStatus, K0sGenEta, K0sGenpT, K0sGenR, K0sGenStatus, K0sPiCandStatus, K0sPiEff, k0sTracksFound, ksCandStatus, ksEffVsEta_denom, ksEffVsEta_num, ksEffVsPt_denom, ksEffVsPt_num, ksEffVsR_denom, ksEffVsR_num, ksFakeDauRadDist, ksFakeVsEta_denom, ksFakeVsEta_num, ksFakeVsPt_denom, ksFakeVsPt_num, ksFakeVsR_denom, ksFakeVsR_num, ksMassAll, ksTkEffVsEta_num, ksTkEffVsPt_num, ksTkEffVsR_num, ksTkFakeVsEta_num, ksTkFakeVsPt_num, ksTkFakeVsR_num, LamCandEta, LamCandpT, LamCandR, LamCandStatus, lamCandStatus, lamEffVsEta_denom, lamEffVsEta_num, lamEffVsPt_denom, lamEffVsPt_num, lamEffVsR_denom, lamEffVsR_num, lamFakeDauRadDist, lamFakeVsEta_denom, lamFakeVsEta_num, lamFakeVsPt_denom, lamFakeVsPt_num, lamFakeVsR_denom, lamFakeVsR_num, LamGenEta, LamGenpT, LamGenR, LamGenStatus, lamMassAll, LamPiCandStatus, LamPiEff, lamTkEffVsEta_num, lamTkEffVsPt_num, lamTkEffVsR_num, lamTkFakeVsEta_num, lamTkFakeVsPt_num, lamTkFakeVsR_num, TrackingParticle::momentum(), nKs, nLam, noTPforK0sCand, noTPforLamCand, TrackingParticle::parentVertex(), TrackingParticle::pdgId(), realK0sFound, realK0sFoundEff, realLamFound, realLamFoundEff, recoRecoToSimCollectionToken_, recoSimToRecoCollectionToken_, recoVertexCompositeCandidateCollection_k0s_Token_, recoVertexCompositeCandidateCollection_lambda_Token_, mathSSE::sqrt(), trackingParticleCollection_Eff_Token_, trackingParticleCollectionToken_, vec_recoVertex_Token_, and TrackingParticle::vertex().

                                                                             {

  using std::cout;
  using std::endl;
  using namespace edm;
  using namespace std;

  //cout << "In analyze(), getting collections..." << endl;
  // Get event setup info, B-field and tracker geometry
  ESHandle<MagneticField> bFieldHandle;
  iSetup.get<IdealMagneticFieldRecord>().get(bFieldHandle);
  ESHandle<GlobalTrackingGeometry> globTkGeomHandle;
  iSetup.get<GlobalTrackingGeometryRecord>().get(globTkGeomHandle);

  // Make matching collections
  //reco::RecoToSimCollection recSimColl;
  //reco::SimToRecoCollection simRecColl;
   
  Handle<reco::RecoToSimCollection > recotosimCollectionH;
  iEvent.getByToken( recoRecoToSimCollectionToken_, recotosimCollectionH );
  //recSimColl= *( recotosimCollectionH.product() ); 
  
  Handle<reco::SimToRecoCollection> simtorecoCollectionH;
  iEvent.getByToken( recoSimToRecoCollectionToken_, simtorecoCollectionH );
  //simRecColl= *( simtorecoCollectionH.product() );

  edm::Handle<TrackingParticleCollection>  TPCollectionEff ;
  iEvent.getByToken( trackingParticleCollection_Eff_Token_, TPCollectionEff );
  const TrackingParticleCollection tPCeff = *( TPCollectionEff.product() );

  edm::ESHandle<TrackAssociatorBase> associatorByHits;
  iSetup.get<TrackAssociatorRecord>().get("TrackAssociatorByHits", associatorByHits);

  //VertexAssociatorBase* associatorByTracks;

  //  edm::ESHandle<VertexAssociatorBase> theTracksAssociator;
  //  iSetup.get<VertexAssociatorRecord>().get("VertexAssociatorByTracks",theTracksAssociator);
  //  associatorByTracks = (VertexAssociatorBase *) theTracksAssociator.product();

  // Get tracks
  Handle< View<reco::Track> > trackCollectionH;
  iEvent.getByToken( edmView_recoTrack_Token_, trackCollectionH );

  Handle<SimTrackContainer> simTrackCollection;
  iEvent.getByToken( edmSimTrackContainerToken_, simTrackCollection );
  const SimTrackContainer simTC = *(simTrackCollection.product());

  Handle<SimVertexContainer> simVertexCollection;
  iEvent.getByToken( edmSimVertexContainerToken_, simVertexCollection );
  const SimVertexContainer simVC = *(simVertexCollection.product());

  //Get tracking particles
  //  -->tracks
  edm::Handle<TrackingParticleCollection>  TPCollectionH ;
  iEvent.getByToken( trackingParticleCollectionToken_, TPCollectionH );
  const View<reco::Track>  tC = *( trackCollectionH.product() );

  // Select the primary vertex, create a new reco::Vertex to hold it
  edm::Handle< std::vector<reco::Vertex> > primaryVtxCollectionH;
  iEvent.getByToken( vec_recoVertex_Token_, primaryVtxCollectionH );
  const reco::VertexCollection primaryVertexCollection   = *(primaryVtxCollectionH.product());

  reco::Vertex* thePrimary = 0;
  std::vector<reco::Vertex>::const_iterator iVtxPH = primaryVtxCollectionH->begin();
  for(std::vector<reco::Vertex>::const_iterator iVtx = primaryVtxCollectionH->begin();
      iVtx < primaryVtxCollectionH->end();
      iVtx++) {
    if(primaryVtxCollectionH->size() > 1) {
      if(iVtx->tracksSize() > iVtxPH->tracksSize()) {
    iVtxPH = iVtx;
      }
    }
    else iVtxPH = iVtx;
  }
  thePrimary = new reco::Vertex(*iVtxPH);

  //cout << "Done with collections, associating reco and sim..." << endl;
 

  //reco::RecoToSimCollection r2s = associatorByHits->associateRecoToSim(trackCollectionH,TPCollectionH,&iEvent );
  //reco::SimToRecoCollection s2r = associatorByHits->associateSimToReco(trackCollectionH,TPCollectionH,&iEvent );

//  reco::VertexRecoToSimCollection vr2s = associatorByTracks->associateRecoToSim(primaryVtxCollectionH, TVCollectionH, iEvent, r2s);
//  reco::VertexSimToRecoCollection vs2r = associatorByTracks->associateSimToReco(primaryVtxCollectionH, TVCollectionH, iEvent, s2r);

  //get the V0s;   
  edm::Handle<reco::VertexCompositeCandidateCollection> k0sCollection;
  edm::Handle<reco::VertexCompositeCandidateCollection> lambdaCollection;
  iEvent.getByToken( recoVertexCompositeCandidateCollection_k0s_Token_, k0sCollection );
  iEvent.getByToken( recoVertexCompositeCandidateCollection_lambda_Token_, lambdaCollection );

  //make vector of pair of trackingParticles to hold good V0 candidates
  std::vector< pair<TrackingParticleRef, TrackingParticleRef> > trueK0s;
  std::vector< pair<TrackingParticleRef, TrackingParticleRef> > trueLams;
  std::vector<double> trueKsMasses;
  std::vector<double> trueLamMasses;

  // Do vertex calculations //
/*
  if( k0sCollection->size() > 0 ) {
    for(reco::VertexCompositeCandidateCollection::const_iterator iK0s = k0sCollection->begin();
    iK0s != k0sCollection->end();
    iK0s++) {
      // Still can't actually associate the V0 vertex with a TrackingVertexCollection.
      //  Is this a problem?  You bet.
      reco::VertexCompositeCandidate::CovarianceMatrix aErr;
      iK0s->fillVertexCovariance(aErr);
      reco::Vertex tVtx(iK0s->vertex(), aErr);
      reco::VertexCollection *tVtxColl = 0;
      tVtxColl->push_back(tVtx);
      reco::VertexRef aVtx(tVtxColl, 0);
      //if(vr2s.find(iK0s->vertex()) != vr2s.end()) {
      if(vr2s.find(aVtx) != vr2s.end()) {
    //cout << "Found it in the collection." << endl;
        std::vector< std::pair<TrackingVertexRef, double> > vVR 
      = (std::vector< std::pair<TrackingVertexRef, double> >) vr2s[aVtx];
      }
    }
  }
*/
  // Do fake rate calculation //

  //cout << "Starting K0s fake rate calculation" << endl;
  // Kshorts
  double numK0sFound = 0.;
  double mass = 0.;
  std::vector<double> radDist;
  //  radDist.clear();
  //cout << "K0s collection size: " << k0sCollection->size() << endl;
  if ( k0sCollection->size() > 0 ) {
    //cout << "In loop" << endl;

    vector<reco::TrackRef> theDaughterTracks;
    for( reco::VertexCompositeCandidateCollection::const_iterator iK0s = k0sCollection->begin();
     iK0s != k0sCollection->end();
     iK0s++) {
      //cout << "In loop 2" << endl;
      // Fill mass of all K0S
      ksMassAll->Fill( iK0s->mass() );
      // Fill values to be histogrammed
      K0sCandpT = (sqrt( iK0s->momentum().perp2() ));
      K0sCandEta = iK0s->momentum().eta();
      K0sCandR = (sqrt( iK0s->vertex().perp2() ));
      K0sCandStatus = 0;
      //cout << "MASS" << endl;
      mass = iK0s->mass();
      //cout << "Pushing back daughters" << endl;

      theDaughterTracks.push_back( (*(dynamic_cast<const reco::RecoChargedCandidate *> (iK0s->daughter(0)) )).track() );
      theDaughterTracks.push_back( (*(dynamic_cast<const reco::RecoChargedCandidate *> (iK0s->daughter(1)) )).track() );
       
      //cout << "1" << endl;
      for (int itrack = 0; itrack < 2; itrack++) {
    K0sPiCandStatus[itrack] = 0;
      }

      std::vector< std::pair<TrackingParticleRef, double> > tp;
      TrackingParticleRef tpref;
      TrackingParticleRef firstDauTP;
      TrackingVertexRef k0sVtx;

      //cout << "2" << endl;
      // Loop through K0s candidate daugher tracks
      for(View<reco::Track>::size_type i=0; i<theDaughterTracks.size(); ++i){
    // Found track from theDaughterTracks
    RefToBase<reco::Track> track( theDaughterTracks.at(i) );
        
    //if(recSimColl.find(track) != recSimColl.end()) {
    if(recotosimCollectionH->find(track) != recotosimCollectionH->end()) {
      //tp = recSimColl[track];
      tp = (*recotosimCollectionH)[track];
      if (tp.size() != 0) {
        K0sPiCandStatus[i] = 1;
        tpref = tp.begin()->first;

        //if( simRecColl.find(tpref) == simRecColl.end() ) {
        if( simtorecoCollectionH->find(tpref) == simtorecoCollectionH->end() ) {
          K0sPiCandStatus[i] = 3;
        }
        //cout << "3" << endl;
        TrackingVertexRef parentVertex = tpref->parentVertex();
        if(parentVertex.isNonnull()) radDist.push_back(parentVertex->position().R());
         
        if( parentVertex.isNonnull() ) {
          if( k0sVtx.isNonnull() ) {
        if( k0sVtx->position() == parentVertex->position() ) {
          if( parentVertex->nDaughterTracks() == 2 ) {
            if( parentVertex->nSourceTracks() == 0 ) {
              // No source tracks found for K0s vertex; shouldn't happen, but does for evtGen events
              K0sCandStatus = 6;
            }
            
            for( TrackingVertex::tp_iterator iTP = parentVertex->sourceTracks_begin();
             iTP != parentVertex->sourceTracks_end(); iTP++) {
              if( (*iTP)->pdgId() == 310 ) {
            //cout << "4" << endl;
            K0sCandStatus = 1;
            realK0sFound++;
            numK0sFound += 1.;
            std::pair<TrackingParticleRef, TrackingParticleRef> pair(firstDauTP, tpref);
            // Pushing back a good V0
            trueK0s.push_back(pair);
            trueKsMasses.push_back(mass);
              }
              else {
            K0sCandStatus = 2;
            if( (*iTP)->pdgId() == 3122 ) {
              K0sCandStatus = 7;
            }
              }
            }
          }
          else {
            // Found a bad match because the mother has too many daughters
            K0sCandStatus = 3;
          }
        }
        else {
          // Found a bad match because the parent vertices from the two tracks are different
          K0sCandStatus = 4;
        }
          }
          else {
        // if k0sVtx is null, fill it with parentVertex to compare to the parentVertex from the second track
        k0sVtx = parentVertex;
        firstDauTP = tpref;
          }
        }//parent vertex is null
      }//tp size zero
    }
    else {
      //cout << "5" << endl;
      K0sPiCandStatus[i] = 2;
      noTPforK0sCand++;
      K0sCandStatus = 5;
      theDaughterTracks.clear();
    }
      }
      //cout << "6" << endl;
      theDaughterTracks.clear();
      // fill the fake rate histograms
      if( K0sCandStatus > 1 ) {
    //cout << "7" << endl;
    ksFakeVsR_num->Fill(K0sCandR);
    ksFakeVsEta_num->Fill(K0sCandEta);
    ksFakeVsPt_num->Fill(K0sCandpT);
    ksCandStatus->Fill((float) K0sCandStatus);
    fakeKsMass->Fill(mass);
    for( unsigned int ndx = 0; ndx < radDist.size(); ndx++ ) {
      ksFakeDauRadDist->Fill(radDist[ndx]);
    }
      }
      if( K0sCandStatus == 5 ) {
    ksTkFakeVsR_num->Fill(K0sCandR);
    ksTkFakeVsEta_num->Fill(K0sCandEta);
    ksTkFakeVsPt_num->Fill(K0sCandpT);
      }
      ksFakeVsR_denom->Fill(K0sCandR);
      ksFakeVsEta_denom->Fill(K0sCandEta);
      ksFakeVsPt_denom->Fill(K0sCandpT);
    }
  }
  //cout << "Outside loop, why would it fail here?" << endl;
  //double numK0sFound = (double) realK0sFound;
  //cout << "numK0sFound: " << numK0sFound << endl;
  nKs->Fill( (float) numK0sFound );
  numK0sFound = 0.;

  //cout << "Starting Lambda fake rate calculation" << endl;

  double numLamFound = 0.;
  mass = 0.;
  radDist.clear();
  // Lambdas
  if ( lambdaCollection->size() > 0 ) {
    //cout << "In lam loop." << endl;
    
    vector<reco::TrackRef> theDaughterTracks;
    for( reco::VertexCompositeCandidateCollection::const_iterator iLam = lambdaCollection->begin();
     iLam != lambdaCollection->end();
     iLam++) {
      // Fill mass plot with ALL lambdas
      lamMassAll->Fill( iLam->mass() );
      // Fill values to be histogrammed
      LamCandpT = (sqrt( iLam->momentum().perp2() ));
      LamCandEta = iLam->momentum().eta();
      LamCandR = (sqrt( iLam->vertex().perp2() ));
      LamCandStatus = 0;
      mass = iLam->mass();
      
      //cout << "Lam daughter tracks" << endl;
      theDaughterTracks.push_back( (*(dynamic_cast<const reco::RecoChargedCandidate *> (iLam->daughter(0)) )).track() );
      theDaughterTracks.push_back( (*(dynamic_cast<const reco::RecoChargedCandidate *> (iLam->daughter(1)) )).track() );
      
      for (int itrack = 0; itrack < 2; itrack++) {
    LamPiCandStatus[itrack] = 0;
      }
      
      std::vector< std::pair<TrackingParticleRef, double> > tp;
      TrackingParticleRef tpref;
      TrackingParticleRef firstDauTP;
      TrackingVertexRef LamVtx;
      // Loop through Lambda candidate daughter tracks
      for(View<reco::Track>::size_type i=0; i<theDaughterTracks.size(); ++i){
    // Found track from theDaughterTracks
    //cout << "Looping over lam daughters" << endl;
    RefToBase<reco::Track> track( theDaughterTracks.at(i) );
    
    //if(recSimColl.find(track) != recSimColl.end()) {
    if(recotosimCollectionH->find(track) != recotosimCollectionH->end()) {
      //tp = recSimColl[track];
      tp = (*recotosimCollectionH)[track];
      if (tp.size() != 0) {
        LamPiCandStatus[i] = 1;
        tpref = tp.begin()->first;

        //if( simRecColl.find(tpref) == simRecColl.end() ) {
        if( simtorecoCollectionH->find(tpref) == simtorecoCollectionH->end() ) {
          LamPiCandStatus[i] = 3;
        }
        TrackingVertexRef parentVertex = tpref->parentVertex();
        if( parentVertex.isNonnull() ) radDist.push_back(parentVertex->position().R());
         
        if( parentVertex.isNonnull() ) {
          if( LamVtx.isNonnull() ) {
        if( LamVtx->position() == parentVertex->position() ) {
          if( parentVertex->nDaughterTracks() == 2 ) {
            if( parentVertex->nSourceTracks() == 0 ) {
              // No source tracks found for K0s vertex; shouldn't happen, but does for evtGen events
              LamCandStatus = 6;
            }

            for( TrackingVertex::tp_iterator iTP = parentVertex->sourceTracks_begin();
             iTP != parentVertex->sourceTracks_end(); ++iTP) {
              if( abs((*iTP)->pdgId()) == 3122 ) {
            LamCandStatus = 1;
            realLamFound++;
            numLamFound += 1.;
            std::pair<TrackingParticleRef, TrackingParticleRef> pair(firstDauTP, tpref);
            // Pushing back a good V0
            trueLams.push_back(pair);
            trueLamMasses.push_back(mass);
              }
              else {
            LamCandStatus = 2;
            if( abs((*iTP)->pdgId() ) == 310 ) {
              LamCandStatus = 7;
            }
              }
              //if(iTP != parentVertex->sourceTracks_end()) {
              //cout << "Bogus check 1" << endl;
              //}
            }
          }
          else {
            // Found a bad match because the mother has too many daughters
            LamCandStatus = 3;
          }
        }
        else {
          // Found a bad match because the parent vertices from the two tracks are different
          LamCandStatus = 4;
        }
          }
          else {
        // if lamVtx is null, fill it with parentVertex to compare to the parentVertex from the second track
        LamVtx = parentVertex;
        firstDauTP = tpref;
          }
        }//parent vertex is null
      }//tp size zero
    }
    else {
      LamPiCandStatus[i] = 2;
      noTPforLamCand++;
      LamCandStatus = 5;
      theDaughterTracks.clear();
    }
      }
      theDaughterTracks.clear();
      // fill the fake rate histograms
      //cout << "Fill lam fake rate histos" << endl;
      if( LamCandStatus > 1 ) {
    //cout << "fake 1" << endl;
    //cout << "fake 1.5" << endl;
    lamFakeVsR_num->Fill(LamCandR);
    //cout << "fake 2" << endl;
    lamFakeVsEta_num->Fill(LamCandEta);
    //cout << "fake 3" << endl;
    lamFakeVsPt_num->Fill(LamCandpT);
    //cout << "fake 4" << endl;
    lamCandStatus->Fill((float) LamCandStatus);
    //cout << "fake 5" << endl;
    fakeLamMass->Fill(mass);
    //cout << "fake 6" << endl;
    for( unsigned int ndx = 0; ndx < radDist.size(); ndx++ ) {
      lamFakeDauRadDist->Fill(radDist[ndx]);
    }
      }
      //cout << "Fill lam Tk fake histos" << endl;
      if( K0sCandStatus == 5 ) {
    lamTkFakeVsR_num->Fill(LamCandR);
    lamTkFakeVsEta_num->Fill(LamCandEta);
    lamTkFakeVsPt_num->Fill(LamCandpT);
      }
      //cout << "Fill denominators" << endl;
      lamFakeVsR_denom->Fill(LamCandR);
      lamFakeVsEta_denom->Fill(LamCandEta);
      lamFakeVsPt_denom->Fill(LamCandpT);
    }
  }
  //cout << "Filling numLamFound" << endl;
  nLam->Fill( (double) numLamFound );
  numLamFound = 0.;


  // Do efficiency calculation //

  //cout << "Starting Lambda efficiency" << endl;
  // Lambdas

  for(TrackingParticleCollection::size_type i = 0; i < tPCeff.size(); i++) {
    TrackingParticleRef tpr1(TPCollectionEff, i);
    const TrackingParticle* itp1 = tpr1.get();
    if( (itp1->pdgId() == 211 || itp1->pdgId() == 2212)
    && itp1->parentVertex().isNonnull()
    && abs(itp1->momentum().eta()) < 2.4
    && sqrt( itp1->momentum().perp2() ) > 0.9) {
      bool isLambda = false;
      if( itp1->pdgId() == 2212 ) isLambda = true;
      if( itp1->parentVertex()->nDaughterTracks() == 2 ) {

    TrackingVertexRef piCand1Vertex = itp1->parentVertex();
    for(TrackingVertex::tp_iterator iTP1 = piCand1Vertex->sourceTracks_begin();
        iTP1 != piCand1Vertex->sourceTracks_end(); iTP1++) {
      if( abs((*iTP1)->pdgId()) == 3122 ) {
        //double motherpT = (*iTP1)->pt();
        //       ----->>>>>>Keep going here
        for(TrackingParticleCollection::size_type j=0;
        j < tPCeff.size();
        j++) {
          TrackingParticleRef tpr2(TPCollectionEff, j);
          const TrackingParticle* itp2 = tpr2.get();
          int particle2pdgId;
          if (isLambda) particle2pdgId = -211;
          else particle2pdgId = -2212;
          if( itp2->pdgId() == particle2pdgId
          && itp2->parentVertex().isNonnull()
          && abs(itp2->momentum().eta()) < 2.4
          && sqrt(itp2->momentum().perp2()) > 0.9) {
        if(itp2->parentVertex() == itp1->parentVertex()) {
          // Found a good pair of Lambda daughters
          TrackingVertexRef piCand2Vertex = itp2->parentVertex();
          for (TrackingVertex::tp_iterator iTP2 = piCand2Vertex->sourceTracks_begin();
               iTP2 != piCand2Vertex->sourceTracks_end(); 
               ++iTP2) {
            LamGenEta = LamGenpT = LamGenR = 0.;
            LamGenStatus = 0;
            for(int ifill = 0;
            ifill < 2;
            ifill++) {
              // do nothing?
            }
            if( abs((*iTP2)->pdgId()) == 3122 ) {
              // found generated Lambda
              
              LamGenpT = sqrt((*iTP2)->momentum().perp2());
              LamGenEta = (*iTP2)->momentum().eta();
              LamGenR = sqrt(itp2->vertex().perp2());
              genLam++;
              if(trueLams.size() > 0) {
            int loop_1 = 0;
            for(std::vector< pair<TrackingParticleRef, TrackingParticleRef> >::const_iterator iEffCheck = trueLams.begin();
                iEffCheck != trueLams.end();
                iEffCheck++) {
              //cout << "In LOOP" << endl;
              if( itp1->parentVertex() == iEffCheck->first->parentVertex()
                  && itp2->parentVertex() == iEffCheck->second->parentVertex() ) {
                realLamFoundEff++;
                //V0Producer found the generated Lambda
                LamGenStatus = 1;
                //cout << "Maybe it's here.." << endl;
                goodLamMass->Fill(trueLamMasses[loop_1]);
                //cout << "Did we make it?" << endl;
                break;
              }
              else {
                //V0Producer didn't find the generated Lambda
                LamGenStatus = 2;
              }
              loop_1++;
            }
              }
              else {
            //No V0 cand found, so V0Producer didn't find the generated Lambda
            LamGenStatus = 2;
              }
              std::vector< std::pair<RefToBase<reco::Track>, double> > rt1;
              std::vector< std::pair<RefToBase<reco::Track>, double> > rt2;
              
              //if( simRecColl.find(tpr1) != simRecColl.end() ) {
              if( simtorecoCollectionH->find(tpr1) != simtorecoCollectionH->end() ) {
            //rt1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) simRecColl[tpr1];
            rt1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) (*simtorecoCollectionH)[tpr1];
            if(rt1.size() != 0) {
              LamPiEff[0] = 1; //Found the first daughter track
              edm::RefToBase<reco::Track> t1 = rt1.begin()->first;
            }
              }
              else {
            LamPiEff[0] = 2;//First daughter not found
              }
              //if( (simRecColl.find(tpr2) != simRecColl.end()) ) {
              if( (simtorecoCollectionH->find(tpr2) != simtorecoCollectionH->end()) ) {
            //rt2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) simRecColl[tpr2];
            rt2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) (*simtorecoCollectionH)[tpr2];
            if(rt2.size() != 0) {
              LamPiEff[1] = 1;//Found the second daughter track
              edm::RefToBase<reco::Track> t2 = rt2.begin()->first;
            }
              }
              else {
            LamPiEff[1] = 2;//Second daughter not found
              }
              
              if( LamGenStatus == 1
              && (LamPiEff[0] == 2 || LamPiEff[1] == 2) ) {
            // Good Lambda found, but recoTrack->trackingParticle->recoTrack didn't work
            LamGenStatus = 4;
            realLamFoundEff--;
              }
              if( LamGenStatus == 2
              && (LamPiEff[0] == 2 || LamPiEff[1] == 2) ) {
            // Lambda not found because we didn't find a daughter track
            LamGenStatus = 3;
              }
              //cout << "LamGenStatus: " << LamGenStatus << ", LamPiEff[i]: " << LamPiEff[0] << ", " << LamPiEff[1] << endl;
              // Fill histograms
              if(LamGenR > 0.) {
            if(LamGenStatus == 1) {
              lamEffVsR_num->Fill(LamGenR);
            }
            if((double) LamGenStatus < 2.5) {
              lamTkEffVsR_num->Fill(LamGenR);
            }
            lamEffVsR_denom->Fill(LamGenR);
              }
              if(abs(LamGenEta) > 0.) {
            if(LamGenStatus == 1) {
              lamEffVsEta_num->Fill(LamGenEta);
            }
            if((double) LamGenStatus < 2.5) {
              lamTkEffVsEta_num->Fill(LamGenEta);
            }
            lamEffVsEta_denom->Fill(LamGenEta);
              }
              if(LamGenpT > 0.) {
            if(LamGenStatus == 1) {
              lamEffVsPt_num->Fill(LamGenpT);
            }
            if((double) LamGenStatus < 2.5) {
              lamTkEffVsPt_num->Fill(LamGenpT);
            }
            lamEffVsPt_denom->Fill(LamGenpT);
              }
            }
          }
        }
          }
        }
      }
    }
      }
    }
  }

  //Kshorts

  //cout << "Starting Kshort efficiency" << endl;
  for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
    TrackingParticleRef tpr1(TPCollectionEff, i);
    const TrackingParticle* itp1 = tpr1.get();
    // only count the efficiency for pions with |eta|<2.4 and pT>0.9 GeV. First search for a suitable pi+
    if ( itp1->pdgId() == 211 
     && itp1->parentVertex().isNonnull() 
     && abs(itp1->momentum().eta()) < 2.4 
     && sqrt(itp1->momentum().perp2()) > 0.9) {
      if ( itp1->parentVertex()->nDaughterTracks() == 2 ) {
    TrackingVertexRef piCand1Vertex = itp1->parentVertex();        
    //check trackingParticle pion for a Ks mother
    for (TrackingVertex::tp_iterator iTP1 = piCand1Vertex->sourceTracks_begin();
         iTP1 != piCand1Vertex->sourceTracks_end(); ++iTP1) {
      //iTP1 is a TrackingParticle
      if ( (*iTP1)->pdgId()==310 ) {
        //with a Ks mother found for the pi+, loop through trackingParticles again to find a pi-
        for (TrackingParticleCollection::size_type j=0; j<tPCeff.size(); j++){
          TrackingParticleRef tpr2(TPCollectionEff, j);
          const TrackingParticle* itp2 = tpr2.get();
          
          if ( itp2->pdgId() == -211 && itp2->parentVertex().isNonnull()  
           && abs(itp2->momentum().eta()) < 2.4 
           && sqrt(itp2->momentum().perp2()) > 0.9) {
        //check the pi+ and pi- have the same vertex
        if ( itp2->parentVertex() == itp1->parentVertex() ) {
          TrackingVertexRef piCand2Vertex = itp2->parentVertex();          
          for (TrackingVertex::tp_iterator iTP2 = piCand2Vertex->sourceTracks_begin();
               iTP2 != piCand2Vertex->sourceTracks_end(); ++iTP2) {
            //iTP2 is a TrackingParticle
            K0sGenEta = K0sGenpT = K0sGenR = 0.;
            K0sGenStatus = 0;
            if( (*iTP2)->pdgId() == 310 ) {
              K0sGenpT = sqrt( (*iTP2)->momentum().perp2() );
              K0sGenEta = (*iTP2)->momentum().eta();
              K0sGenR = sqrt(itp2->vertex().perp2());
              genK0s++;
              int loop_2 = 0;
              if( trueK0s.size() > 0 ) {
            for( std::vector< pair<TrackingParticleRef, TrackingParticleRef> >::const_iterator iEffCheck = trueK0s.begin();
                 iEffCheck != trueK0s.end();
                 iEffCheck++) {
              //if the parent vertices for the tracks are the same, then the generated Ks was found
              if (itp1->parentVertex()==iEffCheck->first->parentVertex() &&
                  itp2->parentVertex()==iEffCheck->second->parentVertex())  {
                realK0sFoundEff++;
                K0sGenStatus = 1;
                //cout << "Maybe here?" << endl;
                goodKsMass->Fill(trueKsMasses[loop_2]);
                //cout << "We made it...." << endl;
                break;
              }
              else {
                K0sGenStatus = 2;
              }
            }
              }
              else {
            K0sGenStatus = 2;
              }

              // Check if the generated Ks tracks were found or not
              // by searching the recoTracks list for a match to the trackingParticles

              std::vector<std::pair<RefToBase<reco::Track>, double> > rt1;
              std::vector<std::pair<RefToBase<reco::Track>, double> > rt2;
              
              //if( simRecColl.find(tpr1) != simRecColl.end() ) {
              if( simtorecoCollectionH->find(tpr1) != simtorecoCollectionH->end() ) {
            rt1 = (std::vector< std::pair<RefToBase<reco::Track>, double> >) (*simtorecoCollectionH)[tpr1];
//simRecColl[tpr1];
            if(rt1.size() != 0) {
              //First pion found
              K0sPiEff[0] = 1;
              edm::RefToBase<reco::Track> t1 = rt1.begin()->first;
            }
              }
              else {
            //First pion not found
            K0sPiEff[0] = 2;
              }
              
              //if( simRecColl.find(tpr2) != simRecColl.end() ) {
              if( simtorecoCollectionH->find(tpr2) != simtorecoCollectionH->end() ) {
            rt2 = (std::vector< std::pair<RefToBase<reco::Track>, double> >) (*simtorecoCollectionH)[tpr2];
//simRecColl[tpr2];
            if(rt2.size() != 0) {
              //Second pion found
              K0sPiEff[1] = 1;
              edm::RefToBase<reco::Track> t2 = rt2.begin()->first;
            }
              }
              else {
            K0sPiEff[1] = 2;
              }
              //cout << "Status: " << K0sGenStatus << ", K0sPiEff[i]: " << K0sPiEff[0] << ", " << K0sPiEff[1] << endl;
              if(K0sGenStatus == 1
             && (K0sPiEff[0] == 2 || K0sPiEff[1] == 2)) {
            K0sGenStatus = 4;
            realK0sFoundEff--;
              }
              if(K0sGenStatus == 2
             && (K0sPiEff[0] == 2 || K0sPiEff[1] == 2)) {
            K0sGenStatus = 3;
              }
              if(K0sPiEff[0] == 1 && K0sPiEff[1] == 1) {
            k0sTracksFound++;
              }
              //Fill Histograms
              if(K0sGenR > 0.) {
            if(K0sGenStatus == 1) {
              ksEffVsR_num->Fill(K0sGenR);
            }
            if((double) K0sGenStatus < 2.5) {             
              ksTkEffVsR_num->Fill(K0sGenR);
            }
            ksEffVsR_denom->Fill(K0sGenR);
              }
              if(abs(K0sGenEta) > 0.) {
            if(K0sGenStatus == 1) {
              ksEffVsEta_num->Fill(K0sGenEta);
            }
            if((double) K0sGenStatus < 2.5) {
              ksTkEffVsEta_num->Fill(K0sGenEta);
            }
            ksEffVsEta_denom->Fill(K0sGenEta);
              }
              if(K0sGenpT > 0.) {
            if(K0sGenStatus == 1) {
              ksEffVsPt_num->Fill(K0sGenpT);
            }
            if((double) K0sGenStatus < 2.5) {
              ksTkEffVsPt_num->Fill(K0sGenpT);
            }
            ksEffVsPt_denom->Fill(K0sGenpT);
              }
            }
          }
        }
          }
        }
      }
    }
      }
    }
  }

  delete thePrimary;
}

void V0Validator::beginRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 67 of file V0Validator.cc.

References DQMStore::book1D(), DQMStore::cd(), dirName, fakeKsMass, fakeLamMass, goodKsMass, goodLamMass, ksAbsoluteDistResolution, ksCandStatus, ksEffVsEta, ksEffVsEta_denom, ksEffVsEta_num, ksEffVsPt, ksEffVsPt_denom, ksEffVsPt_num, ksEffVsR, ksEffVsR_denom, ksEffVsR_num, ksFakeDauRadDist, ksFakeVsEta, ksFakeVsEta_denom, ksFakeVsEta_num, ksFakeVsPt, ksFakeVsPt_denom, ksFakeVsPt_num, ksFakeVsR, ksFakeVsR_denom, ksFakeVsR_num, ksMassAll, ksTkEffVsEta, ksTkEffVsEta_num, ksTkEffVsPt, ksTkEffVsPt_num, ksTkEffVsR, ksTkEffVsR_num, ksTkFakeVsEta, ksTkFakeVsEta_num, ksTkFakeVsPt, ksTkFakeVsPt_num, ksTkFakeVsR, ksTkFakeVsR_num, ksXResolution, ksYResolution, ksZResolution, lamAbsoluteDistResolution, lamCandStatus, lamEffVsEta, lamEffVsEta_denom, lamEffVsEta_num, lamEffVsPt, lamEffVsPt_denom, lamEffVsPt_num, lamEffVsR, lamEffVsR_denom, lamEffVsR_num, lamFakeDauRadDist, lamFakeVsEta, lamFakeVsEta_denom, lamFakeVsEta_num, lamFakeVsPt, lamFakeVsPt_denom, lamFakeVsPt_num, lamFakeVsR, lamFakeVsR_denom, lamFakeVsR_num, lamMassAll, lamTkEffVsEta, lamTkEffVsEta_num, lamTkEffVsPt, lamTkEffVsPt_num, lamTkEffVsR, lamTkEffVsR_num, lamTkFakeVsEta, lamTkFakeVsEta_num, lamTkFakeVsPt, lamTkFakeVsPt_num, lamTkFakeVsR, lamTkFakeVsR_num, lamXResolution, lamYResolution, lamZResolution, nKs, nLam, DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and theDQMstore.

                                                                          {
  //std::cout << "Running V0Validator" << std::endl;
  //theDQMstore = edm::Service<DQMStore>().operator->();
  //std::cout << "In beginJob() at line 1" << std::endl;
  //edm::Service<TFileService> fs;

  theDQMstore->cd();
  std::string subDirName = dirName + "/EffFakes";
  theDQMstore->setCurrentFolder(subDirName.c_str());

  ksEffVsR = theDQMstore->book1D("K0sEffVsR", 
              "K^{0}_{S} Efficiency vs #rho", 40, 0., 40.);
  ksEffVsEta = theDQMstore->book1D("K0sEffVsEta",
                "K^{0}_{S} Efficiency vs #eta", 40, -2.5, 2.5);
  ksEffVsPt = theDQMstore->book1D("K0sEffVsPt",
               "K^{0}_{S} Efficiency vs p_{T}", 70, 0., 20.);;

  ksTkEffVsR = theDQMstore->book1D("K0sTkEffVsR", 
              "K^{0}_{S} Tracking Efficiency vs #rho", 40, 0., 40.);
  ksTkEffVsEta = theDQMstore->book1D("K0sTkEffVsEta",
                "K^{0}_{S} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
  ksTkEffVsPt = theDQMstore->book1D("K0sTkEffVsPt",
               "K^{0}_{S} Tracking Efficiency vs p_{T}", 70, 0., 20.);

  ksEffVsR_num = theDQMstore->book1D("K0sEffVsR_num", 
              "K^{0}_{S} Efficiency vs #rho", 40, 0., 40.);
  ksEffVsEta_num = theDQMstore->book1D("K0sEffVsEta_num",
                "K^{0}_{S} Efficiency vs #eta", 40, -2.5, 2.5);
  ksEffVsPt_num = theDQMstore->book1D("K0sEffVsPt_num",
               "K^{0}_{S} Efficiency vs p_{T}", 70, 0., 20.);;

  ksTkEffVsR_num = theDQMstore->book1D("K0sTkEffVsR_num", 
              "K^{0}_{S} Tracking Efficiency vs #rho", 40, 0., 40.);
  ksTkEffVsEta_num = theDQMstore->book1D("K0sTkEffVsEta_num",
                "K^{0}_{S} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
  ksTkEffVsPt_num = theDQMstore->book1D("K0sTkEffVsPt_num",
               "K^{0}_{S} Tracking Efficiency vs p_{T}", 70, 0., 20.);;


  ksEffVsR_denom = theDQMstore->book1D("K0sEffVsR_denom", 
              "K^{0}_{S} Efficiency vs #rho", 40, 0., 40.);
  ksEffVsEta_denom = theDQMstore->book1D("K0sEffVsEta_denom",
                "K^{0}_{S} Efficiency vs #eta", 40, -2.5, 2.5);
  ksEffVsPt_denom = theDQMstore->book1D("K0sEffVsPt_denom",
               "K^{0}_{S} Efficiency vs p_{T}", 70, 0., 20.);;


  lamEffVsR = theDQMstore->book1D("LamEffVsR",
               "#Lambda^{0} Efficiency vs #rho", 40, 0., 40.);
  lamEffVsEta = theDQMstore->book1D("LamEffVsEta",
                 "#Lambda^{0} Efficiency vs #eta", 40, -2.5, 2.5);
  lamEffVsPt = theDQMstore->book1D("LamEffVsPt",
                "#Lambda^{0} Efficiency vs p_{T}", 70, 0., 20.);


  lamTkEffVsR = theDQMstore->book1D("LamTkEffVsR",
               "#Lambda^{0} TrackingEfficiency vs #rho", 40, 0., 40.);
  lamTkEffVsEta = theDQMstore->book1D("LamTkEffVsEta",
                 "#Lambda^{0} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
  lamTkEffVsPt = theDQMstore->book1D("LamTkEffVsPt",
                "#Lambda^{0} Tracking Efficiency vs p_{T}", 70, 0., 20.);

  lamEffVsR_num = theDQMstore->book1D("LamEffVsR_num",
               "#Lambda^{0} Efficiency vs #rho", 40, 0., 40.);
  lamEffVsEta_num = theDQMstore->book1D("LamEffVsEta_num",
                 "#Lambda^{0} Efficiency vs #eta", 40, -2.5, 2.5);
  lamEffVsPt_num = theDQMstore->book1D("LamEffVsPt_num",
                "#Lambda^{0} Efficiency vs p_{T}", 70, 0., 20.);


  lamTkEffVsR_num = theDQMstore->book1D("LamTkEffVsR_num",
               "#Lambda^{0} TrackingEfficiency vs #rho", 40, 0., 40.);
  lamTkEffVsEta_num = theDQMstore->book1D("LamTkEffVsEta_num",
                 "#Lambda^{0} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
  lamTkEffVsPt_num = theDQMstore->book1D("LamTkEffVsPt_num",
                "#Lambda^{0} Tracking Efficiency vs p_{T}", 70, 0., 20.);


  lamEffVsR_denom = theDQMstore->book1D("LamEffVsR_denom",
               "#Lambda^{0} Efficiency vs #rho", 40, 0., 40.);
  lamEffVsEta_denom = theDQMstore->book1D("LamEffVsEta_denom",
                 "#Lambda^{0} Efficiency vs #eta", 40, -2.5, 2.5);
  lamEffVsPt_denom = theDQMstore->book1D("LamEffVsPt_denom",
                "#Lambda^{0} Efficiency vs p_{T}", 70, 0., 20.);

  //theDQMstore->cd();
  //subDirName = dirName + "/Fake";
  //theDQMstore->setCurrentFolder(subDirName.c_str());


  ksFakeVsR = theDQMstore->book1D("K0sFakeVsR",
               "K^{0}_{S} Fake Rate vs #rho", 40, 0., 40.);
  ksFakeVsEta = theDQMstore->book1D("K0sFakeVsEta",
                 "K^{0}_{S} Fake Rate vs #eta", 40, -2.5, 2.5);
  ksFakeVsPt = theDQMstore->book1D("K0sFakeVsPt",
                "K^{0}_{S} Fake Rate vs p_{T}", 70, 0., 20.);
  ksTkFakeVsR = theDQMstore->book1D("K0sTkFakeVsR",
               "K^{0}_{S} Tracking Fake Rate vs #rho", 40, 0., 40.);
  ksTkFakeVsEta = theDQMstore->book1D("K0sTkFakeVsEta",
                 "K^{0}_{S} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
  ksTkFakeVsPt = theDQMstore->book1D("K0sTkFakeVsPt",
                "K^{0}_{S} Tracking Fake Rate vs p_{T}", 70, 0., 20.);

  ksFakeVsR_num = theDQMstore->book1D("K0sFakeVsR_num",
               "K^{0}_{S} Fake Rate vs #rho", 40, 0., 40.);
  ksFakeVsEta_num = theDQMstore->book1D("K0sFakeVsEta_num",
                 "K^{0}_{S} Fake Rate vs #eta", 40, -2.5, 2.5);
  ksFakeVsPt_num = theDQMstore->book1D("K0sFakeVsPt_num",
                "K^{0}_{S} Fake Rate vs p_{T}", 70, 0., 20.);
  ksTkFakeVsR_num = theDQMstore->book1D("K0sTkFakeVsR_num",
               "K^{0}_{S} Tracking Fake Rate vs #rho", 40, 0., 40.);
  ksTkFakeVsEta_num = theDQMstore->book1D("K0sTkFakeVsEta_num",
                 "K^{0}_{S} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
  ksTkFakeVsPt_num = theDQMstore->book1D("K0sTkFakeVsPt_num",
                "K^{0}_{S} Tracking Fake Rate vs p_{T}", 70, 0., 20.);

  ksFakeVsR_denom = theDQMstore->book1D("K0sFakeVsR_denom",
               "K^{0}_{S} Fake Rate vs #rho", 40, 0., 40.);
  ksFakeVsEta_denom = theDQMstore->book1D("K0sFakeVsEta_denom",
                 "K^{0}_{S} Fake Rate vs #eta", 40, -2.5, 2.5);
  ksFakeVsPt_denom = theDQMstore->book1D("K0sFakeVsPt_denom",
                "K^{0}_{S} Fake Rate vs p_{T}", 70, 0., 20.);

  lamFakeVsR = theDQMstore->book1D("LamFakeVsR",
                "#Lambda^{0} Fake Rate vs #rho", 40, 0., 40.);
  lamFakeVsEta = theDQMstore->book1D("LamFakeVsEta",
                  "#Lambda^{0} Fake Rate vs #eta", 40, -2.5, 2.5);
  lamFakeVsPt = theDQMstore->book1D("LamFakeVsPt",
                 "#Lambda^{0} Fake Rate vs p_{T}", 70, 0., 20.);
  lamTkFakeVsR = theDQMstore->book1D("LamTkFakeVsR",
                "#Lambda^{0} Tracking Fake Rate vs #rho", 40, 0., 40.);
  lamTkFakeVsEta = theDQMstore->book1D("LamTkFakeVsEta",
                  "#Lambda^{0} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
  lamTkFakeVsPt = theDQMstore->book1D("LamTkFakeVsPt",
                 "#Lambda^{0} Tracking Fake Rate vs p_{T}", 70, 0., 20.);

  lamFakeVsR_num = theDQMstore->book1D("LamFakeVsR_num",
                "#Lambda^{0} Fake Rate vs #rho", 40, 0., 40.);
  lamFakeVsEta_num = theDQMstore->book1D("LamFakeVsEta_num",
                  "#Lambda^{0} Fake Rate vs #eta", 40, -2.5, 2.5);
  lamFakeVsPt_num = theDQMstore->book1D("LamFakeVsPt_num",
                 "#Lambda^{0} Fake Rate vs p_{T}", 70, 0., 20.);
  lamTkFakeVsR_num = theDQMstore->book1D("LamTkFakeVsR_num",
                "#Lambda^{0} Tracking Fake Rate vs #rho", 40, 0., 40.);
  lamTkFakeVsEta_num = theDQMstore->book1D("LamTkFakeVsEta_num",
                  "#Lambda^{0} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
  lamTkFakeVsPt_num = theDQMstore->book1D("LamTkFakeVsPt_num",
                 "#Lambda^{0} Tracking Fake Rate vs p_{T}", 70, 0., 20.);

  lamFakeVsR_denom = theDQMstore->book1D("LamFakeVsR_denom",
                "#Lambda^{0} Fake Rate vs #rho", 40, 0., 40.);
  lamFakeVsEta_denom = theDQMstore->book1D("LamFakeVsEta_denom",
                  "#Lambda^{0} Fake Rate vs #eta", 40, -2.5, 2.5);
  lamFakeVsPt_denom = theDQMstore->book1D("LamFakeVsPt_denom",
                 "#Lambda^{0} Fake Rate vs p_{T}", 70, 0., 20.);

  theDQMstore->cd();
  subDirName = dirName + "/Other";
  theDQMstore->setCurrentFolder(subDirName.c_str());

  nKs = theDQMstore->book1D("nK0s",
             "Number of K^{0}_{S} found per event", 60, 0., 60.);
  nLam = theDQMstore->book1D("nLam",
              "Number of #Lambda^{0} found per event", 60, 0., 60.);

  ksXResolution = theDQMstore->book1D("ksXResolution",
                   "Resolution of V0 decay vertex X coordinate", 50, 0., 50.);
  ksYResolution = theDQMstore->book1D("ksYResolution",
                   "Resolution of V0 decay vertex Y coordinate", 50, 0., 50.);
  ksZResolution = theDQMstore->book1D("ksZResolution",
                   "Resolution of V0 decay vertex Z coordinate", 50, 0., 50.);
  lamXResolution = theDQMstore->book1D("lamXResolution",
                "Resolution of V0 decay vertex X coordinate", 50, 0., 50.);
  lamYResolution = theDQMstore->book1D("lamYResolution",
                "Resolution of V0 decay vertex Y coordinate", 50, 0., 50.);
  lamZResolution = theDQMstore->book1D("lamZResolution",
                "Resolution of V0 decay vertex Z coordinate", 50, 0., 50.);
  ksAbsoluteDistResolution = theDQMstore->book1D("ksRResolution",
                      "Resolution of absolute distance from primary vertex to V0 vertex",
                      100, 0., 50.);
  lamAbsoluteDistResolution = theDQMstore->book1D("lamRResolution",
                       "Resolution of absolute distance from primary vertex to V0 vertex",
                       100, 0., 50.);

  ksCandStatus = theDQMstore->book1D("ksCandStatus",
              "Fake type by cand status",
              10, 0., 10.);
  lamCandStatus = theDQMstore->book1D("ksCandStatus",
              "Fake type by cand status",
              10, 0., 10.);

  double minKsMass = 0.49767 - 0.07;
  double maxKsMass = 0.49767 + 0.07;
  double minLamMass = 1.1156 - 0.05;
  double maxLamMass = 1.1156 + 0.05;
  int ksMassNbins = 100;
  double ksMassXmin = minKsMass;
  double ksMassXmax = maxKsMass;
  int lamMassNbins = 100;
  double lamMassXmin = minLamMass;
  double lamMassXmax = maxLamMass;

  fakeKsMass = theDQMstore->book1D("ksMassFake",
                 "Mass of fake K0S",
                 ksMassNbins, minKsMass, maxKsMass);
  goodKsMass = theDQMstore->book1D("ksMassGood",
                 "Mass of good reco K0S",
                 ksMassNbins, minKsMass, maxKsMass);
  fakeLamMass = theDQMstore->book1D("lamMassFake",
                  "Mass of fake Lambda",
                  lamMassNbins, minLamMass, maxLamMass);
  goodLamMass = theDQMstore->book1D("lamMassGood",
                  "Mass of good Lambda",
                  lamMassNbins, minLamMass, maxLamMass);

  ksMassAll = theDQMstore->book1D("ksMassAll",
                  "Invariant mass of all K0S",
                  ksMassNbins, ksMassXmin, ksMassXmax);
  lamMassAll = theDQMstore->book1D("lamMassAll",
                   "Invariant mass of all #Lambda^{0}",
                   lamMassNbins, lamMassXmin, lamMassXmax);

  ksFakeDauRadDist = theDQMstore->book1D("radDistFakeKs",
                   "Production radius of daughter particle of Ks fake",
                   100, 0., 15.);
  lamFakeDauRadDist = theDQMstore->book1D("radDistFakeLam",
                    "Production radius of daughter particle of Lam fake",
                    100, 0., 15.);

  /*

  ksEffVsRHist = new TH1F("K0sEffVsR", 
              "K^{0}_{S} Efficiency vs #rho", 40, 0., 40.);
  ksEffVsEtaHist = new TH1F("K0sEffVsEta",
                "K^{0}_{S} Efficiency vs #eta", 40, -2.5, 2.5);
  ksEffVsPtHist = new TH1F("K0sEffVsPt",
               "K^{0}_{S} Efficiency vs p_{T}", 70, 0., 20.);;
  ksFakeVsRHist = new TH1F("K0sFakeVsR",
               "K^{0}_{S} Fake Rate vs #rho", 40, 0., 40.);
  ksFakeVsEtaHist = new TH1F("K0sFakeVsEta",
                 "K^{0}_{S} Fake Rate vs #eta", 40, -2.5, 2.5);
  ksFakeVsPtHist = new TH1F("K0sFakeVsPt",
                "K^{0}_{S} Fake Rate vs p_{T}", 70, 0., 20.);

  ksTkEffVsRHist = new TH1F("K0sTkEffVsR", 
              "K^{0}_{S} Tracking Efficiency vs #rho", 40, 0., 40.);
  ksTkEffVsEtaHist = new TH1F("K0sTkEffVsEta",
                "K^{0}_{S} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
  ksTkEffVsPtHist = new TH1F("K0sTkEffVsPt",
               "K^{0}_{S} Tracking Efficiency vs p_{T}", 70, 0., 20.);;
  ksTkFakeVsRHist = new TH1F("K0sTkFakeVsR",
               "K^{0}_{S} Tracking Fake Rate vs #rho", 40, 0., 40.);
  ksTkFakeVsEtaHist = new TH1F("K0sTkFakeVsEta",
                 "K^{0}_{S} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
  ksTkFakeVsPtHist = new TH1F("K0sTkFakeVsPt",
                "K^{0}_{S} Tracking Fake Rate vs p_{T}", 70, 0., 20.);

  ksEffVsRHist_denom = new TH1F("K0sEffVsR_denom", 
              "K^{0}_{S} Efficiency vs #rho", 40, 0., 40.);
  ksEffVsEtaHist_denom = new TH1F("K0sEffVsEta_denom",
                "K^{0}_{S} Efficiency vs #eta", 40, -2.5, 2.5);
  ksEffVsPtHist_denom = new TH1F("K0sEffVsPt_denom",
               "K^{0}_{S} Efficiency vs p_{T}", 70, 0., 20.);;
  ksFakeVsRHist_denom = new TH1F("K0sFakeVsR_denom",
               "K^{0}_{S} Fake Rate vs #rho", 40, 0., 40.);
  ksFakeVsEtaHist_denom = new TH1F("K0sFakeVsEta_denom",
                 "K^{0}_{S} Fake Rate vs #eta", 40, -2.5, 2.5);
  ksFakeVsPtHist_denom = new TH1F("K0sFakeVsPt_denom",
                "K^{0}_{S} Fake Rate vs p_{T}", 70, 0., 20.);

  lamEffVsRHist = new TH1F("LamEffVsR",
               "#Lambda^{0} Efficiency vs #rho", 40, 0., 40.);
  lamEffVsEtaHist = new TH1F("LamEffVsEta",
                 "#Lambda^{0} Efficiency vs #eta", 40, -2.5, 2.5);
  lamEffVsPtHist = new TH1F("LamEffVsPt",
                "#Lambda^{0} Efficiency vs p_{T}", 70, 0., 20.);
  lamFakeVsRHist = new TH1F("LamFakeVsR",
                "#Lambda^{0} Fake Rate vs #rho", 40, 0., 40.);
  lamFakeVsEtaHist = new TH1F("LamFakeVsEta",
                  "#Lambda^{0} Fake Rate vs #eta", 40, -2.5, 2.5);
  lamFakeVsPtHist = new TH1F("LamFakeVsPt",
                 "#Lambda^{0} Fake Rate vs p_{T}", 70, 0., 20.);

  lamTkEffVsRHist = new TH1F("LamTkEffVsR",
               "#Lambda^{0} TrackingEfficiency vs #rho", 40, 0., 40.);
  lamTkEffVsEtaHist = new TH1F("LamTkEffVsEta",
                 "#Lambda^{0} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
  lamTkEffVsPtHist = new TH1F("LamTkEffVsPt",
                "#Lambda^{0} Tracking Efficiency vs p_{T}", 70, 0., 20.);
  lamTkFakeVsRHist = new TH1F("LamTkFakeVsR",
                "#Lambda^{0} Tracking Fake Rate vs #rho", 40, 0., 40.);
  lamTkFakeVsEtaHist = new TH1F("LamTkFakeVsEta",
                  "#Lambda^{0} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
  lamTkFakeVsPtHist = new TH1F("LamTkFakeVsPt",
                 "#Lambda^{0} Tracking Fake Rate vs p_{T}", 70, 0., 20.);

  lamEffVsRHist_denom = new TH1F("LamEffVsR_denom",
               "#Lambda^{0} Efficiency vs #rho", 40, 0., 40.);
  lamEffVsEtaHist_denom = new TH1F("LamEffVsEta_denom",
                 "#Lambda^{0} Efficiency vs #eta", 40, -2.5, 2.5);
  lamEffVsPtHist_denom = new TH1F("LamEffVsPt_denom",
                "#Lambda^{0} Efficiency vs p_{T}", 70, 0., 20.);
  lamFakeVsRHist_denom = new TH1F("LamFakeVsR_denom",
                "#Lambda^{0} Fake Rate vs #rho", 40, 0., 40.);
  lamFakeVsEtaHist_denom = new TH1F("LamFakeVsEta_denom",
                  "#Lambda^{0} Fake Rate vs #eta", 40, -2.5, 2.5);
  lamFakeVsPtHist_denom = new TH1F("LamFakeVsPt_denom",
                 "#Lambda^{0} Fake Rate vs p_{T}", 70, 0., 20.);

  nKsHist = new TH1F("nK0s",
             "Number of K^{0}_{S} found per event", 60, 0., 60.);
  nLamHist = new TH1F("nLam",
              "Number of #Lambda^{0} found per event", 60, 0., 60.);

  ksXResolutionHist = new TH1F("ksXResolution",
                   "Resolution of V0 decay vertex X coordinate", 50, 0., 50.);
  ksYResolutionHist = new TH1F("ksYResolution",
                   "Resolution of V0 decay vertex Y coordinate", 50, 0., 50.);
  ksZResolutionHist = new TH1F("ksZResolution",
                   "Resolution of V0 decay vertex Z coordinate", 50, 0., 50.);
  lamXResolutionHist = new TH1F("lamXResolution",
                "Resolution of V0 decay vertex X coordinate", 50, 0., 50.);
  lamYResolutionHist = new TH1F("lamYResolution",
                "Resolution of V0 decay vertex Y coordinate", 50, 0., 50.);
  lamZResolutionHist = new TH1F("lamZResolution",
                "Resolution of V0 decay vertex Z coordinate", 50, 0., 50.);
  ksAbsoluteDistResolutionHist = new TH1F("ksRResolution",
                      "Resolution of absolute distance from primary vertex to V0 vertex",
                      100, 0., 50.);
  lamAbsoluteDistResolutionHist = new TH1F("lamRResolution",
                       "Resolution of absolute distance from primary vertex to V0 vertex",
                       100, 0., 50.);

  ksCandStatusHist = new TH1F("ksCandStatus",
              "Fake type by cand status",
              10, 0., 10.);
  lamCandStatusHist = new TH1F("ksCandStatus",
              "Fake type by cand status",
              10, 0., 10.);

  double minKsMass = 0.49767 - 0.07;
  double maxKsMass = 0.49767 + 0.07;
  double minLamMass = 1.1156 - 0.05;
  double maxLamMass = 1.1156 + 0.05;
  fakeKsMassHisto = new TH1F("ksMassFake",
                 "Mass of fake K0s",
                 100, minKsMass, maxKsMass);
  goodKsMassHisto = new TH1F("ksMassGood",
                 "Mass of good reco K0s",
                 100, minKsMass, maxKsMass);
  fakeLamMassHisto = new TH1F("lamMassFake",
                  "Mass of fake Lambda",
                  100, minLamMass, maxLamMass);
  goodLamMassHisto = new TH1F("lamMassGood",
                  "Mass of good Lambda",
                  100, minLamMass, maxLamMass);

  ksFakeDauRadDistHisto = new TH1F("radDistFakeKs",
                   "Production radius of daughter particle of Ks fake",
                   100, 0., 15.);
  lamFakeDauRadDistHisto = new TH1F("radDistFakeLam",
                    "Production radius of daughter particle of Lam fake",
                    100, 0., 15.);*/


  //std::cout << "Histograms booked" << std::endl;

  /*ksEffVsRHist->Sumw2();
  ksEffVsEtaHist->Sumw2();
  ksEffVsPtHist->Sumw2();
  ksTkEffVsRHist->Sumw2();
  ksTkEffVsEtaHist->Sumw2();
  ksTkEffVsPtHist->Sumw2();
  ksFakeVsRHist->Sumw2();
  ksFakeVsEtaHist->Sumw2();
  ksFakeVsPtHist->Sumw2();
  ksTkFakeVsRHist->Sumw2();
  ksTkFakeVsEtaHist->Sumw2();
  ksTkFakeVsPtHist->Sumw2();

  lamEffVsRHist->Sumw2();
  lamEffVsEtaHist->Sumw2();
  lamEffVsPtHist->Sumw2();
  lamTkEffVsRHist->Sumw2();
  lamTkEffVsEtaHist->Sumw2();
  lamTkEffVsPtHist->Sumw2();
  lamFakeVsRHist->Sumw2();
  lamFakeVsEtaHist->Sumw2();
  lamFakeVsPtHist->Sumw2();
  lamTkFakeVsRHist->Sumw2();
  lamTkFakeVsEtaHist->Sumw2();
  lamTkFakeVsPtHist->Sumw2();*/

}

void V0Validator::endRun ( const edm::Run &  iRun, const edm::EventSetup &  iSetup  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 1196 of file V0Validator.cc.

References DQMStore::save(), theDQMRootFileName, and theDQMstore.

                                                                        {
  //theDQMstore->showDirStructure();
  if(theDQMRootFileName.size() && theDQMstore) {
    theDQMstore->save(theDQMRootFileName);
  }
}

Member Data Documentation

std::string V0Validator::dirName

private

Definition at line 275 of file V0Validator.h.

Referenced by beginRun().

edm::EDGetTokenT<edm::SimTrackContainer> V0Validator::edmSimTrackContainerToken_

private

Definition at line 280 of file V0Validator.h.

Referenced by analyze().

edm::EDGetTokenT<edm::SimVertexContainer> V0Validator::edmSimVertexContainerToken_

private

Definition at line 281 of file V0Validator.h.

Referenced by analyze().

edm::EDGetTokenT< edm::View<reco::Track> > V0Validator::edmView_recoTrack_Token_

private

Definition at line 279 of file V0Validator.h.

Referenced by analyze().

MonitorElement* V0Validator::fakeKsMass

private

Definition at line 261 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::fakeLamMass

private

Definition at line 263 of file V0Validator.h.

Referenced by analyze(), and beginRun().

int V0Validator::genK0s

private

Definition at line 109 of file V0Validator.h.

Referenced by analyze(), and V0Validator().

int V0Validator::genLam

private

Definition at line 109 of file V0Validator.h.

Referenced by analyze(), and V0Validator().

MonitorElement* V0Validator::goodKsMass

private

Definition at line 262 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::goodLamMass

private

Definition at line 264 of file V0Validator.h.

Referenced by analyze(), and beginRun().

float V0Validator::K0sCandEta

private

Definition at line 104 of file V0Validator.h.

Referenced by analyze().

int V0Validator::k0sCandFound

private

Definition at line 110 of file V0Validator.h.

Referenced by V0Validator().

float V0Validator::K0sCandpT

private

Definition at line 104 of file V0Validator.h.

Referenced by analyze().

float V0Validator::K0sCandR

private

Definition at line 104 of file V0Validator.h.

Referenced by analyze().

unsigned int V0Validator::K0sCandStatus

private

Definition at line 105 of file V0Validator.h.

Referenced by analyze().

float V0Validator::K0sGenEta

private

Definition at line 102 of file V0Validator.h.

Referenced by analyze().

float V0Validator::K0sGenpT

private

Definition at line 102 of file V0Validator.h.

Referenced by analyze().

float V0Validator::K0sGenR

private

Definition at line 102 of file V0Validator.h.

Referenced by analyze().

unsigned int V0Validator::K0sGenStatus

private

Definition at line 105 of file V0Validator.h.

Referenced by analyze().

unsigned int V0Validator::K0sPiCandStatus[2]

private

Definition at line 106 of file V0Validator.h.

Referenced by analyze().

unsigned int V0Validator::K0sPiEff[2]

private

Definition at line 106 of file V0Validator.h.

Referenced by analyze().

int V0Validator::k0sTracksFound

private

Definition at line 110 of file V0Validator.h.

Referenced by analyze().

MonitorElement* V0Validator::ksAbsoluteDistResolution

private

Definition at line 249 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksCandStatus

private

Definition at line 258 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksEffVsEta

private

Definition at line 181 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksEffVsEta_denom

private

Definition at line 210 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksEffVsEta_num

private

Definition at line 194 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksEffVsPt

private

Definition at line 182 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksEffVsPt_denom

private

Definition at line 211 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksEffVsPt_num

private

Definition at line 195 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksEffVsR

private

Definition at line 180 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksEffVsR_denom

private

Definition at line 209 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksEffVsR_num

private

Definition at line 193 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksFakeDauRadDist

private

Definition at line 270 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksFakeVsEta

private

Definition at line 187 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksFakeVsEta_denom

private

Definition at line 207 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksFakeVsEta_num

private

Definition at line 200 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksFakeVsPt

private

Definition at line 188 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksFakeVsPt_denom

private

Definition at line 208 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksFakeVsPt_num

private

Definition at line 201 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksFakeVsR

private

Definition at line 186 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksFakeVsR_denom

private

Definition at line 206 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksFakeVsR_num

private

Definition at line 199 of file V0Validator.h.

Referenced by analyze(), and beginRun().

float V0Validator::KsGenX

private

Definition at line 103 of file V0Validator.h.

float V0Validator::KsGenY

private

Definition at line 103 of file V0Validator.h.

float V0Validator::KsGenZ

private

Definition at line 103 of file V0Validator.h.

MonitorElement* V0Validator::ksMassAll

private

Definition at line 266 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksTkEffVsEta

private

Definition at line 184 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksTkEffVsEta_num

private

Definition at line 197 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksTkEffVsPt

private

Definition at line 185 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksTkEffVsPt_num

private

Definition at line 198 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksTkEffVsR

private

Definition at line 183 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksTkEffVsR_num

private

Definition at line 196 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksTkFakeVsEta

private

Definition at line 190 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksTkFakeVsEta_num

private

Definition at line 203 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksTkFakeVsPt

private

Definition at line 191 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksTkFakeVsPt_num

private

Definition at line 204 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksTkFakeVsR

private

Definition at line 189 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksTkFakeVsR_num

private

Definition at line 202 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::ksXResolution

private

Definition at line 246 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksYResolution

private

Definition at line 247 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::ksZResolution

private

Definition at line 248 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamAbsoluteDistResolution

private

Definition at line 253 of file V0Validator.h.

Referenced by beginRun().

float V0Validator::LamCandEta

private

Definition at line 104 of file V0Validator.h.

Referenced by analyze().

int V0Validator::lamCandFound

private

Definition at line 110 of file V0Validator.h.

Referenced by V0Validator().

float V0Validator::LamCandpT

private

Definition at line 104 of file V0Validator.h.

Referenced by analyze().

float V0Validator::LamCandR

private

Definition at line 104 of file V0Validator.h.

Referenced by analyze().

unsigned int V0Validator::LamCandStatus

private

Definition at line 105 of file V0Validator.h.

Referenced by analyze().

MonitorElement* V0Validator::lamCandStatus

private

Definition at line 259 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamEffVsEta

private

Definition at line 221 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamEffVsEta_denom

private

Definition at line 217 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamEffVsEta_num

private

Definition at line 234 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamEffVsPt

private

Definition at line 222 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamEffVsPt_denom

private

Definition at line 218 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamEffVsPt_num

private

Definition at line 235 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamEffVsR

private

Definition at line 220 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamEffVsR_denom

private

Definition at line 216 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamEffVsR_num

private

Definition at line 233 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamFakeDauRadDist

private

Definition at line 271 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamFakeVsEta

private

Definition at line 227 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamFakeVsEta_denom

private

Definition at line 214 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamFakeVsEta_num

private

Definition at line 240 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamFakeVsPt

private

Definition at line 228 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamFakeVsPt_denom

private

Definition at line 215 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamFakeVsPt_num

private

Definition at line 241 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamFakeVsR

private

Definition at line 226 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamFakeVsR_denom

private

Definition at line 213 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamFakeVsR_num

private

Definition at line 239 of file V0Validator.h.

Referenced by analyze(), and beginRun().

float V0Validator::LamGenEta

private

Definition at line 102 of file V0Validator.h.

Referenced by analyze().

float V0Validator::LamGenpT

private

Definition at line 102 of file V0Validator.h.

Referenced by analyze().

float V0Validator::LamGenR

private

Definition at line 102 of file V0Validator.h.

Referenced by analyze().

unsigned int V0Validator::LamGenStatus

private

Definition at line 105 of file V0Validator.h.

Referenced by analyze().

float V0Validator::LamGenX

private

Definition at line 103 of file V0Validator.h.

float V0Validator::LamGenY

private

Definition at line 103 of file V0Validator.h.

float V0Validator::LamGenZ

private

Definition at line 103 of file V0Validator.h.

MonitorElement* V0Validator::lamMassAll

private

Definition at line 267 of file V0Validator.h.

Referenced by analyze(), and beginRun().

unsigned int V0Validator::LamPiCandStatus[2]

private

Definition at line 106 of file V0Validator.h.

Referenced by analyze().

unsigned int V0Validator::LamPiEff[2]

private

Definition at line 106 of file V0Validator.h.

Referenced by analyze().

MonitorElement* V0Validator::lamTkEffVsEta

private

Definition at line 224 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamTkEffVsEta_num

private

Definition at line 237 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamTkEffVsPt

private

Definition at line 225 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamTkEffVsPt_num

private

Definition at line 238 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamTkEffVsR

private

Definition at line 223 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamTkEffVsR_num

private

Definition at line 236 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamTkFakeVsEta

private

Definition at line 230 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamTkFakeVsEta_num

private

Definition at line 243 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamTkFakeVsPt

private

Definition at line 231 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamTkFakeVsPt_num

private

Definition at line 244 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::lamTkFakeVsR

private

Definition at line 229 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamTkFakeVsR_num

private

Definition at line 242 of file V0Validator.h.

Referenced by analyze(), and beginRun().

int V0Validator::lamTracksFound

private

Definition at line 110 of file V0Validator.h.

MonitorElement* V0Validator::lamXResolution

private

Definition at line 250 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamYResolution

private

Definition at line 251 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::lamZResolution

private

Definition at line 252 of file V0Validator.h.

Referenced by beginRun().

MonitorElement* V0Validator::nKs

private

Definition at line 255 of file V0Validator.h.

Referenced by analyze(), and beginRun().

MonitorElement* V0Validator::nLam

private

Definition at line 256 of file V0Validator.h.

Referenced by analyze(), and beginRun().

int V0Validator::noTPforK0sCand

private

Definition at line 110 of file V0Validator.h.

Referenced by analyze(), and V0Validator().

int V0Validator::noTPforLamCand

private

Definition at line 110 of file V0Validator.h.

Referenced by analyze(), and V0Validator().

int V0Validator::realK0sFound

private

Definition at line 109 of file V0Validator.h.

Referenced by analyze(), and V0Validator().

int V0Validator::realK0sFoundEff

private

Definition at line 109 of file V0Validator.h.

Referenced by analyze(), and V0Validator().

int V0Validator::realLamFound

private

Definition at line 109 of file V0Validator.h.

Referenced by analyze(), and V0Validator().

int V0Validator::realLamFoundEff

private

Definition at line 109 of file V0Validator.h.

Referenced by analyze(), and V0Validator().

edm::EDGetTokenT<reco::RecoToSimCollection> V0Validator::recoRecoToSimCollectionToken_

private

Definition at line 276 of file V0Validator.h.

Referenced by analyze().

edm::EDGetTokenT<reco::SimToRecoCollection> V0Validator::recoSimToRecoCollectionToken_

private

Definition at line 277 of file V0Validator.h.

Referenced by analyze().

edm::EDGetTokenT<reco::VertexCompositeCandidateCollection> V0Validator::recoVertexCompositeCandidateCollection_k0s_Token_

private

Definition at line 283 of file V0Validator.h.

Referenced by analyze().

edm::EDGetTokenT<reco::VertexCompositeCandidateCollection> V0Validator::recoVertexCompositeCandidateCollection_lambda_Token_

private

Definition at line 283 of file V0Validator.h.

Referenced by analyze().

std::string V0Validator::theDQMRootFileName

private

Definition at line 274 of file V0Validator.h.

Referenced by endRun().

DQMStore* V0Validator::theDQMstore

private

Definition at line 178 of file V0Validator.h.

Referenced by beginRun(), endRun(), and V0Validator().

edm::EDGetTokenT<TrackingParticleCollection> V0Validator::trackingParticleCollection_Eff_Token_

private

Definition at line 278 of file V0Validator.h.

Referenced by analyze().

edm::EDGetTokenT<TrackingParticleCollection> V0Validator::trackingParticleCollectionToken_

private

Definition at line 278 of file V0Validator.h.

Referenced by analyze().

edm::EDGetTokenT< std::vector<reco::Vertex> > V0Validator::vec_recoVertex_Token_

private

Definition at line 282 of file V0Validator.h.

Referenced by analyze().